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.

Symplectic system based analytical solution for bending of rectangular orthotropic plates on Winkler elastic foundation

This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.

An analytical symplectic approach to the vibration analysis of orthotropic graphene sheets

A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach. A Hamiltonian system is established by introducing a total unknown vector consisting of the displacement amplitude, rotation angle, shear force, and bending moment. The high-order governing differential equation of the vibration of SLGSs is transformed into a set of ordinary differential equations in symplectic space. Exact solutions for free vibration are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary conditions. Vibration modes are expressed in terms of the symplectic eigenfunctions. In the numerical examples, comparison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given. A parametric study of the natural frequency is also included.

A THREE-DIMENSIONAL SOLUTION FOR LAMINATED ORTHOTROPIC RECTANGULAR PLATES WITH VISCOELASTIC INTERFACES

When a body consists completely or even partly of viscoelastic materials, its response under static loading will be time-dependent. The adhesives used to glue together single plies in laminates usually exhibit a certain viscoelastic characteristic in a high temperature environment. In this paper, a laminated orthotropic rectangular plate with viscoelastic interfaces, described by the Kelvin-Voigt model, is considered. A power series expansion technique is adopted to approximate the time-variation of various field quantities. Results indicate that the response of the laminated plate with viscoelastic interfaces changes remarkably with time, and is much different from that of a plate with spring-like or viscous interfaces.

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

Time-reverse location of microseismic sources in viscoelastic orthotropic anisotropic medium based on attenuation compensation

Time reversal is a key component of time-reverse migration and source location using wavefield extrapolation.The implementation of time reversal depends on the time symmetry of wave equations in acoustic and elastic media.This symmetry in time is no longer valid in attenuative medium.Not only the velocity is anisotropic in shale oil and gas reservoirs,but also the attenuation is usually anisotropic,which can be characterized by viscoelastic orthotropic media.In this paper,the fractional order viscoelastic anisotropic wave equation is used to decouple the energy dissipation and the velocity dispersion.By changing the sign of the dissipation term during backpropagation,the anisotropic attenuation is compensated and the time symmetry is restored.The attenuation compensation time-reverse location algorithm can eff ectively locate the source in viscoelastic orthotropic media.Compared to cases without attenuation compensation or using isotropic attenuation compensation,this method can remove location error caused by anisotropic attenuation and improve the imaging eff ect of the source.This paper verifi es the eff ectiveness of the method through theoretical analysis and model testing.

THE EVALUATION OF STRESS INTENSITY FACTORS OF PLANE CRACK FOR ORTHOTROPIC PLATE WITH EQUAL PARAMETER BY F2LFEM

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

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.

Stress field near interface crack tip of double dissimilar orthotropic composite materials

In this paper, double dissimilar orthotropic composite materials interfacial crack is studied by constructing new stress functions and employing the method of composite material complex. When the characteristic equations＇ discriminants △1 〉 0 and △2 〉0, the theoretical formula of the stress field and the displacement field near the mode I interface crack tip are derived, indicating that there is no oscillation and interembedding between the interfaces of the crack.

Measuring Stress Distributions of Orthotropic Composite Material in Plane Stress State by the Lock-in Infrared Thermography Technique

Feasibility of measuring stress distributions of orthotropic composite materials and structures in plane stress state by the lock-in infrared thermography technique is analyzed and stress distributions of a lap joint structure made of a kind of glass reinforced plastic composite lamination plates under tensile loadings are obtained by the lock-in infrared thermography technique.Feasibility and credibility of using this technique to measure stress distributions of orthotropic composite materials and structures in plane stress state are proved by comparing the results with the data given by the digital speckle correlation method.

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.

THE PROBLEMS OF THE NONLINEAR UNSYMMETRICAL BENDING FOR CYLINDRICALLY ORTHOTROPIC CIRCULA RPLATE（I）

In this paper, ihe probleins of nonlinear unsymmeirucal bending for cylindricallyorthotropic circular plale are sludied by using “ the method of two-variabie”￣[1], and theuniformly valid asympiotic soluiions of Nth-order .lor ε_1 and Mth-order for ε_2 areobiained

Effect of Non-Homogeneity on Thermally Induced Vibration of Orthotropic Visco-Elastic Rectangular Plate of Linearly Varying Thickness

The analysis presented here is to study the effect of non-homogeneity on thermally induced vibration of orthotropic visco-elastic rectangular plate of linearly varying thickness. Thermal vibrational behavior of non-homogeneous rectangular plates of variable thickness having clamped boundary conditions on all the four edges is studied. For non–homogeneity of the plate material, density is assumed to vary linearly in one direction. Using the method of separation of variables, the governing differential equation is solved. An approximate but quite convenient frequency equation is derived by using Rayleigh-Ritz technique with a two-term deflection function. Time period and deflection at different points for the first two modes of vibration are calculated for various values of temperature gradients, non- homogeneity constant, taper constant and aspect ratio. Comparison studies have been carried out with non-homogeneous visco-elastic rectangular plate to establish the accuracy and versatility.

Stresses of orthotropic laminated beams subjected to high temperature and mechanical load

Thermo-elastic analysis of simply-supported orthotropic laminated beams subjected to high temperature and mechanical load is presented on the basis of the exact two-dimensional thermoelasticity theory.The beam is composed of several orthotropic layers,each with temperaturedependent material properties.The governing equation for each layer is analytically solved using the state space method.The displacement and stress solutions of the beam are obtained using the transfer-matrix method.A numerical example is included to study the effects of temperature on the mechanical responses of a sandwich beam.The results reveal two main effects of temperature:(i)inducing deformations and stresses by itself;(ii)affecting the deformations and stresses induced by the mechanical load.

ELASTODYNAMIC SOLUTION OF A NON-HOMOGENEOUS ORTHOTROPIC HOLLOW CYLINDER

A theoretical method for analyzing the axisymmetric plane strain elastodynamic problem of a non-homogeneous orthotropic hollow cylinder is developed. Firstly, a new dependent variable is introduced to rewrite the governing equation, the boundary conditions and the initial conditions. Secondly, a special function is introduced to transform the inhomogeneous boundary conditions to homogeneous ones. By virtue of the orthogonal expansion technique, the equation with respect to the time variable is derived, of which the solution can be obtained. The displacement solution is finally obtained, which can be degenerated in a rather straightforward way into the solution for a homogeneous orthotropic hollow cylinder and isotropic solid cylinder as well as that for a non-homogeneous isotropic hollow cylinder. Using the present method, integral transform can be avoided and it can be used for hollow cylinders with arbitrary thickness and subjected to arbitrary dynamic loads. Numerical results are presented for a non-homogeneous orthotropic hollow cylinder subjected to dynamic internal pressure.

NATURAL FREQUENCY FOR RECTANGULAR ORTHOTROPIC CORRUGATED-CORE SANDWICH PLATES WITH ALL EDGES SIMPLY-SUPPORTED

A simple approach to reduce the governing equations for orthotropic corrugated core sandwich plates to a single equation containing only one displacement function is presented, and the exact solution of the natural frequencies for rectangular corrugated-core sandwich plates with all edges simply-supported is obtained. Furthermore, two special cases of practical interests are discussed in details.

TRANSIENT THERMAL STRESSES IN AN ORTHOTROPIC HOLLOW CYLINDER FOR AXISYMMETRIC PROBLEMS

For the thermoelastic dynamic axisymmetric problem of a finite orthotropic hollow cylinder,one comes closer to reality to involve the effect of axial strain than to consider the plane strain case only.However,additional mathematical difficulties should be encountered and a different solution procedure should be developed.By the separation of variables,the thermoelastic axisymmetric dynamic problem of an orthotropic hollow cylinder taking account of the axial strain is transformed to a Volterra integral equation of the second kind for a function of time,which can be solved efficiently and quickly by the interpolation method.The solutions of displacements and stresses are obtained. It is noted that the present method is suitable for an orthotropic hollow cylinder with an arbitrary thickness subjected to arbitrary axisymmetric thermal loads.Numerical comparison is made to show the effect of the axial strain on the displacements and stresses.

Analysis of stress intensity factor in orthotropic bi-material mixed interface crack

Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With con- sideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial dif- ferential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit, undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship be- tween the bi-material stress intensity factor and the ratio of bi-material Young＇s modulus are given in the numerical analysis.

NONLINEAR FREE VIBRATION OF ORTHOTROPIC SHALLOW SHELLS OF REVOLUTION UNDER THE STATIC LOADS

In this paper, the axisymmetric nonlinear free vibration problems of cylindrically orthotropic shallow thin spherical and conical shells under uniformly distributed static loads are studied by using MWR and Lindstedt-Poincare perturbation method, from which, the characteristic relation between frequency ratio and amplitude is obtained. The effects of static loads, geometric and material parameters on vibrational behavior of shells are also discussed.