纤维增强复合材料层合板屈曲性态分析的边界元法

本文用边界元法分析了纤维增强复合材料正交各向异性层合板的屈曲性态。为了克服在用边界元法求解正交各向异性层合板屈曲时寻求相应的基本解的困维,本文采用了双重傅立叶级数和引用等效荷载的概念,建立了层合板屈曲临界荷载的特征方程。算例说明了本文方法的可行性和有效性。      

复合材料边界元素法研究

边界元法是一种迅速发展的数值解法.本文研究了在复合材料层合板中运用平面二次边界元法的有关问题.首先,讨论了将层合板等效为一般正交各向异性单层板的方法;其次,给出了偏轴情况下的基本解,研究了直线型二次边界元偏轴影响系数的解析计算以及域内应力表达式;最后对算例进行了计算,所得结果与相应的解析解和有限元相符甚好.     

正交各向异性板动力响应的边界元方法

本文讨论了正交各向异性板动力响应的边界元方法,发展了一种新的求近似基本解的方法,由该方法得到的基本解,给出了弹性薄板诸问题近似基本解的统一形式。应用这个基本解,得到了正交各向异性板和弹性地基板稳态强迫振动的边界元方程。文中的算例表明,本文的方法具有相当高的精度。     

热荷载作用下正交各向异性层合板振动响应的精确解

本文从三维弹性力学出发,导出了层合板第 i 层在任意变温 T(x,y,z,t)下的状态方程,解此状态方程,可得正交各向异性层合板在热荷载作用下的振动响应。     

正交异性矩形薄板自由振动的一般解析解

本文建立了正交异性矩形薄板自由振动横向位移函数微分方程的一般解。可用来求解任意边界矩形板的振动问题。作为例题,计算了复合材料悬臂层合板的频率。其结果与试验值基本吻合。     

Winkler-Pasternak地基上四边固支层合板的振动

本文在文[1]的基础上,利用Bolotin方法,在考虑横向剪切变形的情况下,首次获得了Winkler-Pasternak地基上四边固支对称正交铺设层合板自由振动问题的渐近解析解.对于不同层数、不同长厚比的层合板以及在不同地基反力系数的情况下,文中进行了具体的数值计算,给出了相应的数值结果.作为本文一个简化算例,计算了弹性地基上正交各向异性厚板自由振动的固有频率,与其他文献所给结果非常一致.      

用样条有限点法分析正交各向异性层合板的振动与稳定问题

本文采用样条有限点法,对多种边界条件下正交各向异性层合板的振动与稳定问题,进行了系统的分析,编制了通用程度。计算结果表明:样条有限点法在分析计算正交各向异性层合板的振动与稳定问题时,具有自由度少、收敛快、精度高、计算格式简便,易于在微机上计算等优点。     

正交各向异性厚板振动分析的边界积分方程法

本文采用正交各向异性厚板静力问题的基本解作为边界积分方程的核函数,利用加权残数法建立了正交各向异性厚板振动分析的边界积分方程。文中详细地讨论了边界积分方程的数值处理过程并给出了若干数值算例以论证本文方法的正确性。      

复合材料层合板智能结构主动振动控制的边界元法

利用边界元法模拟智能结构的振动控制，推导出具有压电传感器及致动器的复合材料层合板的边界积分方程，应用负速度反馈控制律，研究了复合材料层合板智能结构主动振动控制问题，算例分析证明该方程的正确性。     

基于新修正偶应力理论的Mindlin层合板自由振动分析

该文基于各向异性修正偶应力理论建立一个Mindlin层合板（跨厚比10~20的中厚板）自由振动模型。该理论偶应力曲率张量不对称,但偶应力弯矩对称。利用Hamilton原理推导振动微分方程和边界条件。新模型可退化为修正偶应力层合薄板振动模型和经典Mindlin层合板振动模型。以正交铺设简支方板为例计算了偶应力模型的自振频率,分析偶应力Mindlin层合板的自由振动尺度效应。算例表明,该文建立的新修正偶应力层合板模型能够用于分析细观尺度下Mindlin层合板的自由振动及尺度效应。     

Harmonic analysis of shear deformable orthotropic cracked plates using the Boundary Element Method

In this work, the modal and harmonic analysis of orthotropic shear deformable cracked plates using a direct time-domain Boundary Element Method formulation based on the elastostatic fundamental solution of the problem is presented. The Radial Integration Method was used for the treatment of domain integrals involving distributed domain applied loads and those related with inertial mass forces. Numerical examples are presented to demonstrate the efficiency and accuracy of the proposed formulation.     

Elasto-plastic Analysis of Two-dimensional Orthotropic Bodies with the Boundary Element Method

The Boundary Element Method (BEM) is introduced to analyze the elasto-plastic problems of 2-D orthotropic bodies. With the help of known boundary integral equations and fundamental solutions, a numerical scheme for elasto-plastic analysis of 2-D orthotropic problems with the BEM is developed. The Hill orthotropic yield criterion is adopted in the plastic analysis. The initial stress method and tangent predictor-radial return algorithm are used to determine the stress state in solving the nonlinear equation with the incremental iteration method. Finally, numerical examples show that the BEM is effective and reliable in analyzing elasto-plastic problems of orthotropic bodies.     

A BEM-based meshless solution to buckling and vibration problems of orthotropicplates

In this paper a BEM-based meshless solution is presented to buckling and vibration problems of Kirchhoff orthotropic plates with arbitrary shape. The plate is subjected to compressive centrally applied load together with arbitrarily transverse distributed or concentrated loading, while its edges are restrained by the most general linear boundary conditions. The resulting buckling and vibration problems are described by partial differential equations in terms of the deflection. Both problems are solved employing the Analog Equation Method (AEM). According to this method the fourth-order partial differential equation describing the response of the orthotropic plate is converted to an equivalent linear problem for an isotropic plate subjected only to a fictitious load under the same boundary conditions. The AEM is applied to the problem at hand as a boundary-only method by approximating the fictitious load with a radial basis function series. Thus, the method retains all the advantages of the pure BEM using a known simple fundamental solution. Example problems are presented for orthotropic plates, subjected to compressive or vibratory loading, to illustrate the method and demonstrate its efficiency and its accuracy.     

Quasi-static poroelastic boundary element formulation based on the convolution quadrature method

Convolution Quadrature Method (CQM)-based Boundary Element formulations are up to now used only in dynamic formulations. The main difference to usual time-stepping BE formulations is the way to solve the convolution integral appearing in most time-dependent integral equations. In the CQM formulation, this convolution integral is approximated by a quadrature rule whose weights are determined by the Laplace transformed fundamental solutions and a linear multi-step method. In principle, for quasi-static poroelasticity there is no need to apply the CQM because time-dependent fundamental solutions are available. However, these fundamental solutions are highly complicated yielding very sensitive algorithms. On the contrary, the CQM based BE formulation proposed here is very robust and yields comparable results to other methodologies. This formulation is tested in 2-d in comparison with a Finite Element Method and analytical results.     

A general 3D BEM/FEM coupling applied to elastodynamic continua/frame structures interaction analysis

This paper presents a procedure for coupling general finite element models with three‐dimensional bodies modelled by the Boundary Element Method (BEM). Shells, plates and frames are modelled by the Finite Element Method (FEM) and coupled to the BEM domain directly or by means of rigid blocks. The coupling is used for the analysis of buildings connected to half‐space by means of rigid footings, piles or plates in bending and other problems where combinations of different types of sub‐domains are required, composite domains for instance. Several numerical examples are analysed to demonstrate the robustness and accuracy of the proposed scheme. Copyright © 1999 John Wiley & Sons, Ltd.     

A new method for the determination of the flexural rigidity of an orthotropic plate

In this paper a new method for the determination of flexural rigidities in orthotropic plate bending problems is presented. Boundary integral equations are established for the curvatures and the deflections inside the domain. By a simple discretization of the boundary and the inside plate, the elimination of curvatures is possible. If the fundamental solution of isotropic plates is chosen, then a linear system of n equations with three unknowns is obtained. These equations are provided by the knowledge of the deflections inside the plates, and the unknowns are the flexural rigidities. By using the least square method, the computation of these rigidities becomes easy.     

Dynamic analysis of Timoshenko beams by the boundary element method

A Boundary Element Method formulation is developed for the dynamic analysis of Timoshenko beams. Based on the use of not time dependent fundamental solutions a formulation of the type called as Domain Boundary Element Method arises. Beside the typical domain integrals containing the second order time derivatives of the transverse displacement and of the rotation of the cross-section due to bending, additional domain integrals appear: one due to the loading and the other two due to the coupled differential equations that govern the problem. The time-marching employs the Houbolt method. The four usual kinds of beams that are pinned–pinned, fixed–fixed, fixed–pinned and fixed–free, under uniformly distributed, concentrated, harmonic concentrated and impulsive loading, are analyzed. The results are compared with the available analytical solutions and with those furnished by the Finite Difference Method.     

Dynamic behavior of frame structures by boundary integral procedures

The Boundary Element Method is applied to synthesize a set of Boundary Integral Equations representing the uncoupled axial and flexural dynamic behavior of rectilinear Bernoulli–Euler beam elements in the frequency domain. In the sequence, these structural elements are coupled by the sub-region technique to model two-dimensional frame structures, in which the axial and flexural behaviors are coupled. This methodology is used to accurately recover modal data, eigenfrequencies and eigenmodes, of two frame structures. The usual Boundary Element procedure is recast to deliver simultaneously the values of variables at the element boundaries and at an arbitrary number of internal nodes. The inclusion of internal nodes allow to recover the structure eigenmodes and makes feasible the coupling of the assembled systems with a surrounding environment, for instance, an acoustic field. The results obtained are compared with a standard Finite Element eigenvalue analysis. It is shown that for increasing response frequencies, the Boundary Element scheme delivers modal data within a degree of accuracy, which is only obtained by the conventional Finite Element Method with considerable finer meshes.     

Displacement discontinuity formulation for modeling cracks in orthotropic shear deformable plates

A displacement discontinuity formulation is presented for modeling cracks in orthotropic Reisnner plates. Fundamental solutions for displacement discontinuity are derived for the first time using a Fourier transform method. Boundary integral equations are presented in terms of discontinuity rotations on the crack surfaces for opening mode problems. As the fundamental solutions have singularity of O (1/r ²), Chebyshev polynomials of the second kind are used to evaluate the integral equations. By solving for coefficients of the Chebyshev polynomials, the stress intensity factors at the crack tips are obtained directly. Comparisons are made with solutions using the finite element method to demonstrate that the displacement discontinuity method is an efficient and accurate method for solving crack problems in orthotropic Reissner plates.