拟协调等腰梯形壳元显式几何刚度阵及屈曲分析

采用拟协调元方法推导了等腰梯形薄壳元的显式几何刚度阵，用于组合结构屈曲分析的计算，结果表明，这种单元的几何刚度阵收敛快、精度好。     

一类梁系结构线性稳定性分析的简化处理

本文给出了旋转壳与支承支柱系统组合结构线性稳定性分析的简化处理方案，即将空间梁系结构每根梁的单元刚度了四及几何阵按离散富氏系数变量叠加得等效刚度阵及几何阵。从所给算例智 ，这一简化处理是成功的，也说明稳定性分析反映的是结构的总体效应，局部区域的简化处理对整体结构临界载荷影响不大，但可使计算量大大减少。本文的计算方案及程序已应用于实际的冷却塔设计。     

一类梁系结构线性稳定性分析的简化处理

本文给出了旋转壳与支承支柱系统组合结构线性稳定性分析的简化处理方案,即将空间梁系结构每根梁的单元刚度阵及几何阵按离散富氏系数变量叠加得等效刚度阵及几何阵。从所给算例知,这一简化处理是成功的,也说明稳定性分析反映的是结构的总体效应,局部区域的简化处理对整体结构临界载荷影响不大,但可使计算量大大减少。本文的计算方案及程序已应用于实际的冷却塔设计。     

对“杂交/混合有限元增量大位移分析”一文的讨论

原文提出了一个杂交/混合几何非线性有限元列式,根据这个列式,在满足位移边界条件后,单元的切线刚度阵为 K=G~TH~(-1)G(1)单元的节点力向量为 Q=(Q_1+Q_2)-G~TH~-1L(2)两式的符号均请参阅原文。(1)式的切线刚度阵K与初应力无关,也就是说K中不包     

采用离散Kirchhoff假定的三角形板单元的稳定性分析

1.满足离散Kirchhoff假定三角形板单元的几何刚度阵用有限元法和特征值问题的数值解相结合来求解板的稳定问题,是解决实际问题的重要和有效的方法.     

等几何壳体分析与形状优化

首先基于Reissner-Mindlin理论进行了三维壳体等几何分析,而后基于此对三维壳体进行形状优化,提出了形状优化中灵敏度的全解析计算方法,包括位移应变阵、雅克比阵和刚度阵等相对控制顶点位置的灵敏度解析计算公式;通过实例验证了壳体等几何分析和灵敏度全解析计算方法的有效性。与传统的基于网格的灵敏度半解析计算方法相比,基于NURBS的灵敏度全解析计算具有精确、计算效率高的特点,且可以避免优化迭代中的网格畸变。     

等几何壳体分析与形状优化

首先基于Reissner-Mindlin理论进行了三维壳体等几何分析,而后基于此对三维壳体进行形状优化,提出了形状优化中灵敏度的全解析计算方法,包括位移应变阵、雅克比阵和刚度阵等相对控制顶点位置的灵敏度解析计算公式;通过实例验证了壳体等几何分析和灵敏度全解析计算方法的有效性。与传统的基于网格的灵敏度半解析计算方法相比,基于NURBS的灵敏度全解析计算具有精确、计算效率高的特点,且可以避免优化迭代中的网格畸变。     

二次特征矩阵表示的特征值有界性分析

采用二次特征矩阵近似表示精确有限元和精确动态子结构分析中得出的超越或非线性动态刚度阵，并证明了在满足一定条件的前提下，二次特征阵给出的特征值是精确刚度阵特征值的上界或下界。     

弹塑性大变形分析的一致性高阶无单元伽辽金法

采用无单元伽辽金法求解弹塑性大变形问题。充分利用无单元法易于建立高阶近似函数的优点，位移采用二阶移动最小二乘近似。在更新拉格朗日方法的框架下，通过对控制方程弱形式的线性化建立了内力率的表达式，并区分为材料和几何两部分。采用Hughes-Winget算法更新应力，建立了Newton-Raphson迭代求解所需的一致切线刚度阵。刚度阵的数值积分采用近来针对小变形分析建立的二阶一致三点积分格式QC3（Quadratically Consistent 3-point integration scheme）。数值结果证明了本文方法分析弹塑性大变形问题的有效性和优越性。     

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

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

近似平衡多项式加速度动力显式算法

利用已知初始时刻的信息，建立一种可以取到任意阶高精度的多项式加速度单步隐式算法。在该隐式方法中，待采解方程纽系数矩阵中质量阵的系数远远大于阻尼阵和剐度阵的系数，略去非对角阻尼阵和非对角刚度阵对方程组的影响，得到一种近似平衡多项式加速度动力显式计算方法。此方法的精度主要由加速度多项式插值的项数、步长、质量阵的每件数、质量刚度比（质量阵和刚度阵的范数之比）决定。在此基础上给出了这种算法的通式，进行了精度分析，结果表明：如果时间步长h足够短，n次加速度近似平衡动力显式算法的精度可以达到O（hn＋1）。算例采用5次加速度近似平衡显式算法，计算结果的精确性证明了本算法的可行性。     

动态有限元的阻抗阵

对于动态有限元,宜引用阻抗阵来代替文献中使用的刚度及质量矩阵。基于虚位移原理,本文建立了阻抗阵各阶层式的计算公式.最后以薄膜为例,计算了前九阶频率。计算结果表明,动态有限元比静态元精度高。     

复合材料箱型结构的屈曲性态研究

采用带旋转自由度参数C^0类板（壳）单元，研究复合材料箱型结构的初始屈曲性态，在推导单元几何刚度阵时，同时考虑出平面位移和面内位移对膜应变的影响；着重讨论了不同材料正方形截面梁的长细比与欧拉临界力的关系，复合材料典型盒段的截面形式和铺设角，铺设次序对初始屈曲性态的影响。     

大变形薄板多体系统的动力学建模

基于非线性弹性理论,考虑剪切应变和横向应变,用绝对节点坐标法建立了大变形矩形薄板的动力学变分方程;为了提高非线性刚度阵的计算效率,根据非线性刚度阵与广义坐标阵的函数关系式,在非线性刚度阵中分离出广义坐标阵,从而避免了每个时间步长的单元刚度阵的积分运算。在此基础上,引入运动学约束关系,建立了大变形薄板系统第一类拉格朗日方程,对重力作用下大变形二连板进行数值仿真。计算结果表明:随着薄板的柔度增大,低频的弯曲变形与高频拉伸变形的耦合愈加显著;此外,系统机械能守恒验证了该模型正确性。     

弹塑性大变形分析的一致性高阶无单元伽辽金法

采用无单元伽辽金法求解弹塑性大变形问题。充分利用无单元法易于建立高阶近似函数的优点,位移采用二阶移动最小二乘近似。在更新拉格朗日方法的框架下,通过对控制方程弱形式的线性化建立了内力率的表达式,并区分为材料和几何两部分。采用Hughes-Winget算法更新应力,建立了Newton-Raphson迭代求解所需的一致切线刚度阵。刚度阵的数值积分采用近来针对小变形分析建立的二阶一致三点积分格式QC3(Quadratically Consistent 3-point integration scheme)。数值结果证明了本文方法分析弹塑性大变形问题的有效性和优越性。     

对实心弹性薄板优化的进一步讨论

回顾了实心弹性薄板优化研究的发展历程。给出了由无限密和无限细的肋骨加强的环板在常肋骨密度时的解析解；采用精确刚度阵的有限元法计算了具有有限根肋骨加强的板的柔顺性，它低于光滑优化解的柔顺性。这些结果进一步说明近期文献中出现的枕头形的光滑解并不是几何受约束实心弹性薄板优化问题的全局最优解。     

考虑局部切向约束接触问题的直接刚度法

针对工程中大量存在的切向滑移受到约束的接触问题,提出了基于Lagrange乘子的点-面及点-点接触直接刚度法,该方法生成的接触协调条件可以直接组装到结构刚度阵中,从而可直接用于考虑存在结构几何非线性及材料非线性的接触问题分析中。采用这一算法进行了考虑拱坝横缝张合效应的地震响应分析,结果表明考虑切向滑移约束对接触中的张开度有明显影响。     

计算具有区间参数结构的固有频率的优化方法

基于区间函数的单向包含性质，把具有区间非确定参数结构的固有频率所在区间范围问题 转化成两个全局优化问题，并采用一种实数编码遗传算法求取问题的全局解. 用一种能够求 得剪切型结构和弹簧质量系统特征值范围精确解的单调分析方法进行检验. 在一 些文献中，直接采用区间数运算法则和有限元法得到结构区间刚度阵和区间质量阵，并把关 于该区间刚度阵和区间质量阵的广义区间特征值问题的特征值区间作为待求的非确定性结构 的特征值所在的区间范围，该方法易于扩大问题的解域. 算例表明，可望得到结构 固有频率区间范围的准确解.     

拟协调SemiLoof和Loof扁壳元

1 前言目前工程上对任意壳结构分析常用的平板壳元和扁壳元存在着一个问题:如果一个结点周围的单元共面或近似共面,则对应垂直这个面的法向转角θ_N 的刚度是零或近似为零,这时若θ_N 没有被约束住,则结构刚度阵是奇异的,造成了求解困难.现有的一些解     

一种高阶非线性钢筋混凝土平板单元

在文献［１］高阶钢筋混凝土膜元的基础上，提出了一种具有１４个节点的高阶钢筋混凝土平板单元．钢筋混凝土材料模型仍采用Ｖｅｃｃｈｉｏ的抹平旋转裂缝模型，几何非线性仍采用所谓总体Ｌａｇｒａｎｇｅ列式法，非线性方程采用割线刚度位移增量迭代法．数值算例表明本文的方法是可靠的，高阶平板单元虽然列式复杂，但与低阶元相比，其计算量要少且精度要高