含任意椭圆核各向异性板杂交应力有限元

基于含椭圆核有限大各向异性板弹性问题的复变函数级数解,应用杂交变分原理建立了一种与常规有限元相协调的含任意椭圆核各向异性板杂交应力有限元.单元内的应力场和位移场采用满足平衡方程、几何方程与物理方程的复变函数级数解,假设的复变函数级数解精确满足椭圆核边界处的位移协调条件和应力连续条件,单元外边界上的位移场按常规有限元位移场假设,单元内椭圆核的长轴可以与材料主轴不重合.单元刚度矩阵采用Gauss积分求得,并给出了建立刚度矩阵的主要公式和推倒过程.数值计算结果表明该单元具有计算精度高、计算工作量小等优点.     

基于面积坐标和解析试函数的厚薄通用板单元

求出了Mindlin-Reissner厚板理论控制方程的基本解析解,以其作为试探函数构造了两个基于第二类四边形面积坐标的广义协调厚薄通用板单元AATF-PQ4a和AATF-PQ4b;在此基础上计算了承受均布荷载方板的中心挠度和中心弯矩。结果表明:两种单元随网格加密而趋近于精确解,表现出较高的精度和较好的收敛性,而且AATF-PQ4b的收敛速度比AATF-PQ4a快;对于薄板极限情况,这两种厚薄通用板单元自动退化为薄板单元,完全无剪切闭锁现象发生;两种板单元均对网格畸变不敏感,具有较好的稳定性,适用于实际工程计算。     

用通用单元计算板的固有振动和弹性稳定

本文用依据二类变量广义变分原理构造的一种新型的考虑横向剪切影响的任意四边形四节点通用单元[8],研究板的固有振动和弹性稳定问题,计算了各种类型板的固有频率、振型和临界载荷,得到了满意的结果,此单元具有通用性好、公式推导和程序编制简单等特点,并有显式的刚度矩阵和质量矩阵,可节省计算时间。     

厚薄板通用三角形广义协调元TCGC-T9的固有振动分析

研究了厚薄板通用三角形广义协调单元TCGC-T9的固有振动特性.通过编制程序,采用子空间迭代法求得不同板厚的前5阶固有频率值,并将计算结果与理论解和有限元软件ABAQUS的结果进行比较.结果表明:TCGC-T9单元在固有振动分析方面收敛速度快,数值稳定性好、实用性强(无需区分厚板和薄板).     

厚薄板通用三角形位移元

本文构造出了一种具有九个自由度的厚薄板通用三角形位移单元，并给出了单元刚度矩阵显式。这种单元以其简洁的常规位移元列式可在相当宽的板厚变化范围内（包括板厚为零）获得很高的计算精度，其结果可与相应的矩形单元相比。而且不会出现剪切自锁。     

中厚板弯曲问题的自然单元法

自然单元法是一种新兴的无网格数值计算方法,基于Reissner-Mindlin板弯曲理论,将自然单元法应用于平板弯曲问题的计算中,给出了相关的公式,推导了总体刚度矩阵和荷载列阵的计算列式.算例分析表明,自然单元法应用于中厚板的弯曲问题具有较高的计算精度,并可用于Winkler地基上基础板的计算.同时指出,对于厚跨比较小的薄板,由于对挠度和中面法线转角采用相同的插值形式,当板厚变薄时夸大了虚假的剪切变形影响,因而表现出剪切自锁现象.对进一步开发厚薄板通用的计算程序作了初步探讨.     

薄膜横向振动的有限元分析

构造了一种3节点三角形膜单元,以适用于平面薄膜横向振动的有限元分析.在给出单元形函数的基础上,根据最小势能原理建立了薄膜自由振动方程,并推导了单元刚度矩阵和单元质量矩阵.研究结果表明,单元刚度矩阵和单元质量矩阵形式简单,且自由度少;通过两个典型算例,证明3节点三角形膜单元的计算结果非常接近理论解,同时可以达到很高的精度...     

约束层阻尼板动力学问题的半解析解

利用条形传递函数方法(SDTFM)得到了约束层阻尼(CLD)板动力学问题的半解析解.首先对CLD板沿纵向离散成多个条形单元,基于Hamilton原理推导出条形单元的刚度矩阵和质量矩阵,仿照有限元法组集得到系统的总刚度矩阵和总质量矩阵.经Laplace变换后引入状态向量,采用分布参数传递函数方法在状态空间内建立CLD板的控制方程并进行求解.最后以对边固支和悬臂CLD板为例,得到了板的动力学特性和频响曲线,并与NASTRAN或相关文献结果进行了比较,吻合良好,验证了该方法的有效性.从推导过程和算例可以看出,该方法所需的单元数目少,获得的是半解析解,计算效率高且准确可靠.     

有限元软件中通用梁单元的统一形成方法

本文根据作者参与多个有限元软件开发的经验,讨论并给出了满足多种工程需求的通用梁单元的刚度矩阵及其单元载荷向量形成方法,特别是基于柔度法的Euler梁单元、Timoshenko梁单元和变截面梁单元的统一形成方式,以及不同梁端条件下的单元刚度矩阵与载荷向量的统一形成方式,便于工程有限元软件开发人员参考借鉴,并希望能对有限元方法的教学实践有所帮助。     

基于开放式结构有限元系统 SiPESC.FEMS平板壳单元研发

基于开放式结构有限元系统SiPESC.FEMS的单元计算模块的设计模式,研发设计一种通用的平板壳单元计算框架。考虑板壳单元的组合关系和程序编制过程中的重用性及灵活性等特点,采用了软件设计中的构造器(Builder)模式实现不同的组合单元。本框架具有很好的通用性和可扩展性,为有限元程序研发提供了一个新的方式;同时,系统能够处理复杂荷载和边界条件,并可灵活实现不同类型单元的组合分析。本文利用此方法构造五种平板壳单元,通过数值算例分析对比讨论其性能,为选取合适的平板壳单元类型进行结构数值分析提供参考。     

Flexural vibration analysis of an eccentric annular Mindlin plate

A weak-form quadrature element method is presented to study the flexural vibrations of an eccentric annular Mindlin plate. Typical combinations of boundary conditions are considered and the natural frequencies are obtained for both thin and moderately thick plates. All results are verified using the commercial computer code ANSYS. Excellent agreement is reached in all cases. Comparison of the present predictions with other available results for thin plates is also made.     

曲线加筋Kirchhoff-Mindlin板自由振动分析

相比传统加筋板,曲线加筋板能够更充分地发挥材料力学性能.在加筋板力学分析中,厚板通常采用Reissner-Mindlin理论,然而当板厚较薄时易出现剪切自锁,离散的Kirchhoff-Mindlin理论采用假设剪切应变场可避免该问题.针对曲线加筋Kirchhoff-Mindlin板自由振动分析,采用离散的Kirchhoff-Mindlin三角形单元和Timoshenko曲梁单元分别模拟板和加强筋,根据板的位移插值函数及筋板交界面的位移协调条件,建立基于板单元位移自由度的有限元方程.为了验证方法的有效性和准确性,采用直线加筋薄板、曲线加筋薄板和厚板3种模型进行算例研究,通过收敛性和精度分析来选择合理的有限元网格密度.直线加筋薄板前20阶固有频率均与文献结果吻合良好;曲线加筋板算例中,本文方法满足收敛条件的板单元数目为2469,Nastran模型板单元数目为6243;本文所得曲线加筋板固有频率与Nastran计算结果最大误差为3.4%.研究结果表明,本文方法无需筋板单元共节点,可使用较少的有限元网格数量,并能够保证计算精度;在离散Kirchhoff-Mindlin三角形板单元基础上构造Timoshenko梁单元可同时适用于曲线加筋薄板与厚板自由振动分析.     

Stability and vibration of shear deformable plates––first order and higher order analyses

This work presents the highly accurate numerical calculation of the natural frequencies and buckling loads for thick elastic rectangular plates with various combinations of boundary conditions. The Reissener–Mindlin first order shear deformation plate theory and the higher order shear deformation plate theory of Reddy have been applied to the plate’s analysis. The governing equations and the boundary conditions are derived using the dynamic version of the principle of minimum of the total energy. The solution is obtained by the extended Kantorovich method. This approach is combined with the exact element method for the vibration and stability analysis of compressed members, which provides for the derivation of the exact dynamic stiffness matrix including the effect of in-plane and inertia forces. The large number of numerical examples demonstrates the applicability and versatility of the present method. The results obtained by both shear deformation theories are compared with those obtained by the classical thin plate’s theory and with published results. Many new results are given too.     

基于平面弹性—板弯曲模拟关系的薄板有限元列式——从膜单元到薄板单元

本文全面讨论了基于平面弹性－－板弯曲模拟关系的薄板有限单元的理论和方法,由于直接对弯矩函数进行插值,c1连续性的要求得以自然避免,薄板单元可以直接在c0连续的层面上加以构造,无需借用Reissner-Mindlin的中厚板理论,由之引发的闭锁问题也得以避免,本文系统地阐明了平面弹性膜单元与薄板弯曲单元的对应关系,及由平面弹性膜单元的向薄板弯曲单元转换的一整套方法。为薄板单元的构造提供了一条新的有余的途径,文中给出了对应于平面弹性膜单元CST,LST,Q4,Q8的薄板单元,我们称之为MPS板单元,MPS板元以挠度和转角为自由度,便于实际应用,和其它板单元相比具有非常高的精度。     

The influences of both offset and mass moment of inertia of a tip mass on the dynamics of a centrifugally stiffened visco-elastic beam

The present study is concerned with the out-of-plane vibrations of a rotating, internally damped (Kelvin-Voigt model) Bernoulli-Euler beam carrying a tip mass. The centroid of the tip mass, possessing also a mass moment of inertia is offset from the free end of the beam and is located along its extended axis. This system can be thought of as an extremely simplified model of a helicopter rotor blade or a blade of an auto-cooling fan. The differential eigenvalue problem is solved by using Frobenius method of solution in power series. The characteristic equation is then solved numerically. The simulation results are tabulated for a variety of the nondimensional rotational speeds, tip mass, tip mass offset, mass moment of inertia and internal damping parameters. These are compared with the results of a conventional finite element modeling as well, and excellent agreement is obtained. Some numerical results are given in graphical form. The numerical results obtained, indicate clearly that the tip mass offset and mass moment of inertia are important parameters on the eigencharacteristics of rotating beams so that they have to be included in the modeling process.     

Effects of rotary inertia on sub- and super-critical free vibration of an axially moving beam

The most important issue in the vibration study of an engineering system is dynamics modeling. Axially moving continua is often discussed without the inertia produced by the rotation of the continua section. The main goal of this paper is to discover the effects of rotary inertia on the free vibration characteristics of an axially moving beam in the sub-critical and super-critical regime. Specifically, an integro-partial-differential nonlinear equation is modeled for the transverse vibration of the moving beam based on the generalized Hamilton principle. Then the effects of rotary inertia on the natural frequencies, the critical speed, post-buckling vibration frequencies are presented. Two kinds of boundary conditions are also compared. In super-critical speed range, the straight configuration of the axially moving beam loses its stability. The buckling configurations are derived from the corresponding nonlinear static equilibrium equation. Then the natural frequencies of the post-buckling vibration of the super-critical moving beam are calculated by using local linearization theory. By comparing the critical speed and the vibration frequencies in the sub-critical and super-critical regime, the effects of the inertia moment due to beam section rotation are investigated. Several interesting phenomena are disclosed. For examples, without rotary inertia, the study overestimates the stability of the axially moving beam. Moreover, the relative differences between the super-critical fundamental frequencies of the two theories may increase with an increasing beam length.     

Vibration band gap analysis of a new periodic beam model using GDQR method

In this study, a new periodic beam model is introduced. This beam consists of the concentrated rigid masses and tapered beam elements with linearly variable width. The theoretical equations are derived by employing the Euler-Bernoulli beam and the Bloch–Floquet theorem and then solved using the generalized differential quadrature rule method to calculate the first two band gaps. The effects of the mass, mass moment of inertia and taper ratio on the widths and central frequencies of the first two band gaps are investigated. Results show that the wide band gaps at low frequency ranges can be obtained by changing the geometrical parameters. This is of interest for different applications of the band gap phenomenon such as broadband piezoelectric energy harvesting. Finally, the finite element simulation (ANSYS software) is used to validate the analytical method and good agreement is found.     

An improved HSFR method for natural vibration analysis of an immersed cylinder pile with a tip mass

Immersed cylinder piles are usually modelled as immersed cantilever cylinder columns carrying a tip mass and rotary moment of inertia. In this paper, the equations of motion of an immersed cylinder pile along transversal modes of vibration are developed. Compressibility of water and structural damping are included in the formulation. Natural frequencies of the immersed pile are obtained from the developed equations using harmonic sweep frequency response analyses. The proposed method is applied to numerical examples, and the results obtained are shown satisfactory when compared to other numerical solutions in the literature, or to finite element solutions and experimental data.     

平稳/非平稳激励下中厚圆柱壳随机振动响应的基准解

圆柱壳结构被广泛应用于航空航天、船舶、汽车工程等领域.由于服役环境复杂,圆柱壳会受到随机激励作用,从而产生随机振动响应.本文针对考虑横向剪切变形和转动惯量的中厚圆柱壳,将虚拟激励法拓展到连续体结构,高效获得了各类随机激励下响应均方根的基准解.首先,开展了简支条件下中厚圆柱壳的自由振动分析,精确求得各阶自振频率和解析振型...