曲线加筋Kirchhoff-Mindlin板自由振动分析

相比传统加筋板,曲线加筋板能够更充分地发挥材料力学性能.在加筋板力学分析中,厚板通常采用Reissner-Mindlin理论,然而当板厚较薄时易出现剪切自锁,离散的Kirchhoff-Mindlin理论采用假设剪切应变场可避免该问题.针对曲线加筋Kirchhoff-Mindlin板自由振动分析,采用离散的Kirchhoff-Mindlin三角形单元和Timoshenko曲梁单元分别模拟板和加强筋,根据板的位移插值函数及筋板交界面的位移协调条件,建立基于板单元位移自由度的有限元方程.为了验证方法的有效性和准确性,采用直线加筋薄板、曲线加筋薄板和厚板3种模型进行算例研究,通过收敛性和精度分析来选择合理的有限元网格密度.直线加筋薄板前20阶固有频率均与文献结果吻合良好;曲线加筋板算例中,本文方法满足收敛条件的板单元数目为2469,Nastran模型板单元数目为6243;本文所得曲线加筋板固有频率与Nastran计算结果最大误差为3.4%.研究结果表明,本文方法无需筋板单元共节点,可使用较少的有限元网格数量,并能够保证计算精度;在离散Kirchhoff-Mindlin三角形板单元基础上构造Timoshenko梁单元可同时适用于曲线加筋薄板与厚板自由振动分析.     

A parametric study for thick plates resting on elastic foundation with variable soil depth

The soil depth is generally considered to be constant for the analysis of plates resting on elastic foundation in the literature. However, it is most reasonable to have a variable subsoil depth as the plate dimensions get larger. In present study, linearly varying subsoil depth is considered as well as constant, linear and quadratic variation of modulus of elasticity with subsoil depth. Also, a parametric study is performed to demonstrate the behavior of thick plates on elastic foundations with variable soil depth. Modified Vlasov Model is used for the analysis of the plate foundation system, and 8-noded Mindlin plate element incorporating shear strain throughout plate thickness is used for the finite element model. Numerical examples are obtained from the literature to compare results and to show the influence of variable soil stratum depth on the behavior of plates. Displacements, bending moments, and shear forces are presented in tabular and graphical formats. As far as results are compared, it can be concluded that variable soil depth significantly affects the variation of the displacements and therefore the internal forces of the plate while keeping it constant ends up with unrealistic results.     

304L/Q235B大面积金属板爆炸焊接物质点法模拟分析

大面积金属板材304L/Q235B的爆炸焊接过程涉及炸药爆轰、金属板材的高速碰撞和塑性变形等。采用有限元法计算模拟这个问题时,网格单元会发生扭曲畸变现象,导致计算精度下降,甚至出现单元负体积而使计算终止,并且炸药爆轰形成气体产物飞散过程也很难模拟。为了能模拟大面积金属板材的爆炸焊接整个过程并获得合理的技术工艺参数,采用物质点法进行三维数值模拟分析。物质点法作为一种无网格法,在模拟冲击动力学问题中主要采用显式积分算法。通过将拉格朗日质点单元与固定的欧拉背景网格相结合,可以实现爆炸焊接的复板与基板的高速碰撞、炸药滑移爆轰、金属板面的塑性变形过程的数值模拟,并给出爆炸复合板材的形变、有效塑性应变和复板与基板的碰撞速度的计算结果。采用物质点法模拟的复合板材变形与爆炸焊接实验结果基本一致。计算复板与基板的碰撞速度这个重要的物理参数时,物质点法与Richter理论公式的相对误差不超过13%。数值计算和实验结果表明,物质点法在数值精度和计算效率方面具有优势,物质点法是研究金属焊接爆炸问题的一种有效数值方法。     

Boundary Element Analysis of Raft Foundations on Piles

The boundary element method is used in the formulation of models for the analysis of raft foundations on piles. Two models are considered: a Kirchhoff plate on a layered elastic half-space and a Kirchhoff plate on a Winkler soil. The plates are modelled using conforming boundary elements and the piles by using linear finite elements. Mindlin's solution is used as influence function within the half-space while Boussinesq's solution, a precursor although a particular case of Mindlin's solution, is used to derive the deflections of the soil surface. The models are used in the analysis of some raft foundations on piles and the results and relative merits are discussed.     

An investigation of the static and dynamic behavior of electrically actuated rectangular microplates

We present an investigation of the static and dynamic behavior of the nonlinear von-Karman plates when actuated by the nonlinear electrostatic forces. The investigation is based on a reduced order model developed using the Galerkin method, which rely on modeshapes and in-plane shape functions extracted using a finite element method. In this study, a fully clamped microplate is considered. We investigate the static behavior and the effect of different non-dimensional design parameters. The static results are validated by comparison with the results calculated by a finite element model. The forced-vibration response of the plate is then investigated when the plate is excited by a harmonic AC load superimposed to a DC load. The dynamic behavior is examined near the primary and secondary (superharmonic and subharmonic) resonances. The microplate shows a strong hardening behavior due to the cubic nonlinearity of mid-plane stretching. However, the behavior switches to softening as the DC load is increased. Finally, near-square plates are studied to understand the effect of geometric imperfections of microplates.     

Large deflection analysis of plates on elastic foundation by the boundary element method

A boundary element method is developed for the large deflection analysis of thin elastic plates resting on elastic foundation. The subgrade reaction may depend linearly (Winkler-type) or nonlinearly on the deflection as well as on the point coordinates (nonhomogeneous subgrade). Moderately large deflections are examined as described by the von Karman equations. The plate may have arbitrary shape and its boundary may be subjected to any type of boundary condition. The proposed method uses the fundamental solution of the linear plate theory and treats the nonlinearities as well as the subgrade reaction as unknown domain forces. Numerical results are presented to illustrate the method and demonstrate its effectiveness and accuracy.     

基于高阶剪切变形理论的四边形求积元板单元及其应用

近年来由各类新型复合材料或功能梯度材料构成的板结构在工程领域得到了广泛应用,其显著特点是材料性能沿板厚变化.为合理考虑横向剪切应变,许多学者基于Reddy高阶剪切变形理论,构建了不同的有限元单元对该类板结构进行分析,但其中满足$C^{1}$连续条件的单元相对较少.本文基于Reddy高阶剪切变形理论,采用求积元方法,建立了$C^{1}$连续的四边形板单元.利用该单元对均质材料、复合材料、功能梯度材料构成的等厚度矩形板、变厚度矩形板及等厚度斜板的线弹性弯曲和自由振动问题进行了计算分析,并与现有文献中的相应计算结果进行了对比.研究表明:基于高阶剪切变形理论的四边形求积元板单元具有较高的计算效率和良好的适应性,文中各类材料构成的等变厚度矩形板及等厚度斜板均只需1个单元即可得到理想的计算结果.对于等/变厚度矩形板,可仅使用9$\times$9个积分点,而对于等厚度斜板,随着斜角的增大,所需积分点的数目逐渐增多至15$\times $15.该四边形求积元板单元可进一步用于新型复合材料板的非线性分析.      

功能梯度板动力分析的光滑有限元法研究

为提高功能梯度板动响应问题求解精度,基于一阶剪切变形板理论,提出了求解功能梯度板自由振动问题的Cell-Based 光滑有限元格式．功能梯度板材料属性沿厚度方向呈梯度连续变化,计算系统刚度矩阵时在光滑域内进行光滑梯度操作,可提高求解精度．采用Cell-Based 光滑有限元法,讨论了长厚比、形状因子和边界条件对两种典型功能梯度方板无量纲频率参数的影响,并与FEM法的解和文献中的解做了对比．数值算例的结果表明,光滑梯度操作可改善有限元系统的刚度,Cell-Based 光滑有限元法的计算结果更加逼近真实解,从而为功能梯度材料的进一步应用打下基础．     

Thermal buckling analysis of truss-core sandwich plates

Truss-core sandwich plates have received much attention in virtue of the high values of strength-to-weight and stifness-to-weight as well as the great ability of impulseresistance recently. It is necessary to study the stability of sandwich panels under the influence of the thermal load. However, the sandwich plates are such complex threedimensional(3D) systems that direct analytical solutions do not exist, and the finite element method(FEM) cannot represent the relationship between structural parameters and mechanical properties well. In this paper, an equivalent homogeneous continuous plate is idealized by obtaining the efective bending and transverse shear stifness based on the characteristics of periodically distributed unit cells. The first order shear deformation theory for plates is used to derive the stability equation. The buckling temperature of a simply supported sandwich plate is given and verified by the FEM. The efect of related parameters on mechanical properties is investigated. The geometric parameters of the unit cell are optimized to attain the maximum buckling temperature. It is shown that the optimized sandwich plate can improve the resistance to thermal buckling significantly.     

无单元法在薄板稳定问题中的应用

用无单元法研究了薄板的弹性稳定问题，从滑动最小二乘法和变分原理出发导出了薄板的无单元法几何刚度矩阵，编制了相应的计算程序，并给出了算例，结果表明，方法合理可行，且精度高于有限元。     

Numerical simulation of the ground effect using the boundary element method

As is well known, the lift of a wing passing over the ground becomes larger than that of a wing in a finite air field because of the ground effect. Owing to its special aerodynamic characteristics and applications, the problem of the ground effect has become increasingly common. In this paper some investigations were conducted to calculate the unsteady aerodynamic forces for long and short ground plates by means of boundary element techniques. In order to calculate the pressure variation on a long ground plate, the steady boundary element method was used. However, when using a short ground plate, the boundary element method was modified to treat the unsteady aerodynamic phenomena. Experimental studies were also made for both ground plates to confirm the validity of the numerical results. At low angles of attack the qualitative behaviour of the unsteady aerodynamic pressure on both ground plates was well predicted by the boundary element methods and qualitative agreement is found between the calculated and measured results. © 1997 John Wiley & Sons, Ltd.     

复合材料正交加筋板在横向载荷下挠度计算的半解析法

本文采用一种新的半解析法,即独特利用Heaviside函数建立与加筋板等效的变刚度模型来开展复合材料双向正交加筋板在横向载荷下的弯曲挠度分析．此模型可以准确地描述筋条在板面上的分布,以及由于筋条的存在而导致的板面刚度不均匀分布．使用Galerkin加权残值法求解该模型的控制方程,得到不同边界条件和载荷情况下的级数解．对于双向正交加筋板,将此半解析法的结果与传统均匀化方法和使用商业有限元软件ABAQUS建立的有限元模型所得到的弯曲挠度结果比较,验证了此方法的准确性和优越性．不同于传统均匀化方法,本双向正交加筋板的弯曲挠度半解析法可精确、有效地获取加筋间的局部弯曲挠度,可以促进复合材料结构的设计分析与优化的研究进展．     

Secondary buckling of a flat plate under uniaxial compression: Part 2. Analysis of clamped plate by F.E.M. and comparison with experiments

In Part 1, a theoretical analysis was used to study the secondary buckling of a simply supported plate. In Part 2, a clamped plate is analysed by the finite element method. The stability criterion for a non-linear post-buckling equilibrium state is evaluated by the sign of the determinant of the stiffness matrix. It should be noted that the secondary buckling loads of clamped plates are unexpectedly smaller than those of simply supported plates and are only one and a half times the primary buckling loads. In previous analyses, only a quarter segment of the plate was considered by assuming a stable equilibrium state with symmetrical mode. However, instability can also be predicted by considering the unsymmetrical mode over the whole plate. Results of experimental analysis of secondary instability of clamped square plates under uniaxial compression agreed with the numerical results.     

基于应变梯度有限元的单层石墨烯振动研究

建立了单层石墨烯等效非局部薄板的一种新的有限元模型,并运用有限元法分析不同边界条件下单层石墨烯振动的小尺度效应。给出了基于弹性应变梯度理论下Kirchhoff板的振动方程。发展了一种4节点24自由度的板单元,用于离散化求解考虑微纳结构尺度效应的高阶微分方程。在研究四边简支板振动时,考虑应变梯度的非局部弹性有限元数值计算结果与理论分析结果相一致。用有限元方法研究了不同尺寸、振动波长、振动模态阶数、边界条件类型以及非局部参数的单层石墨烯振动。     

Nonlinear transient response of strain rate dependent composite laminated plates using multiscale simulation

The effects of strain rate dependency and inelasticity on the transient responses of composite laminated plates are investigated. A micromechanics model which accounts for the transverse shear stress effect, the effect of strain rate dependency and the effect of inelasticity is used for analyzing the mechanical responses of the fiber and matrix constituents. The accuracy of the micromechanics model under transverse shear loading is verified by comparing the results with those obtained using a general purpose finite element code. A higher order laminated plate theory is extended to capture the inelastic deformations of the composite plate and is implemented using the finite element technique. A complete micro–macro numerical procedure is developed to model the strain rate dependent behavior of inelastic composite laminates by implementing the micromechanics model into the finite element model. Parametric studies of the transient responses of composite plates are conduced. The effects of geometry, ply stacking sequence, material models, boundary conditions and loadings are investigated. The results show that the strain rate dependency and inelasticity influence the transient responses of composite plates via two significantly different mechanisms.     

弯矩作用下薄板振动的仿真分析

针对工程中广泛存在的弯矩作用下板结构的振动问题,分别建立四边固定、三边固定一边自由、两对边固定两对边自由的3种不同边界条件下板的振动模型;基于有限元法计算振动板的频率响应,计算了各节点振动速度的平方和;运用部分追加法正交试验方案分别对板厚、边界条件、板的损耗因子三因素三水平、激励点位置四水平进行了正交试验。实例结果表明,以速度平方和的大小为目标,影响薄板振动的主次因素顺序为:损耗因子、板厚、边界、激励点位置;当振动板两对边固定另两对边自由支承、板厚为0.014m、在(0.5m,0.4m)点处激励且板的损耗因子为0.0008时,此组合为薄板结构振动最小的最优组合。     

模态分析和有限元分析相结合识别材料结构刚度

本文论述了利用模态分析和有限元分析相结合来识别复合材料板的刚度系数的方法。该方法取决于（１）正确的有限元模型；（２）可靠的实验模态分析数据和正确的相关准则；（３）快速而又稳妥的估算方法，该方法对有限元模型动力修正也是有用的     

Elasto-plastic buckling analysis for perforated steel plates subject to uniform compression

This study investigates the elasto-plastic buckling behaviour of simply supported square and rectangular thin steel plates having elliptic cut-outs by means of finite element method. Plates with simply supported in the out-of-plane direction are applied uniform compression in long-edge direction. A50 steel was used in the analysis and the focus was on the effect of plate aspect ratio, elliptical hole size, elliptical hole angle, elliptical hole location and slenderness ratio on buckling behaviour. It was found in the study that as the plate slenderness ratio increases, the critical buckling stress decreases for all the perforated plates.     

Transient dynamic response of debonded sandwich plates predicted with finite element analysis

Dynamic transient response of a composite sandwich plate with a penny-shaped debonded zone has been studied by using the finite element analysis within the ABAQUS/Explicit code in this paper. In order to accurately predict the response of the debonded sandwich plate to impulsive loading, contact–impact and sliding conditions along the damaged skin-to-core interface were imposed in the model through a kinematic predictor/corrector contact algorithm. The accuracy of the finite element (FE) model used was verified by comparing between numerical predictions and experimental data known in literature for the frequency spectrum of a cracked polycarbonate laminated beam containing a delamination. By analyzing nonlinear aspects of the transient dynamics of the sandwich plate, it is shown that the presence of the debond significantly alters its short-term response. In this respect, a considerable influence of contact events within the debonded region on the plate’s global dynamic response was found out. These results were presented in both time and frequency domains. The predictions performed also showed that the FE model applied would be useful for nondestructive evaluation of defects in composite sandwich plates, and for studying dynamic response of such plates to impact.     

Static analysis of simply supported functionally graded and layered magneto-electro-elastic plates

In this article, static analysis of functionally graded, anisotropic and linear magneto-electro-elastic plates have been carried out by semi-analytical finite element method. A series solution is assumed in the plane of the plate and finite element procedure is adopted across the thickness of the plate such a way that the three-dimensional character of the solution is preserved. The finite element model is derived based on constitutive equation of piezomagnetic material accounting for coupling between elasticity, electric and magnetic effect. The present finite element is modeled with displacement components, electric potential and magnetic potential as nodal degree of freedom. The other fields are calculated by post-computation through constitutive equation. The functionally graded material is assumed to be exponential in the thickness direction. The numerical results obtained by the present model are in good agreement with available functionally graded three-dimensional exact benchmark solutions given by Pan and Han [Pan, E., Han, F., in press. Green’s function for transversely isotropic piezoelectric functionally graded multilayered half spaces. Int. J. Solids Struct.]. Numerical study includes the influence of the different exponential factor, magneto-electro-elastic properties and effect of mechanical and electric type of loading on induced magneto-electro-elastic fields. In addition further study has been carried out on non-homogeneous transversely isotropic FGM magneto-electro-elastic plate available in the literature [Chen, W.Q., Lee, K.Y., Ding, H.J., 2005. On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic plates].