Mathematical construction of a Reissner–Mindlin plate theory for composite laminates

A Reissner–Mindlin theory for composite laminates without invoking ad hoc kinematic assumptions is constructed using the variational-asymptotic method. Instead of assuming a priori the distribution of three-dimensional displacements in terms of two-dimensional plate displacements as what is usually done in typical plate theories, an exact intrinsic formulation has been achieved by introducing unknown three-dimensional warping functions. Then the variational-asymptotic method is applied to systematically decouple the original three-dimensional problem into a one-dimensional through-the-thickness analysis and a two-dimensional plate analysis. The resulting theory is an equivalent single-layer Reissner–Mindlin theory with an excellent accuracy comparable to that of higher-order, layer-wise theories. The present work is extended from the previous theory developed by the writer and his co-workers with several sizable contributions: (a) six more constants (33 in total) are introduced to allow maximum freedom to transform the asymptotically correct energy into a Reissner–Mindlin model; (b) the semi-definite programming technique is used to seek the optimum Reissner–Mindlin model. Furthermore, it is proved the first time that the recovered three-dimensional quantities exactly satisfy the continuity conditions on the interface between different layers and traction boundary conditions on the bottom and top surfaces. It is also shown that two of the equilibrium equations of three-dimensional elasticity can be satisfied asymptotically, and the third one can be satisfied approximately so that the difference between the Reissner–Mindlin model and the second-order asymptotical model can be minimized. Numerical examples are presented to compare with the exact solution as well as the classical lamination theory and the first-order shear-deformation theory, demonstrating that the present theory has an excellent agreement with the exact solution.     

考虑电场梯度的挠曲电纳米板弯曲性能分析

具有尺度依赖的挠曲电效应在器件的设计中扮演着越来越关键的角色, 研究人员在微纳米尺度多物理场分析中进行了大量工作. 基于考虑挠曲电和电场梯度效应的弹性介电材料非经典理论, 以二维纳米板为例, 通过理论建模, 分析纳米板在弯曲问题中的力?电耦合行为. 根据Mindlin假设给出板的位移场和电势场的一阶截断, 选取板的材料为立方晶体(m3m点群), 将广义三维本构方程代入到高阶应力、高阶偶应力、高阶电位移和高阶电四极矩的表达式中得到相应的二维本构方程, 利用弹性电介质变分原理得到板的控制方程和边界上的线积分等式, 分别将二维本构方程和边界上外法线的方向余弦代入, 得到板的高阶弯曲方程、高阶电势方程以及对应的四边简支边界条件. 利用四边简支矩形板的高阶弯曲方程、高阶电势方程和相应的边界条件, 根据Navier解理论, 求解纳米板的电势场, 重点分析电场梯度对板内一阶电势的影响. 数值计算结果表明: 电场梯度对纳米板中由挠曲电效应产生的一阶电势有削弱作用, 且材料参数g₁₁越大, 一阶电势受到的削弱越大; 同时电场梯度的存在消除了纳米板在受横向集中载荷作用时一阶电势的奇异性. 本文是对具有挠曲电效应和电场梯度效应的纳米板结构分析理论的一个扩展, 为微纳米尺度器件的结构设计提供参考.      

A new simple third-order shear deformation theory of plates

An improved simple third-order shear deformation theory for the analysis of shear flexible plates is presented in this paper. This new plate theory is composed of three parts: the simple third-order kinematics of displacements reduced from the higher-order displacement field derived previously by the author; a system of 10th-order differential equilibrium equations in terms of the three generalized displacements of bending plates; five boundary conditions at each edge of plate boundaries. Although the resulting displacement field is the same as that proposed by Murthy, the variational consistent governing equations and the associated proper boundary conditions are derived and identified in this work for the first time in the literature. The applications and accuracy of the present shear deformation theory of plates are demonstrated by analytically solving the differential governing equations of a twisting plate, a bending beam and two bending plates to which the 3-D elasticity solutions are available, and excellent agreements are achieved even for the torsion of a plate with square cross-section as well the local effects of stresses at plate boundaries can be characterized accurately. These analytical solutions clearly show that the simple third-order shear deformation theory developed in this work indeed gives better results than the first-order shear deformation theories and other simple higher-order shear deformation theories, since the present third-order shear flexible theory is based on a more rigorous kinematics of displacements and consists of not only a system of variational consistent differential equations, but also a group of consistent boundary conditions associated with the differential equations. The present simple third-order shear deformation theory can easily be applied to the static and dynamic finite element analysis of laminated plates just like the applications of other popular shear flexible plate theories, and improved results could be obtained from the present simple third-order shear deformable theories of plates.     

Scattering of flexural wave in a thin plate with multiple circular holes by using the multipole Trefftz method

The scattering of flexural wave by multiple circular holes in an infinite thin plate is analytically solved by using the multipole Trefftz method. The dynamic moment concentration factor (DMCF) along the edge of circular holes is determined. Based on the addition theorem, the solution of the field represented by multiple coordinate systems centered at each circle can be transformed into one coordinate system centered at one circle, where the boundary conditions are given. In this way, a coupled infinite system of simultaneous linear algebraic equations is derived as an analytical model for the scattering of flexural wave by multiple holes in an infinite plate subject to the incident flexural wave. The formulation is general and is easily applicable to dealing with the problem containing multiple circular holes. Although the number of hole is not limited in our proposed method, the numerical results of an infinite plate with three circular holes are presented in the truncated finite system. The effects of both incident wave number and the central distance among circular holes on the DMCF are investigated. Numerical results show that the DMCF of three holes is larger than that of one, when the space among holes is small and meanwhile the specified direction of incident wave is subjected to the plate.     

An approximate analysis of thickness-stretch waves in an elastic plate

Two-dimensional equations for coupled extensional and thickness-stretch waves in an elastic plate are simplified by eliminating the extensional displacements in a systematic manner; the result is a single equation governing thickness-stretch motions. A similar reduction is also performed for coupled extensional, thickness-stretch, and symmetric thickness-shear waves. The procedure is similar to that used in the thickness-shear approximation, wherein the flexural displacement is eliminated from the equations for coupled flexural and thickness-shear motions. The resulting equations are used to discuss the energy-trapped vibration of plates in thickness-stretch modes.     

Mindlin 矩形微板的热弹性阻尼解析解

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.      

ANLYTICAL SOLUTION OF THERMOELASTIC DAMPING IN RECTANGULAR MINDLIN MICRO PLATES 1)

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.     

Free vibration analysis of moderately thick functionally graded plates on elastic foundation using the extended Kantorovich method

Free vibration analysis of moderately thick rectangular FG plates on elastic foundation with various combinations of simply supported and clamped boundary conditions are studied. Winkler model is considered to describe the reaction of elastic foundation on the plate. Governing equations of motion are obtained based on the Mindlin plate theory. A semi-analytical solution is presented for the governing equations using the extended Kantorovich method together with infinite power series solution. Results are compared and validated with available results in the literature. Effects of elastic foundation, boundary conditions, material, and geometrical parameters on natural frequencies of the FG plates are investigated.     

应变梯度Mindlin板边值问题的讨论

实验和分子动力学计算结果表明, 当材料/结构的特征尺寸降为微纳米量级时, 他们将表现出明显的尺度效应, 因此能否建立精确表征其力学行为的连续介质力学模型具有重要的理论和现实意义. 尽管现有文献对非经典Mindlin板的力学行为进行了大量研究, 但该模型的变分自洽的边值问题是近年来未攻克的科学问题之一. 基于简化的应变梯度理论给出了各向同性Mindlin板应变能的表达式, 通过变分原理和张量分析, 得到了Mindlin板变分自洽的边值问题及其对应角点条件的位移微分表达式. 本文Mindlin板模型的边值问题可退化为相应的Timoshenko梁和Kirchhoff板模型的边值问题, 验证了本文结果的有效性. 研究结果发现, 该Mindlin板模型的控制方程是一个解耦后横向振动具有12阶的偏微分方程, 因此需要每个板边提供6个边界条件. 角点条件由双应力(double stress)产生, 并与经典的剪力、弯矩和扭矩沿截面的法向梯度有关. 本文首次澄清了应变梯度Mindlin板存在角点条件这一事实, 所得的变分结果有望为其有限元法和伽辽金法等数值方法提供理论依据.      

Early short time transient response of finite L-shaped Mindlin plate

The early short time transient response of an L-shaped plate is investigated by the method of reverberation ray matrix (MRRM) and the generalized ray method (GRM). The MRRM and GRM have been successfully used to study the transient elastic wave transmitting in the beams, planar trusses, space frames and infinite layered media. In the present paper, these two methods are extended to investigate the flexural wave propagation in a finite L-shaped Mindlin plate. The reverberation-ray matrix representing the multi-reflected and scattered waves in the L-shaped plate is presented. The rays in the GRM propagating in the L-shaped plate are determined by the shock source and the boundary conditions. From the numerical results, it indicates that the early short time transient responses of L-shaped plate calculated by the MRRM are the same as those by the GRM. The early short time transient acceleration responses are very large, while the early short time transient displacement responses are very small for the L-shaped plate subjected to the unit impulse. Furthermore, the effects of the different shock signals and the plate thickness on the early short time transient responses of the finite L-shaped Mindlin plates are also analyzed.     

Higher-order shear deformable theories for flexure of sandwich plates—Finite element evaluations

A simple isoparametric finite element formulation based on a higher-order displacement model for flexure analysis of multilayer symmetric sandwich plates is presented. The assumed displacement model accounts for non-linear variation of inplane displacements and constant variation of transverse displacement through the plate thickness. Further, the present formulation does not require the fictitious shear correction coefficient(s) generally associated with the first-order shear deformable theories. Two sandwich plate theories are developed: one in which the free shear stress conditions on the top and bottom bounding planes are imposed and another, in which such conditions are not imposed. The validity of the present development(s) is established through, numerical evaluations for deflections/stresses/stress-resultants and their comparisons with the available three-dimensional analyses/closed-form/other finite element solutions. Comparison of results from thin plate. Mindlin and present analyses with the exact three-dimensional analyses yields some important conclusions regarding the effects of the assumptions made in the CPT and Mindlin type theories. The comparative study further establishes the necessity of a higher-order shear deformable theory incorporating warping of the cross-section particularly for sandwich plates.     

基于应变梯度中厚板单元的石墨烯振动研究

基于应变梯度理论建立了单层石墨烯等效明德林(Mindlin) 板动力学方程,推导了四边简支明德林中厚板自由振动固有频率的解析解. 提出了一种考虑应变梯度的4 节点36 自由度明德林板单元,利用虚功原理建立了单层石墨烯的等效非局部板有限元模型. 通过对石墨烯振动问题的研究,验证了应变梯度有限元计算结果的收敛性. 运用该有限元法研究了尺寸、振动模态阶数以及非局部参数对石墨烯振动特性的影响. 研究表明,这种单元能够较好地适用于研究考虑复杂边界条件石墨烯的尺度效应问题. 基于应变梯度理论的明德林板所获得石墨烯的固有频率小于基于经典明德林板理论得到的结果. 尺寸较小、模态阶数较高的石墨烯振动尺度效应更加明显. 无论采用应变梯度理论还是经典弹性本构关系,考虑一阶剪切变形的明德林板模型预测的固有频率低于基尔霍夫(Kirchho) 板所预测的固有频率.      

Stability and vibration of shear deformable plates––first order and higher order analyses

This work presents the highly accurate numerical calculation of the natural frequencies and buckling loads for thick elastic rectangular plates with various combinations of boundary conditions. The Reissener–Mindlin first order shear deformation plate theory and the higher order shear deformation plate theory of Reddy have been applied to the plate’s analysis. The governing equations and the boundary conditions are derived using the dynamic version of the principle of minimum of the total energy. The solution is obtained by the extended Kantorovich method. This approach is combined with the exact element method for the vibration and stability analysis of compressed members, which provides for the derivation of the exact dynamic stiffness matrix including the effect of in-plane and inertia forces. The large number of numerical examples demonstrates the applicability and versatility of the present method. The results obtained by both shear deformation theories are compared with those obtained by the classical thin plate’s theory and with published results. Many new results are given too.     

On consistent plate theories

Summary  Applying the uniform-approximation technique, consistent plate theories of different orders are derived from the basic equations of the three-dimensional linear theory of elasticity. The zeroth-order approximation allows only for rigid-body motions of the plate. The first-order approximation is identical to the classical Poisson-Kirchhoff plate theory, whereas the second-order approximation leads to a Reissner-type theory. The proposed analysis does not require any a priori assumptions regarding the distribution of either displacements or stresses in thickness direction. Received 10 January 2002; accepted for publication 16 April 2002     

Analytical solution for the cylindrical bending vibration of piezoelectric composite plates

An analytical solution for the cylindrical bending vibrations of linear piezoelectric laminated plates is obtained by extending the Stroh formalism to the generalized plane strain vibrations of piezoelectric materials. The laminated plate consists of homogeneous elastic or piezoelectric laminae of arbitrary thickness and width. Fourier basis functions for the mechanical displacements and electric potential that identically satisfy the equations of motion and the charge equation of electrostatics are used to solve boundary value problems via the superposition principle. The coefficients in the infinite series solution are determined from the boundary conditions at the edges and continuity conditions at the interfaces between laminae, which are satisfied in the sense of Fourier series. The formulation admits different boundary conditions at the edges of the laminated piezoelectric composite plate. Results for laminated elastic plates with either distributed or segmented piezoelectric actuators are presented for different sets of boundary conditions at the edges.     

两相材料V形切口应力强度因子边界元分析

建立了边界元法计算两相材料粘结V形切口奇异应力场的新途径。在V形切口尖端挖出一小扇形,将该扇形弧线边界的位移和面力表示为有限项奇性指数和特征角函数的线性组合,其组合系数即为广义应力强度因子,将该组合回代到在被挖去小扇形后的剩余结构内建立的边界积分方程,离散后可求解出组合系数,获得两相材料粘结V形切口尖端的应力强度因子。算例证明了本文方法的有效性。     

功能梯度中厚圆/环板轴对称弯曲问题的解析解

基于一阶剪切变形板理论，导出了热/机载荷作用下，位移形式的功能梯度 中厚圆/环板轴对称弯曲问题的控制方程，获得了问题的位移和内力的一般解析解. 作为特 例，分别研究了边缘径向固定和可动的夹紧和简支的4种实心功能梯度圆板，给出了它们的 解，并分析了热/机载荷作用下解的形态，讨论了横向剪切变形、材料梯度常数和边界条件， 对板的轴对称弯曲行为的影响.     

Diffraction of surface waves by floating elastic plates

Based on the dynamical theories of water waves and Mindlin thick plates, the diffraction of surface waves by a floating elastic plate is presented by using the Wiener–Hopf technique. Firstly, the problem is related to a wave guide in water of finite depth, which is analysed to determine the poles. The resulting hybrid boundary value problem is reduced to solving an infinite system of linear algebraic equations. The results obtained are compared with those calculated by an alternative analysis, and with experimental data. Finally, the effects of the geometric and physical parameters on the distribution of deflection and bending moments in plates are analysed and discussed.     

ANALYSIS AND EXPERIMENT ON TRANSIENT DYNAMIC RESPONSE IN FINITE MINDLIN PLATE

The transient wave propagation in the finite rectangular Mindlin plate is investigated by the analytical and experimental methods. The generalized ray method (GRM) which has been successfully applied to study the transient responses of beams, planar trusses, space frames and infinite layered media is extended to investigate the transient wave propagation and early short time transient response in finite Mindlin plate. Combining the wave solution, the shock source and the boundary conditions, the ray groups transmitted in the finite rectangular plate can be determined. Numerical simulations and experiments are performed and compared with each other. The results show that the transient wave propagation and early short time transient responses in the finite plate can be studied using the GRM. The early short time transient accelerations are very large for the finite plate subjected to the unit impulse, while the early short time transient displacements are very small. The early short time transient accelerations under the unit impulse are much larger than those under the unit step impulse. The thickness and material characteristics have remarkable effects on the early short time transient responses.     

Zum Randwertproblem der schwingenden schubelastischen Platte

Übersicht Über das Prinzip von Hamilton werden mittels verschiedener Ritzansätze für die Verschiebungen bzw. Spannungen Bewegungsgleichungen für die schwingende Platte hergeleitet. Abhängig von der Phasengeschwindigkeit ebener Biegewellen wird ein Schubbeiwert definiert, der im Vergleich mit der exakten Lösung der unendlichen Platte als Gütekriterium dient.

---

On the boundary value problem of vibrating plates

Summary Following Hamilton's principle the equations of motion of a vibrating plate are derived by means of an expansion of the displacements respectively stresses using Ritz method. Depending upon the wave velocity of straight-crested flexural waves a shear coefficient is defined to check the above solutions with the known exact solution for the infinite plate.