Dynamic geometrically nonlinear analysis of transversely compressible sandwich plates

In order to construct a plate theory for a thick transversely compressible sandwich plate with composite laminated face sheets, the authors make simplifying assumptions regarding distribution of transverse strain components in the thickness direction. The in-plane stresses and σ_yy (Fig. 1) are computed from the constitutive equations, and the improved values of transverse stress components and σ_zz need to be computed by integration of pointwise equations of motion in a post-process stage of the finite element analysis. The improved values of the transverse strains can also be computed in the post-process stage by substituting the improved transverse stresses into the constitutive relations. A problem of cylindrical bending of a simply supported plate under a uniform load on the upper surface is considered, and comparison is made between the displacements, the in-plane stress and the improved transverse stresses (obtained by integration of the pointwise equations of motion), computed from the plate theory, with the corresponding values of exact elasticity solutions. In this comparison, a good agreement of both solutions is achieved. In the finite element analysis of sandwich plates in cylindrical bending with small thickness-to-length ratios, the shear locking phenomenon does not occur. The model of a sandwich plate in cylindrical bending, presented in this paper, has a wider range of applicability than the models presented in literature so far: it can be applied to the sandwich plates with a wide range of ratios of thickness to the in-plane dimensions, with both thin and thick face sheets (as compared to the thickness of the core) and to the sandwich plates with both transversely rigid and transversely compressible face sheets and cores.     

FRP strengthened walls with delaminated regions – Geometrically nonlinear response

The paper studies the geometrically nonlinear behavior of walls that are strengthened with fiber reinforced polymer (FRP) composite materials but include pre-existing delaminated regions. The paper uses an analytical–numerical methodology. Three specially tailored finite elements that correspond to perfectly bonded regions, to delaminated regions where the debonded layers are in contact, and to delaminated regions where the debonded layers are not in contact are presented. All finite elements are based on a high order multi layered plate theory. The geometrical nonlinearity is introduced by means of the Von Karman nonlinear strains whereas the contact nonlinearity is handled iteratively. The validity and convergence of the finite element models is demonstrated for each type of element through comparison with closed form analytical solutions available for specific cases. The unified model that combines the three types of finite element is then used for studying the nonlinear behavior of a locally delaminated FRP strengthened wall under in-plane normal and in-plane shear loads. Finally, conclusions regarding the effect of the delamination on the response of the strengthening system, on the conditions that evolve in the bonded region that surrounds the delamination, and on the global response of the multi-layered structure are drawn. Additional conclusions regarding the application of the modeling approach to other delamination sensitive layered structural systems close the paper.     

Thermal buckling of functionally graded circular plates based on higher order shear deformation plate theory

In this research, thermal buckling of circular plates compose of functionally graded material (FGM) is considered. Equilibrium and stability equations of a FGM circular plate under thermal loads are derived, based on the higher order shear deformation plate theory (3rd order plate theory). Assuming that the material properties vary as a power form of the thickness coordinate variable z and using the variational method, the system of fundamental partial differential equations is established. A buckling analysis of a functionally graded circular plate (FGCP) under various types of thermal loads is carried out and the result are given in closed-form solutions. The results are compared with the critical buckling temperature obtained for FGCP based on first order (1st order plate theory) and classical plate theory (0 order plate theory) given in the literature. The study concludes that higher order shear deformation theory accurately predicts the behavior of FGCP, whereas the first order and classical plate theory overestimates buckling temperature.     

复合材料层合板精化高阶理论及其精化三角形板单元

提出一种新的精化高阶理论，该理论满足层间位移、应力连续条件，由此建立了三角形精化板单元。该单元满足单元间Ｃ１类弱连续条件，其收敛性得到保证，且具有简单、高效率的优点。     

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.     

基于自然单元法的功能梯度中厚板自由振动分析

自然单元法是一种基于自然邻近插值的无网格数值方法.相对于移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,也不需要任何人为的参数.基于一阶剪切变形板理论,利用自然单元法对功能梯度中厚板的自由振动进行了数值分析.功能梯度板材料属性沿厚度方向呈梯度连续变化.由于自然邻近插值函数具有Kronecker delta函数性质,可以直接施加本质边界条件.通过本文给出的方法,对不同梯度指数和不同边界条件的功能梯度中厚板的振动频率进行了计算.通过与文献结果的对比验证了自然单元法求解的有效性.     

Some results for approximate strain and rotation tensor formulations in geometrically non-linear Reissner-Mindlin plate theory

Finite element deflection and stress results are presented for four flat plate configurations and are computed using kinematically approximate (rotation tensor, strain tensor or both) non-linear Reissner-Mindlin plate models. The finite element model is based on a mixed variational principle and has both displacement and force field variables. High order interpolation of the field variables is possible through p-type discretization. Results for some of the higher order approximate models are given for what appears to be the first time. It is found that for the class of example problems examined, exact strain tensor but approximate rotation tensor theories can significantly improve the solution over approximate strain tensor models such as the von Kármán and moderate rotation models when moderate deflections/rotations are present. However, for each of the problems examined (with the exception of a postbuckling problem) the von Kármán and moderate rotation model results compared favorably with the higher order models for deflection magnitudes which could be reasonably expected in typical aeroelastic configurations.     

A comparative analysis of numerical approaches to the mechanics of elastic sheets

Numerically simulating deformations in thin elastic sheets is a challenging problem in computational mechanics due to destabilizing compressive stresses that result in wrinkling. Determining the location, structure, and evolution of wrinkles in these problems has important implications in design and is an area of increasing interest in the fields of physics and engineering. In this work, several numerical approaches previously proposed to model equilibrium deformations in thin elastic sheets are compared. These include standard finite element-based static post-buckling approaches as well as a recently proposed method based on dynamic relaxation, which are applied to the problem of an annular sheet with opposed tractions where wrinkling is a key feature. Numerical solutions are compared to analytic predictions of the ground state, enabling a quantitative evaluation of the predictive power of the various methods. Results indicate that static finite element approaches produce local minima that are highly sensitive to initial imperfections, relying on a priori knowledge of the equilibrium wrinkling pattern to generate optimal results. In contrast, dynamic relaxation is much less sensitive to initial imperfections and can generate low-energy solutions for a wide variety of loading conditions without requiring knowledge of the equilibrium solution beforehand.     

基于Layerwise离散层理论的复合材料夹层板屈曲分析

复合材料夹层结构由于面板和芯层力学特性差异较大,屈曲分析时要分层考虑各层的剪切变形.基于Reddy的Layerwise离散层理论,假设每一层变形服从一阶剪切变形理论,在统一的位移场描述下,推导建立了一种用于复合材料夹层结构屈曲分析的四节点四边形板单元,并采用混合插值方法对单元的剪切锁定进行了修正.分别对三种典型的夹层板结构进行线性屈曲有限元分析,并将计算结果与文献中已有结果进行了对比.结果表明:论文的分析方法能离散考虑各层的力学特性,将结构离散为多层时,计算结果与三维弹性理论或高阶板理论吻合;将结构等效为单层时,计算结果与基于一阶剪切变形理论的文献结果吻合,验证了单元的有效性.     

矩形板屈曲问题的一个小波解

利用Wavelet-Galerkin法分析了四边固支与四边简支矩形板的屈曲问题.以小波作为基函数表示板的挠度,推导出屈曲系数及屈曲模态的计算过程.数值计算给出了不同边长比的矩形板的屈曲系数及屈曲半波数.与传统的三角函数作基函数的Galerkin法及有限元法结果比较,结果表明在一定条件下小波可以作为试函数解决结构力学的屈...     

新的圆柱曲板剪切屈曲实验的研究

提出一个用两块曲板以单向加载进行的Ｒ１０００ｍｍ的圆柱曲板剪切屈曲新实验方案。其原理清楚、装置简单紧凑、节省试材、实验进行方便和不需扭转设备等优点,而实验结果和成熟的工作符合良好．     

ANALYTICAL RELATIONS BETWEEN EIGENVALUES OF CIRCULAR PLATE BASED ON VARIOUS PLATE THEORIES

Based on the mathematical similarity of the axisymmetric eigenvalue problems of a circular plate between the classical plate theory(CPT), the first-order shear deformation plate theory(FPT) and the Reddy's third-order shear deformation plate theory (RPT), analytical relations between the eigenvalues of circular plate based on various plate theories are investigated. In the present paper, the eigenvalue problem is transformed to solve an algebra equation. Analytical relationships that are expressed explicitly between various theories are presented. Therefore, from these relationships one can easily obtain the exact RPT and FPT solutions of critical buckling load and natural frequency for a circular plate with CPT solutions. The relationships are useful for engineering application, and can be used to check the validity, convergence and accuracy of numerical results for the eigenvalue problem of plates.     

An adaptive remeshing strategy for viscoelastic fluid flow simulations

In the last few years, we have developed a fairly general adaptive finite element solution procedure which can be applied to a large variety of problems. In this paper, this strategy is briefly recalled and applied to the solution of two-dimensional viscoelastic fluid flow problems. A log-conformation formulation recently introduced by Fattal and Kupferman [R. Fattal, R. Kupferman, Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation, J. Non-Newtonian Fluid Mech. 126 (2005) 23-37] was implemented in order to improve the convergence properties of the numerical scheme. We confirm some results obtained in Hulsen, Fattal and Kupferman [M. Hulsen, R. Fattal, R. Kupferman, Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithm, J. Non-Newtonian Fluid Mech. 127 (2005) 27-39] and in some instances, we show that mesh adaptation allows to almost automatically reproduce accurate results obtained on very fine structured meshes.     

板桁结构考虑局部屈曲空间计算的横向条带板段单元法

用有限个横向条带法构造了板桁组合结构板段考虑局部屈曲的空间位移模式.基于三维连续体的虚功增量方程,导出了横向条带板段单元的UL列式,并考虑了板段单元位形变化的影响.此计算方法自由度少,计算精度高,能用于大型板桁结构的几何非线性分析.文末计算了广东西江桥板桁组合结构模型梁,计算结果与实测结果吻合较好.     

Variational formulation of the static Levinson beam theory

In this communication, we provide a consistent variational formulation for the static Levinson beam theory. First, the beam equations according to the vectorial formulation by Levinson are reviewed briefly. By applying the Clapeyron's theorem, it is found that the stresses on the lateral end surfaces of the beam are an integral part of the theory. The variational formulation is carried out by employing the principle of virtual displacements. As a novel contribution, the formulation includes the external virtual work done by the stresses on the end surfaces of the beam. This external virtual work contributes to the boundary conditions in such a way that artificial end effects do not appear in the theory. The obtained beam equations are the same as the vectorially derived Levinson equations. Finally, the exact Levinson beam finite element is developed.     

Coarse-graining atomistic dynamics of brittle fracture by finite element method

In this work we present the finite element (FE) implementation of an atomistic formulation of balance equations and its application to coarse-grained (CG) simulation of dynamic fracture. First, we simulate a notched specimen that contains about 1.8 million atoms by the CG-FE method, and we compare the CG-FE results with that by all-atom molecular dynamics (MD) simulations. We find that CG-FE simulations with about 5% degrees of freedom of the MD simulation can capture the essential dynamic features, not in exact correspondence, but qualitatively and quantitatively similar to that obtained by MD simulations. We then proceed to simulate a series of micron-sized specimens by the CF-FE method. We find that it is the interaction of the forward propagating crack with the stress waves being reflected back by the boundaries of the specimen that triggers the dynamic instability and hence the branching of cracks in micron-sized specimens. The potential application of the method and future work for improvements are discussed.     

Geometric nonlinear formulation and discretization method for a rectangular plate undergoing large overall motions

In the present paper, the geometric nonlinear formulation is developed for dynamic stiffening of a rectangular plate undergoing large overall motions. The dynamic equations, which take into account the stiffening terms, are derived based on the virtual power principle. Finite element method is employed for discretization of the plate. The simulation results of a rotating rectangular plate obtained by using such geometric nonlinear formulation are compared with those obtained by the conventional linear method without consideration of the stiffening effects. The application limit of the conventional linear method is clarified according to the frequency error. Furthermore, the accuracy of the assumed mode method is investigated by comparison of the results obtained by using the present finite element method and those obtained by using the assumed mode method.     

一种有效的厚薄板壳单元

对于厚薄板壳分析提出了一种有效的拟协调位移型三角型单元。采用Ｍｉｎｄｌｉｎ变形假定，在横向剪切刚度项中引入特殊的罚函数，有效地解除了剪切闭锁和抑制了零能模式，并扩大了对跨厚比的适用范围，计算精度也明显提高。文中给出了各种材料和形状板壳的线性和非线性分析。