扩展有限元与遗传算法相结合的结构缺陷反演分析

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

一种XFEM断裂分析的裂尖单元新型改进函数

提出了一种适用于裂尖改进单元的新型改进函数, 基于三角变换的方法, 保留裂纹尖端场的应力奇异性和裂纹上、下表面的位移不连续性, 将常规扩展有限元法裂尖改进单元的4 项改进函数缩减为2 项, 裂尖改进单元的结点由常规的8 个改进自由度减少为4 个. 采用2 个正交的水平集函数表征材料内部裂纹面, 详细阐述了改进单元类型的判别方法, 给出一种改进单元的分区域积分方案. 最后, 若干断裂力学问题经典算例的数值计算结果表明:建议的裂尖改进函数具有较高的数值精度, 该方法是十分有效的.     

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

压电材料动态裂纹问题的奇异积分方程法

利用广义Betti-Rayleigh 互易公式给出了二维压电材料非渗透裂纹问题的一般解和奇异积分方程,其中未知函数为裂纹上的位移间断和电势间断的导数. 在理论分析的基础上,使用高斯-切比雪夫求积公式及Lubich 卷积积分方法建立了问题的数值求解方法,并给出典型算例的广义动应力强度因子随时间变化的规律.     

含缺陷结构的塑性极限分析

结合极限分析中的数学规划理论和有限元技术,提出了三维含缺陷结构极限分析的数学规划方法,并采用罚函数法引入塑性不可压条件.对于考虑多组独立变化载荷联合作用的情况,提出了加载路径射线辐射求解方案,并基于这种射线辐射状的加载路径,推导了多组载荷联合作用下结构塑性极限上限分析的数崐学规划格式,编制了相应的有限元程序.文中的数值结果表明了该方法的正确性与有效性.     

ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD 1)

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.     

NUMERICAL SIMULATION OF CRACK GROWTH IN BRITTLE MATRIX OF PARTICLE REINFORCED COMPOSITES USING THE XFEM TECHNIQUE

Crack growth in particle reinforced composites is significantly influenced by the character of the reinforcement particles.To accurately deal with such problems,the extended finite element method (XFEM...     

用于板料成形反向模拟的特殊应力更新算法

提出了一种改进的反向模拟法,以最终构型为研究对象,采用Euler坐标系,基于虚功原理获得有限元列式. 改进的反向模拟法采用了一种基于塑性流动理论的本构方程,可以充分考虑应变历史对塑性变形的影响. 为了避免流动理论应力更新算法过程中关于未知量\Delta\lambda 的非线性方程的求解,引入等效应力思想,无需Newton-Raphson迭代直接计算未知量\Delta \lambda . 盒形件的拉深实例中,传统的基于塑性形变本构方程的反向模拟法和改进的基于塑性流动本构方程的反向模拟法计算结果,分别与基于增量有限元法的正向数值模拟求解器LS-DYNA计算结果进行对比. 通过获得的坯料轮廓、成形极限图、等效应变分布、计算效率等的比较,验证了所提出的基于塑性流动理论本构模型的应力更新算法的有效性.      

功能梯度材料结构的热应力边界元分析

导出了一种新的可对功能梯度材料结构进行二维和三维热应力分析的积分方程, 利用该方程并结合多区域边界元三步求解技术, 可对由任意多种介质组成的变物性参数复合结构进行热应力分析; 采用与弹性模量无关的Kelvin解作为问题的基本解, 导出的积分方程含有由材料的非均质性以及温度变化引起的域积分; 使用径向积分法将所有的域积分转换成等价的边界积分, 从而建立起只需要边界离散的无内部网格边界元算法; 最后对两种典型的飞行器防热结构进行二维和三维热应力算例分析, 并通过与有限元计算结构对比验证结果的正确性.      

关于动力分析精细积分算法精度的讨论

对动力问题分析的精细积分算法的精度问题进行深入研究,并在此基础上提出对原有的算法的改进策略,改进后的算法可以较好地克服算法精度对积分时间步长的依赖性问题。     

空间杆系结构非线性动力稳定性分析的一种实用算法

本文对任意空间杆系结构进行非线性动力稳定性分析。     

结构地震反应分析的逐步微分积分方法

将结构地震时程反应分析问题离散化为一系列初值问题的求解, 在每个离散时步内应用微分积分(DQ)法则, 建立了一种结构地震反应的逐步DQ数值求解方法. 根据地震地面运动及结构体系的动态特征和物理本质, 假定离散时步内地面加速度为线性分布. 推导并给出了该分布模式假定下利用DQ法逐步求解体系地震反应的数值格式, 通过一个双自由度体系的数值实验阐述了结构地震反应DQ分析方法的应用. 算例表明, 体系地震反应DQ解的稳定性和精度比较好, 在较大的离散步距条件下仍然可以获得较好的计算结果.      

轴对称组合弹性体初应力分析的边界元传递矩阵法

基于边界元与传递矩阵法耦合基础上，提出了一种分析轴对称组合弹性体初应力的新方法。利用边界元和传递矩阵进行分域计算轴对称组合弹性体的初应力，无需在整个结构上进行矩阵组装，只需在轴对称子午面的边界上进行离散。故该方法具有输入数据少，计算精度高及所需计算机内存小等优点，适合在微机上求解大型复杂轴对称组合结构问题。文中给出了一个组合厚壁圆筒的计算实例，分别与有限元解和理论解进行比较；结果表明：本文所提出的     

DYNAMIC ANALYSIS OF BOUNDED DOMAINS BY SBFE AND THE IMPROVED CONTINUED-FRACTION EXPANSION

基于比例边界有限元理论框架,通过采用连分式展开和引入辅助变量,将有限域的动力刚度矩阵和质量矩阵采用高阶的矩阵表示. 采用改进的连分式法求解比例边界有限元方程中的动力刚度矩阵. 通过增加连分式展开的阶数,该求解方法能包含动力分析的主要频率范围. 针对结构自由度较多的系统当连分式阶数逐渐增大时,原连分式算法可能会造成矩阵运算病态的问题,提出采用改进的连分式算法能有效地提高数值计算稳定性.通过对一正八边形的自由振动分析及矩形平面的时域分析,算例结果表明改进算法的鲁棒性更强,适合大规模系统的动力分析.     

用不同模量有限元分析坝体应力和变形

探讨拉,压不同模量有限元法在重力坝应力,变形分析中的应用,利用该方法计算坝体裂缝,接触,稳定问题明显地提高了分析精度,本指出,给定拉伸时的弹性模量和泊松比是重要的可能的。在岩土力学,冰力学,陶瓷和塑料的力学分析中,使不同模量有限元法必定会有实际意义。     

任意的拉格朗日欧拉边界元－有限元混合法分析物体撞水响应

计算物体的撞水响应目前已有了一些专用的算法．本文在分析和比较这些算法的基础上，提出了一个解撞水问题的任意的拉格朗日欧拉边界元－有限元混合法（ＡＬＥ－ＢＥ－ＦＥＭ），这个方法不仅充分发挥了边界元法计算半空间流场的优越性，而且还能计及液面大晃动的非线性边界条件和物体变形所造成的影响．文中给出圆柱刚体和楔形刚柱体两个撞水算例，结果有力表明该方法的可靠性和有效性。     

内压和自重作用下浅埋管道围岩的应力解1)

对于浅埋管道而言,土体自重对其围岩应力场影响不可忽略。基于Verruijt提供浅埋洞室问题求解基本方法,同时考虑自重和内压作用,研究了半无限域中圆形孔洞的围岩应力问题。利用复变函数、保角映射和傅里叶变换等方法,推导了考虑围岩自重和内压作用下半无限域中孔洞应力的解析解,分析了内压和孔洞大小及径深比对考虑自重作用下的围岩应力的影响,为管道设计和施工安全提供了理论依据。     

薄板统计分析的随机边界元法

本文建立了分析含随机材料参数的薄板弯曲问题的随机边界元法。基于Ｔａｙｌｏｒ级数展开技术，分别得到了位移的均值和一阶偏差的边界积分方程，发现材料参数的随机性可作为一个等效的随机荷载处理，从而得到边界位移或边界力的均值和协方差，并进一步求出内点位移和力矩的均值和协方差，最后用本文方法计算了两个算例，并对结果进行了必要的分析。     

用无网格局部径向点插值法分析非均质中厚板的弯曲问题

用无网格局部径向点插值法分析了非均质中厚板的弯曲问题.利用虚位移原理推导了中厚板的离散系统方程.采用径向基函数耦合多项式基函数来近似试函数,用四次样条函数作为加权残值公式中的权函数.所构造成的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件.此方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化.算例结果表明这种无网格方法具有效率高、精度高和易于实现等优点.