张拉结构非线性分析两节点曲线单元有限元法

提出适合张拉结构几何非线性分析的两节点曲线单元有限元方法．假定索元的初始形状呈二次抛物线,根据单索的平衡条件、几何和物理关系建立了索元的位移函数;由拉格朗日应变的定义建立了可以考虑任意次高阶位移影响的索元轴向应变的精确表达式,并基于拉格明日描述方法和虚功原理得到了索元的非线性平衡方程与切线刚度矩阵．采用荷载增量法与Newton－Raphson法相结合的混合法进行了实例计算,结果表明：本文方法的精度明显优于两节点直线索单元,适合于大跨度索阿、索穹顶等张拉结构的几何非线性分析．     

有限元广义应力影响函数定理及算法

当一个定向的单位集中力在结构上移动时,结构某处某应力分量的值是该力位置的函数,称为应力影响函数.本文作者曾经提出并证明了连续体和有限元模型的平均应力影响函数定理,设计了算法,编制了程序,可用于求各类单元组成的复杂结构的应力影响函数,并已应用于工程中。本文发展上述理论,提出有限元模型广义应力影响函数定理,据此可用间接法一次求得任意单元中任一点应力分量或由应力分量的线性组合表示的物理量的影响函数。讨论了单元基准节点位移模式的单元特征向量求法。     

薄板弯曲问题的P型杂交解析有限元方法

本文采用两套变量构造有限元试函数空间，在单元内部要求试函数精确满足平衡微分方程，在单元边界上对位移和转角分别用Ｐｅａｎｏ升阶函数插值，然后利用广义变分原理建立了一种薄板弯曲问题的Ｐ型杂交解析有限方法，与常规有限元法相比，该方法不心进行过细的网格剖分，通过增加单元插值多项式的阶数Ｐ来提高精度，此外，该方法还具有积分计算只需在单元边界上进行、单元钢度矩阵和载荷向量具有嵌入结构、协调程度可以自动控制等优     

划线法构建Serendipity单元插值函数的几点思考

插值函数的构建在有限元分析中起着至关重要的作用,本文从一个变节点四边形等参单元的插值函数的构建出发,分析了划线法构建Serendipity单元插值函数时应注意的几点问题,提出了划线法构建插值函数时选择直线或曲线的原则,并通过3个例子验证了该结论.     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 2D ELASTICITY

针对边界元法中高阶单元中几乎奇异积分计算难题，解剖了二维边界元法高阶单元的几何特征，定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近，从而分离出积分核中主导的奇异函数部分，其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分，再对奇异核函数积分导出解析公式，从而建立了一种新的半解析法，用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例，采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力，并与线性边界元正则化算法结果作了比较，结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构，采用正则化算法的线性边界元分析比有限元有显著优势，而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。     

爆炸荷载作用下地下拱形结构动力响应样条小波有限元研究

采用小波有限元方法研究爆炸荷载作用下地下结构的动力响应,克服在模拟过程中由于材料的奇异性、地质条件的复杂性和加载的快速性出现的应力集中和计算效率低下,根据样条有限点法,构造了单向和双向区间B样条圆环扇形小波单元,用尺度函数作为插值函数;结合工程实例,通过Matlab软件编程对爆炸荷载作用地下拱形结构的动力响应进行了数值模拟,并与应用传统有限元程序模拟结果进行对比。结果表明,小波有限元用很少单元取得了较高的精度,计算效率比传统有限元提高一倍。     

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

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

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.     

结构动态应力计算方法研究

复杂结构动态应力的准确计算是一个没有圆满解决的问题。本文以最小余能原理为基础，提出了计算结构动态应力的最小伤痛有法。该方法采用二次分析思想，首先采用常规有限元对结构进行适当离散，计算输出结构所需应力区域的有限元结点位移和加速度动力时程反应，再应用最小余能法计算所求部位的动态应力值。这种方法的优点是它可以与现有的有限元程序有机结合，方便使用；动应力在区域内的分布规律可以由计算者根据具体情况而确定，一般情况下，可以选用二次曲线来逼近动应力在区域内的实际分布，避免了常规有限元法计算结构动应力时必须对单元形函数求导的做法，从而提高了动应力计算精度。计算结果表明：本文方法计算结构动态应力结果较常规有限元法的计算结果有明显改进，特别是当结构变化剧烈时，改进效果更为明显。     

区间B样条小波平面弹性及Mindlin板单元构造研究

基于二维张量积区间B样条小波及小波有限元理论,构造了一类用于分析弹性力学平面问题和中厚板问题的C0型区间B样条小波板单元。在二维小波单元的构造过程中,传统多项式插值被二维区间B样条小波尺度函数取代,进而构造形状函数和单元。与小波Galerkin方法不同,本文构造的区间B样条小波单元通过转换矩阵将无明确物理意义的小波插值系数转换到物理空间。区间B样条小波单元同时具有传统有限元和B样条函数数值逼近精度高及多种用于结构分析的基函数的优点。数值算例表明:与传统有限元和解析解相比,本文构造的二维小波单元具有求解精度高,单元数量和自由度少等优点。     

一种空间缆索结构静力分析的解析元法

将空间缆索结构简化为具有拉伸刚度的质点系,给出了缆索结构空间解析元法的基本方程和求解方法,单元间的作用力与坐标变化的关系可以用解析法得到,对所得到的反映结构特性的质点系方程组进行力的平衡迭代,求解方程组.采用自动的动态可变步长的迭代方法,能够提高计算效率,保证收敛.这种方法既考虑了几何非线性,又适用于材料非线性的计算,比有限元法优越之处还在于,它不用求解线性方程组,所以适用范围广,允许求解多自由度的几何可变体系,而有限元法在求解此类问题时经常不收敛.     

悬索结构非线性分析的滑移索单元法

本文提出了一种悬索结构存在在滑移单元时的非线性计算方法,利用文献（１）提出的五节点曲线索单元有限元基本方程,建立了滑移单元非线性分析的计算理论、方法及公式,并进行了实例计算,结果表明,本文方法是行之有效的,可供悬索结构设计、施工时参考。     

降落伞膜索非线性有限元模拟

用膜单元和索单元模拟降落伞织物绳索系统,基于完全拉格朗日格式的非线性有限元方法编程计算降落伞的结构动力学特性。采用增量与迭代混合方法改善非线性计算的收敛特性并结合HHT隐式时间推进方法减小整体迭代计算量。使用修正应力应变张量导数的方法模拟膜单元单向应力状态并针对膜单元和索单元分别进行了非线性有限元计算验证。最后针对C-9型降落伞建立三维有限元模型,根据设定流速对伞衣施加均匀压强载荷,将模拟展开的结果与使用相同模型、不同方法商业软件的文献进行对比,显示了隐式非线性有限元方法模拟降落伞膜索系统大变形动力学的能力。     

预应力钢桁架结构分析中的摩擦滑移索单元

为精确模拟预应力钢桁架中连续长索在支撑点的滑动,本文创建了一种考虑摩擦力影响的新单元。被称为摩擦滑移索单元的新单元有三个节点,中间节点为支撑点。本文首先利用弹性悬链线的解析解,建立了弹性悬链线单元,并推导了单元两端点的张拉刚度。摩擦滑移索单元由两个弹性悬链线单元组合而成,根据支撑点处索的滑动方向、索力差、滑动摩擦力和滑移刚度调整两索段的原长,使支撑点两侧的索力满足给定的摩擦关系。算例验证了新单元算法的正确性和高效性,对设计的预应力钢桁架的分析,显示了张拉过程中及使用荷载作用下新单元在结构分析中的应用。新单元可直接用于常规的有限元分析,研究在工作状态或施工中存在连续长索滑移的索结构。     

悬链线索单元算法的改进

H．B．Jayaraman在20世纪80年代推导的悬链线索元有限元法计算精度高，特别适用于精度要求比较高的大型索结构。但是，当索原长Lu的取值与悬索两节点之间的直线长度相近时，迭代不易收敛，甚至发散。提出了当该迭代不收敛时，应采用的迭代策略。计算结果表明，该方法准确，计算精度高，可供悬索结构设计、施工时参考。     

索托桥结构分析的滑移索单元法

为在索托桥的结构分析中精确模拟连续长索的滑动,本文创建了一种新的单元。被称为“滑移索单元”的新单元有三个节点,以点接触的形式模拟索从下方绕过滑轮,它可以通过自动调整两侧索段的长度而使单元处于平衡状态,从而简化了计算。新单元算法的推导基于有限元分析的基本原理和弹性悬链线的解析解,并利用了平衡状态时单元内力之间的关系。本文介绍滑移索单元的推导过程,用设计的算例验证了它的正确性,分析了连续长索的滑移对索托桥桥面竖向变形的影响。新单元可以直接用于常规的有限元分析中,研究处于工作状态或在施工中的索结构。     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

基于欧拉描述的两节点索单元非线性有限元法

本文针对柔性悬索结构几何非线性分析的特点,提出了一种用欧拉描述来表示的两节点索单元非线性有限元模型,在索元变形后的位置上由虚功能建立了非线有限元基本方程及切线刚度矩阵。这样建立的非线性有限元分析方法可充分考虑拉索的几何非线性特性的影响并给悬索结构的初始平衡分析带来方便,算例结果表明,本文方法是精确有效的。     

Super element approach to cable passing through multiple pulleys

A parametric super element model for cable passing through multiple pulleys is presented in this study for the static analysis of structures. The amounts of cable passages over pulleys are introduced as additional degrees-of-freedoms in the finite element model and the relationship between cable tensions at the two sides of each pulley is imposed based on the friction law or empirical data. The proposed finite element model is firstly verified by a simple pulley cable system and then applied to the analysis of real complex engineering structures. The verification results satisfy the static equilibrium and deformation compatibility conditions of the structural system and basic engineering principles. With the application of the proposed super element model, the global deformation and stress distribution for structures with multiple-pulley cable systems can be effectively and accurately computed. Numerical results for structural analysis show that the effect of friction of pulleys on the cable tensions is significant and the friction-free and fixed models both give unrealistic and incorrect results in cable tensions in some cases.