一个高精度八结点六面体单元

本文提出一个用于分析弹性空间问题的八结点六面体单元Q_(c11)用等参拟协调元方法,建立了该单元的单元函数显式。Q_(c11)单元具有24个结点参数和9个内部自由度,对长方体单元,其结果与Q_6单元相同,能精确描述纯弯曲,对任意六面体单元能保证通过分片试验。该单元具有构造简单,计算方便,精度高等优点。实际上,这里提供了一种构造非协调等参元单元函数的方法。     

关于计算结构力学的一些看法

评述了我国计算结构力学发展的方向性问题,强调了学科的应用性质和发展计算结构力学应用软件的迫切性,尤其应注意微机的工程应用,以适合我国国情。     

有旋转自由度的高精度四边形单元

从具有三次位移模式的平面等参元出发 ,通过对单元自由度进行线性变换并引入连续介质力学中关于旋转度的定义 ,构造出角结点有旋转自由度的高精度四边形单元。该单元与平面等参元有相同的单元特征及精度     

计算杆件结构的自由度

为了简化和统一结构力学教材中杆件体系计算自由度W的算法,本文提出了"直接法"。该方法以相邻节点间的杆件和结点为分析对象,对于含m根杆件、g个单刚结点和h个单铰结点的平面杆件体系,直接法给出计算自由度公式W=3m-(3g+2h)。相比于已有方法,直接法研究对象清晰,数学运算统一,适合于编程计算,便于初学者准确而快速计算复杂体系的W;也为结构超静定次数的确定提供了一种简单有效的数值算法。     

16自由度的中厚板弯曲矩形单元

基于有限条带思想,引入结点扭率自由度,利用深梁单元的位移模式建立了一个4结点16自由度中厚板弯曲高阶单元,此单元是薄板单元BFS-16的推广形式,其特点是单元的横向位移、转角位移、剪应变位移模式直接构造,在边界上位移模式与深梁单元一致,方便与梁单元叠加,适应于带加劲肋的板弯曲问题分析,用于薄壁结构时可考虑翘曲。实例计算显示,此单元精度高,计算稳定,收敛快,无剪切闭锁现象,能较好地反映中厚板的边界效应。     

有旋转自由度的高精度三角形单元

本文从三次平面等参元出发，构造出角结点有旋转自由度的三角形单元。该单元与三次平面等参元有相同的单元特征及精度，同时具有连续介质力学中关于旋转度的定义。     

两个自由度的单质点体系自振频率及振型计算

对于具有两个自由度的单质点体系自由振动,先假定一个振动方向,求出该振动方向的柔度系数,对柔度系数求极值,满足柔度极大值的方向即是第1振型方向,由振型正交性原理可确定第2振型. 由此将两个自由度体系的计算变成了单自由度体系的计算,根据求出的柔度系数,可方便地求出自振频率.     

一种精确积分的四结点板弯曲单元

本文从Ｍｉｎｄｌｉｎ／Ｒｅｉｓｓｎｅｒ理论出发,采用一种新的平行四边形母单元和相应的形函数推导四结点板弯曲单元刚度矩阵的精确积分解。弯曲应变和横向剪切应变分别采用不同的插值公式构成单元刚度矩阵。理论和算例分析表明本文方法克服了“闭锁”现象并能应用于很薄的板,单元刚度矩阵计算速度比采用数值积分计算的同类单元的快四倍。     

结点位移计算的一种简单方法

推导了一种结点位移计算的简单方法,在平面桁架结点位移的计算和静不定杆系结构的变形几何关系的确定中应用方便.     

关于壳体的极限承载能力

一、前言 薄壳结构在近代各种工程技术部门中到处得到广泛的应用和重视。普通的实用弹性理论只能给出材料在弹性限度以内的应力和变形状态,这时壳体在强度方面的储备,还远没有竭尽,因此壳体极限承载能力的研究是近代塑性力学中最有实际意义的课题之一。在这方面,伊留兴()、然尼村()、拉包脱诺夫(     

12结点三维等参奇异单元的构造和应用

通过改变三维8结点六面体等参单元的结点位置、结点数目和形函数,构造了一种12结点三维等参奇异单元,该单元的应力场具有1/(√r)奇异性,可以模拟裂缝前沿的奇异应力场;该单元的位移模式在其中两个坐标方向是线性变化的,因此,该单元与线性单元连接时不需要过渡单元,仍能保证交界面位移协调,克服了20结点三维等参奇异单元不能与线性单元协调连接的缺陷;文章最后将该奇异单元布置在裂缝前沿,应用有限元法计算了三点弯曲梁预制裂缝前沿的应力强度因子,该结果与规范公式计算值基本一致.     

壳体计算的位势——有限元方法

本文提出了壳体计算中的一种位势——有限元方法。这种方法的特点是首先利用位势理论将壳体理论的边值问题化成一组积分微分方程,然后再将之离散而求解。与一般有限元方法相比,位势——有限元法的形状函数只需保证位移本身的连续性;另外,离散化过程并不需改变壳体原来的几何形状,因而可以方便地计算各种变曲率壳体。     

结点应力修正计算中边值问题的简便处理

在有限元分析中,对于利用数值积分的曲面(线)单元,我们发现在积分点上的应力精度远高于单元边界结点上的应力精度。本文以积分点的应力计算为基础,用带线性的约束的最小二乘法,简便地处理了曲面(线)单元结点应力修正计算中的边值问题。     

改进的八至二十结点三维非协调等参元

本文对R,L,Taylor等人提出的平面非协调元QM6满足分片检验修正条件,提出了进一步的修正建议,并将其应用于八至二十结点等参元,从而得到了一类适应空间不规则网格应力分析的非协调元QMM6,该单元与QM6单元相比,具有计算量小,在畸变网格条件下精度高等优点。     

薄膜单元和六面体单元的耦合计算模型

本文提出薄膜单元和六面体二次等参单元的耦合计算模型.由所编制的计算程序进行实例计算,计算结果与解析解一致.     

一个简洁有效的带旋转自由度的平面三角形单元

从原始的Allman单元[1]出发,构造出一个简洁有效的带有旋转自由度的平面三角形单元。该单元基本保持了Allman单元的精度,同时消除了多余的零能模式,而且使转动自由度等同于连续介质力学中关于旋转分量的定义,并对旋转自由度自身的正确性进行了考核。     

非匹配结点的有限单元等参插值方法及其应用研究

针对两零件的异构网格单元结点在接触界面不能相互匹配导致结点属性不能连续过渡和传递的问题,提出非匹配结点的有限单元等参插值方法,通过构建非匹配结点的形函数和修正原结点的形函数,将结点属性值的影响范围限制在可控的局部区域,从而实现两异构网格结点属性在接触界面的连续过渡和传递。通过两个异构的四边形单元网格的结点属性在接触界面的过渡实例和啮合齿轮的接触分析应用,验证了该方法的正确性和有效性。     

两个新型的八结点块体单元—三维离散算子差分法

给出了弹性力学三维问题的离散算子差分法 ,讨论离散算子差分法在三维问题中的特点 ,意在为该方法的进一步发展提供依据 ,为应用弱形式进行数值求解的研究提供参考。本文从弹性力学平衡方程更为一般的弱形式出发 ,给出了含边界参数的弱形式方程。由该方程不仅可以得到有限元法 ,还可得到离散算子差分法。给出了两个八结点块体单元 ,虽然单元中位移函数是非协调的 ,不需特殊处理便可保证离散格式收敛 ,并对单元位移有十分好的反映能力。     

计算曲面壳体等效节点热载荷的新方法

提出了一种利用平面壳单元计算曲面壳体热应力问题等效节点热载荷的新方法, 具有较高的计算精度. 首先在平面壳元的理论基础上, 在壳单元切平面建立局部坐标系; 然后根据提出的理论, 利用单元节点整体坐标直接计算壳单元等效节点热载荷积分方程中所需的未知量, 如: 形函数对局部坐标的导数、从对局部坐标积分转换到自然坐标积分时的雅可比行列式等; 最后, 根据提出的算法求出从局部坐标转换到整体坐标的转换矩阵, 进而求出整体坐标系下壳单元等效节点热载荷. 通过与商用软件ANSYS 的计算结果进行对比分析, 证明提出的方法是正确而且精确的.