Reliability assessment of different plate theories for elastic wave propagation analysis in functionally graded plates

The importance of elastic wave propagation problem in plates arises from the application of ultrasonic elastic waves in non-destructive evaluation of plate-like structures. However, precise study and analysis of acoustic guided waves especially in non-homogeneous waveguides such as functionally graded plates are so complicated that exact elastodynamic methods are rarely employed in practical applications. Thus, the simple approximate plate theories have attracted much interest for the calculation of wave fields in FGM plates. Therefore, in the current research, the classical plate theory (CPT), first-order shear deformation theory (FSDT) and third-order shear deformation theory (TSDT) are used to obtain the transient responses of flexural waves in FGM plates subjected to transverse impulsive loadings. Moreover, comparing the results with those based on a well recognized hybrid numerical method (HNM), we examine the accuracy of the plate theories for several plates of various thicknesses under excitations of different frequencies. The material properties of the plate are assumed to vary across the plate thickness according to a simple power-law distribution in terms of volume fractions of constituents. In all analyses, spatial Fourier transform together with modal analysis are applied to compute displacement responses of the plates. A comparison of the results demonstrates the reliability ranges of the approximate plate theories for elastic wave propagation analysis in FGM plates. Furthermore, based on various examples, it is shown that whenever the plate theories are used within the appropriate ranges of plate thickness and frequency content, solution process in wave number-time domain based on modal analysis approach is not only sufficient but also efficient for finding the transient waveforms in FGM plates.     

各向异性立方晶体的弹性格林函数及其应用

本文计算了各向异性立方晶体的弹性格林函数的级数展开式,给出了直到二级近似的展开式的系数。将所得结果应用于弹性偶极子模型,给出了对称中心所产生的位移场及两对称中心间的弹性相互作用的表示式。应用于Cu,K等强各向异性立方晶体,虽然级数的收敛较慢,但所得关于对称中心的位移场,及二对称中心间的互作用能的数值结果,竟与基于点阵的不连续性作出的点阵静力学计算所得的结果基本一致。从而表明,本文给出的直到二级近似的弹性格林函数的解析表示提供了一个可以普遍应用的简便的方法。它可以较准确地描述立方晶体的某些力学行为。 关键词：     

Free vibration of antisymmetric,angle-ply laminated plates including transverse shear deformation by the finite element method

A finite element formulation of the equations governing the laminated anisotropic plate theory of Yang, Norris and Stavsky, is presented. The theory is a generalization of Mindlin's theory for isotropic plates to laminated anisotropic plates and includes shear deformation and rotary inertia effects. Finite element solutions are presented for rectangular plates of antisymmetric angle-ply laminates whose material properties are typical of a highly anisotropic composite material. Two sets of material properties that are typical of high modulus fiber-reinforced composites are used to show the parametric effects of plate aspect ratio, length-to-thickness ratio, number of layers and lamination angle. The numerical results are compared with the closed form results of Bert and Chen. As a special case, numerical results are presented for thick isotropic plates, and are compared with those for 3-D linear elasticity theory and Mindlin's thick plate theory.     

二相介质中一种位错构型的弹性场及所受的力

本文建立了二相介质中与相界斜交位错的一种理想构型。根据位错各向异性弹性理论的普遍方法,以及关于界面的Gebbia-Eshelby定理,计算得到此位错在两相中的弹性场,以及位错线段所受“像力”的数值解。本文所提出的模型及算法具有一般适用性。结果可用于进一步计算此位错与其它缺陷的相互作用。对于考虑这种斜交位错在无扩散相变中的作用,本工作也是一个有意义的起点。 关键词：     

不同铺设角度下层合板结构参数对声功率影响*

本文研究了不同纤维铺设角度的层合板结构参数对声功率的影响,从而为层合板的低噪声设计提供理论依据。通过分层有限元理论获得层合板结构动力学响应,基于声辐射模态概念分析不同铺设角度下层合板不同铺设方式、宽厚比和弹性模量比对其声功率影响。结果表明,层合板的铺设角度和宽厚比对复合材料层合板结构的声辐射功率影响较大。首先相同铺设角度的层合板,改变弹性模量比,声功率变化不明显;其次改变不同的铺设角度,宽厚比较小的层合板声功率下降的空间更大,更易于声功率的降低。层合板作为结构件时,从降低声功率角度而言总体上对称铺设结构比单向铺设层合板结构有优势,并且相同铺设角度下,反对称铺设层合板可获得更小的辐射声功率。     

ON GROUP VELOCITY OF ELASTIC WAVES, IN AN ANISOTROPIC PLATE

A proof is offered for the equivalence of the group velocity and the energy velocity for elastic waves in a free anisotropic homogenous plate. The proof is valid in the case where the group velocity may be negative while the wave number is positive.     

Elastic wave field computation in multilayered nonplanar solid structures: a mesh-free semianalytical approach

Multilayered solid structures made of isotropic, transversely isotropic, or general anisotropic materials are frequently used in aerospace, mechanical, and civil structures. Ultrasonic fields developed in such structures by finite size transducers simulating actual experiments in laboratories or in the field have not been rigorously studied. Several attempts to compute the ultrasonic field inside solid media have been made based on approximate paraxial methods like the classical ray tracing and multi-Gaussian beam models. These approximate methods have several limitations. A new semianalytical method is adopted in this article to model elastic wave field in multilayered solid structures with planar or nonplanar interfaces generated by finite size transducers. A general formulation good for both isotropic and anisotropic solids is presented in this article. A variety of conditions have been incorporated in the formulation including irregularities at the interfaces. The method presented here requires frequency domain displacement and stress Green's functions. Due to the presence of different materials in the problem geometry various elastodynamic Green's functions for different materials are used in the formulation. Expressions of displacement and stress Green's functions for isotropic and anisotropic solids as well as for the fluid media are presented. Computed results are verified by checking the stress and displacement continuity conditions across the interface of two different solids of a bimetal plate and investigating if the results for a corrugated plate with very small corrugation match with the flat plate results.     

EXPERIMENTAL AND THEORETICAL INVESTIGATION OF THE LINEAR AND NON-LINEAR DYNAMIC BEHAVIOUR OF A GLARE 3 HYBRID COMPOSITE PANEL

A theoretical model based on Hamilton's principle and spectral analysis is used to study the non-linear free vibration of hybrid composite plates made of Glare 3, a new aircraft structural material. It consists of alternating layers of metal- and fibre-reinforced composites. In previous work, the theoretical model has been used to calculate the first non-linear mode of fully clamped rectangular composite fibre-reinforced plastic (CFRP) laminated plates. This study concerns determination of the linear dynamic properties of the Glare 3 hybrid composite rectangular panel (G3HCRP) such as natural frequencies and mode shapes. The theoretical model is used to calculate the fundamental non-linear mode shape and associated flexural behaviour of the fully clamped G3HCRP. A series of experimental investigations have been conducted using a G3HCRP in order to determine linear dynamic properties. The response due to random excitation was investigated and the experimental measurements are analyzed and discussed. Comparisons are made with finite element predictions and response estimates given by the ESDU method, the latter being a “design guide” approach used by industry. Concerning the non-linear analysis, the results are given for various plate aspect ratios and vibration amplitudes, showing a higher increase of the induced bending stress near the clamps at large deflections. Comparisons between the dynamic behaviour of an isotropic plate and G3HCRP at large vibration amplitudes are presented and good results are obtained.     

Normal mode expansion method for laser-generated ultrasonic Lamb waves in orthotropic thin plates

A quantitative theory for modeling the laser-generated transient ultrasonic Lamb waves, which propagates along arbitrary directions in orthotropic plates, is presented by employing an expansion method of generalized Lamb wave modes. The displacement field is expressed by a summation of the symmetric and antisymmetric modes in the surface stress-free orthotropic plate, and therefore the theory is particularly appropriate for waveform analyses of Lamb waves in thin plates because one needs only to evaluate several lower modes. The transient waveforms excited by the thermoelastic expansion and the oil-coating evaporation are analyzed for a transversely isotropic thin plate. The results show that the theory provides a quantitative analysis to characterize anisotropic elastic stiffness properties of orthotropic plates by laser-generated Lamb wave detection.     

Elasticity theory of a twisted stack of plates

We present an elastic model of B-form DNA as a stack of thin, rigid plates or base pairs that are not permitted to deform. The symmetry of DNA and the constraint of plate rigidity limit the number of bulk elastic constants contributing to a macroscopic elasticity theory of DNA to four. We derive an effective twist-stretch energy in terms of the macroscopic stretch along and relative excess twist about the DNA molecular axis. In addition to the bulk stretch and twist moduli found previously, we obtain a twist-stretch modulus with the following remarkable properties: 1) it vanishes when the radius of the helical curve following the geometric center of each plate is zero, 2) it vanishes with the elastic constant K₂₃ that couples compression normal to the plates to a shear strain, if the plates are perpendicular to the molecular axis, and 3) it is nonzero if the plates are tilted relative to the molecular axis. This implies that a laminated helical structure carved out of an isotropic elastic medium will not twist in response to a stretching force, but an isotropic material will twist if it is bent into the shape of a helix. Received: 4 July 1997 / Received in final form: 16 October 1997 / Accepted: 21 October 1997     

Detection and localization of inclusions in plates using inversion of point actuated surface displacements

A numerical simulation is carried out demonstrating the use of plate surface vibration measurements for detecting and locating inclusions within the structure. A finite element code is used to calculate normal surface displacement for both steel and mortar plates subjected to a monochromatic point force. The data is generated for the homogeneous plate and the identical plate within which exists a small rectangular inclusion. It is observed that when the elastic modulus of the inclusion is orders of magnitude lower than the base material, resonances of the inclusion produce large local displacements that are readily observed in the raw displacement data. For more modest moduli differences, there are no such directly observable effects. In this case, three inverse algorithms are used to process the displacement data. The first two are local inversion techniques that each yield a spatial map of the elastic modulus normalized by density. These algorithms successfully detect and localize the inclusion based on its modulus difference from that of the base plate. The third technique uses a form of the inhomogeneous equation of motion to obtain the induced force distribution connected with the inclusion. The spatial mapping of this force also successfully detects and localizes the inclusion.     

On the non-linear vibrations of rectangular plates

A solution, based on a one-term mode shape, for the large amplitude vibrations of a rectangular orthotropic plate, simply supported on all edges or clamped on all edges for movable and immovable in-plane conditions, is found by using an averaging technique that helps to satisfy the in-plane boundary conditions. This averaging technique for satisfying the immovable in-plane conditions can be used to resolve many anisotropic and skew plate problems where otherwise, when a stress function is used, the integration of the u and v equations becomes difficult, if not impossible. The results obtained herein are compared with those available in the literature for the isotropic case and excellent agreement is found. Results available for the one-term mode shape solutions of these problems are compared and the non-linear effect is presented as functions of aspect ratio and of the orthotropic elastic constants of the plate. The results are further compared with those based on the Berger method and the detailed comparative studies show that the use of the Berger approximation for large deflection static and dynamic problems and its extension to anisotropic plates, skew plates, etc., can lead to quite inaccurate results.     

ESTIMATION OF 3D FIELD QUANTITIES AND ENERGY FLUX ASSOCIATED WITH ELASTIC WAVES IN A BAR

A method to estimate dispersion relations and warping associated with elastic wave propagation in a bar is presented. The method is based on Hamilton's principle. It is shown how the theoretical model together with strain measurements can be used to evaluate three dimensional (3D) field quantities like displacements and stresses at an arbitrary position in the bar, as well as energy flux through an arbitrary cross-section of the bar. It is also shown how redundant measurements can be used to increase the accuracy. The method is general and can be applied to any mode of wave propagation, isotropic or anisotropic linearly elastic material, and any cross-sectional geometry. Here, it is applied to longitudinal waves in a split Hopkinson pressure bar with linear elastic isotropic material behaviour and square cross-section. In particular, axial displacement, axial stress and energy flux are evaluated at a free end of the bar in order to test the method. The method is also used to estimate the Poisson ratio of the bar material, by measuring axial and transverse strains at the same axial position.     

Non-linear flexural waves in thin layers

Non-linear flexural waves in thin plates or layers have been analyzed in this paper. The equation of motion of the plate is derived assuming that the motion is antisymmetric about the mid-plane of the plate and that the plate is thin. The plate is considered to be elastic. The Von Karman non-linear strains and Landau elastic constants have been used to model geometric and material non-linearities, respectively. An asymptotic analysis of wave motion is presented using the method of multiple scales. Evolution equations are derived for small amplitude traveling flexural elastic waves. Numerical results show waveform distortion, amplitude amplification, and harmonic generation.     

Sound radiation from shear deformable stiffened laminated plates

Sound radiation from shear deformable stiffened laminated plates is studied theoretically. The equations of motion for the composite laminated plate are derived on the basis of the first-order shear deformation plate theory. Two sets of parallel stiffeners interact with the laminated plate only through the normal line forces. By using the Fourier wavenumber transform and the stationary phase method, the far-field sound pressure is described analytically. Sound pressure given by the first-order shear deformation plate theory and the classical thin plate theory is compared, and the differences of sound pressure are shown in the high frequency range for an isotropic plate. Sound pressure and the transverse displacement spectra are presented to illustrate the effects of force location, stiffeners and angle-ply layers. Sound radiation from symmetric and antisymmetric composite plates with multiple loadings is also investigated.     

Dynamic analysis of circular and sector thick,layered plates

The finite element method is extended to the free vibration analysis of laminated thick plates with curved boundaries. Two elements are developed on the basis of Mindlin's thick plate theory in which the effects of thickness-shear deformation and rotary inertia are included. Both elements are derived in polar co-ordinates and can be joined together to handle annular as well as circular laminated anisotropic plate problems. Since axisymmetry has not been assumed, variations in material properties in the tangential direction can be dealt with. Numerical results are presented to demonstrate the influence of geometrical shape as well as that of thickness-shear deformation on the free vibrations of both homogeneous and layered plates. Comparisons between the numerical results obtained and those presented by other investigators confirm the accuracy of the new elements. The elements also can be used in the analysis of rectangular plates by assuming very large radii and very small subtended angle values.     

Approximate closed form solutions for free vibration of polar orthotropic circular plates

For free vibrations of polar orthotropic plate, simple approximate closed form solutions for mode shapes and its natural frequencies were obtained using the Rayleigh-Ritz method. Coordinate function satisfying the natural boundary conditions and the predetermined coefficients was adapted, which results in compact expressions and enables to readily calculate symmetric and nonsymmetric natural frequencies for arbitrary values of the elastic constants. The derived formulation can be used in designing of circular plates such as wood disk, which are naturally endowed with material orthotropy as well as fiber reinforced composite materials. The model can easily be used for the evaluation of parametric studies on dynamic behaviors and nondestructive methods during the initial design process.     

水下多层均匀材料的声特性

本文研究水和空气中由不同厚度的两层弹性板及消声橡胶共同组成的分层介质在斜入射时的声特性,并采用数值方法求解多层结构声场的线性方程组。得到反射,透射系数,当板的厚度远小于波长时,采用弹性薄壳振动模型和弹性波模型进行计算,它们的结果基本相同,但是薄壳振动模型更便于计算和理论分析;对系统反射系数和透射系数的计算结果表明,系统增加消声橡胶层后,极大地改变后射系数和透射系数的频率响应;频率升高,反射系数减小     

Complete mode identification for resonance ultrasound spectroscopy

This study is devoted to deducing exact elastic constants of an anisotropic solid material without using any advance information on the elastic constants by incorporating a displacement-distribution measurement into resonant ultrasound spectroscopy (RUS). The usual RUS method measures free-vibration resonance frequencies of a solid and compares them with calculations to find the most suitable set of elastic constants by an inverse calculation. This comparison requires mode identification for the measured resonance frequencies, which has been difficult and never been free from ambiguity. This study then adopts a laser-Doppler interferometer to measure the displacement-distribution patterns on a surface of the vibrating specimen mounted on pinducers; comparison of the measured displacement distributions with those computed permits us to correctly identify the measured resonance frequencies, leading to unmistakable determination of elastic constants. Because the displacement patterns are hardly affected by the elastic constants, an exact answer is surely obtained even when unreasonable elastic constants are used as initial guesses at the beginning of the inverse calculation. The usefulness of the present technique is demonstrated with an aluminum alloy and a langasite crystal.     

NON-DESTRUCTIVE EVALUATION OF HORIZONTAL CRACK DETECTION IN BEAMS USING TRANSVERSE IMPACT

Numerical and experimental studies for crack detection in beam employing transverse impact are presented. In the numerical study, a beam model of wave propagation is developed to calculate the time history of beam response before, over and after the crack region. It is expected that the resulting wave in the beam will be scattered by the crack and will carry information on the location and geometry of the crack. Experiments using a scanning laser vibrometer on specimens containing simulated crack are then conducted to verify the numerical results. Comparison study between the numerical results and experimental observations are conducted; good correlation between theory and experiment is observed. The beam model of wave propagation and adaptive multilayer perceptron networks (MLP) are then used for inverse identification of crack parameters (i.e., crack location, depth and length) in the beams. Time-domain displacement responses calculated using the present beam model containing predetermined crack parameters are used as training data for the MLP. Once the MLP is trained, the MLP networks are then employed for inverse determination of an unknown crack in a beam using experimental displacement responses measured with a scanning laser vibrometer. Examples show that the procedure performs well for the determination of a wide range of values for the location, depth and length of the crack.