In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method （F2LFEM）. The general solutio...In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method （F2LFEM）. The general solution of an orthotropic crack problem is obtained by assimilating the problem with isotropic crack problem, and is employed as the global interpolation function in F2LFEM. In the neighborhood of crack tip of the crack plate, the fractal geometry concept is introduced to achieve the similar meshes having similarity ratio less than one and generate an infinitesimal mesh so that the relationship between the stiffness matrices of two adjacent layers is equal. A large number of degrees of freedom around the crack tip are transformed to a small set of generalized coordinates. Numerical examples show that this method is efficient and accurate in evaluating the stress intensity factor （SIF）.展开更多
For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions （or between the derivatives of eig...For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions （or between the derivatives of eigenfunctions） of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors （SIFs）, so any order SIF can be extracted based on a certain known solution of displacement （an analytic result or a numerical result）. Many numerical examples based on the finite element method of lines （FEMOL） show that the present method is very powerful and efficient.展开更多
To forecast the modification and evolution of a residual stress field in a butt-welded plate, the residual stress field was numerically evaluated by a finite element method code. The plate was originated by welding an...To forecast the modification and evolution of a residual stress field in a butt-welded plate, the residual stress field was numerically evaluated by a finite element method code. The plate was originated by welding and suffered after chip- forming machining. The residual stress field was taken as a pre-stress condition for the next cutting simulations. Residual stress distribution after machining was then compared to that in the initial residual stress field. Numerical results have been compared to experimental measurements. That shows the potentiality as well as the limitations of numerical techniques. Three major contributions are summarized as follows ： Longitudinal residual stress distribution in welded plates is deeply changed by mechanical tooling; Planing and cutting increases peak values of transversal residual stress and in general introduces some unevenness in distribution along X-direction; Comparison of experimental and numerical values is generally satisfactory for longitudinal residual stress.展开更多
Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are o...Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.展开更多
To predict the behavior of geogrids embedded in sand under pullout loading conditions, the two dimensional plane-stress finite element model was presented. The interactions between soil and geogrid were simulated as ...To predict the behavior of geogrids embedded in sand under pullout loading conditions, the two dimensional plane-stress finite element model was presented. The interactions between soil and geogrid were simulated as non-linear springs, and the stiffness of the springs was determined from simple tests in the specially designed pullout box. The predicted behavior of the geogrid under pullout load agrees well with the observed data including the load-displacement properties, the displacement distribution along the longitudinal direction and the mobilization of the frictional and bearing resistance. (Edited author abstract) 8 Refs.展开更多
文摘In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method （F2LFEM）. The general solution of an orthotropic crack problem is obtained by assimilating the problem with isotropic crack problem, and is employed as the global interpolation function in F2LFEM. In the neighborhood of crack tip of the crack plate, the fractal geometry concept is introduced to achieve the similar meshes having similarity ratio less than one and generate an infinitesimal mesh so that the relationship between the stiffness matrices of two adjacent layers is equal. A large number of degrees of freedom around the crack tip are transformed to a small set of generalized coordinates. Numerical examples show that this method is efficient and accurate in evaluating the stress intensity factor （SIF）.
基金Project supported by the National Natural Science Foundation of China (Nos. 59525813 and 19872066).
文摘For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions （or between the derivatives of eigenfunctions） of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors （SIFs）, so any order SIF can be extracted based on a certain known solution of displacement （an analytic result or a numerical result）. Many numerical examples based on the finite element method of lines （FEMOL） show that the present method is very powerful and efficient.
文摘To forecast the modification and evolution of a residual stress field in a butt-welded plate, the residual stress field was numerically evaluated by a finite element method code. The plate was originated by welding and suffered after chip- forming machining. The residual stress field was taken as a pre-stress condition for the next cutting simulations. Residual stress distribution after machining was then compared to that in the initial residual stress field. Numerical results have been compared to experimental measurements. That shows the potentiality as well as the limitations of numerical techniques. Three major contributions are summarized as follows ： Longitudinal residual stress distribution in welded plates is deeply changed by mechanical tooling; Planing and cutting increases peak values of transversal residual stress and in general introduces some unevenness in distribution along X-direction; Comparison of experimental and numerical values is generally satisfactory for longitudinal residual stress.
基金The project supported by the National Natural Science Foundation of China （50275073）
文摘Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.
文摘To predict the behavior of geogrids embedded in sand under pullout loading conditions, the two dimensional plane-stress finite element model was presented. The interactions between soil and geogrid were simulated as non-linear springs, and the stiffness of the springs was determined from simple tests in the specially designed pullout box. The predicted behavior of the geogrid under pullout load agrees well with the observed data including the load-displacement properties, the displacement distribution along the longitudinal direction and the mobilization of the frictional and bearing resistance. (Edited author abstract) 8 Refs.
