多层复合壳体三维振动分析的谱-微分求积混合法

对于较厚的多层复合壳体,其振动位移沿厚度方向呈锯齿形变化且层间剪切和拉、压应力呈三维耦合状态,采用传统的等效单层理论分析已不能满足精度要求.建立不受结构厚度、铺层材料性质和铺层方式限制的三维分析方法具有重要的研究价值.本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析新方法——谱-微分求积混合法.该方法应用三维弹性理论对独立铺层进行精确建模,有效克服了二维简化理论对横向变形以及层间应力估计不确切的缺点;引入微分求积技术对铺层进行数值离散,将三维偏微分问题转化为二维偏微分问题,降低了求解维度和难度;应用广义谱方法近似地表述离散计算面上的场变量,将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组,避免了对超越方程的求解.数值验证结果表明该方法收敛性好,计算精度高.     

集中载荷作用下梯度复合材料梁的小波-微分求积法分析

将梯度复合材料梁作为平面应力问题处理,采用小波和微分求积混合法,对集中荷载作用下结构的响应进行了分析.考虑材料特性参数沿高度方向呈梯度分布,在该方向上采用广义微分求积法进行离散;鉴于广义微分求积法求解集中荷载问题精度不高的缺点,在梁的长度方向上引入对突变信号敏感的小波插值函数.数值计算表明,小波-微分求积混合法不仅保留了广义微分求积法高效的优点,而且能够很好地模拟结构局部化特征.     

弹性地基上变厚度矩形板自由振动的GDQ法求解

基于广义微分求积法(GDQ法),对弹性地基上变厚度矩形板横向自由振动的控制微分方程及其不同边界条件进行离散,研究了其自由振动的频率特性。数值计算得到了不同长宽比?、不同厚度变化参数?、不同地基参数K条件下以及简支或固定边界条件下弹性地基上变厚度矩形板的量纲为一的振动频率,并与已有文献进行了比较。结果表明:运用广义微分求积法对弹性地基上变厚度矩形板的频率求解结果在退化到K=0时与幂级数解的结果非常吻合;在条件相同的情况下,采用广义微分求积法仅需较少的节点(N=M=13)就能达到满意的求解精度。本文的研究为求解此类问题的低阶、高阶振动频率提供了一种简便有效的数值方法。     

基于边界值方法的微分动力系统数值计算方法

对高维非线性初值问题,微分求积法在每一步的积分过程中需要求解一个更高维的非线性方程组,因而计算量巨大。基于微分求积法与边界值方法两者之间的关系,可以将广义向后差分方法和扩展的隐式梯形积分方法看作是经典微分求积法的稀疏表达形式。将广义向后差分方法以及扩展的隐式梯形积分方法这两类边界值方法应用于微分动力系统的数值计算,提出了一类新的数值计算方法。理论分析及算例结果表明,对高维非线性微分初值问题的数值计算,本文方法相对于经典的微分求积法具有更高的计算效率。     

用微分求积方法计算二维薄板在超音速流中的非线性颤振

采用了一种微分求积方法将二维薄板在超音速气流作用下的非线性动力学方程离散为常微分方程,并用Runge-Kutta数值方法进行了计算.为验证微分求积方法的结果,与伽辽金方法计算结果进行了比较,取得了一致的结果.微分求积法的计算结果用分叉图、相平面、时域曲线以及功率谱进行了描述,结果表明在特定的参数区间存在混沌运动,而通向混沌的道路是经过一系列周期倍化分叉产生的.     

有纤维增强的双层圆柱壳体的层间应力解析解

本文对一个两端简支具有任意厚度的双层圆柱壳休,在均布外压或内压作用下而引起的层间应力的基本性质进行了研究,为了简单起见,双层壳体采用一层带[O °]铺向的纤维增强层和另一层各向同性壳组成,文中直接应用三维弹性理论精确求解这一轴对称问题。每层壳体的位移和应力场都表达成富里叶和富氏-贝塞尔函数两个无穷级数之和的形式,文中还作了一些实例计算,数值计算结果绘成的曲线图分別表明其几何参数、材料常数、载荷形式以及铺层次序对层间应力的影响,其中显示出在某些情况下壳体两层之间将可能形成一个拉应力区,这一现象似应值得引起设计工作者的重视,此外,由于本文求解方法是严格的,因此所得分析结果可以作为比较或评价其它近似方法的依据。     

求解具有弹性接头桩基动力学大变形的微分—代数方法

本文采用弧坐标首先建立了求解具有弹性接头的桩基大变形分析的非线性动力学微分方程,其中, 广义Winkler模型用来模拟土对桩基的抗力.其次,在空间域内应用微分求积单元法来离散非线性微分方程组,并给出了处理弹性接头处连接条件的微分求积单元公式,得到了时间域内的一组微分-代数方程,采用二阶向后差分来代替二阶时间导数离散微分-代数方程组,得到一组离散化的非线性代数方程,应用Newton-Raphson方法求解了该非线性代数方程组.最后给出了数值算例,得到了桩基在顶部处受到组合动载荷作用时的响应,考察了弹性接头的刚度、位置对桩基动力学行为的影响.     

薄板的多变量小波有限元分析

基于区间B样条小波和广义变分原理,提出了多变量小波有限元法,构造了一种新的薄板多变量小波有限单元.由广义变分原理推导结构的多变量有限元列式,区间B样条小波尺度函数作为插值函数构造的多变量小波有限元法中,广义应力和应变被作为独立变量进行插值,避免了传统方法中应力应变求解的微分运算,减小了计算误差.区间B样条小波良好的数值...     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

混合微分求积法在功能梯度材料平板柱形弯曲问题中的应用

利用混合微分求积法,对任意荷载作用下不同材料梯度分布的功能梯度材料平板柱形弯曲问题进行了分析。针对广义微分求积法求解集中荷载问题精度不高的缺点,本文利用小波微分求积法进行了改进。由于小波对突变信号具有良好的自适应描述能力,因此在平板宽度方向上,利用小波微分求积法可以有效地处理集中荷载;而在材料梯度变化的板厚方向上,则利用广义微分求积法计算量小且精度高的特点进行离散计算。计算表明,混合微分求积法不仅保留了广义微分求积法高效的特点,而且能有效地求解任意荷载作用的问题。通过算例,分析了在机械荷载作用下,材料不同梯度形式、平板上下表面材料性质差异对功能梯度平板结构响应的影响。     

Three-dimensional free vibration analysis of rotating laminated conical shells: layerwise differential quadrature (LW-DQ) method

This paper focuses on the free vibration analysis of thick, rotating laminated composite conical shells with different boundary conditions based on the three-dimensional theory, using the layerwise differential quadrature method (LW-DQM). The equations of motion are derived applying the Hamilton’s principle. In order to accurately account for the thickness effects, the layerwise theory is used to discretize the equations of motion and the related boundary conditions through the thickness of the shells. Then, the equations of motion as well as the boundary condition equations are transformed into a set of algebraic equation applying the DQM in the meridional direction. This study demonstrates the applicability, accuracy, stability and the fast rate of convergence of the present method, for free vibration analyses of rotating thick laminated conical shells. The presented results are compared with those of other shell theories obtained using conventional methods and a special case where the angle of the conical shell approaches zero, that is, a cylindrical shell and excellent agreements are achieved.     

Free vibration of functionally graded truncated conical shells under internal pressure

The influence of internal pressure on the free vibration behavior of functionally graded (FG) truncated conical shells are investigated based on the first-order shear deformation theory (FSDT) of shells. The initial mechanical stresses are obtained by solving the static equilibrium equations. Using Hamilton’s principle and by including the influences of initial stresses, the free vibration equations of motion around this equilibrium state together with the related boundary conditions are derived. The material properties are assumed to be graded in the thickness direction. The differential quadrature method (DQM) as an efficient and accurate numerical tool is adopted to discretize the governing equations and the related boundary conditions. The convergence behavior of the method is numerically investigated and its accuracy is demonstrated by comparing the results in the limit cases with existing solutions in literature. Finally, the effects of internal pressure together with the material property graded index, the semi-vertex angle and the other geometrical parameters on the frequency parameters of the FG truncated conical shells subjected to different boundary conditions are studied.     

横观各向同性层状场地对环形简谐荷载的位移响应

本文采用横观各同性层状粘弹性模型拟半空间上之的层状场地,用阻尼器取代其下部的半空间以吸收振动能量。     

A simple and accurate mixed FE-DQ formulation for free vibration of rectangular and skew Mindlin plates with general boundary conditions

A simple and accurate mixed finite element-differential quadrature formulation is proposed to study the free vibration of rectangular and skew Mindlin plates with general boundary conditions. In this technique, the original plate problem is reduced to two simple bar (or beam) problems. One bar problem is discretized by the finite element method (FEM) while the other by the differential quadrature method (DQM). The mixed method, in general, combines the geometry flexibility of the FEM and high accuracy and efficiency of the DQM and its implementation is more easier and simpler than the case where the FEM or DQM is fully applied to the problem. Moreover, the proposed formulation is free of the shear locking phenomenon that may be encountered in the conventional shear deformable finite elements. A simple scheme is also presented to exactly implement the mixed natural boundary conditions of the plate problem. The versatility, accuracy and efficiency of the proposed method for free vibration analysis of rectangular and skew Mindlin plates are tested against other solution procedures. It is revealed that the proposed method can produce highly accurate solutions for the natural frequencies of rectangular and skew Mindlin plates with general boundary conditions.     

Circular shallow spherical shells with central circular hole under simultaneous actions of arbitrary unsteady temperature field and arbitrary dynamic normal load

In this paper, under assumption that tempeature is linearly distributed along the thickness of the shell, we deal with problems as indicated in the title and obtain general solutions of them which are expressed in analytic form.In the first part, we investigate free vibration of circular shallow spherical shells with circular holes at the center under usual arbitrary boundary conditions. As an example, we calculate fundamental natural frequency of a circular shallow spherical shell whose edge is fixed (m=0). Results we get are expressed in analytic form and check well with E. Reissner’s [1]. Method for calculating frequency equation is recently suggested by Chien Wei-zang and is to be introduced in appendix 3.In the second part, we investigate forced vibration of shells as indicated in the title under arbitrary harmonic temperature field and arbitrary harmonic dynamic normal load.In the third part, we investigate forced vibration of the above mentioned shells with initial conditions under arbitrary unsteady temperature field and arbitrary normal load.In appendix 1 and 2, we discuss how to express displacement boundary conditions with stress function and boundary conditions in the case m=1.     

FREE VIBRATION OF PIEZOELECTRIC CYLINDRICAL SHELLS

Three displacement functions are introduced to represent each mechanical displacementaccording to the 3-D theory in this paper,By expanding the displacement functions and the electric po-tential in orthogonal series,the free vibration equation of piezoelectric cylindrical shells satisfying SS3boundary conditions can be obtained.The equation was solved by utilizing Bessel functions with com-plex arguments.Results are presented graphically as well as in table,and compared with those of otherreferences.Some frequencies that were missing in Ref.[9]are discovered.     

ELASTIC DYNAMIC ANALYSIS OF MODERATELY THICK PLATE USING MESHLESS LRPIM

A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function （RBF） coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.     

A layerwise boundary integral equation model for layers and layered media

A hybrid method is presented for the analysis of layers, plates, and multilayered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multilayered system using a total potential energy formulation. The layerwise laminate theory of Reddy is employed to develop a layerwise, two-dimensional, displacement-based, hybrid boundary element model that assumes piecewise continuous distribution of the displacement components through the system's thickness. A one-dimensional finite element model is used for the analysis of the multilayered system through its thickness, and integral Fourier transforms are used to obtain the exact solution for the in-plane problem. Explicit expressions are obtained for the fundamental solution of a typical infinite layer (element) assuming linear displacement distribution through its thickness. This fundamental solution is given in a closed form in the cartesian space, and it can be applied in the two-dimensional boundary integral equation model to analyze layered structures with finite dimensions. The proposed method provides a simple, efficient, and versatile model for a three-dimensional analysis of thick plates or multilayered systems.     

非均质中厚板的无网格LRPIM动力学分析

用局部加权残值法建立了非均质中厚板的局部径向点插值离散系统方程,采用无网格局部径向点插值法分析了非均质中厚板的自由振动和强迫振动问题。用径向基函数耦合多项式基函数来近似试函数,用四次样条函数做为加权残值法中的权函数。所构造的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件。该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行。在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化。计算结果表明,利用该方法计算非均质中厚板的自由振动和强迫振动问题可以得到具有较高精度的解。     

人行荷载下梁式结构振动分析的DQ-IQ混合法

为了评估人行荷载作用下梁式结构的振动舒适度，利用微分求积-积分求积，即DQ-IQ混合法求解移动荷载作用下梁的振动响应。人行荷载作用下梁式结构的振动控制方程是含Dirac函数的偏微分方程，首先利用IQ法离散与时间相关的Dirac函数，再利用DQ法把控制方程转化为二阶常系数微分方程，最后利用Newmark算法求解微分方程。以某钢结构连廊为例，利用DQ法计算结构自振频率并与解析解进行对比，结果验证了节点选取和边界条件施加的合理性，再利用DQ-IQ混合法和振型叠加法分别计算了不同行走步频下连廊的响应，计算结果表明，DQ-IQ混合法具有较高的可靠性和精确性。DQ-IQ混合法也可以推广到诸如车辆荷载作用下路面或桥梁的动力响应等其他移动荷载下结构的振动分析。