On the assumption that the yield criterion of orthotropic materials is isomorphic with Huber-Mises criterion of isotropic materials, we put forward a dimensionless stress yield criterion, and obtained the associated p...On the assumption that the yield criterion of orthotropic materials is isomorphic with Huber-Mises criterion of isotropic materials, we put forward a dimensionless stress yield criterion, and obtained the associated plastic flow law. Using experimental stress-strain curves in various simple stress states, generalized effective stress-strain formulae may be derived correspondingly in various forms.展开更多
Orthotropic materials weakened by a doubly periodic array of cracks under far-field antiplane shear are investigated, where the fundamental cell contains four cracks of unequal size. By applying the mapping technique,...Orthotropic materials weakened by a doubly periodic array of cracks under far-field antiplane shear are investigated, where the fundamental cell contains four cracks of unequal size. By applying the mapping technique, the elliptical function theory and the theory of analytical function boundary value problems, a closed form solution of the whole-field stress is obtained. The exact formulae for the stress intensity factor at the crack tip and the effective antiplane shear modulus of the cracked orthotropic material are derived. A comparison with the finite element method shows the efficiency and accuracy of the present method. Several illustrative examples are provided, and an interesting phenomenon is observed, that is, the stress intensity factor and the dimensionless effective modulus are independent of the material property for a doubly periodic cracked isotropic material, but depend strongly on the material property for the doubly periodic cracked orthotropic material. Such a phenomenon for antiplane problems is similar to that for in-plane problems. The present solution can provide benchmark results for other numerical and approximate methods.展开更多
On the basis of the existing fundamental solutions ofdisplacements, further improvement is made, and then the generalfundamental solutions of both plane elastic and plane plasticproblems for ortho- tropic materials ar...On the basis of the existing fundamental solutions ofdisplacements, further improvement is made, and then the generalfundamental solutions of both plane elastic and plane plasticproblems for ortho- tropic materials are obtained. Two parametersbased on material constants a_1, a_2 are used to derive the rele-vant expressions in a real variable form. Additionally, an analyticalmethod of solving the singular integral for the internal stresses isintroduced, and the corresponding result are given. If a_1=a_2=1, allthe expres- sions obtained for orthotropy can be reduced to thecorresponding ones for isotropy. Because all these expres- sions andresults can be directly used for both isotropic problems andorthotropic problems, it is convenient to use them in engineeringwith the boundary element method (BEM).展开更多
In this work,we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials(FGMs).A recovery type error indicator combined with quadtree decomposition is emp...In this work,we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials(FGMs).A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement.The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency.The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations.The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle.The study is then extended to the analysis of orthotropic FGMs.It is observed that,if the gradation in fracture properties is neglected,the material gradient plays a secondary role,with the fracture behaviour being dominated by the orthotropy of the material.However,when the toughness increases along the crack propagation path,a substantial gain in fracture resistance is observed.展开更多
The purpose of the present work is to study the buckling problem with plate/shell topology optimization of orthotropic material. A model of buckling topology optimization is established based on the independent, conti...The purpose of the present work is to study the buckling problem with plate/shell topology optimization of orthotropic material. A model of buckling topology optimization is established based on the independent, continuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function (PF) as filter functions are introduced to recognize the element mass, the element stiffness matrix, and the element geometric stiffness matrix. The filter functions of the orthotropic material stiffness are deduced. Then these filter functions are put into buckling topology optimization of a differential equation to analyze the design sensitivity. Furthermore, the buckling constraints are approximately expressed as explicit functions with respect to the design variables based on the first-order Taylor expansion. The objective function is standardized based on the second-order Taylor expansion. Therefore, the optimization model is translated into a quadratic program. Finally, the dual sequence quadratic programming (DSQP) algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions (CEF and PF) are applied to solve the optimal model. Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.展开更多
This paper deals with the mathematical modelling and 3D FEM study of the energy release rate(ERR)in the band crack’s front contained in the orthotropic thick rectangular plate which is stretched or compressed initial...This paper deals with the mathematical modelling and 3D FEM study of the energy release rate(ERR)in the band crack’s front contained in the orthotropic thick rectangular plate which is stretched or compressed initially before the loading of the crack's edge planes.The initial stretching or compressing of the plate causes uniformly distributed normal stress to appear acting in the direction which is parallel to the plane on which the band crack is located.After the appearance of the initial stress in the plate it is assumed that the crack's edge planes are loaded with additional uniformly distributed normal forces and the ERR caused with this additional loading is studied.The corresponding boundary value problem is formulated within the scope of the so-called 3D linearized theory of elasticity which allows the initial stress on the values of the ERR to be taken into consideration.Numerical results on the influence of the initial stress,anisotropy properties of the plate material,the crack’s length and its distance from the face planes of the plate on the values of the ERR,are presented and discussed.In particular,it is established that for the relatively greater length of the crack’s band,the initial stretching of the plate causes a decrease,but the initial compression causes an increase in the values of the ERR.展开更多
The behaviors of an interface crack between dissimilar orthotropic elastic half-planes subjected to uniform tension was reworked by use of the Schmidt method.By use of the Fourier transform,the problem can be solved w...The behaviors of an interface crack between dissimilar orthotropic elastic half-planes subjected to uniform tension was reworked by use of the Schmidt method.By use of the Fourier transform,the problem can be solved with the help of two pairs of dual integral equations,of which the unknown variables are the jumps of the displacements across the crack surfaces.Numerical examples are provided for the stress intensity factors of the cracks.Contrary to the previous solution of the interface crack,it is found that the stress singularity of the present interface crack solution is of the same nature as that for the ordinary crack in homogeneous materials.When the materials from the two half planes are the same,an exact solution can be otained.展开更多
Forced and natural vibrations of a rectangular pre-stressed orthotropic compositeplate containing two neighboring cylindrical cavities whose cross sections are rectangularwith rounded-off corners are investigated nume...Forced and natural vibrations of a rectangular pre-stressed orthotropic compositeplate containing two neighboring cylindrical cavities whose cross sections are rectangularwith rounded-off corners are investigated numerically. It is assumed that all the end surfacesof the rectangular pre-stressed composite plate are simply supported and subjected to auniformly distributed normal time-harmonic force on the upper face plane. The consideredproblem is formulated within the Three-Dimensional Linearized Theory of Elastic Waves inInitially Stressed Bodies (TDLTEWISB). The influence of mechanical and geometricalparameters as well as the initial stresses and the effect of cylindrical cavities on the dynamicalcharacteristics of the rectangular orthotropic composite plate are analyzed and discussed.展开更多
The relation between J-integral near model I crack tip in the orthotropic plateand displacement derivative is derived in this paper. Meanwhile,Ihe relation betweenstress intensity factor K_1 and displacement is also ...The relation between J-integral near model I crack tip in the orthotropic plateand displacement derivative is derived in this paper. Meanwhile,Ihe relation betweenstress intensity factor K_1 and displacement is also given in this paper.With stickingfilm moire interferometry method, the three-point bending beam is tested,thus thevalues of J-integral and K_1 can be obtained.from the displacement field,and then thetruth of relation formula between J-integral and K_1 in the orthotrtopic compositematerials is experimentally verified.展开更多
The plane crack problem of an orthotropic functionally graded strip under concentrated loads is studied. The edge crack is perpendicular to the boundary and the elastic property of the material is assumed to vary depe...The plane crack problem of an orthotropic functionally graded strip under concentrated loads is studied. The edge crack is perpendicular to the boundary and the elastic property of the material is assumed to vary depending on thickness. By using an integral transform method, the present problem can be reduced to a single integral equation which is solved numerically. The influences of parameters such as the nonhomogeneity constant and the geometry parameters on the stress intensity factors (SIFs) are studied. It is found that the nonhomogeneity constant has important influences on the SIFs.展开更多
Advanced design based on the concept of orthotropic structure includes better use of materials, less weight compared to the equivalent isotropic construction and controlled effectively reserve resistance in all its se...Advanced design based on the concept of orthotropic structure includes better use of materials, less weight compared to the equivalent isotropic construction and controlled effectively reserve resistance in all its segments. In this case a calculation of critical load is exposed using the FDM (Finite Difference Method) concept of thin plates subjected to complex loads due to forces in the middle-plane. Results of calculation model, discussed in this paper, are given in graphic form. Presented results should serve as an indicator of the expansion of theoretical base of similar models, which can be reasonably use by researchers and engineers in their practices, and by students for educational purposes.展开更多
Stress intensity factors (SIFs) for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method (SBFEM). In the SBFEM, the analytical adva...Stress intensity factors (SIFs) for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method (SBFEM). In the SBFEM, the analytical advantage of the solution in the radial direction allows SIFs to be directly determined from its definition, therefore no special crack-tip treatment is necessary. Furthermore anisotropic material behavior can be treated easily. Different distributions of surface tractions are considered for the center and double-edge-cracked disks. The benchmark examples are modeled and an excellent agreement between the results in the present study and those in published literature is found. It shows that SBFEM is effective and possesses high accuracy. The SIFs of the cracked orthotropic material circular disks subjected to different surface tractions are also evaluated. The technique of substructure is applied to handle the multiple cracks problem.展开更多
Orthotropic nonlinear elastic materials are common in nature and widely used by various industries.However,there are only limited constitutive models available in today's commercial software(e.g.,ABAQUS,ANSYS,etc....Orthotropic nonlinear elastic materials are common in nature and widely used by various industries.However,there are only limited constitutive models available in today's commercial software(e.g.,ABAQUS,ANSYS,etc.)that adequately describe their mechanical behavior.Moreover,the material parameters in these constitutive models are also difficult to calibrate through low-cost,widely available experimental setups.Therefore,it is paramount to develop new ways to model orthotropic nonlinear elastic materials.In this work,a data-driven orthotropic nonlinear elastic(DDONE)approach is proposed,which builds the constitutive response from stress–strain data sets obtained from three designed uniaxial tensile experiments.The DDONE approach is then embedded into a finite element(FE)analysis framework to solve boundary-value problems(BVPs).Illustrative examples(e.g.,structures with an orthotropic nonlinear elastic material)are presented,which agree well with the simulation results based on the reference material model.The DDONE approach generally makes accurate predictions,but it may lose accuracy when certain stress–strain states that appear in the engineering structure depart significantly from those covered in the data sets.Our DDONE approach is thus further strengthened by a mapping function,which is verified by additional numerical examples that demonstrate the effectiveness of our modified approach.Moreover,artificial neural networks(ANNs)are employed to further improve the computational efficiency and stability of the proposed DDONE approach.展开更多
In this paper, the propagation of guided thermoelastic waves in laminated orthotropic plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green-Naghdi (GN) generaliz...In this paper, the propagation of guided thermoelastic waves in laminated orthotropic plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green-Naghdi (GN) generalized thermoelastic theory (without energy dissipation). The coupled wave equations and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The validity of the method is confirmed through a comparison. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding pure elastic plate are also shown to analyze the influence of the thermoelasticity on elastic modes. The displacement and temperature distributions are shown to discuss the differences between the elastic modes and thermal modes.展开更多
文摘On the assumption that the yield criterion of orthotropic materials is isomorphic with Huber-Mises criterion of isotropic materials, we put forward a dimensionless stress yield criterion, and obtained the associated plastic flow law. Using experimental stress-strain curves in various simple stress states, generalized effective stress-strain formulae may be derived correspondingly in various forms.
基金supported by the National Natural Science Foundation of China (No.10672008).
文摘Orthotropic materials weakened by a doubly periodic array of cracks under far-field antiplane shear are investigated, where the fundamental cell contains four cracks of unequal size. By applying the mapping technique, the elliptical function theory and the theory of analytical function boundary value problems, a closed form solution of the whole-field stress is obtained. The exact formulae for the stress intensity factor at the crack tip and the effective antiplane shear modulus of the cracked orthotropic material are derived. A comparison with the finite element method shows the efficiency and accuracy of the present method. Several illustrative examples are provided, and an interesting phenomenon is observed, that is, the stress intensity factor and the dimensionless effective modulus are independent of the material property for a doubly periodic cracked isotropic material, but depend strongly on the material property for the doubly periodic cracked orthotropic material. Such a phenomenon for antiplane problems is similar to that for in-plane problems. The present solution can provide benchmark results for other numerical and approximate methods.
文摘On the basis of the existing fundamental solutions ofdisplacements, further improvement is made, and then the generalfundamental solutions of both plane elastic and plane plasticproblems for ortho- tropic materials are obtained. Two parametersbased on material constants a_1, a_2 are used to derive the rele-vant expressions in a real variable form. Additionally, an analyticalmethod of solving the singular integral for the internal stresses isintroduced, and the corresponding result are given. If a_1=a_2=1, allthe expres- sions obtained for orthotropy can be reduced to thecorresponding ones for isotropy. Because all these expres- sions andresults can be directly used for both isotropic problems andorthotropic problems, it is convenient to use them in engineeringwith the boundary element method (BEM).
基金E.Martínez-Paneda acknowledges financial support from the Royal Commission for the 1851 Exhibition through their Research Fellowship programme(RF496/2018).
文摘In this work,we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials(FGMs).A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement.The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency.The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations.The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle.The study is then extended to the analysis of orthotropic FGMs.It is observed that,if the gradation in fracture properties is neglected,the material gradient plays a secondary role,with the fracture behaviour being dominated by the orthotropy of the material.However,when the toughness increases along the crack propagation path,a substantial gain in fracture resistance is observed.
基金supported by the National Natural Science Foundation of China (Grants 11072009, 11172013)the Beijing Education Committee Development Project (Grant SQKM2016100 05001)the Beijing University of Technology Basic Research Fund (Grant 001000514313003)
文摘The purpose of the present work is to study the buckling problem with plate/shell topology optimization of orthotropic material. A model of buckling topology optimization is established based on the independent, continuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function (PF) as filter functions are introduced to recognize the element mass, the element stiffness matrix, and the element geometric stiffness matrix. The filter functions of the orthotropic material stiffness are deduced. Then these filter functions are put into buckling topology optimization of a differential equation to analyze the design sensitivity. Furthermore, the buckling constraints are approximately expressed as explicit functions with respect to the design variables based on the first-order Taylor expansion. The objective function is standardized based on the second-order Taylor expansion. Therefore, the optimization model is translated into a quadratic program. Finally, the dual sequence quadratic programming (DSQP) algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions (CEF and PF) are applied to solve the optimal model. Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.
文摘This paper deals with the mathematical modelling and 3D FEM study of the energy release rate(ERR)in the band crack’s front contained in the orthotropic thick rectangular plate which is stretched or compressed initially before the loading of the crack's edge planes.The initial stretching or compressing of the plate causes uniformly distributed normal stress to appear acting in the direction which is parallel to the plane on which the band crack is located.After the appearance of the initial stress in the plate it is assumed that the crack's edge planes are loaded with additional uniformly distributed normal forces and the ERR caused with this additional loading is studied.The corresponding boundary value problem is formulated within the scope of the so-called 3D linearized theory of elasticity which allows the initial stress on the values of the ERR to be taken into consideration.Numerical results on the influence of the initial stress,anisotropy properties of the plate material,the crack’s length and its distance from the face planes of the plate on the values of the ERR,are presented and discussed.In particular,it is established that for the relatively greater length of the crack’s band,the initial stretching of the plate causes a decrease,but the initial compression causes an increase in the values of the ERR.
文摘The behaviors of an interface crack between dissimilar orthotropic elastic half-planes subjected to uniform tension was reworked by use of the Schmidt method.By use of the Fourier transform,the problem can be solved with the help of two pairs of dual integral equations,of which the unknown variables are the jumps of the displacements across the crack surfaces.Numerical examples are provided for the stress intensity factors of the cracks.Contrary to the previous solution of the interface crack,it is found that the stress singularity of the present interface crack solution is of the same nature as that for the ordinary crack in homogeneous materials.When the materials from the two half planes are the same,an exact solution can be otained.
文摘Forced and natural vibrations of a rectangular pre-stressed orthotropic compositeplate containing two neighboring cylindrical cavities whose cross sections are rectangularwith rounded-off corners are investigated numerically. It is assumed that all the end surfacesof the rectangular pre-stressed composite plate are simply supported and subjected to auniformly distributed normal time-harmonic force on the upper face plane. The consideredproblem is formulated within the Three-Dimensional Linearized Theory of Elastic Waves inInitially Stressed Bodies (TDLTEWISB). The influence of mechanical and geometricalparameters as well as the initial stresses and the effect of cylindrical cavities on the dynamicalcharacteristics of the rectangular orthotropic composite plate are analyzed and discussed.
文摘The relation between J-integral near model I crack tip in the orthotropic plateand displacement derivative is derived in this paper. Meanwhile,Ihe relation betweenstress intensity factor K_1 and displacement is also given in this paper.With stickingfilm moire interferometry method, the three-point bending beam is tested,thus thevalues of J-integral and K_1 can be obtained.from the displacement field,and then thetruth of relation formula between J-integral and K_1 in the orthotrtopic compositematerials is experimentally verified.
文摘The plane crack problem of an orthotropic functionally graded strip under concentrated loads is studied. The edge crack is perpendicular to the boundary and the elastic property of the material is assumed to vary depending on thickness. By using an integral transform method, the present problem can be reduced to a single integral equation which is solved numerically. The influences of parameters such as the nonhomogeneity constant and the geometry parameters on the stress intensity factors (SIFs) are studied. It is found that the nonhomogeneity constant has important influences on the SIFs.
文摘Advanced design based on the concept of orthotropic structure includes better use of materials, less weight compared to the equivalent isotropic construction and controlled effectively reserve resistance in all its segments. In this case a calculation of critical load is exposed using the FDM (Finite Difference Method) concept of thin plates subjected to complex loads due to forces in the middle-plane. Results of calculation model, discussed in this paper, are given in graphic form. Presented results should serve as an indicator of the expansion of theoretical base of similar models, which can be reasonably use by researchers and engineers in their practices, and by students for educational purposes.
基金financially supported by the National Natural Youth Foundation of China (Grant Nos. 51109134,51009019, 11102118 and 51208310)the Liaoning Province Education Administration Foundation (Grant No. L2010413)+1 种基金the China Postdoctoral Science Foundation (Grant No. 2011M500557)the Natural Science Foundation of Liaoning Province (Grant No.20102164)
文摘Stress intensity factors (SIFs) for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method (SBFEM). In the SBFEM, the analytical advantage of the solution in the radial direction allows SIFs to be directly determined from its definition, therefore no special crack-tip treatment is necessary. Furthermore anisotropic material behavior can be treated easily. Different distributions of surface tractions are considered for the center and double-edge-cracked disks. The benchmark examples are modeled and an excellent agreement between the results in the present study and those in published literature is found. It shows that SBFEM is effective and possesses high accuracy. The SIFs of the cracked orthotropic material circular disks subjected to different surface tractions are also evaluated. The technique of substructure is applied to handle the multiple cracks problem.
基金The support of Project MKF20210033 is acknowledged.
文摘Orthotropic nonlinear elastic materials are common in nature and widely used by various industries.However,there are only limited constitutive models available in today's commercial software(e.g.,ABAQUS,ANSYS,etc.)that adequately describe their mechanical behavior.Moreover,the material parameters in these constitutive models are also difficult to calibrate through low-cost,widely available experimental setups.Therefore,it is paramount to develop new ways to model orthotropic nonlinear elastic materials.In this work,a data-driven orthotropic nonlinear elastic(DDONE)approach is proposed,which builds the constitutive response from stress–strain data sets obtained from three designed uniaxial tensile experiments.The DDONE approach is then embedded into a finite element(FE)analysis framework to solve boundary-value problems(BVPs).Illustrative examples(e.g.,structures with an orthotropic nonlinear elastic material)are presented,which agree well with the simulation results based on the reference material model.The DDONE approach generally makes accurate predictions,but it may lose accuracy when certain stress–strain states that appear in the engineering structure depart significantly from those covered in the data sets.Our DDONE approach is thus further strengthened by a mapping function,which is verified by additional numerical examples that demonstrate the effectiveness of our modified approach.Moreover,artificial neural networks(ANNs)are employed to further improve the computational efficiency and stability of the proposed DDONE approach.
基金supported by the National Natural Science Foundation of China (No. 10802027)
文摘In this paper, the propagation of guided thermoelastic waves in laminated orthotropic plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green-Naghdi (GN) generalized thermoelastic theory (without energy dissipation). The coupled wave equations and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The validity of the method is confirmed through a comparison. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding pure elastic plate are also shown to analyze the influence of the thermoelasticity on elastic modes. The displacement and temperature distributions are shown to discuss the differences between the elastic modes and thermal modes.