组合杂交有限元法对Wilson元的改进

利用组合杂交有限元法能够在几乎不增加计算量的情况下增强低阶位移格式数值精度的特点,作者讨论了采用最简常应力模式时组合杂交格式对非协调Wilson元的改进.数值试验表明改进后的格式能达到比Wilson元更高的数值精度.     

采用常应力模式的8节点组合杂交六面体有限元

作者讨论了采用常应力模式的组合杂交有限元方法对协调三线性H8-六面体元的改进.位移逼近使用两种方式:连续等参三线性插值和非协调Wilson位移模式.数值试验表明这两种新单元能显著改进H8-元的粗网格精度.由于应力参数可在单元水平消去,新方法的计算量与H8-元相当.     

试探函数不满足divu=0的定常Stokes方程的类Wilson元

利用研究二阶问题的非协调类Wilson任意四边形单元的特殊收敛性,通过新的技巧,并基于Falk等人的基础,讨论了类Wilson元对定常Stokes方程逼近的一种不需要试探函数满足divu=0条件的有限元方法.同时给出了与传统协调元及非协调元完全相同的最优误差估计,从而大大扩展了其应用范围.     

关于奇应变拟协调元的一点注记

1 引言用奇应变元与常规元结合的有限元法来计算裂纹尖端的应力强度因子,由于它们各自取不同的位移模式,所以单元与单元的接触面上是不协调的.通常不考虑这种不协调性.这种做法必然会产生两个问题:(a)对不协调的单元不能直接使用虚功原理来列式;(b)奇异单元必须取得很小,从而要补充过渡单元.     

考虑变厚度的杂交应力元

提出了用于变厚度板平面应力问题分析的杂交元列式,其中以位移场和膜力场为独立未知场量.基于Pian-Sumihara四边形杂交单元的五参数应力场和常规双线性位移场构造一个变厚度平面四边形杂交元.通过若干数值算例验证了新单元的正确性和优越性.     

二维时间分数阶扩散方程Wilson元的收敛性分析

利用Wilson元提出了一类二维时间分数阶扩散方程的新的全离散逼近格式.基于单元的性质,在不需要外推和插值后处理技术的前提下,得到了u的比传统的H1-范数更大的模意义下相应的O(h~2+τ~(2-α/2))阶的误差分析结果,正好比通常的关于Wilson元的误差估计高出一阶.这里,h,τ分别表示空间剖分参数和时间步长.     

具有转角自由度的曲面矩形扁壳元

扁壳单元中引入结点转角自由度可以在不增加结点的情况下,增加位移场的阶次,提高计算精度,从而显著地提高单元性能。利用广义协调薄板单元RGC-12的位移函数作为扁壳元的法向位移,利用广义协调矩形膜元的位移函数作为扁壳面的切向位移,构造了一个具有转角自由度的4结点广义协调曲面矩形扁壳元。并通过实例分析对该单元的收敛性和精度进行了验证,数值算例的结果表明,该单元收敛稳定迅速,用较少的自由度就获得了很高的精度,是一个性能良好的壳单元。     

试探函数不满足divu^→=0的定常Stokes方程的类Wilson元

利用研究二阶问题的非协调类Wilson任意四边形单元的特殊收敛性，通过新的技巧，并基于Falk等人的基础，讨论了类Wilson元对定常Stokes方程逼近的一种不需要试探函数满足divu^→=0条件的有限元方法。同时给出了与传统协调元及非协调元完全相同的最优误差估计，从而大大扩展了其应用范围。     

杂交宏元方法的收敛性分析

对一个4节点的低阶杂交应力四边形宏元方法进行了理论分析.该宏元采用连续分片线性位移插值逼近和分片独立设计的5参数自平衡应力模式.分析表明，单元上采用连续分片线性位移与采用等参双线性位移具有等价性，从而证明了有限元解的存在唯一性，并导出相应的误差估计.     

无剪切闭锁的厚薄板矩形单元的构造

在板弯曲的有限元计算中 ,因为剪切闭锁现象导致厚板元无法自动退化为薄板元 ,工程中很难构造既适用于厚板又适用于薄板的计算单元 ,因此采用由龙驭球教授提出的广义协调元法 ,构造了矩形厚薄板弯曲通用计算单元 .通过铁木辛柯梁理论 ,建立了板单元边界和板单元内部的位移插值场函数 .计算结果表明该单元可顺利解决剪切闭锁现象 ,可在实际工程中推广应用     

Poisson方程的组合杂交有限元方法

将能量优化思想应用到Poisson方程的一类边值问题数值求解中,得到一种组合杂交有限元计算格式,利用Wilson非协调位移和Taylor非协调位移构造了两组低阶四边形有限元．相应于双二次元Q2,计算量较小且精度接近．     

一种耦合算法的三维弹性空腔声振特性分析

针对三角形单元边界适应性强、位移插值简单、自由度低的特点,利用面积坐标构造了平板三角形单元位移插值函数,以对平板进行了动力学建模分析。此后,为了对三维弹性空腔声振问题进行较高精度的数值分析,弹性板采用三角形单元进行离散处理,而声腔域使用了在声学问题上具有精度高的光滑有限元方法建模。充分考虑结构和声场的耦合作用,建立了该类结构声振分析的一种新模型,并分析了弹性空腔声振特性。通过数值计算对比表明,这种耦合模型计算规模较小、计算结果精确,可进一步应用于形状更复杂的三维弹性空腔声振特性分析中。     

两阶椭圆边值问题的cascadic多重网格方法

证明了两阶椭圆边值问题的ｃａｓｃａｄｉｃ多重网格方法,对于ＰＩ非协调元、Ｃａｒｅｙ非协调元、Ｗｉｌｓｏｎ非协调元均可以达到最优。     

应用功的互等定理推导不连续位移基本解

不连续位移基本解在边界元法中得到了广泛的应用,本文应用功的互等定理推导了该基本解,指出了开尔文解中的某些应力分量就是不连续位移基本解中某些位移分量的一般规律.这种推导方法不但十分简单而且可以推广到各种区域的不连续位移基本解的推求.     

柔性结构非线性分析的杆单元有限元法

本文根据非线性弹性理论，采用拉格朗日坐标描述法，在三维整体坐标系下，直接根据应变的定义建立了精确应变的非线性几何关系．从而可以考虑任意次高阶位移的影响，得到了轴向应力应变的直接关系，根据虚功原理推导了杆单元非线性有限元增量方程及切线刚度矩阵，利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法进行了实例计算．结果表明本法精度很高，适合于柔性结构分析、设计时采用．     

索穹顶结构非线性分析的曲线索单元有限元法

针对目前在国内外引起广泛关注的大跨度索穹顶结构建议的了一种五节点曲线索单元有限元法，利用四次多项式作为位移插值函数，由应变的定义建立了可以考虑任意高阶位移影响的非线性应几何关系，从非线性弹性理论出发并基于Ｌａｇｒａｎｇｉａｎ坐标推导了胡限元基本方程，进行了实例计算。     

悬索结构大位移分析改进的两节点索单元模型

基于非线性弹性理论,提出了适合悬索结构大位移分析的改进的节点索单元有限元模型,在 移函数的表达式中计及了索元自重垂度的影响,从而克服了目前悬索结构普遍采用的现代了点索元无法直接计及自重垂度影响的缺点。     

矩形板动力刚化有限元分析

在有限元方法中首次引入了单元耦合形函数（阵）,以此将单元弹性位移表示成为单元结点位移的二阶小量形式,利用几何非线性的应变－位移关系式,在小变形假设条件下确定了单元耦合形函数,在此基础上,根据Ｋａｎｅ方程,运用模态坐标压缩,并采用适当的线性化处理,得到了包含动力刚度项的线性动力学方程,针对矩形板编制了动力刚化有限元分析程序,仿真算例证明了理论和算法的正确性。     

悬索结构非线性有限元分析

提出了一种适合于悬索结构空间非线性分析的三节点等参数索单元有限元模型,在三维整体坐标下,由就的定义得到了可以考虑任意高阶位移影响的非线性应变几何关系,采用Ｌａｇｒａｎｇｉａｎ坐标描述法,并由虚功原理建立了非线性有限无限元增量方程及切线刚度和矩阵,算例计算表明,本文方法可用于各种悬索结构的计算分析。     

子域样条有限元与无限元耦合法分析相邻建筑物的相互作用

运用惯性耦合法分析地基上相邻建筑物的相互作用，并将样条有限元法分析直接承受建筑物的桩基和土体，用半解析无限元法分析外围半无限土体。此方法不仅具有输入数据少，计算工作量小，并且避免了人为有限边界的影响。数值计算表明，当两结构相邻较近时，单体与双体在自振特性，尤其在位移、基底总剪力等方面均有差异。