Three-dimensional solutions for the stress fields in functionally graded cylindrical panel with finite length and subjected to thermal/mechanical loads

Three-dimensional thermo-elastic analysis of a functionally graded cylindrical panel with finite length and subjected to nonuniform mechanical and steady-state thermal loads are carried out in this paper. Thermal and mechanical properties of the functionally graded material are assumed to be temperature independent and continuously vary in the radial direction of the panel. Analytical solutions for the temperature and stress fields expressed in terms of trigonometric and power series for the simply supported boundary conditions are derived and graphically presented.     

Mixed-mode I/II fields around a crack with a cohesive zone ahead of the crack tip

Considering a cohesive zone ahead of the crack tip, the linear elastic crack problem under Mode I, Mode II or mixed-mode conditions is formulated in an elliptic coordinate system, so that the cohesive surfaces are conveniently represented by straight line segments. It is shown that the displacement and stress fields around the crack tip and the cohesive zone, expressed in terms of elliptic coordinates, have a simple mathematical form, which does not contain a stress singularity at the crack tip due to the existence of the cohesive zone.     

Extension of the virtual fields method to elasto-plastic material identification with cyclic loads and kinematic hardening

The virtual fields method (VFM) has been specifically developed for solving inverse problems from dense full-field data. This paper explores recent improvements regarding the identification of elasto-plastic models. The procedure has been extended to cyclic loads and combined kinematic/isotropic hardening. A specific attention has also been given to the effect of noise in the data. Indeed, noise in experimental data may significantly affect the robustness of the VFM for solving such inverse problems. The concept of optimized virtual fields that minimize the noise effects, previously developed for linear elasticity, is extended to plasticity in this study. Numerical examples with models combining isotropic and kinematic hardening have been considered for the validation. Different load paths (tension, compression, notched specimen) have shown that this new procedure is robust when applied to elasto-plastic material identification. Finally, the procedure is validated on experimental data.     

Optimal design of bars under nonuniform tension with respect to mixed creep rupture time

The first attempt of finding of optimal shape for bars in presence of body forces with respect to mixed creep rupture is made. For given volume of the bar, distribution of initial cross-section, ensuring the longest life-time to mixed rupture is sought. The finite strain theory and physical law in form of Norton's law generalized for true stresses and logarithmic strains are applied. Using the method of parametric optimization, the best of linear and quadratic functions describing the initial shape of the bar are found. The shape of initial strength is corrected in a way leading to longer life-time. Results of both approaches are compared.     

The singularity analysis of the stress and strain fields for Mode I fracture in elastic-plastic state

The stress and strain singularities of power hardening material for Mode I fracrure are analysed according to the fundamental equations of elastic-plastic mechanics. It is found that the singularities of all stress and strain components do not change in the thick direction, and neither the six stress components nor the six strain components have the same singularity.     

Local stress distribution and yielding of steel sheets due to concentrated transverse loads

In this paper, a theoretical investigation on local stress distribution and yielding of steel sheets subjected to concentrated transverse loads, caused by bolts, is presented. The stress state produced by the above loading leads to the final extraction of the bolt. For both elastic and elastoplastic cases, it is assumed that only a circular area of the metal sheet is affected by the load of the bolt and hence, one can study a representative circular plate of radius a and thickness t, fixed at its ends with a hole of radius b at its center. In addition, it is assumed that the force acting on the metal sheet by the bolt is applied at the edge of the hole at r =?b.     

Interface damage and its effect on vibrations of slab track under temperature and vehicle dynamic loads

This paper presents a three-dimensional finite element model to investigate the interface damage occurred between prefabricated slab and CA (cement asphalt) mortar layer in the China Railway Track System (CRTS-II) slab track system. In the finite element model, a cohesive zone model with a non-linear constitutive law is introduced and utilized to model the damage, cracking and delamination at the interface. Combining with the temperature field database obtained from the three-dimensional transient heat transfer analysis, the interface damage evolution as a result of temperature change is analyzed. A three-dimensional coupled dynamic model of a vehicle and the slab track is then established to calculate the varying rail-supporting forces which are utilized as the inputs to the finite element model. The non-linearities of the wheel–rail contact geometry, the wheel–rail normal contact force and the wheel–rail tangential creep force are taken into account in the model. Setting the maximum interface damaged state calculated under temperature change as the initial condition, the interface damage evolution and its influence on the dynamic response of the slab track are investigated under the joint action of the temperature change and vehicle dynamic load. The analysis indicates that the proposed model is capable of predicting the initiation and propagation of cracks at the interface. The prefabricated slab presents lateral warping, resulting in severe interface damage on both the sides of the slab track along the longitudinal direction during temperature drop process, while the interface damage level does not change significantly under vehicle dynamic loads. The interface damage has great effects on the dynamic responses of the slab track.     

Age- and direction-related adaptations of lumbar vertebral trabecular bone with respect to apparent stiffness and tissue level stress distribution

The objective of this study was to study the age-related adaptation of lumbar vertebral trabecular bone at the apparent level, as well as the tissue level in three orthogonal directions. Ninety trabecular specimens were obtained from six normal L4 vertebral bodies of six male cadavers in two age groups, three aged 62 years and three aged 69 years, and were scanned using a high-resolution micro-computed tomography （micro-CT） system, then converted to micro- finite element models to do micro-finite element analyses. The relationship between apparent stiffness and bone volume fraction, and the tissue level yon Mises stress distribution for each trabecular specimen when compressed separately in the longitudinal direction, medial-lateral and anterior-posterior directions （transverse directions） were derived and compared between two age groups. The results showed that at the apparent level, trabecular bones from 69-year group had stiffer bone structure relative to their volume fractions in all three directions, and in both age groups, changes in bone volume fraction could explain more variations in apparent stiffness in the longitudinal direction than the transverse directions; at the tissue level, aging had little effect on the tissue von Mises stress distributions for the compressions in all the three directions. The novelty of the present study was that it provided quantitative assessments on the age and direction- related adaptation of Chinese male lumbar vertebral trabecular bone from two different levels： stiffness at the apparent level and stress distribution at the tissue level. It may help to understand the failure mechanisms and fracture risks of vertebral body associated with aging and direction for the prevention of fracture risks in elder individuals.     

Hydride-induced embrittlement and fracture in metals—effect of stress and temperature distribution

A mathematical model for the hydrogen embrittlement of hydride forming metals has been developed. The model takes into account the coupling of the operating physical processes, namely: (i) hydrogen diffusion, (ii) hydride precipitation, (iii) non-mechanical energy flow and (iv) hydride/solid-solution deformation. Material damage and crack growth are also simulated by using de-cohesion model, which takes into account the time variation of energy of de-cohesion, due to the time-dependent process of hydride precipitation. The bulk of the material, outside the de-cohesion layer, is assumed to behave elastically. The hydrogen embrittlement model has been implemented numerically into a finite element framework and tested successfully against experimental data and analytical solutions on hydrogen thermal transport (in: Wunderlich, W. (Ed.), Proceedings of the European Conference on Computational Mechanics, Munich, Germany, 1999, J. Nucl. Mater. (2000a) 279 (2-3) 273). The model has been used for the simulation of Zircaloy-2 hydrogen embrittlement and delayed hydride cracking initiation in (i) a boundary layer problem of a semi-infinite crack, under mode I loading and constant temperature, and (ii) a cracked plate, under tensile stress and temperature gradient. The initial and boundary conditions in case (ii) are those encountered in the fuel cladding of light water reactors, during operation. The effects of near-tip stress intensification as well as of temperature gradient on hydride precipitation and material damage have been studied. The numerical simulation predicts hydride precipitation at a small distance from the crack-tip. When the remote loading is sufficient, the near-tip hydrides fracture. Thus a microcrack is generated, which is separated from the main crack by a ductile ligament, in agreement with experimental observations.     

Transient temperature fields and residual stress fields of metallic materials under welding

Based on the situation of welding thermal conduction and thermo-elasto-plasticityresearch,Ms paper explores some problems in this field.First,the boundary elementmethod for nonlinear problems is improved by linearization of nonlinear problems and usedin welding thermal conduction analysis.Second,the thermo-elasto-plastic finite elementmethod is used for the welding stress calculation,in which the phase transformation isconsidered by the"equivalent linear expansion coefficient method".The comparison of the calculated results with experimental data shows that themethods provided in this paper are available.     

Cohesive zone laws for void growth — II. Numerical field projection of elasto-plastic fracture processes with vapor pressure

Modeling ductile fracture processes using Gurson-type cell elements has achieved considerable success in recent years. However, incorporating the full mechanisms of void growth and coalescence in cohesive zone laws for ductile fracture still remains an open challenge. In this work, a planar field projection method, combined with equilibrium field regularization, is used to extract crack-tip cohesive zone laws of void growth in an elastic-plastic solid. To this end, a single row of void-containing cell elements is deployed directly ahead of a crack in an elastic-plastic medium subjected to a remote K-field loading; the macroscopic behavior of each cell element is governed by the Gurson porous material relation, extended to incorporate vapor pressure effects. A thin elastic strip surrounding this fracture process zone is introduced, from which the cohesive zone variables can be extracted via the planar field projection method. We show that the material's initial porosity induces a highly convex traction-separation relationship — the cohesive traction reaches the peak almost instantaneously and decreases gradually with void growth, before succumbing to rapid softening during coalescence. The profile of this numerically extracted cohesive zone law is consistent with experimentally determined cohesive zone law in Part I for multiple micro-crazing in HIPS. In the presence of vapor pressure, both the cohesive traction and energy are dramatically lowered; the shape of the cohesive zone law, however, remains highly convex, which suggests that diffusive damage is still the governing failure mechanism.     

Interfacial stress analysis and prediction of debonding for a thin plate bonded to a curved substrate

This paper focuses on the analytical and numerical modeling of the interface between a rigid substrate with simple constant curvature and a thin bonded plate. The interfacial behavior is modeled by independent cohesive laws in the normal and tangential directions, coupled with a mixed-mode fracture criterion. The newly developed analytical model determines the interfacial shear and normal stress distributions as functions of the substrate curvature, during the various behavioral stages of the interface prior to the initiation of debonding. The model is also able to predict the debonding load and the effective bond length. In the numerical model the interface is modeled by zero-thickness node-to-segment contact elements, in which both the geometrical relationships between the nodes of the discretized problem and the interface constitutive laws are suitably defined. Numerical results and comparisons between the predictions of the two models are presented.     

First- and second-order sensitivities of beams with respect to cross-sectional cracks

This paper deals with a variational formulation for the sensitivity problem of beam systems in the context of deformable solids with cracks. Natural frequencies are defined as state variables involved in the energy functional of the system, while the cracks length and position are treated as system parameters. The hierarchical equation system is formed and solved for the first and second derivatives of the natural frequency functions with respect to the cracks length and position. An analytical procedure for calculations of the second-order sensitivities of natural frequencies is proposed for the non-symmetrical equation system operator. Numerical algorithms are worked out and implemented computationally. Analytical and numerical aspects of the problem are discussed in detail through a number of illustrative results.The support of this work by the State Committee for Scientific Research (KBN) under Grant No. 4-050-0148/17-98-00 is gratefully acknowledged.     

Non-linear extended thermodynamics of real gases with 6 fields

We establish extended thermodynamics (ET) of real gases with 6 independent fields, i.e., the mass density, the velocity, the temperature and the dynamic pressure, without adopting the near-equilibrium approximation. We prove its compatibility with the universal principles (the entropy principle, the Galilean invariance and the stability), and obtain the symmetric hyperbolic system with respect to the main field. In near-equilibrium we recover the previous results. The correspondence between the ET 6-field theory and Meixner׳s theory of relaxation processes is discovered. The internal variable and the non-equilibrium temperature in Meixner׳s theory are expressed in terms of the quantities of the ET 6-field theory, in particular, the dynamic pressure. As an example, we present the cases of a rarefied polyatomic gas and study the monatomic-gas limit where the system converges to the Euler system of a perfect fluid.     

Dynamic photoelastic study of wave fields in elastic plates with stress concentrators

The method of dynamic photoelasticity is used to study the dynamic failure of structural members in the form of plates with a curvilinear (circular or elliptic) hole and an isolated crack under impulsive loading. The time-dependences of the stress intensity factors and the crack tip velocity are investigated for two types of models __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 12, pp. 84–92, December 2005.