含面内转动自由度的广义协调曲面矩形扁壳元

扁壳单元中引入结点转角自由度可以在不增加结点的情况下,增加位移场的阶次,提高计算精度,从而显著地提高单元性能。同时在单元中引入泡状位移场,能有效地扩大了单元位移场的解空间,所构造的单元具有计算精度高、对计算网格畸变不敏感的优良特性。本文利用广义协调薄板单元RGC-12的位移函数作为扁壳元的法向位移,利广义协调矩形膜元的位移函数作为扁壳面的切向位移,通过附加面内转动自由度构造了一个具有24个自由度的4结点广义协调曲面矩形扁壳元GRC-S24。在此基础上再增加一个广义泡状位移,又构造了一个具有更高计算精度的曲面矩形扁壳元GRC-S24M。并通过实例分析对这两个单元的收敛性和精度进行了验证。     

一种对网格凹凸畸变不敏感的高应力精度4结点带旋转自由度平面杂交应力函数单元

论文给出一种简单的高性能带旋转自由度4结点四边形平面单元.该单元的理论基础与卞学鐄先生首个杂交应力元相似,也从最小余能原理出发,无需单元内部位移场.但是应力场试函数采用Ariy应力函数的基本解析解(基于直角坐标的多项式),并强调其对坐标的对称和完备性.这样假设的应力场可以同时满足平衡和协调方程,因而更加合理.而单元边界位移则采用著名的Allman模式(采用局部坐标),即考虑结点转动自由度的二次协调位移.与其他同类单元相比,本文单元对位移和应力展现出更高的精度,特别是应力解答尤其突出.更有趣的是,单元对网格畸变非常不敏感,即使单元退化为三角形和凹四边形,仍然能保持较高的计算精度.此外,单元没有方向依赖性等缺点.     

基于解析试函数三角形带旋转自由度膜元的显式格式研究

三角形单元是有限元分析中常用的单元.在平面单元内引入结点转动自由度,可以提高单元位移场的阶次,在不增加单元结点的前提下提高单元性能.论文利用问题基本解析解作为试函数来构造带旋转自由度的三角形单元ATF-R3H,采用了杂交应力函数单元模式,确保了单元优良的抗畸变性能和较高应力计算精度.论文利用直角坐标与三角形面积坐标的线性关系,以及面积坐标幂函数在三角形域内和边界上的积分公式,直接给出单元刚度矩阵的显式表达式,从而避免了大量数值积分,提高了计算效率.数值算例表明,显式格式的ATF-R3H单元具有良好的性能.     

基于拟应变法的壳单元大变形分析及在冲击动力问题中的应用

为了简便有效地解决板壳结构的大变形问题，本文针对八节点相对自由度壳单元进行研究。该单元的位移场由壳的中面节点位移和上表面节点的相对位移组成，不带有转动变量。所有的研究都是基于完全的三维位移、应力、应变场。采用拟应变法，对应变场另行假设，能够改善该单元在大变形情况下的计算精度。通过引入Wilson非协调模式，构造了大变形情况下的拟应变场表达式，给出了该单元用于解决非线性动力分析问题的有限元求解方程。通过算例表明，本文针对相对自由度壳单元提出的方法及推导的公式，能够解决冲击动力问题中的大变形问题。     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法，在此基础上建立起了构造等参非协调元的新方法，作为实例，构造出两个可以给出单刚显式的四结点平面非协调新单元。     

含一个无外力圆柱面带转动的三维杂交应力元

根据一种修正的余能原理,建立了一类具有一个无外力圆柱表面及结点含转动自由度的8 结点新型三维杂交应力元. 单元边界位移场选择二次位移插值函数,且与相邻元协调;单元内假定应力场满足以柱坐标表示的平衡方程及圆柱面上无外力边界条件. 数值算例表明,这种特殊杂交应力元在相当粗的网格下即能十分准确地分析圆弧形槽口附近及曲梁的三维（及二维）的孔边应力分布.      

含一个无外力圆柱面带转动的三维杂交应力元简

根据一种修正的余能原理,建立了一类具有一个无外力圆柱表面及结点含转动自由度的8 结点新型三维杂交应力元. 单元边界位移场选择二次位移插值函数,且与相邻元协调;单元内假定应力场满足以柱坐标表示的平衡方程及圆柱面上无外力边界条件. 数值算例表明,这种特殊杂交应力元在相当粗的网格下即能十分准确地分析圆弧形槽口附近及曲梁的三维（及二维）的孔边应力分布.     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法,在此基础上建立起了构造等参非协调元的新方法。作为实例,构造出两个可以给出单刚显式的四结点平面非协调新单元。     

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

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

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

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

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

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

Geometrically exact beam theory without Euler angles

In modeling highly flexible beams undergoing arbitrary rigid–elastic deformations, difficulties exist in describing large rotations using rotational variables, including three Euler angles, two Euler angles, one principal rotation angle plus three direction cosines of the principal rotation axis, four Euler parameters, three Rodrigues parameters, and three modified Rodrigues parameters. The main problem is that such rotational variables are either sequence-dependent and/or spatially discontinuous because they are not mechanics-based variables. Hence, they are not appropriate for use as nodal degrees of freedom in total-Lagrangian finite-element modeling. Moreover, it is difficult to apply boundary conditions on such discontinuous and/or sequence-dependent rotational variables. This paper presents a new geometrically exact beam theory that uses no rotation variables and has no singular points in the spatial domain. The theory fully accounts for geometric nonlinearities and initial curvatures by using Jaumann strains, exact coordinate transformations, and orthogonal virtual rotations. The derivations are presented in detail, fully nonlinear governing equations and boundary conditions are presented, a finite element formulation is included, and the corresponding governing equations for numerically exact analysis using a multiple shooting method is also derived. Numerical examples are used to illustrate the problems of using rotational variables and to demonstrate the accuracy of the proposed geometrically exact displacement-based beam theory.     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

A simple and efficient triangular finite element for plate bending

A simple and efficient triangular finite element is introduced for plate bending application. The element is a three-node triangular one with three basic degrees of freedom per node and two internal rotation degrees of freedom, using selective reduced integration. Numerical examples indicate that, despite its simplicity, the element is not only competitively accurate, but also useful as a thick/thin triangular plate bending element. It is also pointed out that this element using selective reduced integration is in fact a mixed element.     

应用平板壳单元DKQl6分析板壳结构的几何非线性

应用新近开发的四边形十六自由度离Kirchhoff平板壳单元DKQl6，分析了板壳结构的几何非线性问题，采用Total Lagrange格式，在小应交、中等转动的假定下，建立了该单元几何刚度阵和大位移矩阵．非线性方程采用位移引导或弧长引导的牛顿-拉夫森增量迭代法求解．讨论了网格和加载步效对收敛性的影响，通过对典型算例的计算以及与其它单元的比较，说明了DKQl6单元在板壳结构几何非线性分析中也有良好的精度．     

弹性力学中的一种非协调数值流形方法

通过引入数学和物理双重网格，将插值域与 积分域分别定义在不同的覆盖上，即在数学网格上进行插值函数的构造，物理网格上完成 系统能量泛函积分运算，最后通过覆盖权函数将二者联结在一起. 它的优点是单元网格划 分随意，不受复杂边界形状和二相材料界面的限制，单元可以是任意形状，是较之于有限 元方法更一般的数值模拟方法. 在4节点四边形数值流形方法中，由于单元总体位移函数 包含的完全多项式不完全，使得计算精度不够精确，为此，在单元总体位移函数上附 加非协调位移基本项，使之趋于完全，提出了弹性力学问题的一种改进的数值流形 方法------非协调数值流形方法. 通过内部自由度静力凝聚处理，导出了消除内参后的单元应变矩阵 和单元刚度矩阵，使得在不增加广义节点自由度的前提下，大大提高了数值流形方法的计 算精度和计算效率. 同时对非协调项进行了显式处理，可以对工程实践起到更切实的帮助. 数值试验表明，它们能够保证收敛，有较高的精度，对畸变不敏感，从而证明了该方法的 可行性.     

平面Cosserat模型有限元分析的4和8节点单元与分片试验研究

经典连续体理论不包括物质内部尺度,当考虑应变软化问题时,有限元结果对网格具有很强的依赖性。与经典连续介质力学理论不同,Cosserat连续体模型在传统平动自由度的基础上添加了一独立的旋转自由度,在本构模型中引入了内尺度参数。本文研究了基于Cosserat理论的平面4和8节点等参元以及8（4）节点线、角位移混合插值等参单元,给出Cosserat单元分片试验的实施过程。最后将单元运用到小孔应力集中问题的分析当中,通过计算结果与理论解的比较,表明了4和8节点以及8（4）节点等参元的适用性,为问题的非线性分析打下基础。