Elastodynamic behavior of laminated orthotropic plates under initial stress

A finite element method of analysis of the vibrational and wave propagational characteristics is presented for a laminated orthotropic plate under initial stress. The plate may have an arbitrary number of bonded elastic orthotropic layers, each with distinct thickness, density and mechanical properties, and the analysis is capable of treating a completely arbitrary three-dimensional state of initial stress. Biot's theory for incremental elastic deformations of a stressed solid forms the basis for this study. A homogeneous, isotropic plate under two different states of initial stress was analyzed and their numerical results showed excellent correlation with those from an exact solution. Further examples of a three layer composite plate and a sandwich plate are offered to add some general insight to the physical behavior of such plates.     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

The vibroacoustic response and sound absorption performance of multilayer,microperforated rib-stiffened plates

The vibroacoustic response and sound absorption performance of a structure composed of multilayer plates and one rigid back wall are theoretically analyzed. In this structure, all plates are two-dimensional, microperforated, and periodically rib-stiffened. To investigate such a structural system, semianalytical models of one-layer and multilayer plate structures considering the vibration effects are first developed. Then approaches of the space harmonic method and Fourier transforms are applied to a one-layer plate, and finally the cascade connection method is utilized for a multilayer plate structure. Based on fundamental acoustic formulas,the vibroacoustic responses of microperforated stiffened plates are expressed as functions of a series of harmonic amplitudes of plate displacement, which are then solved by employing the numerical truncation method. Applying the inverse Fourier transform, wave propagation, and linear addition properties, the equations of the sound pressures and absorption coefficients for the one-layer and multilayer stiffened plates in physical space are finally derived. Using numerical examples, the effects of the most important physical parameters—for example, the perforation ratio of the plate, sound incident angles, and periodical rib spacing—on sound absorption performance are examined. Numerical results indicate that the sound absorption performance of the studied structure is effectively enhanced by the flexural vibration of the plate in water.Finally,the proposed approaches are validated by comparing the results of stiffened plates of the present work with solutions from previous studies.     

复合材料中矩形夹杂角端部力学行为分析

提出了一种分析矩形夹杂角端部奇异应力场的新型杂交有限元方法,该方法在分析矩形夹杂角端部奇异应力场时,需要在夹杂端部构造一个超级单元。超级单元的刚度矩阵可以通过夹杂端部特征问题数值解建立。我们用这种方法计算了单向载荷作用下无限大均质板中单个矩形夹杂角端部奇异应力场,并与现有的数值解进行了比较。比较结果表明:本文提出的方法是可行的、有效的,而且数值结果精度高。为说明本文方法适用范围更广,文章最后讨论了各向异性弹性材料和横观各向同性压电材料中矩形夹杂角端部电弹性场行为。     

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

Dynamic stress analysis of a functionally graded material plate with a circular hole

This paper is to study the two-dimensional dynamic stress of a functionally graded material (FGM) plate with a circular hole under plane compressional waves at infinity. With using the method of piece-wise homogeneous layers, the dynamic stress distribution of the FGM plate having radial arbitrary material parameters is derived based on the complex variable method. As examples, numerical results are presented for the FGM plate having given radial shear modulus, density and Poisson’s ratio. It is found that the dynamic stress around the circular hole in the FGM plate can be effectively reduced by choosing the proper change ways of the radial material parameters for different frequency incident wave.     

Interaction between cracks and a circular inclusion in a finite plate with the distributed dislocation method

This paper presents a numerical solution of interaction between cracks and a circular inclusion in a finite plate. Both the boundaries and the cracks are modeled by distributed dislocations. This approach will result in a set of singular integral equations with Cauchy kernels, which can be solved by Gauss–Chebyshev quadratures. Several numerical examples are given to assess the performance of the presented method. The solutions obtained by this method have been checked and confirmed by the finite element analysis results.     

尖锐夹杂角端部局部应力场杂交有限元分析

本文首先利用作者曾提出的一维有限元特征分析方法计算所得到的尖锐夹杂角端部应力奇异指数和奇异应力场、位移场角分布函数,并依据Hellinger-Reissner原理,开发出了一个特殊的、能够反映夹杂角端部局部弹性现象的n结点多边形超级角端部单元,然后将该超级单元与标准的4结点杂交应力单元耦合在一起构建了一种分析异形夹杂角端部奇异弹性场的新型特殊杂交应力有限元方法.文中给出了两个应用算例,算例结果表明:本文方法不仅使用单元少、计算结果精度高,而且适用范围广,可拓展应用于分析复合材料微结构组织与力学行为关系.     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

Stress state of a piezoelectric elastic medium with an arbitrarily oriented triaxial ellipsoidal inclusion

We determine the electrostressed state of a piezoceramic medium with an arbitrarily oriented triaxial ellipsoidal inclusion under homogeneous mechanical and electric loads. Use is made of Eshelby’s equivalent inclusion method generalized to the case of a piezoelectric medium. Solving the problem for a spheroidal cavity with the axis of revolution aligned with the polarization axis demonstrates the high efficiency of the approach. A numerical analysis is carried out. The stress distribution along the surface of the arbitrarily oriented triaxial ellipsoidal inclusion is studied     

含椭圆刚性核有限大复合材料层板弯曲应力分析Bending problem of a finite composite laminated plate with elliptical rigid inclusions

基于经典的复合材料层板理论,将有限大复合材料层板等效成各向异性弹性平板。采用复变函数理论中的Faber级数分析方法,使用最小二乘边界配点法,对含多椭圆刚性核有限大各向异性板弯曲问题进行应力分析,得到了该问题的级数解形式,分析了含椭圆刚性核层板在弯曲载荷下的应力分布,并讨论了形状和结构参数对应力分布的影响。结果表明,本文方法对于分析含多个椭圆形刚性核有限大薄板弯曲应力问题非常有效,该方法具有精度高及计算方便等优点。     

Several three-dimensional solutions for transversely isotropic functionally graded material plate welded with circular inclusion

The problem of a transversely isotropic functionally graded material (FGM) plate welded with a circular inclusion is considered. The analysis starts with the generalized England-Spencer plate theory for transversely isotropic FGM plates, which expresses a three-dimensional (3D) general solution in terms of four analytic functions. Several analytical solutions are then obtained for an infinite FGM plate welded with a circular inclusion and subjected to the loads at infinity. Three different cases are considered, i.e., a rigid circular inclusion fixed in the space, a rigid circular inclusion rotating about the x-, y-, and z-axes, and an elastic circular inclusion with different material constants from the plate itself. The static responses of the plate and/or the inclusion are investigated through numerical examples.     

Identification of a moving boundary for a heat conduction problem in a multilayer medium

In this paper, we give a uniqueness theorem for the moving boundary of a heat problem in a composite medium. Through solving the Cauchy problem of heat equation in each subdomain, we finally find an approximation to the moving boundary for one-dimensional heat conduction problem in a multilayer medium. The numerical scheme is based on the use of the method of fundamental solutions and a discrete Tikhonov regularization technique with the generalized cross-validation choice rule for a regularization parameter. Numerical experiments for five examples show that our proposed method is effective and stable.     

Stress-strain state of an anisotropic plate with an elliptic hole and thin rigid inclusions

The stress-strain state of an anisotropic plate containing an elliptic hole and thin, absolutely rigid, curvilinear inclusions is studied. General integral representations of the solution of the problem are constructed that satisfy automatically the boundary conditions on the elliptic-hole contour and at infinity. The unknown density functions appearing in the potential representations of the solution are determined from the boundary conditions at the rigid inclusion contours. The problem is reduced to a system of singular integral equations which is solved by a numerical method. The effects of the material anisotropy, the degree of ellipticity of the elliptic hole, and the geometry of the rigid inclusions on the stress concentration in the plate are studied. The numerical results obtained are compared with existing analytical solutions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 4, pp. 173–180, July–August, 2007.     

Flexural vibrations of a thin circular elastic inclusion in an unbounded body under the action of a plane harmonic wave

The interaction of a plane harmonic longitudinal wave with a thin circular elastic inclusion is considered. The wave front is assumed to be parallel to the inclusion plane. Since the inclusion is thin, the matrix-inclusion interface conditions (perfect bonding) are formulated on the mid-plane of the inclusion. The bending displacements of the inclusion are determined from the bending equation for a thin plate. The problem is solved using discontinuous Lamé solutions for harmonic vibrations. Therefore, the problem can be reduced to the Fredholm equation of the second kind for a function associated with the discontinuity of normal stresses on the inclusion. The equation obtained is solved by the method of mechanical quadratures using Gaussian quadrature formulas. Approximate formulas for the stress intensity factors are derived. Results from a numerical analysis of the dependence of the SIFs on the dimensionless wave number and the stiffness of the inclusion are presented __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 5, pp. 16–21, May 2008.     

Strength calculations and optimal weight design of multilayer shell-shaped composite products under a set of loads

The stress-strain state of axisymmetric multilayer shells is analyzed using kinematic and static hypotheses that allow for the transverse shear stresses satisfying the necessary equations of state, continuity conditions at the boundaries between the layers and given boundary conditions. A numerical solution of the problem of the stress-strain state for a multilayer bar is compared with the Lekhnitskii solution (for a cantilever beam loaded by a concentrated force and moment) to asses the applicability of the employed bending equations of multilayer shells. It is shown that these solutions are in good agreement. The problem of the initial fracture of the shells considered is formulated using phenomenological strength criteria for each layer. A coordinate-wise descent method in the unit interval is proposed to solve weight optimization problems for multilayer shells of composite materials under combined loading. Regions of safe operating loads and the optimal weight distribution of layer thicknesses are determined for a multilayer bar acted upon by a uniformly distributed load and concentrated force.     

板结构计算模型的新发展

现代复合材料层合板具有高强和轻型的突出优点，从而在军工和民用等诸多领域发挥着重要作用。这种板结构的特点是随着纤维走向的不同，层间材料的物理-力学特性发生剧烈变化。沿板厚方向变形的梯度比较陡峭，并在层间结合面处发生强不连续，呈现zig-zag （锯齿状）现象。这导致横向剪应变在板的静态和动态响应中发生重要作用，不计横向变形的经典组合板计算模型CLPT难以适应现代多层板计算分析的需要。考虑横向剪切变形影响的板的计算模型得到重视和发展。需要指出，现有各种考虑剪切变形影响的计算模型虽然有了很大的发展，但在全面和准确性上仍然存在一定的不足，难以适应现代多层组合板横向力和物理性能多变的情况。模型预测的沿板厚方向位移和应力的变化规律难以通过严格的检验。本文提出的以比例边界有限元为基础的正交各向异性板的数值计算模型，同时可适用于各种薄板与厚板的分析，对现代复合材料层合板的分析具有特殊的优越性。所得到的板的位移、正应力和剪应力沿板厚方向的变化，与三维弹性理论的标准解高度吻合。数值算例进一步表明，随着层间纤维走向的变化，板内位移场和应力场沿板厚方向剧烈变化所呈现的锯齿现象均可以精准地进行模拟。据此，本文建议方法对现代板分析的广泛适应性和高度准确性得到了充分论证。     

An exact solution for a multilayered two-dimensional decagonal quasicrystal plate

By extending the pseudo-Stroh formalism to two-dimensional decagonal quasicrystals, an exact closed-form solution for a simply supported and multilayered two-dimensional decagonal quasicrystal plate is derived in this paper. Based on the different relations between the periodic direction and the coordinate system of the plate, three internal structure cases for the two-dimensional quasicrystal layer are considered. The propagator matrix method is also introduced in order to treat efficiently and accurately the multilayered cases. The obtained exact closed-form solution has a concise and elegant expression. Two homogeneous quasicrystal plates and a sandwich plate made of a two-dimensional quasicrystal and a crystal with two stacking sequences are investigated using the derived solution. Numerical results show that the differences of the periodic direction have strong influences on the stress and displacement components in the phonon and phason fields; different coupling constants between the phonon and phason fields will also cause differences in physical quantities; the stacking sequences of the multilayer plates can substantially influence all physical quantities. The exact closed-form solution should be of interest to the design of the two-dimensional quasicrystal homogeneous and laminated plates. The numerical results can also be employed to verify the accuracy of the solution by numerical methods, such as the finite element and difference methods, when analyzing laminated composites made of quasicrystals.     

Interaction of plane nonstationary waves with a thin elastic inclusion under smooth contact conditions

We solve the problem of determining the stress state near a thin elastic inclusion in the form of a strip of finite width in an unbounded elastic body (matrix) with plane nonstationary waves propagating through it and with the forces exerted by the ambient medium taken into account. We assume that the matrix is in the plane strain state, and the smooth contact conditions are realized on both sides of the inclusion. The method for solving this problem consists in using the integral Laplace transform with respect to time and in representing the stress and displacement images in terms of the discontinuous solution of Lamé equations in the case of plane strain. As a result, the initial problem is reduced to a system of singular integral equations for the transforms of the unknown stress and displacement jumps. To invert the Laplace transform, we use a numerical method based on replacing the Mellin integral by the Fourier series. As a result, we obtain approximate formulas for calculating the stress intensity factors (SIF) for the inclusion, which are used to study the SIF time-dependence and its influence on the values of the inclusion rigidity. We also studied the possibility of considering the inclusions of higher rigidity as absolutely rigid inclusions.     

Stress analysis of a functional graded material plate with a circular hole

This paper is to study the two-dimensional stress distribution of a functional graded material plate (FGMP) with a circular hole under arbitrary constant loads. With using the method of piece-wise homogeneous layers, the stress distribution of the functional graded material plate having radial arbitrary elastic properties is derived based on the theory of the complex variable functions. As examples, numerical results are presented for the FGMPs having given radial Young’s modulus or Poisson’s ratio. It is shown that the stress is greatly reduced as the radial Young’s modulus increased, but it is less influenced by the variation of the Poisson’s ratio. Moreover, it is also found that the stress level varies when the radial Young’s modulus increased in different ways. Thus, it can be concluded that the stress around the circular hole in the FGMP can be effectively reduced by choosing the proper change ways of the radial elastic properties.