FEM Simulations and Optimization about Residual Stresses in Coating Structures with Functionally Graded Materials Layer

An elasto-plastic finite element method is developed to predict the residual stresses of thermal spraying coatings with functionally graded material layer. In numerical simulations, temperature sensitivity of various material constants is included and mix

Modeling and mesoscopic damage constitutive relation of brittle short-fiber-reinforced composites

Aimed at brittle composites reinforced by randomly distributed short-fibers with a relatively large aspect ratio, stiffness modulus and strength, a mesoscopic material model was proposed. Based on the statistical description, damage mechanisms, damage-induced anisotropy, damage rate effect and stress redistribution, the constitutive relation were derived. By taking glass fiber reinforced polypropylene polymers as an example, the effect of initial orientation distribution of fibers, damage-induced anisotropy, and damage-rate effect on macro-behaviors of composites were quantitatively analyzed. The theoretical predictions compared favorably with the experimental results.

Material response and failure mechanism of unidirectional metal matrix composites under impulsive shear loading

The material response and failure mechanism of unidirectional metal matrix composite under impulsive shear loading are investigated in this paper. Both experimental and analytical studies were performed. The shear strength of unidirectional C-f/A356.0 composite and A356.0 aluminum alloy at high strain rate were measured with a modified split Hopkinson torsional bar technique. The results indicated that the carbon fibers did not improve the shear strength of aluminum matrix if the fiber orientation aligned with the shear loading axis. The microscopic inspection of the fractured surface showed a multi-scale zigzag feature which implied a complicated shear failure mechanism in the composite. In addition to testing, the micromechanical stress field in the composite was analyzed by the generalized Eshelby equivalent method (GEEM). The influence of cracking in matrix on the micromechanical stress field was investigated as well. The results showed that the stress distribution in the composite is quite nonhomogeneous and very high shear stress concentrations are found in some regions in the matrix. The high shear stress concentration in the matrix induces tensile cracking at 45 degrees to the shear direction. This in turn aggravates the stress concentration at the fiber/matrix interface and finally leads to a catastrophic failure in the composite. From the correlation between the analysis and experimental results, the shear failure mechanism of unidirectional C-f/A356.0 composite can be elucidated qualitatively.

Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

Thermocapillary flow in an annular liquid layer painted on a moving fiber

In the present paper, a theoretical model is studied on the flow in the liquid annular film, which is ejected from a vessel with relatively higher temperature and painted on the moving solid fiber. A temperature gradient, driving a thermocapillary flow, is formed on the free surface because of the heat transfer from the liquid with relatively higher temperature to the environmental gas with relatively lower temperature. The thermocapillary flow may change the radii profile of the liquid film. This process analyzed is based on the approximations of lubrication theory and perturbation theory, and the equation of the liquid layer radii and the process of thermal hydrodynamics in the liquid layer are solved for a temperature distribution on the solid fiber.

3D Anisotropic Elastoplastic-Damage Model and its Application in Simulating the Behavior of Rock Materials

A 3D anisotropic elastoplastic-damage model was presented based on continuum damage mechanics theory. In this model, the tensor decomposition technique is employed. Combined with the plastic yield rule and damage evolution, the stress tensor in incremental format is obtained. The derivate eigenmodes in the proposed model are assumed to be related with the uniaxial behavior of the rock material. Each eigenmode has a corresponding damage variable due to the fact that damage is a function of the magnitude of the eigenstrain. Within an eigenmodes, different damage evolution can be used for tensile and compressive loadings. This model was also developed into finite element code in explicit format, and the code was integrated into the well-known computational environment ABAQUS using the ABAQUS/Explicit Solver. Numerical simulation of an uniaxial compressive test for a rock sample is used to examine the performance of the proposed model, and the progressive failure process of the rock sample is unveiled.

A micromechanical damage theory for brittle materials with small cracks

For brittle solids containing numerous small cracks, a micromechanical damage theory is presented which accounts for the interactions between different small cracks and the effect of the boundary of a finite solid, and includes growth of the pre-existing small cracks. The analysis is based on a superposition scheme and series expansions of the complex potentials. The small crack evolution process is simulated through the use of fracture mechanics incorporating appropriate failure criteria. The stress-strain relations are obtained from the micromechanics analysis. Typical examples are given to illustrate the potential capability of the proposed theory. These results show that the present method provides a direct and efficient approach to deal with brittle finite solids containing multiple small cracks. The stress-strain relation curves are evaluated for a rectangular plate containing small cracks.

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.

Vaporization of heated materials into discharge plasmas

The vaporization of condensed materials in contact with high-current discharge plasmas is considered. A kinetic numerical method named direct simulation Monte Carlo (DSMC) and analytical kinetic approaches based on the bimodal distribution function approximation are employed. The solution of the kinetic layer problem depends upon the velocity at the outer boundary of the kinetic layer which varies from very small, corresponding to the high-density plasma near the evaporated surface, up to the sound speed, corresponding to evaporation into vacuum. The heavy particles density and temperature at the kinetic and hydrodynamic layer interface were obtained by the analytical method while DSMC calculation makes it possible to obtain the evolution of the particle distribution function within the kinetic layer and the layer thickness.