Plate/shell structure topology optimization of orthotropic material for buckling problem based on independent continuous topological variables

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.

EXACT SOLUTION FOR ORTHOTROPIC MATERIALS WEAKENED BY DOUBLY PERIODIC CRACKS OF UNEQUAL SIZE UNDER ANTIPLANE SHEAR

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.

FURTHER IMPROVEMENT ON FUNDAMENTAL SOLUTIONS OF PLANE PROBLEMS FOR ORTHOTROPIC MATERIALS

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).

A STUDY ON YIELD AND FLOW OF ORTHOTROPIC MATERIALS

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.

Adaptive phase field modelling of crack propagation in orthotropic functionally graded materials

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.

ON J-INTEGRALS NEAR MODELS Ⅰ,Ⅱ CRACK TIPS IN THE PLATE OF ORTHOTROPIC COMPOSITE MATERIAL

In this paper, we prove the independence of the path of the imegrals near model i II track tips in the plate of the orthotropic composite material. Then we derive the computing formulae of the J-integrals at the cases of and by using a complex variable method and redueing J-imegrals to compley farmThe J-integral computational formulae derived this paper have certain reference value for the theoretical researches and the experimental verifications in the plane fracture for composite material.

A STUDY OF J-INTEGRAL OF THE ORTHOTROPIC COMPOSITE MATERIAL

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.

Mathematical Modelling and 3D FEM Analysis of the Influence of Initial Stresses on the ERR in a Band Crack’s Front in the Rectangular Orthotropic Thick Plate

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.

INVESTIGATION OF THE BEHAVIOR OF A GRIFFITH CRACK AT THE INTERFACE BETWEEN TWO DISSIMILAR ORTHOTROPIC ELASTIC HALF-PLANES FOR THE OPENING CRACK MODE

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 an Orthotropic Pre-Stressed Rectangular Plate with Neighboring Two Cylindrical Cavities

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.

Evaluation of Stress Intensity Factors for Multiple Cracked Circular Disks Under Crack Surface Tractions with SBFEM

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.

Derivation of the Orthotropic Nonlinear Elastic Material Law Driven by Low-Cost Data(DDONE)

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.

GUIDED THERMOELASTIC WAVE PROPAGATION IN LAYERED PLATES WITHOUT ENERGY DISSIPATION

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.