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

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

组合结构总体节点参数自动标定的基本定理

在有限元位移法中,单元矩阵(刚度矩阵、几何矩阵、一致质量矩阵、载荷列阵等)往往是在单元局部坐标下计算的(以下简称局部方位);结构整体矩阵的形成,单元应力和变形的计算,都要求确定总体节点参数(自由度)和单 ...     

海洋平台导管架结构分析中管节点局部柔度的等效单元

本文建立了一种在海洋平台导管架有限元分析中代表管节点局部柔度的等效单元,可有效地反映管节点在轴力和弯矩作用下的局部变形性质,并导出了T型、Y型以及K型管节点的局部柔度等效单元的刚度矩阵。本文建立的等效单元可方便地与通用结构分析有限元程序(如SAP5)配合,对海洋平台导管架进行计及管节点局部柔度的精确静动力分析。作者编制了与SAP5配合使用的前处理程序,可自动生成符合SAP5输入格式的等效单元读入刚度矩阵的数据。文中给出了说明等效单元应用方法的简单算例。与有限元法计算的结果比较表明,用等效单元模拟管节点的局部柔度是必要的和准确可行的。     

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

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

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

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

局部有限元模型边界条件的直接传递方法

提出一种将整体分析得到的节点力或节点位移直接传递到精细化局部有限元模型的方法,即部分混合单元法。沿精细化局部有限元模型周边建立一组过渡单元,该组过渡单元采用与整体模型一致的单元类型和模拟方式,其外侧边界上的节点与整体模型节点的相对坐标对应,内侧边界与精细化局部有限元模型采用基于面约束的方式连接。在外侧边界上根据节点坐标对应施加整体分析获得的节点力或节点位移,过渡单元就可直接将边界条件传递到精细化局部有限元模型。通过贵州红水河特大桥钢-混结合段的精细化有限元分析,验证了本文方法的实用性和有效性。     

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

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

焊接管结构的非线性推覆简化模型

有限元分析方法中,采用三维实体单元对焊接管结构进行离散是最为精确的建模方法,但对三维实体模型分析时会划分大量单元,导致计算时间冗长。工程中常把焊接管结构简化为刚架模型,即刚性连接的杆系系统。这种建模方式忽略了节点的局部柔度,过高地估计了整体结构的刚度,使计算结果偏于危险。为解决上述问题,在节点处引入虚拟梁单元,使其等效模拟节点的局部变形,还原支管真实长度。为了验证节点柔度对焊接管结构的影响,对典型的焊接管结构缩尺,进行了推覆试验。分别采用壳单元模型、刚架模型及引入虚拟梁单元的简化模型进行有限元模拟,并与试验结果对比分析。结果表明,刚架模型会过高的估计焊接管结构的刚度,而简化模型计算结果与壳模型、试验结果较为吻合。     

一种提高四边形四节点平面壳单元计算精度的新方法

平面壳单元是由平面应力单元和平板弯曲单元叠加组合而成,具有简单的理论表达,但是它在计算曲面壳体结构时误差较大。为了进一步提高平面壳单元的计算精度,本文提出了一种计算平面壳单元刚度矩阵的新方法。通过该方法在高斯积分点建立多个单元局部坐标系,并保证每个局部坐标系都位于单元在高斯点处的切平面上,从而可以有效适应曲面壳体形状,达到进一步提高平面壳单元计算精度的目的。为了在这种新坐标系下计算单元刚度矩阵,给出了求解形函数对局部坐标的导数、局部到自然坐标系积分转换的雅可比、以及局部到整体坐标系的转换矩阵的新型计算方法。通过将这些新坐标系以及新计算方法运用到平面壳单元DKQ24中,可以有效提高平面壳单元尤其是在计算曲面壳体时的精度。计算结果表明,本文方法和平面壳单元相结合,不仅具有平面壳单元简单的理论表达式,还能得到满意的精度。另外,本文方法还可以应用到其他类型的平面壳单元,为提高其他类型平面壳单元的计算精度提供了一种新的途径和思路,具有广阔的应用前景。     

应变局部化分析的嵌入强间断多尺度有限元法

固体材料的应变局部化行为是导致结构破坏失效的重要因素之一,开展相关数值模拟分析对于结构安全性评估具有重要意义.然而由于材料的非均质和多尺度特性,采用传统数值方法进行求解时通常需要从最小特征尺度离散求解的结构,这将大幅度增加计算规模和成本.针对这一问题,本文提出了一种基于嵌入强间断模型的多尺度有限元方法.该方法从粗细两个尺度离散求解模型,首先在细尺度单元上引入嵌入强间断模型来描述单元间断特性,所附加的跳跃位移自由度则通过凝聚技术进行消除,从而保持细尺度单元刚度阵维度不变.其次,提出了一种增强多节点粗单元技术,其可根据局部化带与粗单元边界相交情况自适应动态地增加粗节点,新构造的增强数值基函数可以捕捉细尺度间断特性,完成物理信息从细单元到粗单元的准确传递以及宏观响应的快速分析;再次,在细尺度解的计算中,将细尺度解分解为降尺度解与单胞局部摄动解,从而消除弹塑性分析时单胞内部的不平衡力.最后,通过两个典型算例分析,并与完全采用细单元的嵌入有限元结果进行对比,验证了所提出算法的正确性与有效性.     

结构动力分析自适应有限元方法综述

结构动力分析自适应有限元方法主要研究有限元动力分析的误差估计理论,建立适用于复杂结构动力分析的有限元网格自适应过程．介绍了结构动力问题自适应有限元方法的重要发展,包括固有振动和动响应分析的误差估计及相应的自适应策略;且简要介绍了几种现有的网格生成技术及其特点．最后指出这种方法存在的问题和今后的研究方向．     

THE INVERSE ANALYSIS OF INTERNAL DEFECTS (INCLUSIONS) IN STRUCTURES USING A COMBINED XFEM AND GA METHOD

建立扩展有限元法与遗传算法相结合的结构缺陷反演分析模型.扩展有限元法通过引入不连续位移模式使得网格剖分无需依赖结构内部的不连续界面.通过改变水平集函数表征结构缺陷（夹杂）的位置和大小，遗传算法在每次迭代过程中具有全局和局部搜索能力，通过评估响应测点的响应量适应度值决定是否进一步迭代.对带有单个圆形缺陷（夹杂）和多个缺陷（夹杂）的结构进行了反演分析，并就响应测点的布置进行了讨论.结果表明，建立的反演分析模型能准确地探测结构存在的单个甚至多个缺陷（夹杂）.     

空间几何构造分析的有限单元法

提出空间杆系几何构造分析的有限单元法，构造了两种单元(链杆单元和准梁单元)的几何约束矩阵，集成为整体矩阵并引入承条件后，通过对其阶数与秩的比较分析确定体系的几何可变性及静定性．本法原理简单，便于计算机实施，结果完备：对于几何不变体系，可指出多余约束的数目；对于几何可变体系，可给出体系的自由度数及相应的运动模态，并确定自由度的常变瞬变性质．     

基于附加虚拟质量的结构损伤识别方法

针对大型土木结构的模态信息往往对局部损伤灵敏度较低的特点,提出通过在各个子结构上附加虚拟质量的方法提高局部灵敏度,实现整体结构的准确损伤识别。该方法无需在结构上布置真实质量,它首先利用虚拟变形法（VDM）可进行结构快速重分析的思想,由实测结构的激励和加速度响应,构造结构附加虚拟质量后虚拟结构的频率响应;然后结合灵敏度分析和附加质量与频率关系,确定所需附加质量和对应具有较高灵敏度的频率;最后分别在每个子结构上附加虚拟质量,联合所有虚拟结构和对应的频率即可准确快速地识别出整体结构各个子结构的损伤。本文通过两层平面框架有限元模型验证了附加质量损伤识别方法的有效性。     

高超声速舵面热气动弹性不确定性及全局灵敏度分析

鉴于高超声速中气动热预测的不确定性影响热气动弹性分析的可靠性,提出一种温度分布参数化模型,基于此模型,对高超声速舵面热气动弹性中气动热的不确定性及全局灵敏度进行分析,分析方法:求解N-S方程得到物面的温度分布,对此温度分布进行参数化,分别采用蒙特卡罗模拟(Monte Carlo simulation,MCS)方法和稀疏网格数值积分(spare grid numerical integration,SGNI)方法生成不确定性及全局灵敏度分析所需样本,对所有样本都进行热气动弹性分析,热气动弹性分析过程为:由样本得到温度分布,基于此温度分布,考虑热应力和材料属性的影响,对结构进行模态分析,将结构模态插值到气动网格,采用基于CFD的当地流活塞理论进行了气动弹性分析.分别在两种飞行状态下进行分析,计算结果表明:(1) M=5,H=15 km,结构固有频率和颤振分析结果的变异系数约为5.83%;(2) M=6,H=15 km,结构和颤振分析结果的变异系数约为8.84%.两种状态下,两个不确定参数的全局灵敏度都在50%左右,两者耦合作用很小,约为0.与MCS方法相比,SGNI方法显著的提高了不确定性分析效率.      

基于变形修正的局部刚体化动力模型简化方法

提出一种基于变形修正的动力模型简化方法. 以复杂工程结构的有限元模型为基础,根 据运动同步性假设将结构分成若干同步性区域. 采用适当的位移模式去逼近位移真值,再通 过总体应变能最小来优化,从而得到整体位移与局部位移模式间的转换关系. 算例结果 表明了简化方法的有效性.     

结构可靠度分析FORM迭代算法的混沌控制

利用混沌控制原理对FORM收敛失败进行控制. 理清了全局性和局部性两类混沌反馈 控制各种方法的内在联系，说明稳定转换法和自适应调节法属于全局混沌反馈控制 方法，自适应调节法可视为稳定转换法的特例. 参 数调节混合法不过是松弛牛顿法的另一种表达形式，它们都属于局部混沌反馈控制方法. 阐 明了混沌反馈控制表达式与工程力学收敛控制迭代算法的对应关系. 也揭示了这些迭代算法 收敛控制措施的功效和局限性. 提出了一个以稳定转换法为主联合松弛牛顿法的混 沌反馈控制方法，对可靠度分析FORM迭代算法实现了周期振荡、分岔和混沌控制.     

沉管隧道横向动力响应分析的多体动力学方法

沉管隧道横向地震响应分析是沉管隧道结构抗震设计的重点.根据沉管隧道的结构特点和接头剪力键的节点构造,建立了沉管隧道的多刚体-剪切铰-黏性阻尼铰多体动力学模型,并应用多体动力学中的离散时间传递矩阵法(multibody system discrete time transfer matriX method,MS—DT-TMM)建立了数学模型及表达式,同时编写Matlab程序进行求解实现.设定算例分别采用MS—DT—TMM和传统有限元法(finite element method,FEM)对沉管隧道横向地震响应进行对比分析,两者吻合较好,表明MS—DT-TMM应用于沉管隧道横向动力响应分析的可行性和有效性.     

能量有限元预示复合材料结构动响应研究进展

复合材料结构高频动响应预示是飞行器等结构设计中的重要研究内容之一.为了探讨精确预示复合材料结构高频动响应方法,分析比较了目前较为通用的3种动力学响应预示方法,指出能量有限元法最适合求解具有各向异性特征的复合材料结构高频动响应问题.紧接着概述了国内外关于能量有限元方法和该方法在复合材料结构高频动响应预示方面的研究进展.在此基础上分析了能量有限元法在预示复合材料结构高频动响应问题中尚待深入研究的问题.