Methods of finite element analysis of arbitrary buckling forms of sandwich plates and shells

Two statements of the problem of arbitrary buckling forms (BFs) (including synphasic, antiphasic, mixed flexural, flexural-shear, and shear forms in the tangential directions) of general-form sandwich shells and two schemes of its solution by the FEM are given. The first of the schemes is based on the use of refined linear equations for determination of the precritical stress-strain state and linearized equations of neutral equilibrium with all parametric addends necessary to determine the critical loads and reveal the possible BFs. The second one uses the general geometrically nonlinear relations of elasticity theory for investigation of the whole deformation process up to buckling in terms of a modified incremental (stepwise) statement of the problem. Examples of solution of particular problems are given.Center for Study of Dynamics and Stability, Tupolev Kazan' State Technical University, Kazan', Tatarstan, Russia. Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 4, pp. 473–486, July–August, 2000.     

Refined geometric nonlinear theory of sandwich shells with a transversely soft core of medium thickness for investigation of mixed buckling forms

A variant of the refined geometric nonlinear theory is suggested for nonshallow shells with a transversely soft core of medium thickness with regard to modifications of metric characteristics across the core thickness. The kinematic relations for the core are derived by sequential integration of the initial three-dimensional equations of elasticity theory along the transverse coordinate. The equations are preliminarily simplified by the assumption that the tangential stress components are equal to zero. With the example of sandwich plates, it is shown that these equations allow us to investigate synphasic, antiphasic, mixed flexural, and mixed flexural-shear buckling forms of load-bearing layers and the core depending on the precritical stress-strain state. Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 1, pp. 95–108, January–February, 2000.     

Calculation of rigidly clamped circular plates in the statement of the problem on contact interaction in a two-layer package

The use of a mechanicomathematical model of bending of thick transversely isotropic plates is illustrated, where the plate is divided in an arbitrary number of equally thick conditional layers. This model allows one to approximately reduce the problem of determination of stresses and displacements in the thick plate to a corresponding contact problem for a bent pack age of layers. The axisymmetric bending of a rigidly clamped package consisting of two plates rigidly fastened together is considered. The results of numerical calculations are presented, which are compared with those obtained within the framework of a refined bending model of plates (with account of transverse compression and shear) and of the Timoshenko model, as well as in the statement of the three-dimensional theory of elasticity. The accuracy of satisfying the boundary conditions in each model is analyzed. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 1, pp. 93–108, January–February, 2007.     

Shear Buckling Form of a Circular Sandwich Ring under Uniform External Pressure

The problem on the stability of a circular sandwich ring under uniform external pressure is considered. It is shown that, along with a mixed flexural-shear buckling form (BF), a pure shear BF can be realized in the core. This form is accompanied by rotation of the load-carrying layers at the cost of the transverse shear strain (constant in the circumferential direction) in the core. It is found that the simplified equations of the theory of shallow shells cannot describe this nonclassical BF. The critical loads corresponding to the shear BF may prove to be smaller than the critical load of the classical mixed flexural BF for a circular sandwich ring of medium thickness at a low shear modulus of the core. The results obtained contribute greatly to the understanding of buckling mechanisms of sandwich structures and supplement the existing classification of BFs.     

A continuous model for investigation of physically nonlinear deformation of orthotropic sandwich plates

A phenomenological yield condition for quasi-brittle and plastic orthotropic materials with initial stresses is suggested. All components of the yield tensor are determined from experiments on uniaxial loading. The reliability estimates of the criterion suggested is discussed. For a plastic material without initial stresses, the given condition transforms into the Marin—Hu criterion. The defining equations of the deformation theory of plasticity with isotropic and “anisotropic” hardening, associated with the yield condition suggested, are obtained. These equations are used as the basis for a highly accurate nonclassical continuous model for nonlinear deformation of thick sandwich plates. The approximations with respect to the transverse coordinate take into account the flexural and nonflexural deformations in transverse shear and compression. The high-order approximations allow us to model the occurrence of layer delamination cracks by introducing thin nonrigid interlayers without violating the continuity concept of the theory. Submitted to the 11th International Conference on Mechanics of Composite Materials (Riga, June 11–15, 2000). Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. pp. 329–340, May–June, 2000.     

Buckling forms of homogeneous and sandwich plates in pure shear in tangential directions

The mathematical problem of the plane shear buckling form (BF) of sandwich plates and plates homogeneous across the thickness in pure shear is considered. The solution to this problem is compared with the solution to the problem of a flexural BF which is realized in these plates with the formation of oblique waves. It is established that, in the case of plates homogeneous across the thickness, the critical loads corresponding to the plane shear BF are maximum for isotropic ones. In real one-layer structural elements manufactured both of isotropic homogeneous and orthotropic composite materials, these critical loads cannot be reached since they exceed considerably the critical loads for the flexural BFs with oblique waves. The critical loads corresponding to the two BFs are comparable only for relatively thick plates. However, the plane shear BF can be realized in sandwich plates earlier than the flexural one even if the plates are thin. Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 2, pp. 215–228, 2000.     

Stress State of Multilayer Thin-Walled Elements with Interfacial Defects of Structure

Based on the discrete-structural theory of thin plates and shells, a calculation model for thin-walled elements consisting of a number of rigid anisotropic layers is put forward. It is assumed that the transverse shear and compression stresses are equal on the interfaces. Elastic slippage is allowed over the interfaces between adjacent layers. The solution to the problem is obtained in a geometrically nonlinear statement with account of the influence of transverse shear and compression strains. The stress-strain state of circular two-layer transversely isotropic plates, both without defects and with a local area of adhesion failure at their center, is investigated numerically and experimentally. It is found that the kinematic and static contact conditions on the interfaces of layered thin-walled structural members greatly affect the magnitude of stresses and strains. With the use of three variants of calculation models, in the cases of perfect and weakened contact conditions between layers, the calculation results for circular plates are compared. It is revealed that the variant suggested in this paper adequately reflects the behavior of layered thin-walled structural elements under large deformations. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 41, No. 6, pp. 761–772, November–December, 2005.     

An unflexural refined model of deformation of multilayered plates on an undeformable foundation

The stress-strain state of multilayered plates on an undeformable foundation is investigated. The design diagram of a transversely loaded plate is formed by supplementing it with a symmetric one about the contact surface of the foundation. The double-thickness plate obtained becomes loaded bilaterally and symmetrically about its midsurface. In such a way, only unflexural deformations can be modeled, which reduces the number of unknowns and the general order of differentiation of the resolving system of equations. The refined continual model developed takes into account the deformations of transverse shear and transverse compression in high iterative approximations. The cases of a rigid contact between the foundation and the plate and frictionless slip of the plate over the foundation are considered. Calculations confirm the efficiency of such an approach, which allows one to obtain solutions qualitatively and quantitatively close to three-dimensional ones. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 42, No. 5, pp. 643–654, September–October, 2006.     

High-precision modeling of deformation of laminated structures

Precise three-dimensional solutions for homogeneous, two- and three-layer plates of symmetric and nonsymmetric structure over their thickness with orthotropic layers in transverse symmetric and antisymmetric loading normal to the surface of the plates are given. It is shown that the character of stressed states under flexural (antisymmetric) and nonflexural (symmetric transverse compression) loading differs greatly. It is noted that the known refined continual models, which take into account the transverse shear and compression, are all essentially flexural and therefore cannot describe the nonflexural deformations well. In particular, continual shear models in symmetric pressure loading lead to zero solutions. A refined nonflexural continual model of deformation of sandwich plates in bilateral symmetric compression is constructed. The general order of resolving differential equations for continual models does not depend on the number of layers. Approximation functions of the transverse coordinate are obtained with the help of well-founded hypotheses. A high-accuracy variant of the flexural continual model is proposed for antisymmetric loading with account of shear and transverse normal strains, as well as a version combining both models mentioned. A method of precise satisfation of all the constitutive relations for the layers, including the conditions of their contact, is proposed, whereas in the known continual models the dependence between the transverse normal stress and strain is satisfied only integrally, or else the Poisson effect is neglected.Ukrainian State Academy of Water Management, Rovno, the Ukraine. Translated from Mekhanika Kompozitnykh Materialov, Vol. 35, No. 1, pp. 13–28, January–February, 1999.     

Deformation of sandwich panels in cylindrical bending by point forces. 1. Analytical construction

A refined model for bending of a three-layered panel with a soft filler is proposed. The modified model permits us to consider the asymmetry of elastic properties and thickness of the outer layer relative to the middle plane of the panel in a composite sandwich structure. In constructing the deformation mechanism, a heterogeneous kinematic model was adopted, which, in contrast to the assumptions for the deformation of the whole stack of layers, features four degrees of displacement freedom permitting consideration of the separate nature of the deformation of the outer layers in bending and of the intermediate layer in transverse compression combined with shear. This approach is postulated according to an energetic evaluation of the deformation of the layers [2]. The specific features of the stress from point forces in cylindrical bending are considered using the operational Laplace method, which avoids the additional difficulties in analyzing the solution convergence arising when it is represented by a series of eigenfunctions of the boundary value problem. The fundamental functions of a twelfth-order set of equations are used to construct the boundary problem reduced to a Cauchy problem. Various boundary effects of the point stress are described using a generalized Dirac function. Variants are examined for the limiting transformation of the model parameters leading to a qualitative change in its kinematics and the corresponding simplified bending models. Institute of Polymer Mechanics, Latvian Academy of Sciences, LV-1006 Riga, Latvia. Translated from Mekhanika Kompozitnykh Materialov, No. 5, pp. 588–611, September–October, 1996.     

Calculation of simply supported circular plates in the statement of the problem of contact interaction in a package of two plates

The possibility of applying the mechanicomathematical model of bending of a package of transversely isotropic plates to approximate calculations of thick plates is investigated. Within the frame work of this model, the problem of axisymmetric bending of a package of two identical plates with simply supported edges is considered. The conditions of rigid contact are given between the plates. Based on the analytical result obtained, the unknown distribution of stresses and displacements across the thickness of the plate is approximated by the distribution of corresponding parameters in the package of two plates. For an isotropic body, the results of numerical calculations are compared with those given by the 3D theory of elasticity. In the case of a transversely isotropic body, a comparison with the results found by the refined bending model of plates (taking into account the transverse compression and shear) and the Timoshenko model is carried out. The accuracy to which the boundary conditions in every model are satisfied is analyzed. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 5, pp. 603–616, September–October, 2007.     

Applied model of a round cylinder reinforced with systems of yarns at large tensile,inflation, and torsional deformations

Equations for a round cylinder weakly reinforced with systems of yarns and subjected to large tensile, inflation, and torsional deformations are presented. Since the degree of filling is small, the model of uniaxial stress state is assumed. The fibers are aligned with spirals on cylindrical surfaces and with radii in the transverse and meridional sections of the cylinder. The equations are obtained in the macroscopically unidimensional statement for the case of cylindrically symmetric strains. Numerical results are given for twisted hollow rubber cylinders reinforced with polymer yarns in the axial and radial directions. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 2, pp. 237–256, March–April, 2007.     

Identifying the elastic moduli of composite plates by using high-order theories. 2. Theoretical-experimental approach

The identification of elastic properties of laminated composite plates from measured eigenfrequencies has been performed. The elastic moduli of the laminates were determined by using a multilevel modeling and a two-step identification procedure. At the first step, based on a genetic algorithm, the Young’s and shear moduli were found, but at the second one, by minimizing an error function, the values of transverse moduli were refined. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 44, No. 2, pp. 207–216, March–April, 2008.     

Non-uniformly scaled buckling modes of reinforcing elements in fiber reinforced plastic

We consider a problem about non-uniformly scaled buckling modes of isolated fiber (without accounting of interaction with the surrounding epoxy) or bundle of fibers, which are structural elements of fiber reinforced plastics under the transverse tension (compression) and shear stresses in prebuckling state. Such initial state is formed in fibers and bundles of fibers at tension-compression tests of flat specimens from cross ply composites with unidirectional fibers. For problem statement we use equations recently constructed by reduction of consistent version of geometrically nonlinear equations of theory of elasticity to one dimensional equations of rectilinear beams. Equations are based on refined shear S. P. Timoshenko model with accounting of tension-compression stresses in transverse directions. We give theoretical explanation of developed phenomenon as reducing shear modulus of elasticity of fiber reinforced plastic during the increasing of shear strains. We show that under the loading process of specimens under review uninterruptedly structure reconstruction of composite trough implementation and uninterruptedly changing of internal buckling modes at changing wave parameter is feasible.     

Refined stability theory for sandwich shells with transversally stiff core. 2. Linearized equations of neutral equilibrium

Neutral equilibrium equations of the refined theory of stability for sandwich shells with a transversally stiff core are constructed and used for studying local mixed forms of stability loss (FSL), as well as admitting different variants of simplification, depending on the type of precritical state and realized FSL. The generalized Reissner variational principle used for deriving the stability equations allows us to refine transverse shear stresses in the core as compared to [1]. A method for a highly accurate definition of these stresses is proposed. Namely, after the integration of three-dimensional equilibrium equations over the transverse coordinate, the number of free constants and the number of static conditions to be satisfied are equalized according to the actual stress distribution across the thickness.Science and Technology Center for Study of Dynamics and Strength. Tupolev Kazan State Technical University, Kazan, Tatarstan, Russia. Translated from Mekhanika Kompozitnykh Materialov, Vol. 33, No. 6, pp. 786–795, November–December, 1997.     

Dynamic properties of composite laminates and panels with highly contrasting directions of anisotropy

An asymptotic analysis of the dynamic stress-strain state of a thin laminar packet of anisotropic layers is presented. The statement is nonclassical, since in layer materials the ratio of elastic moduli in the longitudinal and transverse directions can generate small parameters comparable to the relative half-thickness of the packet, as, for example, in high-strength unidirectional composites. Alternation of strong load-carrying layers and a relatively soft filler with a similar difference in the elastic moduli between the layers is also allowed. The averaged two-dimensional equations and the total stress tensor in the layers are determined. The results are classified with respect to the types of anisotropy and the indices of differences in the elastic moduli. It is shown that first-approximation models lead to kinematic relations similar to those of the theories of high-order shear strains. Institute of Problems of Mechanics, Russian Academy of Sciences, Moscow, Russia. Translated from Mekhanika Kompozitnykh Materialov, Vol. 35, No. 5, pp. 605–614, September–October, 1999.     

Nonlinear free vibrations of multilayer shallow shells with a symmetric structure and with a complicated form of the plan

We propose a method to study free nonlinear vibrations of multilayer shallow shells with a complicated form of the plan. The mathematical statement of the problem is realized in the frame of a refined firstorder theory of the Tymoshenko type. A distinctive specific feature of the work is the application of the theory of R -functions and variational methods to the determination of eigenfunctions used as the basis in the construction of the required solution of a nonlinear problem. The proposed method is tested, and some new problems are solved. As a result, we constructed the amplitude–frequency dependences for spherical shells with a complicated form of the plan.     

A limit surface formulation for plastically compressible polymers

The tightening of industrial safety standards for structures generates a need for refined computational methods, which, among other things, must be able to describe the yield surface and the deformation behaviour of non-reinforced thermoplastics. To describe the plastic behaviour of materials, a potential formulation is suggested. This formulation contains a number of known potentials as special cases. The parameters of the model, which are obtained from test data, are restricted by the convexity condition for the potential. The new model allows one to take into account effects of the second order, for instance, the unequal behaviour under tension and compression, the plastic compressibility, and the Poynting-Swift effect. For each particular choice of the parameters, the Poisson ratio in tension is computed. If the restrictions imposed on the Poisson ratio do not hold, the application of a non-associated flow rule is necessary. A simple non-associated flow rule with different values of Poisson ratio intension and compression is proposed. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 3, pp. 367–384, May–June, 2007.     

Discussion: “On interactive buckling in a sandwich structure” by C. D. Coman

A recent article by Coman (Z Angew Math Phys 2009) on the response of compression sandwich struts made some claims on the quality of the simplified version of the interactive buckling model presented in Hunt and Wadee (Proc R Soc A 454(1972):1197–1216, 1998). Some of these claims are examined in detail herein; it is concluded that great care must be exercised when performing parametric studies with equations that have been derived from simplifying a mechanical model. This is because the resulting system of equations does not necessarily describe the original mechanical system in full, since the key assumptions necessarily change.     

The large time behavior of strong solutions to a P1-approximation arising in radiating hydrodynamics

This paper is concerned with the study of the Cauchy problem to a multi-dimensional P1-approximation model. Based on a known global well-posedness (Danchin and Ducomet in J Evol Equ 14:155–195, 2014), in $$L^{2}$$-critical regularity framework the time decay rates of the constructed global strong solutions are obtained if the low frequencies of the data under a suitable additional condition. The proof mainly relies on an application of Fourier analysis to a mixed parabolic-hyperbolic system, and on a refined time-weighted energy functional. As a by-product, those time-decay rates of $$L^{q}$$–$$L^{r}$$ type are also captured in the critical framework.