Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Deformation behavior and microstructure effect in 2124Al/SiCp composite

Dynamic compression tests were performed by means of a Split Hopkinson Pressure Bar (SHPB). Test materials were 2124Al alloys reinforced with 17% volume fraction of 3, 13 and 37 μm SiC particles, respectively. Under strain rate ε = 2100 l/s, SiC particles have a strong effect on σ_0.2 of the composites and the σ_0.2 increases with different SiC particle size in the following order: 2124Al-alloy → 124Al/SiC_p (37 μm) → 2124Al/SiC_p (13 μm) → 2124Al/SiC_p (3 μm), and the strain hardening of the composites depends mainly on the strain hardening of matrix, 2124A1 alloy. The results of dimensional analysis present that the flow stress of these composites not only depends on the property of reinforcement and matrix but also relates to the microstructure scale, matrix grain size, reinforcement size, the distance between reinforcements and dislocations in matrix. The normalized flow stress here is a function of inverse power of the edge-edge particle spacing, dislocation density and matrix grain size. Close-up observation shows that, in the composite containing SiC particles (3 μm), localized deformation formed readily comparing with other materials under the same loading condition. Microscopic observations indicate that different plastic flow patterns occur within the matrix due to the presence of hard particles with different sizes.

Influence of Indenter Tip Roundness on Hardness Behavior in Nanoindentation

In this paper. the effect of indenter tip roundness on hardness behavior for two typical elastic perfectly plastic materials is studied by means of finite element simulation. A rigid conical indenter of semi apex angle 70.3 degrees fitted smoothly with a spherical tip is employed. It is shown that as the indentation depth increases hardness first rises from zero, reaches a maximum and then decreases slowly approaching asymptotically the limiting value equal to that due to a conical indenter of ideally sharp tip. The range within which hardness varies appreciably is comparable to the radius of the indenter tip. The difference between the maximum value and the limiting value depends on the yield stress over the Young's modulus ratio. The smaller this ratio the greater the difference is. Numerical simulation also provides an opportunity for checking the accuracy and limitations of the widely used Oliver-Pharr method.

Instability in threshold voltage and subthreshold behavior in Hf-In-Zn-O thin film transistors induced by bias-and light-stress

Electrical bias and light stressing followed by natural recovery of amorphous hafnium-indium-zinc-oxide (HIZO) thin film transistors with a silicon oxide/nitride dielectric stack reveals defect density changes, charge trapping and persistent photoconductivity (PPC). In the absence of light, the polarity of bias stress controls the magnitude and direction of the threshold voltage shift (Δ VT), while under light stress, VT consistently shifts negatively. In all cases, there was no significant change in field-effect mobility. Light stress gives rise to a PPC with wavelength-dependent recovery on time scale of days. We observe that the PPC becomes more pronounced at shorter wavelengths. © 2010 American Institute of Physics.

A time domain computation method for dynamic behavior of mooring system

A quasi-steady time domain method is developed for the prediction of dynamic behavior of a mooring system under the environmental disturbances, such as regular or irregular waves, winds and currents. The mooring forces are obtained in a static sense at each instant. The dynamic feature of the mooring cables can be obtained by incorporating the extended 3-D lumped-mass method with the known ship motion history. Some nonlinear effects, such as the influence of the instantaneous change of the wetted hull surface on the hydrostatic restoring forces and Froude-Krylov forces, are included. The computational results show a satisfactory agreement with the experimental ones.

Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

Numerical Investigation On Evolution Of Cylindrical Cellular Detonation

Cylindrical cellular detonation is numerically investigated by solving two-dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.

Simulations of Mechanical Behavior of Polycrystalline Copper with Nano-Twins

Mechanical behavior and microstructure evolution of polycrystalline copper with nano-twins were investigated in the present work by finite element simulations. The fracture of grain boundaries are described by a cohesive interface constitutive model based on the strain gradient plasticity theory. A systematic study of the strength and ductility for different grain sizes and twin lamellae distributions is performed. The results show that the material strength and ductility strongly depend on the grain size and the distribution of twin lamellae microstructures in the polycrystalline copper.

Novel Mechanical Behavior Of Zno Nanorods

A novel stress-strain relation with two stages of linear elastic deformation is observed in [0 0 0 1]-oriented ZnO nanorods under uniaxial tensile loading. This phenomenon results from a phase transformation from wurtzite (WZ, P6(3)mc space group) to a body-centered tetragonal structure with four-atom rings (denoted as BCT-4) belonging to the P4(2)/mnm space group. The analysis here focuses on the effects of nanorod size and temperature on the phase transformation and the associated mechanical behavior. It is found that as size is increased from 19.5 to 45.5 angstrom, the critical stress for nucleation of the transformation decreases by 25% from 21.90 to 16.50 GPa and the elastic moduli of the WZ- and BCT-4-structured nanorods decrease by 24% (from 299.49 to 227.51 GPa) and 38% (from 269.29 to 166.86 GPa), respectively. A significant temperature effect is also observed, with the critical stress for transformation initiation decreasing 87.8% from 17.89 to 2.19 GPa as temperature increases from 300 to 1500 K. (c) 2007 Elsevier B.V. All rights reserved.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Growth Mechanism And Corrosion Behavior Of Ceramic Coatings On Aluminum Produced By Autocontrol Ac Pulse Peo

Ceramic coatings are produced on aluminum alloy by autocontrol AC pulse Plasma Electrolytic Oxidation (PEO) with stabilized average current. Transient signal gathering system is used to study the current, voltage, and the transient wave during the PEO process. SEM, OM, XRD and EDS are used to study the coatings evolution of morphologies, composition and structure. TEM is used to study the micro profile of the outer looser layer and inner compact layer. Polarization test is used to study the corrosion property of PEO coatings in NaCl solution. According to the test results, AC pulse PEO process can be divided into four stages with different aspects of discharge phenomena, voltage and current. The growth mechanism of AC PEO coating is characterized as anodic reaction and discharge sintering effect. PEO coating can increase the corrosion resistance of aluminum alloy by one order or two; however, too long process time is not necessarily needed to increase the corrosion resistance. In condition of this paper, PEO coating at 60 min is the most protective coating for aluminum alloy substrate. (C) 2008 Elsevier B.V. All fights reserved.