Asymptotic analysis for a crack on interface of damaged materials

An asymptotic analysis for a crack lying on the interface of a damaged plastic material and a linear elastic material is presented in this paper. The present results show that the stress distributions along the crack tip are quite similar to those with HRR singularity field and the crack faces open obviously. Material constants n, mu and mo are varied to examine their effects on the resulting stress distributions and displacement distributions in the damaged plastic region. It is found that the stress components sigma(rr), sigma(theta theta), sigma(r theta) and sigma(e) are slightly affected by the changes of material constants n, mu and m(0), but the damaged plastic region are greatly disturbed by these material parameters.

Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

Prevalence Of Shear Banding In Compression Of Zr41Ti14Cu12.5Ni10Be22.5 Pillars As Small As 150 Nm In Diameter

Recently, the size dependence of mechanical behaviors, particularly the yield strength and plastic deformation mode, of bulk metallic glasses (BMG) has created a great deal of interest. Contradicting conclusions have been drawn by different research groups, based on various experiments on different BMG systems. Based on in situ compression transmission electron microscopy (TEM) experiments on Zr41Ti14Cu12.5Ni10Be22.5 (Vit 1) nanopillars, this paper provides strong evidence that shear banding still prevails at specimen length scales as small as 150 nm in diameter. This is supported by in situ and ex situ images of shear bands, and by the carefully recorded displacement bursts under load control its well as load drops under displacement control. Finite element modeling of the stress state within the pillar shows that the unavoidable geometry constraints accompanying such experiments impart a strong effect on the experimental results, including non-uniform stress distributions and high level hydrostatic pressures. The seemingly improved compressive ductility is believed to be due to such geometry constraints. Observations underscore the notion that the mechanical behavior of metallic glasses, including strength and plastic deformation mode, is size independent at least in Vit 1. (C) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Finite-Element Simulations of Earthquakes

This thesis discusses simulations of earthquake ground motions using prescribed ruptures and dynamic failure. Introducing sliding degrees of freedom led to an innovative technique for numerical modeling of earthquake sources. This technique allows efficient implementation of both prescribed ruptures and dynamic failure on an arbitrarily oriented fault surface. Off the fault surface the solution of the three-dimensional, dynamic elasticity equation uses well known finite-element techniques. We employ parallel processing to efficiently compute the ground motions in domains containing millions of degrees of freedom.

Using prescribed ruptures we study the sensitivity of long-period near-source ground motions to five earthquake source parameters for hypothetical events on a strike-slip fault (M_w 7.0 to 7.1) and a thrust fault (M_w 6.6 to 7.0). The directivity of the ruptures creates large displacement and velocity pulses in the ground motions in the forward direction. We found a good match between the severity of the shaking and the shape of the near-source factor from the 1997 Uniform Building Code for strike-slip faults and thrust faults with surface rupture. However, for blind thrust faults the peak displacement and velocities occur up-dip from the region with the peak near-source factor. We assert that a simple modification to the formulation of the near-source factor improves the match between the severity of the ground motion and the shape of the near-source factor.

For simulations with dynamic failure on a strike-slip fault or a thrust fault, we examine what constraints must be imposed on the coefficient of friction to produce realistic ruptures under the application of reasonable shear and normal stress distributions with depth. We found that variation of the coefficient of friction with the shear modulus and the depth produces realistic rupture behavior in both homogeneous and layered half-spaces. Furthermore, we observed a dependence of the rupture speed on the direction of propagation and fluctuations in the rupture speed and slip rate as the rupture encountered changes in the stress field. Including such behavior in prescribed ruptures would yield more realistic ground motions.

板条热容激光器的二维热特性

固体热容激光器（SSHCL）作为高功率固体激光器的一个重要发展方向,引起人们广泛关注。数值模拟激光介质板条在热容方式下工作的温度和应力分布是了解该类激光器工作特性的一种有效手段,采用平面应力近似法导出了半导体激光器抽运热容激光介质板的二维温度和应力分布公式,同时也对二维抽运光吸收密度、介质板温度分布和折射率变化进行了分析与讨论。数值计算的结果表明二维效应的温度分布和应力分布要比一维效应给出的分布更均匀。

Avaliação de ligações entre perfis tubulares em aço submetidas a esforços de flexão.

Os perfis tubulares sem costura são largamente utilizados em diversos países, porém, no Brasil, o uso desses perfis na construção civil era bastante limitado, restringindo-se praticamente a coberturas espaciais. Considerando uma nova realidade para o uso de perfis tubulares, este trabalho apresenta uma análise de ligações tipo T com perfis tubulares quadrados (SHS) para o banzo e para o montante efetuada com base na norma europeia, Eurocode 3, no CIDECT, na NBR 16239:2013 e ISO 14346 através de um modelo em elementos finitos desenvolvido no programa Ansys. Verificou-se a influência do momento atuante no montante no comportamento global das ligações. As não-linearidades físicas e geométricas foram incorporadas aos modelos, a fim de se mobilizar totalmente a capacidade resistente desta ligação. A não-linearidade do material foi considerada através do critério de plastificação de von Mises através da lei constitutiva tensão versus deformação trilinear de forma a exibir um comportamento elasto-plástico com encruamento. A nãolinearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizado. Os resultados numéricos são avaliados para a ligação em estudo quanto aos modos de falha e a distribuição de tensões. A análise dos momentos resistentes obtidos em comparação com os resultados do modelo numérico, apresentou valores excessivamente a favor da segurança no cálculo utilizando as equações de dimensionamento. Um estudo para fatores de correção das equações de dimensionamento é proposto.

Numerical investigation of the active role of the actin cytoskeleton in the compression resistance of cells.

Numerous in-vitro studies have established that cells react to their physical environment and to applied mechanical loading. However, the mechanisms underlying such phenomena are poorly understood. Previous modelling of cell compression considered the cell as a passive homogenous material, requiring an artificial increase in the stiffness of spread cells to replicate experimentally measured forces. In this study, we implement a fully 3D active constitutive formulation that predicts the distribution, remodelling, and contractile behaviour of the cytoskeleton. Simulations reveal that polarised and axisymmetric spread cells contain stress fibres which form dominant bundles that are stretched during compression. These dominant fibres exert tension; causing an increase in computed compression forces compared to round cells. In contrast, fewer stress fibres are computed for round cells and a lower resistance to compression is predicted. The effect of different levels of cellular contractility associated with different cell phenotypes is also investigated. Highly contractile cells form more dominant circumferential stress fibres and hence provide greater resistance to compression. Computed predictions correlate strongly with published experimentally observed trends of compression resistance as a function of cellular contractility and offer an insight into the link between cell geometry, stress fibre distribution and contractility, and cell deformability. Importantly, it is possible to capture the behaviour of both round and spread cells using a given, unchanged set of material parameters for each cell type. Finally, it is demonstrated that stress distributions in the cell cytoplasm and nucleus computed using the active formulation differ significantly from those computed using passive material models.

Fracture behavior of a Zr-metallic glass under dynamic compression

Dynamic planar compressive experiments on a typical tough Zr-BMG (Bulk Metallic Glass) were carried out under impact velocity of 500-600 m/sec and strain rate of 10(6)/s. The fracture surface of samples exhibits different fracture patterns at different parts of the sample. At a corner close to the front loading boundary, fracture patterns from the free edge toward the centre changed from equiaxial veins in microscale to periodic corrugations in nanoscale; in the middle of the sample, the fracture surface contains glazed zones laid out orderly along the same boundary. FEM simulation was performed to investigate the stress distributions in the impacted sample under a short duration impact loading. It has revealed that the fracture patterns changing from the free edge toward the centre were resulted from the fracture modes' changing from the tensile dominant fracture to the shear dominant fracture. Whereas at the middle part of the sample, fracture initiated from several parallel shear bands propagating close to the same boundary is due to a large strain or much higher shear stress in this area.

INVESTIGATING A NEAR-BED SUBMARINE PIPELINE IN A CURRENT

This paper considers the lift forces acting on a pipeline with a small gap between the pipeline and the plane bottom or scoring bottom. A more reasonable fluid force on the pipeline has been obtained by applying the knowledge of modified potential theory (MPT), which includes the influences of the downstream wake. By finite element method, an iteration procedure is used to solve problems of the nonlinear fluid-structure interaction. Comparing the deflection and the stress distributions with the difference sea bottoms, the failure patterns of a spanning pipeline have been discussed. The results are essential for engineers to assess pipeline stability.

Computer simulation of the damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects

The damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects was studied via a Monte Carlo technique combined with a finite element method. A finite element model was constructed to predict the effects of various matrix defect shapes on the stress distributions. The results indicated that a small matrix defect had almost no effect on fiber stress distributions other than interfacial shear stress distributions. Then, a finite element model with a statistical distribution of the fiber strength was constructed to investigate the influences of the spatial distribution and the volume fraction of matrix defects on composite failure. The results showed that it was accurate to use the shear-lag models and Green's function methods to predict the tensile strength of composites even though the axial stresses in the matrix were neglected.

Numerical investigation of the effects of drill geometry on drilling induced delamination of carbon fiber reinforced composites

Drilling is a major process in the manufacturing of holes required for the assemblies of composite laminates in aerospace industry. Simulation of drilling process is an effective method in optimizing the drill geometry and process parameters in order to improve hole quality and to reduce the drill wear. In this research we have developed three-dimensional (3D) FE model for drilling CFRP. A 3D progressive intra-laminar failure model based on the Hashin's theory is considered. Also an inter-laminar delamination model which includes the onset and growth of delamination by using cohesive contact zone is developed. The developed model with inclusion of the improved delamination model and real drill geometry is used to make comparison between the step drill of different stage ratio and twist drill. Thrust force, torque and work piece stress distributions are estimated to decrease by the use of step drill with high stage ratio. The model indicates that delamination and other workpiece defects could be controlled by selection of suitable step drill geometry. Hence the 3D model could be used as a design tool for drill geometry for minimization of delamination in CFRP drilling. © 2013 Elsevier Ltd.

Evaluation of sources of uncertainties in microtensile bond strength of dental adhesive system for different specimen geometries

Objective: The aim of this research is to use finite element analysis (FEA) to quantify the effect of the sample shape and the imperfections induced during the manufacturing process of samples on the bond strength and modes of failure of dental adhesive systems through microtensile test. Using the FEA prediction for individual parameters effect, estimation of expected variation and spread of the microtensile bond strength results for different sample geometries is made. Methods: The estimated stress distributions for three different sample shapes, hourglass, stick and dumbbell predicted by FEA are used to predict the strength for different fracture modes. Parameters such as the adhesive thickness, uneven interface of the adhesive and composite and dentin, misalignment of axis of loading, the existence of flaws such as induced cracks during shaping the samples or bubbles created during application of the adhesive are considered. Microtensile experiments are performed simultaneously to measure bond strength and modes of failure. These are compared with the FEA results. Results: The relative bonding strength and its standard deviation for the specimens with different geometries measured through the microtensile tests confirm the findings of the FEA. The hourglass shape samples show lower tensile bond strength and standard deviation compared to the stick and dumbbell shape samples. ANOVA analysis confirms no significant difference between dumbbell and stick geometry results, and major differences of these two geometries compared to hourglass shape measured values. Induced flaws in the adhesive and misalignment of the angle of application of load have significant effect on the microtensile bond strength. Using adhesive with higher modulus the differences between the bond strength of the three sample geometries increase. Significance: The result of the research clarifies the importance of the sample geometry chosen in measuring the bond strength. It quantifies the effect of the imperfections on the bond strength for each of the sample geometries through a systematic and all embracing study. The results explain the reasons of the large spread of the microtensile test results reported by various researchers working in different labs and the need for standardization of the test method and sample shape used in evaluation of the dentin-adhesive bonding system. © 2007 Academy of Dental Materials.

Tool life and hole surface integrity studies for hole-making of Ti6Al4V Alloy

With a significant growth in the use of titanium alloys in the aviation manufacturing industry, the key challenge of making high-quality holes in the aircraft assembly process needs to be addressed. In this work, case studies deploying traditional drilling and helical milling technologies are carried out to investigate the tool life and hole surface integrity for hole-making of titanium alloy. Results show that the helical milling process leads to much longer tool life, generally lower hole surface roughness, and higher hole subsurface microhardness. In addition, no plastically deformed layer or white layer has been observed in holes produced by helical milling. In contrast, a slightly softened region was always present on the drilled surface. The residual stress distributions within the hole surface, including compressive and tensile residual stress, have also been investigated in detail.