采用有限单元法计算梁任意形状截面特性

采用将梁截面离散化的方式,用数值积分计算截面的几何特性,并根据梁剪切变形和扭转理论,利用变分原理建立截面的有限元法方程,求解任意形状截面的扭转常数、剪切中心以及剪切面积修正系数等特性.本方法适用于各种形式的截面,具有计算精度高及适应性强的特点.根据上述理论编制了相应程序,按照不同的单元划分方式,分别计算出矩形截面截面特性,与理论解进行比较;又对舟山市定海长峙至岙山预应力混凝土连续箱梁截面进行了计算,并与Ansys结果进行比较,均证明采用本文的计算方法能得到满意的结果,且该方法适用于各种形状的截面形式.     

计算叶片力学特性的三维8节点非协调有限单元法

针对叶片结构几何形状复杂的特点，建立了能进行叶片力学特性分析的三维8节点非协调有限单元模型。该模型采用了几何非线性和线弹性模式，来考虑叶片径向刚度远大于弯曲刚度而引起的在不同转速下初应力对叶片的作用，能够正确反映叶片的弯曲和扭转耦合振动模态、叶片组的切向和轴向振动模态以及它的应力状态。根据质量等效和力等效的原理，导出了在非协调单元中计算单元变形能时，记入附加的内部自由度，而计算单元的动能、体积力、表面力、以及阻尼力所做的虚功时，不记入附加的内部自由度，这是不同于目前一些文献中的提法。最后，通过实例计算验证了该模型的正确性。     

混凝土裂缝开展宽度计算的有限单元法

本文采用“最大拉应力单元开裂”准则，改进了非线性分析的迭代流程，使其最大可能地模拟实际结构的开裂情况，进而可以基于分布裂缝模式计算裂缝发生位置，裂缝间距，裂缝开展深度和宽度。算例表明，该法所计算的裂缝分布规律是正确的，计算裂缝宽度与试验实测值及规范公式计算值基本吻合，可以推广应用。     

有限单元质量质心集中方法的探讨

提出利用单元质量质心集中方法来优化单元惯性效应表征、改善单元力学效应匹配的思想: 将单元质量质心集中能够建立精确描述单元质点系惯性效应的动力学方程, 进而达到优化单元惯性效应表征以及改善单元惯性和变形两种力学效应匹配的目的. 本文通过一维问题的频散效应分析等对该方法的有效性进行了详细的分析和验证.     

基于能量平衡原理的有限单元耦合方法

结构非线性数值计算分析应真实反映局部损伤破坏细节,以作为损伤演化全过程分析的依据。对同类构件,有限单元耦合方法可以解决破坏细节与整体模拟的空间尺度差异问题。基于能量平衡原理,建立了梁与实体单元、梁与壳单元以及壳与实体单元的耦合方程,适用于结构的损伤数值计算。对某RC框架结构原位推覆试验的损伤数值分析表明,有限单元耦合模型能正确反映整体结构的承载力和变形性能,并且能准确反映局部损伤破坏细节。     

蜂窝密封动力特性计算方法的改进及应用

针对蜂窝密封的转子动力学特性分析,本文由双控制体模型入手,建立相应的控制方程,将其无量纲化,并通过摄动分析分解为零阶及一阶无量纲方程。通过零阶方程求得定常流场,确定流场是否发生阻塞现象,继而求解一阶方程得到密封的动力特性系数。本文的双控制体模型在蜂窝密封静子和转子表面采用Blasius模型的摩擦系数,且舍弃转子轴心传统的圆轨迹小扰动假设,直接设定为椭圆轨迹小扰动,加之矩阵形式的表达,得到了更为清晰的模型及更为简洁的算法。随后与实验结果相比较,验证模型及计算结果的准确性,并比较一般气体密封的计算结果,体现蜂窝密封在稳定性方面的优越性。最后分析了几种密封参数对其动力特性的影响,揭示某些规律并给出一些建议。     

一个统一单元矩阵计算的方法

1.引言有限元应用软件系统中,普遍采用有限元位移法。但单元矩阵(刚度矩阵、质量矩阵、阻尼矩阵、以至载荷列阵)等的计算,多数是逐一给出,部分是以元素簇给出的。本文将讨论不同类型元素的单元矩阵的统一算法和程序实施技巧,给出计算精密单元和几何非线性单元方便、有效的途径,达到以积累插值函数积累单元矩阵库的目的。     

有限单元法解弹性理论问题的一般方法

本文将位移函数取如座标变量的任意次多项式形式,讨论了用三角形单元族来求解弹性理论平面问题与用四面体单元族来求解弹性理论空间问题的一般方法,获得了通解。在此基础上,我们又编制了通用计算程序,以解决实际计算问题.     

非线性复合材料的一种杂交应力有限单元方法

建立了非线性复合材料模型的杂交应力有限元方法,并在材料主坐标系下提出直接方法计算单元非线性应力场,然后由此计算单元切线刚度矩阵和剩余载荷并转换到整体坐标系下,利用Newton-Raphson方法进行结构的位移迭代。在Hahn-Tsai非线性复合材料杂交元分析中,由位移和应力方程所导出求解单元非线性应力场的简单迭代法是条件收敛的,对较大载荷当迭代位移增加到一定程度以后无法得到应力收敛解。但是,利用本文提出的直接法由于完全避免了非线性应力场迭代,不仅很好地解决了这一问题,而且极大地提高了计算效率。数值算例说明该方法是确实有效的。     

一种改进的无单元方法

使用 1阶或 1阶以上最小滑动二乘法 ( MLS)形函数的无网格伽辽金法 ( EFGM) ,它们的主要缺点是形函数构造复杂、计算费用十分昂贵。本文提出了一种改进的无单元方法 ( IEFM) ,它通过采用 Shepard形函数 ( 0阶 MLS形函数 )对结点的覆盖位移函数加权求和来简化整体近似位移函数的构造 ,且能够避免 EFGM里求解结点形函数时矩阵的求逆及相乘计算。文中的数值算例表明 ,这种改进的 IEFM法收敛快、精度高 ,与标准的EFGM相比其计算时间得到了大幅度的减少     

可变网格的有限单元分析在边坡稳定计算中的应用

本文在非线性有限单元法的基础上，引进八节点的狭长单元来描述边坡中的滑动面，同时对于有限元单元网格可调节点坐标进行变更形成新的网格，于是利用原有基本网格架，通过多次有限元分析求出边坡中可能滑动面的稳定系数曲线Ｋ＝ｆ（α），并进而确定最危险的滑动面，在此基础上开展一系列稳定分析工作。     

谱体积方法单元分割与问题单元标记方法的改进研究

谱体积方法是一种本质上解决网格依赖性的高精度CFD计算方法,本文研究了二维Euler方程的谱体积方法,提出一种基于切比雪夫多项式的单元分割方法,建立了基于WENO的变量限制器方法,并发展了结合谱体积和控制体的问题单元标记方法.采用15°超声速压缩拐角和NACA0012跨声速流动两个典型算例进行验证,结果表明,该分区方法具有更好的计算精度,标记方法可有效识别不连续区域,在较少的网格下即可获得与密网格传统有限体积法相当的计算精度.     

夹层板的有限单元法

本文根据Hoff型夹层板弯曲理论,采用有限单元法导出了夹层板单元的刚度矩阵、质量矩阵、一致的载荷矢量和温度载荷矢量;用以计算一般夹层板的变形、应力或自由振动。其结果与理论解相比,具有较满意的精度。     

计算平面问题应力强度因子的任意高阶奇应变有限单元

W.K.Wilson提出的高阶奇应变圆单元(SSC)把有限单元法和裂纹尖端附近的线弹性解析解结合起来,能有效地计算应力强度因子K_Ⅰ或K_Ⅱ。A.Holston Jr.进一步把对称和反对称位移叠加起来,用以计算复合型应力强度因子K_Ⅰ和K_Ⅱ。但他们在对称和反对称的位移函数中最多只取到4项。胡海昌建议取更多的项,以便把圆单元的半径     

侵蚀滑移计算方法的改进

将二维拉格朗日有限元程序LTZ-2D中侵蚀滑移的计算方法进行了改进.将单元侵蚀判断标准由单一的累积塑性应变标准改为累积塑性应变、单元的硬化因子以及单元相对于滑移线的变形角度的综合判断标准.同时,利用本文的改进算法对侵彻问题进行了数值模拟.计算结果显示,采用改进后的侵蚀滑移计算方法不仅很好地显示了侵彻过程中弹靶的接触面,...     

张力等参有限单元的构造及其在分离流动计算中的应用

本文利用张力样条函数，构造了一种十二节点张力等参有限单元。它不仅具有等参元的一切特点，而且具有张力可调的独特性质，从而在一定程度上可调整插值函数以达到较好的收敛收和计算精度。在分离流动计算中，尤其可以显示出它独特的优越性。本文的两个分离流动算例说明了这一点。     

层合板的非线性有限单元分析

对层合板引入小应变中转动几何非线性假设，位移模式考虑了横向剪切应力及其连续性，有采用罚单元的方式构造了简单，协调的八节点等参单元，并应用于层合板的应力分析，取得了很好的数值结果。     

一维有限单元的节点位移精度

研究高次杆单元和梁单元的节点位移精度问题.首先求出一端固支均匀杆和悬臂梁在任意次多项式形式分布载荷作用下的位移精确解,然后用二次杆单元、五次欧拉梁单元和三次铁木辛柯梁单元求得了节点位移.通过比较有限元解与精确解以及利用静力凝聚方法,发现一次以上杆单元、三次以上欧拉梁单元以及三次以上铁木辛柯梁单元都可以给出精确的端点位移.     

斜交板梁桥的有限单元算法

借助斜坐标与直角坐标的关系,构造了斜交板梁单元的位移函数,由此提出了斜交板梁桥的实用有限单元算法。     

基于mbS模式的有限单元法

mbS模式及其有限元法是在固体和结构分析模型中引入薄膜、弯曲和剪切理论,且采用纯拉压、纯弯和纯剪单元进行分析的数值方法。在时空系中剖分物质单元和时间单元上构造以指数函数和贝塞尔函数为插入函数且按Lagrange插值条件的薄膜、弯曲和剪切等基本位移函数,由此得到更加完备和耦合的固体和结构实体单元的变形模式,根据能量泛函变分原理得到静动力有限元基本方程的一致格式。研究表明,mbS模式及其有限元法可用于梁柱和板壳等结构的静动力分析及屈曲分析。