Scattering of waves by flaws in anisotropic laminated plates

The strip element method (SEM) is used to investigate wave scattering by rectangular flaws in anisotropic laminated plates. The plates containing flaws are divided by junctions into domains in which the SEM is applied. For each domain, SEM equations are obtained which give a relationship between the traction and displacement vector on the vertical boundaries. A set of equations which gives a relationship between the traction and displacement vector on the junctions is then obtained by assembling the SEM equations for all the domains. This set of equations is solved by using the conditions on the junctions. Scattered wave fields in the frequency domain for isotropic and composite plates with rectangular flaws are computed and discussed in comparison with results for corresponding plates without flaws. A technique for determining the length of a rectangular flaw in a plate is also presented. The results presented in this paper are of importance and could be used in the characterization of flaws in anisotropic laminated plates.     

Advanced variable kinematics Ritz and Galerkin formulations for accurate buckling and vibration analysis of anisotropic laminated composite plates

Accurate free-vibrations and linearized buckling analysis of anisotropic laminated plates with different lamination schemes and simply supported boundary condition are addressed in this paper. Approximation methods, such as Rayleigh-Ritz, Galerkin and Generalized Galerkin, based on Principle of Virtual Displacement are derived in the framework of Carrera’s Unified Formulation (CUF). CUF widely used in the analysis of composite laminate beams, plates and shells, have been here formulated both for the same and different expansion orders, for the displacement components, in the thickness layer-plate direction. An extensive assessment of advanced and refined plate theories, which include Equivalent single Layer (ESL), Zig-Zag (ZZ) and Layer-wise (LW) models, with increasing number of displacement variables is provided. Accuracy of the results is shown to increase by refining the theories. Convergence studies are made in order to demonstrate that accurate results are obtained examining thin and thick plates using trigonometric approximation functions. The effects of boundary terms, upon frequency parameters and critical loads are evaluated. The effects of the various parameters (material, number of layers, fiber orientation, thickness ratio, orthotropic ratio) upon the frequencies and critical loads are discussed as well. Numerical results are compared with 3D exact solution when available from the open literature.     

Transient waves in plates of functionally graded materials

A numerical method is proposed for analysing transient waves in plates of functionally graded material (FGM) excited by impact loads. The material properties of the FGM plate have a gradient in the thickness direction and are anisotropic in the plane of the plate. In the present method, the FGM plate is divided into layer elements in the thickness direction. For an accurate modelling of the variation of the material property of FGM plates, it is expressed by second‐order polynomials in the thickness direction within an element. This can further reduce the number of elements to obtain more accurate results effectively. The principle of virtual work is used to develop approximate dynamic equilibrium equations. The displacement response is determined by employing the Fourier transformation and the modal analysis. As examples, the displacement response of FGM plates excited by line, point and distributed loads is calculated. The computations have shown the efficiency of the present method. Copyright © 2001 John Wiley & Sons, Ltd.     

机电耦合载荷下的压电层合板瞬态响应分析

针对压电层合板在机电耦合激振下的瞬态响应问题, 提出一种高效混合数值计算方法。经过位移场、 电势场在厚度方向的离散, 利用机电耦合理论和哈密顿原理, 推导出结构的运动方程。引入傅里叶变换, 得到波数域内运动控制方程。应用模态分析方法求解波数域内的位移场和电势场, 对结果进行傅里叶逆变换, 得到空间域内的瞬态响应。以PZT-5A/0° PVDF铺层两相材料复合压电层合板为算例, 分析了力、 电耦合线载荷激励下, 位移场和电势场的瞬态响应历程与分布规律, 计算结果给出了该结构的动力学基本特征。该方法结合了有限元法、 傅里叶变换和模态分析法, 计算高频载荷激振下的压电层合板瞬态响应较一般有限元法大幅减少了单元的划分。该方法可推广至分析任意机电载荷下的各类铺层材料压电层合板瞬态响应问题。      

水下冲击载荷用下复合材料圆柱壳的弹性响应

研究了平面阶跃冲击波作用下沉浸无限长复合材料圆柱壳的瞬态响应。首先推导了基于Flugge薄壳理论的复合材料圆柱壳的控制运动方程,然后采用反射尾流虚源法建立流体-结构的相互作用,最后采用有限差分法来求解平面冲击波作用下无限长复合材料圆柱壳的控制运动微分方程。详细考察了复合材料纤维方向和壳的半径对复合材料圆柱壳的无因次中面应变、第0阶模态径向位移和第1阶模态径向速度的影响。     

圆锥头弹体冲击复合材料层合板数值模拟

采用ABAQUS软件建立了圆锥头弹体正冲击复合材料层合板的有限元模型,并与已有文献结果进行对比验证了模型的可靠性,进而研究圆锥头弹体以不同的入射角度冲击复合材料层合板时初始速度与剩余速度的关系、复合材料层合板的破坏形态及弹体发生跳弹的规律。结果表明：弹体以90°入射角冲击复合材料层合板,在距离临界速度较大时,弹体的剩余速度与初始速度呈线性关系;不同的初始速度对复合材料层合板的损伤面积和破坏机制也不相同;弹体的入射角度越小、复合材料层合板越厚,越容易产生跳弹现象,并给出了入射角度和铺陈层数对跳弹现象的影响规律。该研究可为各种防护装备的设计和优化提供参考。     

复合材料层合板的声辐射模态幅值分析

低频时, 控制振动结构第一阶声辐射模态伴随系数可有效控制总声功率。通过分层有限元模型可以求解层合板的位移模式。对层合板的固有频率和动态响应进行了理论推导。结合声辐射模态理论, 研究了层合板铺设角度、 弹性模量比、 跨厚比以及阻尼比等结构参数对层合板结构第一阶声辐射模态伴随系数的影响。计算结果表明, 分层理论结合有限元方法可以较准确地计算层合板固有频率, 而且铺设角度和跨厚比对层合板结构声辐射模态影响较大。     

The multiblock method. A new strategy based on domain decomposition for the solution of wave propagation problems

A novel single‐step domain decomposition technique for the elastic wave propagation problem is introduced in this paper, based on the Huygens principle. The method allows an effective and efficient implementation on parallel computers through coarse‐grain multiprocessor computations. The various tests and numerical examples let infer that it is very competitive in comparison with classic substructuring techniques, specially for implicit discretization schemes. Copyright © 2000 John Wiley & Sons, Ltd.     

A closed form solution for linear and nonlinear free vibrations of composite and fiber metal laminated rectangular plates

In this paper, the nonlinear equations of motion for laminated composite rectangular plates based on first order shear deformation theory, which include shear deformation and rotary inertia, have been derived. Then, through introducing a force function, these equations reduced to a set of coupled nonlinear partial differential equations and a compatibility equation. By using the Galerkin method, for the first time, a nonlinear ordinary differential equation is obtained, which includes nonlinear inertia and stiffness terms. By using the multiple time scales method, analytical relations for nonlinear frequency and transverse displacement have been obtained. Results are compared with the literature and good agreement is achieved for both linear and nonlinear frequencies. After proving the validity of our work for isotropic rectangular plates and laminated rectangular plates, linear and nonlinear free vibration of a Fiber Metal Laminate panel have been investigated. Also the effects of some system parameters on the nonlinear frequency have been investigated.     

Dynamic analysis of laminated composite plates subjected to a moving oscillator by FEM

This paper presents a finite element model based on the first order shear deformation theory to investigate the dynamic behavior of laminated composite plates traversed by a moving oscillator. The oscillator model is assumed to be consisting of two nodal masses that are connected by means of a spring-damper unit. The governing equations of motion of two sub-systems are separately integrated by applying the Newmark’s time integration procedure. Then, the obtained equations are coupled and the responses of system components are calculated in each time step. The accuracy of algorithm is verified by comparing the numerical results of static, free vibration and simplified moving force problems analysis with the available exact solutions and other numerical results in the literature. Also, the effects of mass ratio, damping ratio of system components, stiffness of suspension system, velocity and eccentricity of moving oscillator on dynamic responses is parametrically studied. This algorithm can be applied to various boundary conditions, lamination schemes and fiber angels.     

A numerical approach to the analysis of failure modes in anisotropic plates

This study is motivated by the attempt to characterize failure modes of silicon chips commonly used in electronic industries. Previous experimental investigations provided the failure probability of dies made of a single-crystal and produced a large variety of crack patterns, but were not able to elucidate the link between defect distributions and crack initiation and propagation. To get some insight in the fracture activation and propagation mechanisms, we resort to finite element analyses and adopt an explicit methodology for crack tracking, based on the self-adaptive insertion of cohesive elements into a coherent mesh of solid elements. Finite kinematics material models with anisotropic features for both bulk and cohesive surfaces are employed to describe the behavior of single-crystal silicon plates undergoing a particular bending test up to failure. The cohesive model adopted in the calculation is fully anisotropic and newly formulated to accomplish the present study. Numerical simulations considering different material properties were able to ascertain the effects of particular flaws on failure modes of brittle silicon plates.     

Analysis of the integration of cohesive elements in regard to utilization of coarse mesh in laminated composite materials

An analysis of the solver convergence difficulties and erroneous results when large cohesive elements are utilized in delamination propagation simulations in laminated composites is presented. Special focus is given to the numerical integration of the cohesive element force vector and stiffness matrix. The magnitude and variation of the integration error are analyzed, and the results show that contrary to statements found elsewhere in the literature, the 2 × 2 point Newton‐Cotes quadrature, commonly used in commercial software, results in large errors and is one of the limiting factors for using larger elements. By reducing the integration error in the damage process zone of the interface, more accurate results can be obtained, and larger elements can be utilized with less iterations, thereby decreasing the computational cost. Copyright © 2014 John Wiley & Sons, Ltd.     

Transient Response in Cross-Ply Laminated Cylinders and Its Application to Reconstruction of Elastic Constants

An efficient hybrid numerical method is presented for investigating transient response of cross-ply laminated axisymmetric cylinders subjected to an impact load. In this hybrid numerical method, the laminated cylinder is divided into layered cylindrical elements in the thickness direction. The Hamilton principle is used to develop governing equations of the structure. The displacement response is determined by employing the Fourier transformations and the modal analysis. Numerical examples for analyzing transient waves have been provided in axisymmetric laminated cylindrical structures, both for thin cylindrical shells and thick cylinders.
A computational inverse technique is also presented for reconstructing elastic constants of axisymmetric crossply laminated cylinders from the surface displacement response data, using the present hybrid numerical method as the forward solver and neural network as the inverse operator. This technique is utilized to reconstruct the elastic constants of an axisymmetric laminated cylindrical shell.     

An enriched macro finite element for the static analysis of thick general quadrilateral laminated composite plates

This article presents the formulation of an enriched macro finite element based on the trigonometric shear deformation theory for the static analysis of symmetrically laminated composite plates. Shear correction factor is not required because this theory accounts for tangential stress-free boundary conditions on the plate boundary surfaces. The macro element is obtained using the principle of virtual work and Gram-Schmidt orthogonal polynomials as enrichment functions. The implementation of the obtained algorithm is simple and efficient, and allows studying general quadrilateral plates with a single macro element. Several examples are presented to show the capability and applicability of the developed formulation.     

A meshless local radial point collocation method for free vibration analysis of laminated composite plates

In this paper, a meshless local radial point collocation method based on multiquadric radial basis function is proposed to analyze the free vibration of laminated composite plates. This method approximates the governing equations based on first-order shear deformation theory using the nodes in the support domain of any data center. Natural frequencies of the laminated composite plates with various boundary conditions, side-to-thickness ratios, and material properties are computed by present method. The choice of shape parameter, effect of dimensionless sizes of the support domain on accuracy, convergence characteristics are studied by several numerical examples. The results are compared with available published results which demonstrate the accuracy and efficiency of present method.     

A mixed‐enhanced finite‐element for the analysis of laminated composite plates

This paper presents a new 4‐node finite‐element for the analysis of laminated composite plates. The element is based on a first‐order shear deformation theory and is obtained through a mixed‐enhanced approach. In fact, the adopted variational formulation includes as variables the transverse shear as well as enhanced incompatible modes introduced to improve the in‐plane deformation. The problem is then discretized using bubble functions for the rotational degrees of freedom and functions linking the transverse displacement to the rotations. The proposed element is locking free, it does not have zero energy modes and provides accurate in‐plane/out‐of‐plane deformations. Furthermore, a procedure for the computation of the through‐the‐thickness shear stresses is discussed, together with an iterative algorithm for the evaluation of the shear correction factors. Several applications are investigated to assess the features and the performances of the proposed element. Results are compared with analytical solutions and with other finite‐element solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Load increment procedure for post‐buckling analysis of laminated plates under in‐plane loads by finite strip method

A load increment procedure has been presented to integrate with the finite strip method for the post‐buckling analysis of laminated plates when subjected to uniform end shortening. In‐plane loads are introduced to reflect the end shortening effect. The Newton–Raphson procedure is implemented to attend a solution that satisfies the equilibrium condition and at the same time meets the loading requirements. Error associated with loading condition is minimised by adjusting the load factor to preserve the rate of convergence. The enhanced capability can be easily incorporated into the context of both classical and shear deformation plate theories. A range of application has been described. Convergence test and numerical results are presented for isotropic plate and laminates with general lay‐up arrangement. Copyright © 2000 John Wiley & Sons, Ltd.     

基于傅里叶变换的正弦信号源波形失真评价方法

介绍了一种直接利用经典傅里叶变换技术评价正弦信号源波形总失真度的过程和方法。借助于周期的精确测量技术,实现基波参数的精确测量。用时域能量方式计算失真,其内涵全面,包括了谐波、杂波、噪声等全部影响,并对测量系统本身的影响进行了补偿。同时,对不同条件下的失真度测量误差进行了分析,用一组仿真数据的实验结果,验证了正确性及可行性。所述方法可用于正弦信号波形失真度的精确测量和计量校准,尤其适用于大失真度的测量。     

A co-rotational 8-node assumed strain shell element for postbuckling analysis of laminated composite plates and shells

 The formulation of a nonlinear composite shell element is presented for the solution of stability problems of composite plates and shells. The formulation of the geometrical stiffness presented here is exactly defined on the midsurface and is efficient for analyzing stability problems of thin and thick laminated plates and shells by incorporating bending moment and transverse shear resultant forces. The composite element is free of both membrane and shear locking behaviour by using the assumed natural strain method such that the element performs very well as thin shells. The transverse shear stiffness is defined by an equilibrium approach instead of using the shear correction factor. The proposed formulation is computationally efficient and the test results showed good agreement. In addition the effect of the viscoelastic material is investigated on the postbuckling behaviour of laminated composite shells. Received: 6 February 2002 / Accecpted: 6 January 2003 ID=" Present address: School of Civil Engineering, Asian Institute of Technology     

Numerical plate testing for linear two‐scale analyses of composite plates with in‐plane periodicity

A method of numerical plate testing (NPT) for composite plates with in‐plane periodic heterogeneity is proposed. In the two‐scale boundary value problem, a thick plate model is employed at macroscale, while three‐dimensional solids are assumed at microscale. The NPT, which is nothing more or less than the homogenization analysis, is in fact a series of microscopic analyses on a unit cell that evaluates the macroscopic plate stiffnesses. The specific functional forms of microscopic displacements are originally presented so that the relationship between the macroscopic resultant stresses/moments and strains/curvatures to be consistent with the microscopic equilibrated state. In order to perform NPT by using general‐purpose FEM programs, we introduce control nodes to facilitate the multiple‐point constraints for in‐plane periodicity. Numerical examples are presented to verify that the proposed method of NPT reproduces the plate stiffnesses in classical plate and laminate theories. We also perform a series of homogenization, macroscopic, and localization analyses for an in‐plane heterogeneous composite plate to demonstrate the performance of the proposed method. Copyright © 2015 John Wiley & Sons, Ltd.