A method for solving boundary value problems and two-dimensional theories without ad hoc assumptions

A method is proposed with which certain three-dimensional problems may be solved and also approximate two-dimensional theories can be deduced from three-dimensional equations of the theory of elasticity without ad hoc assumptions. To illustrate the method, the axially symmetric problem of hollow cylinders of any length is treated herein. The method, however, can be applied to nonaxisymmetric deformation of cylinders, spheres, plates, and other problems.     

An exact derivation of the thin plate equation

It is shown that, when the traditional assumptions of thin plate theory are taken as exact methematical hypotheses, the desired field and boundary equations can be obtained by mere integration over the thickness of the corresponding equations for a three-dimensional cylindrical body made of a homogeneous, linearly elastictransversely isotropic, constrained material, yet avoiding some inconsistencies usually to be found in textbooks of structural mechanics.     

Optimizing formulation of hybrid model for thin and moderately thick plates

A 20 — DOF hybrid stress element based upon Mindlin plate theory is developed using the optimization design method for thin and moderately thick plates. Numerical tests consist of the convergency and performance to the plates with arbitrary thickness and shape and of the ultimate thin plate problems.Projects Supported by the National Natural Science Foundation of China.     

大变形薄板多体系统的动力学建模

基于非线性弹性理论,考虑剪切应变和横向应变,用绝对节点坐标法建立了大变形矩形薄板的动力学变分方程;为了提高非线性刚度阵的计算效率,根据非线性刚度阵与广义坐标阵的函数关系式,在非线性刚度阵中分离出广义坐标阵,从而避免了每个时间步长的单元刚度阵的积分运算。在此基础上,引入运动学约束关系,建立了大变形薄板系统第一类拉格朗日方程,对重力作用下大变形二连板进行数值仿真。计算结果表明:随着薄板的柔度增大,低频的弯曲变形与高频拉伸变形的耦合愈加显著;此外,系统机械能守恒验证了该模型正确性。     

A thick anisotropic plate element in the framework of an absolute nodal coordinate formulation

In this research, the incorporation of material anisotropy is proposed for the large-deformation analyses of highly flexible dynamical systems. The anisotropic effects are studied in terms of a generalized elastic forces (GEFs) derivation for a continuum-based, thick, and fully parameterized absolute nodal coordinate formulation plate element, of which the membrane and bending deformation effects are coupled. The GEFs are first derived for a fully anisotropic, linearly elastic material, characterized by 21 independent material parameters. Using the same approach, the GEFs are obtained for an orthotropic material, characterized by nine material parameters. Furthermore, the analysis is extended to the case of nonlinear elasticity; the GEFs are introduced for a nonlinear Cauchy-elastic material, characterized by four in-plane orthotropic material parameters. Numerical simulations are performed to validate the theory for statics and dynamics and to observe the anisotropic responses in terms of displacements, stresses, and strains. The presented formulations are suitable for studying the nonlinear dynamical behavior of advanced elastic materials of an arbitrary degree of anisotropy.     

Boundary conditions at the edge of a thin or thick plate bonded to an elastic support

At the clamped edge of a thin plate, the interior transverse deflection ω(x ₁, x₂) of the mid-plane x ₃=0 is required to satisfy the boundary conditions ω=?ω/?n=0. But suppose that the plate is not held fixed at the edge but is supported by being bonded to another elastic body; what now are the boundary conditions which should be applied to the interior solution in the plate? For the case in which the plate and its support are in two-dimensional plane strain, we show that the correct boundary conditions for ω must always have the form % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakqaabeqaaiaabEhacaqGTaWa% aSaaaeaacaGG0aGaae4vamaaCaaaleqabaGaamOqaaaaaOqaaiaaco% dadaqadaqaaiaacgdacqGHsislcaqG2baacaGLOaGaayzkaaaaaiaa% bIgadaahaaWcbeqaaiaackdaaaGcdaWcaaqaaiaabsgadaahaaWcbe% qaaiaackdaaaGccaqG3baabaGaaeizaiaabIhafaqabeGabaaajaaq% baqcLbkacaGGYaaajaaybaqcLbkacaGGXaaaaaaakiabgUcaRmaala% aabaGaaiinaiaabEfadaahaaWcbeqaaiaadAeaaaaakeaacaGGZaWa% aeWaaeaacaGGXaGaeyOeI0IaaeODaaGaayjkaiaawMcaaaaacaqGOb% WaaWbaaSqabeaacaGGZaaaaOWaaSaaaeaacaqGKbWaaWbaaSqabeaa% caGGZaaaaOGaae4DaaqaaiaabsgacaqG4bqcaaubaeqabiqaaaqcaa% saaiaacodaaKaaafaajugGaiaacgdaaaaaaOGaeyypa0Jaaiimaiaa% cYcaaeaadaWcaaqaaiaabsgacaqG3baabaGaaeizaiaabIhaliaacg% daaaGccqGHsisldaWcaaqaaiaacsdacqqHyoqudaahaaWcbeqaaiaa% bkeaaaaakeaacaGGZaWaaeWaaeaacaGGXaGaeyOeI0IaaeODaaGaay% jkaiaawMcaaaaacaqGObWaaSaaaeaacaqGKbWaaWbaaSqabeaacaGG% YaaaaOGaae4DaaqaaiaabsgacaqG4bqbaeqabiqaaaqcaauaaKqzGc% GaaiOmaaqcaawaaKqzGcGaaiymaaaaaaGccqGHRaWkdaWcaaqaaiaa% csdacqqHyoqudaahaaWcbeqaaiaabAeaaaaakeaacaGGZaWaaeWaae% aacaGGXaGaeyOeI0IaaeODaaGaayjkaiaawMcaaaaacaqGObWaaWba% aSqabeaacaGGYaaaaOWaaSaaaeaacaqGKbWaaWbaaSqabeaacaGGZa% aaaOGaae4DaaqaaiaabsgacaqG4bqcaaubaeqabiqaaaqcaasaaiaa% codaaKaaafaajugGaiaacgdaaaaaaOGaeyypa0JaaiimaiaacYcaaa% aa!993A!\[\begin{gathered}{\text{w - }}\frac{{4{\text{W}}^B }}{{3\left( {1 - {\text{v}}} \right)}}{\text{h}}^2 \frac{{{\text{d}}^2 {\text{w}}}}{{{\text{dx}}\begin{array}{*{20}c}2 \\1 \\\end{array} }} + \frac{{4{\text{W}}^F }}{{3\left( {1 - {\text{v}}} \right)}}{\text{h}}^3 \frac{{{\text{d}}^3 {\text{w}}}}{{{\text{dx}}\begin{array}{*{20}c}3 \\1 \\\end{array} }} = 0, \hfill \\\frac{{{\text{dw}}}}{{{\text{dx}}1}} - \frac{{4\Theta ^{\text{B}} }}{{3\left( {1 - {\text{v}}} \right)}}{\text{h}}\frac{{{\text{d}}^2 {\text{w}}}}{{{\text{dx}}\begin{array}{*{20}c}2 \\1 \\\end{array} }} + \frac{{4\Theta ^{\text{F}} }}{{3\left( {1 - {\text{v}}} \right)}}{\text{h}}^2 \frac{{{\text{d}}^3 {\text{w}}}}{{{\text{dx}}\begin{array}{*{20}c}3 \\1 \\\end{array} }} = 0, \hfill \\\end{gathered}\]with exponentially small error as L/h→∞, where 2h is the plate thickness and L is the length scale of ω in the x ₁-direction. The four coefficients W ^B, W^F, Θ ^B , Θ ^F are computable constants which depend upon the geometry of the support and the elastic properties of the support and the plate, but are independent of the length of the plate and the loading applied to it. The leading terms in these boundary conditions as L/h→∞ (with all elastic moduli remaining fixed) are the same as those for a thin plate with a clamped edge. However by obtaining asymptotic formulae and general inequalities for Θ ^B , W ^F, we prove that these constants take large values when the support is ‘soft’ and so may still have a strong influence even when h/L is small. The coefficient W ^F is also shown to become large as the size of the support becomes large but this effect is unlikely to be significant except for very thick plates. When h/L is small, the first order corrected boundary conditions are w=0,% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakeaadaWcaaqaaiaabsgacaqG% 3baabaGaaeizaiaabIhaliaacgdaaaGccqGHsisldaWcaaqaaiaacs% dacqqHyoqudaahaaWcbeqaaiaabkeaaaaakeaacaGGZaWaaeWaaeaa% caGGXaGaeyOeI0IaaeODaaGaayjkaiaawMcaaaaacaqGObWaaSaaae% aacaqGKbWaaWbaaSqabeaacaGGYaaaaOGaae4DaaqaaiaabsgacaqG% 4bqbaeqabiqaaaqcaauaaKqzGcGaaiOmaaqcaawaaKqzGcGaaiymaa% aaaaGccqGH9aqpcaGGWaGaaiilaaaa!5DD4!\[\frac{{{\text{dw}}}}{{{\text{dx}}1}} - \frac{{4\Theta ^{\text{B}} }}{{3\left( {1 - {\text{v}}} \right)}}{\text{h}}\frac{{{\text{d}}^2 {\text{w}}}}{{{\text{dx}}\begin{array}{*{20}c}2 \\1 \\\end{array} }} = 0,\]which correspond to a hinged edge with a restoring couple proportional to the angular deflection of the plate at the edge.     

Stability of a thick elastic plate under thrust

When a rectangular plate of incompressible neo-Hookean elastic material is subjected to a thrust, bifurcations of the flexural or barreling types become possible at certain critical values of the compression ratio. The states of pure homogeneous deformation corresponding to these critical compression ratios are states of neutral equilibrium. Their stability is investigated on the basis of an energy criterion, without restriction on the thickness of the plate.The critical state corresponding to the lowest order flexural mode is found to be stable (unstable) if the aspect ratio (thickness/length) is less (greater) than about 0.2. Agreement with the classical Euler theory is established in the limiting case of small aspect ratio.     

On the equation of deflection of a thick plate

Summary Approximate equations of the deflection of a thick plate are derived from fundamental equations of a three-dimensional elastic body, by expanding components of displacement into power series in platethickness h and then by truncating at appropriate terms of 0(h ⁿ ).The method proposed here enables us to give a systematic treatment to obtain approximate equations with any desired accuracy in the sense of increasing order of h. As an example, a thick plate is treated under distributed pressure acting at its upper surface.

---

Über die Gleichung zur Durchbiegung einer dicken Platte

Übersicht Näherungsgleichungen für die Durchbiegung einer dicken Platte werden aus den drei-dimensionalen Fundamentalgleichungen der Elastizität hergeleitet, wobei die Verschiebungen erst nach Potenzreihen der Plattendicke h entwickelt und danach an passenden Termen von 0(h ⁿ ) abgebrochen werden.Die hier präsentierte Methode ermöglicht uns ein systematisches Verfahren, um die Näherungsgleichungen mit einer beliebig gewünschten Genauigkeit im Sinne zunehmender Ordnung von h zu gewinnen. Als ein Beispiel wird eine dicke Platte unter verteilter Belastung an ihrer Deckfläche behandelt.

     

Propagation of two-dimensional plastic waves in a thick plate

A study is made of the propagation and interaction of two-dimensional waves of high amplitude in a thick plate. A monotonically decreasing pressure is applied to the surface of the plate. Deformations are assumed to be large; the problem is formulated and solved in Lagrangian variables. An approximate method for constructing the fronts of the shock waves is proposed. The pressure and particle velocity at an arbitrary point and at an arbitrary instant of time are determined by the method of characteristics. A numerical example is given.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 100–106, January–Febraury, 1971.     

An equation of motion for a thick viscoelastic plate

In this paper an equation of motion is presented for a general thick viscoelastic plate, including the effects of shear deformation, extrusion deformation and rotatory inertia. This equation is the generalization of equations of motion for the corresponding thick elastic plate, and it can be degenerated into several types of equations for various special cases.     

Symmetrical contact problem of a thick transversely isotropic plate

The symmetrical problem of the contact between a spherical indenter and a thick transversely isotropic plate is solved using the techniques of the Hankel transform. Solutions are written as the sums of the associated half-space solution and plate thickness effect terms. The normalized normal contact stress and the surface radial stress are obtained and calculated numerically for both composite materials and metallic substances. The example materials give both real and complex characteristic roots. A general method of calculation is described to determine the maximum tensile radial stress and the maximum compressive normal contact stress.The plate thickness effects on the contact stress and on the tensile radial stress are shown to be increasingly important with decreasing thickness. The effects that the material anisotropy has on the magnification of the contact stress and the maximum tensile radial stress are clearly revealed in the numerical results presented. The effects of material anisotropy for the composite are compared to those for the metallic substances.     

The existence of buckled states on a perforated thin plate

On the basis of the generalized yon Kàrmàn theory for perforated thin plates established in [1, 2], the existence of buckled states for perforated plates subjected to self-equilibrating inplane forces along each edge systematically is investigated. This work completely generalizes the results in [3, 4].     

A 3D shell-like approach using element-free Galerkin method for analysis of thin and thick plate structures

A new efficient meshless method based on the element-free Galerkin method is proposed to analyze the static deformation of thin and thick plate structures in this paper. Using the new 3D shell-like kinematics in analogy to the solid-shell concept of the finite element method, discretization is carried out by the nodes located on the upper and lower surfaces of the structures. The approximation of all unknown field variables is carried out by using the moving least squares (MLS) approximation scheme in the in-plane directions, while the linear interpolation is applied through the thickness direction. Thus, different boundary conditions are defined only using displacements and penalty method is used to enforce the essential boundary conditions. The constrained Galerkin weak form, which incorporates only displacement degrees of freedom (d.o.f.s), is derived. A modified 3D constitutive relationship is adopted in order to avoid or eliminate some self-locking effects. The numeric efficiency of the proposed meshless formulation is illustrated by the numeric examples.     

Bifurcation conditions for a thick elastic plate under thrust

A thick rectangular plate of incompressible isotropic elastic material is subjected to a pure homogeneous deformation by tensile forces or thrusts applied to a pair of opposite faces. The theory of small deformations superposed on finite deformations is applied to determine the critical conditions under which bifurcation solutions (i.e. adjacent equilibrium positions) can exist. The adjacent equilibrium positions considered are those for which the superposed deformation is two-dimensional and is coplanar with the loading force and the thickness direction of the plate, the faces of the plate normal to its thickness being force-free. A number of theorems relating to the critical conditions for superposed deformations of the flexural and barreling types are derived under conditions on the strain-energy function more general than those employed in earlier work. It is also shown how these results can be applied to the determination of the bifurcation conditions corresponding to any specified strain-energy function.     

On the deflection of a thick plate under a sinusoidal load

Summary The equations of deflection of a thick plate presented in a previous paper are solved for an infinite plate under sinusoidally distributed pressure acting at its upper surface.The deflection and stress of the plate are calculated from the solutions of the approximate equations. The results are compared with those calculated from the exact solution and also with those obtained from other approximate methods.

---

Über die Durchbiegung einer dicken Platte unter einer sinusförmigen Last

Übersicht Die in einer früheren Arbeit gegebenen Gleichungen für die Durchbiegung einer dicken Platte werden für den Fall einer unendlichen Platte gelöst, deren Deckfläche durch einen sinusförmigen verteilten Druck belastet wird.Aus den Lösungen der Näherungsgleichungen werden die Durchbiegung und die Spannung einer dicken Platte berechnet. Die Ergebnisse werden mit dem aus der exakten Lösung erlangten Resultat sowie mit den aus anderen Näherungsverfahren erhaltenen Ergebnissen verglichen.