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

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

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

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

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

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

带旋转自由度拟协调三角形板壳单元

将拟协调三角形罚函数板单元和Allman二次膜位移插值模式相结合,通过在膜内增加一个旋转自由度参数,构造一种新的Mindlin三角形板壳单元。采用这一单元对板壳结构的线性和几何非线性问题作了分析。数值计算表明,该单元不仅克服了用板单元拟合壳体分析中的病态,而且使单元在保持弯曲精度的同时大大提高了膜内变形的精度。     

两个新的空间单元

给出了两种改进协调元的方法，并建立了两个新的八节点空间单元。第一个单元 是直接构造不协调的单元函数，并增加对不协调函数的分片试验约束．该单元的应力 计算精度较作者以前提出的单元Ｑｃ１１有较大的改善。第二个单元是广义杂交模型． 它是以广义变分原理为根据．通过调整单元内的应力、应变参数．实现了用杂交法建 立高精度的八节点空间单元．相对卞学提出的基于Ｒｅｉｓｓｎｅｒ原理的应力杂交模型， 它避免了推导时所需的几何摄动。     

用样条函数求解板梁组合结构

本文根据子结构的原理,用样条函数作为位移试函数,按最小二乘配点法。分析板梁组合结构,在IBM微机上所得的算例表明,这种方法精度令人满意.可用于加权残值法中关于组合体结构的计算问题。图2,表1,参4。     

精化直接刚度法及九参数三角形薄板单元

提出一种可直接用于精化不协调元的精化直接刚度法．其列式与直接刚度法类同，单元间的协调条件平均满足；既能保证收敛又能提高精度．精化直接刚度法的变分根据是广义变分原理．用精化直接刚度法对著名的九参数三角形Ｚｉｅｎｋｉｅｗｉｃｚ薄板单元进行了精化，建立了新的九参数三角形板元ＲＴ9和ＲＴ１０．数值结果表明本文建立的单元收敛且精度高．     

Kirchhoff板弯曲问题的改进的复变量无单元Galerkin方法

基于改进的复变量移动最小二乘法,建立了Kirchhoff板弯曲问题的改进的复变量无单元Galerkin方法.相对于移动最小二乘法,改进的复变量移动最小二乘法采用一维基函数建立二维问题的逼近函数,提高了形函数计算效率.由改进的复变量移动最小二乘法建立Kirchhoff板的挠度逼近函数,根据Kirchhoff板弯曲问题的Galerkin弱形式建立离散方程,并应用罚函数法施加本质边界条件,推导了Kirchhoff板弯曲问题的改进的复变量无单元Galerkin法的公式.通过对4个典型算例进行计算和分析,说明了本文建立的Kirchhoff板弯曲问题的改进的复变量无单元Galerkin方法的有效性,并通过分析数值解的精度对本文方法中如何选取合适的基函数、权函数、影响域比例参数、节点分布和罚因子进行了讨论.数值算例说明了本文方法具有较好的收敛性和较高的计算精度.     

一种混合型厚薄通用四边形板单元

据Reissner假定[1]，运用广义变分原理，以节点的位移和弯矩，扭矩为基本未知量，构造了一种混合型板单元 ，该元既可用于分析需考虑效应的中厚板，又可用于分析不计剪切效应的薄板，同时还可解决带有文克尔地基与变温作用的问题，例题计算表明：该元收敛快且内力与位移均有较高精度。     

正交异性钢桥面板局部振动计算的组合板梁单元法

针对大跨度钢桥中采用梯形加劲肋的正交异性钢桥面板,提出一种计算钢桥面板局部振动的组合板梁单元法.其中顶板用平板壳单元分析,梯形加劲肋视为板梁单元,两者的位移模式根据板与肋的变形协调关系建立.组合板梁单元的刚度矩阵通过能量变分原理得到,利用各子单元的形函数可以求得组合板梁单元的一致质量矩阵与一致荷载列阵,进而得到相应的有...     

复合材料层合板的应力杂交有限元分析

本文根据修正的余能变分原理,构造了一个适合于复合材料纤维层合板特点的矩形应力杂交板弯单元,该单元能够考虑横向剪切变形的影响和局部扭曲效应.在动力分析中,我们选取单元内部位移场与边界位移场相协调.通过计算表明,本文单元具有自由度少,使用方便,计算精度高等特点.     

矩形薄板与无限层土地基的相互作用

用半解析单元法分析了矩形薄板与无限层土地基的静动力相互作用，提出了板单元和地基单元半解析位移函数。采用在水平两方向解析，在竖向进行离散化的方法，使三维相互作用的问题简化为一维数值计算问题。根据半解析单元法编制的程序操作简单，收敛速度快，提高了计算精度，大大减少了计算工作量。数值计算表明，在工程上必须注意基础附近的表层土的弹性模量的确定，而较深处层土的特性可粗略估计。     

任意边界条件下矩形板的面内自由振动特性

采用改进傅里叶级数法(IFSM)对矩形板在任意边界下的面内自由振动特性进行了研究.将结构的位移容许函数表示为包含正弦三角级数的改进傅里叶级数,正弦三角级数的引入能够有效地解决在边界处存在的不连续或者跳跃现象;将位移容许函数的未知傅里叶展开系数看作广义变量,采用能量原理建立结构的能量泛函,结合Rayleigh-Ritz法对未知傅里叶展开系数求极值,将矩形板的面内问题转换为一个标准特征值求解问题.通过大量的数值算例,并与现有文献中解及有限元方法计算结果进行对比,验证了文中方法的正确性,结果还显示文中方法具有良好的收敛速度与计算精度.     

非协调平面等参元NQ_9的改进—NQ_7元

本文对位移函数在直角坐标下二次完备的等参非协调元NQ_9作了改进。在不影响函数完备性的条件下,将非协调位移参数由原来的5个减少到3个。文中给出若干算例及数值比较,新单元NQ_7具有优越的数值性。     

杂交位移模型的实施及其与拟协调元的关系

本文详细地讨论了杂交位移模型与放松连续性条件的胡海昌－■津广义变分原理的联系，并做为广义杂交元的特殊情况，以三角形变厚度簿板单元为例，讨论了杂交位移模型的实施问题． 本文还证明了拟协调元方法的一种形式可等价于杂交位移模型．     

有限元增强型分片检验

针对常应力分片检验理论上的不严格和不能做Mindlin 板非零常剪力及细观应变梯度理论非零常应变梯度曲率分片检验的问题, 基于对应齐次阶微分方程的放松连续条件的不协调元的变分原理, 建立了通过分片检验的单体条件及被检验单元的收敛条件: 除通过分片检验外, 单元函数还应包含刚体位移和常应变模式,无伪零能模式和满足弱连续条件. 建立了对应非齐次阶微分方程的放松连续条件的不协调元的变分原理和增强型分片检验条件及单体条件, 通过增强型分片检验条件的单元的收敛条件是单元函数应包含刚体位移和满足平衡的非零应变模式, 无伪零能模式和新的弱连续条件. 提出的增强型分片检验条件是对齐次和非齐次阶微分方程的分片检验统一提法. 对Mindlin 板问题建立了非零常剪力分片检验, 对细观偶应力-应变梯度理论问题建立了非零常应变梯度曲率C^0-1分片检验.     

广义节点无网格法基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法,在阐述这种方法基本原理的同时,针对线弹性力学问题给出了其计算列式.与传统无网格方法相比,这种方法更具有一般性,当选取0阶广义节点位移插值函数时便可得到传统的无网格法;可以通过提高广义节点位移插值函数的阶数降低完备基函数的次数,从而可减少支持域内节点的数目并保证计算精度.最后通过一端承受剪力悬臂梁和中间开口无穷板算例分析,论证了这种方法的合理性.     

基于位移形函数和力形函数的梁柱单元理论对比研究

基于位移形函数的梁柱单元理论被广泛应用于框架结构梁柱单元的有限元分析,但其计算精度受到位移形函数的限制而难以提高.对比基于位移形函数的相关理论,总结基于力形函数的梁柱单元理论,进行有限元程序编制,计算过程中采用柱面弧长法来进行迭代求解.程序计算结果对比表明,基于力形函数的梁柱单元计算效率和计算精度均高于基于位移形函数的梁柱单元.     

无单元迦辽金法在对称叠层板屈曲问题中的应用

无单元法是一种新的数值算法,它信息处理简单,其基本思想是在计算域上用一些离散的点由最小二乘法来拟合场函数,从而摆脱单元的限制.因为它应用了滑动最小二乘法,所以与有限元不同,它的近似场函数不一定精确通过计算点,这会使计算稍趋复杂,但与这种方法带来的方便性相比,这点缺陷是微不足道的,因此这种方法愈来愈引起人们的关注.将无单元法用于对称叠层板的屈曲问题,通过计算不同边界条件下各向同性板和对称叠层板的屈曲荷载参数,验证了该法的有效性,且精度较高,有着广泛的工程应用前景.     

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

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