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

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

有限元显式算法解大变形冲击响应中的网格稳定性研究

1 引言结构在冲击载荷作用下的瞬态响应问题是许多工程问题提出的,如结构耐撞性、防护工程等.这类问题的数值解法已成为计算结构力学中的一个研究热点,有限元显式算法,     

金属刚粘塑性变形的能量泛函与动力显式算法有限元分析

设计出刚粘塑性动态变形场的能量泛函，采用罚函数法和Ｌａｇｒａｎｇｉａｎ乘子法取消对运动容许速度场满足体积不可压缩的约束条件，运用虚功原理和广义变分原理，推导出刚粘塑性变形的动力分析显式算法有限元方程，以及速率形式的中心差分求解的时间积分显式，给出了两个金属变形的算例。     

基于显式有限元方法的两相介质弹塑性动力反应计算分析

针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的弹塑性动力反应进行计算和分析,将计算结果与相应的弹性动力反应的计算结果进行对比;对本文应用的弹塑性...     

弹塑性损伤本构关系的显式积分算法研究

首先引入弹塑性损伤本构关系,分别从材料软化与残余应变两个方面,描述伪脆性材料的非线性行为.针对结构动力分析中的强非线性问题,给出了弹塑性损伤本构关系的显式积分算法.算法中引入算子分解的思想,将弹塑性本构关系分成塑性与损伤两个模块.首先求解塑性模块,根据有效应空间塑性演化公式,采用前进欧拉算法,直接构造塑性演化的预测值,并且根据屈服函数的漂移构造了误差限公式,作为衡量显式算法精度的指标.将塑性模块求解的结果代入损伤模块,可以方便地求得损伤变量的演化,并最终得到更新后的应力.整个求解过程不需要迭代,可最大程度的算法稳定性.将论文建立的本构关系显式算法与结构分析显式算法结合,构造了结构显式分析方法,并模拟了两个经典算例,算例结果验证了论文方法的有效性.     

显式方法精确模拟形状记忆合金循环荷载下的变形行为

本文提出全新的有限弹塑性J2流方程,用来显式、精确地模拟SMAs（形状记忆合金）材料在循环加载-卸载条件下从塑性逐渐转变为伪弹性的变形行为．首先,改进流动法,使得本构方程耦合屈服中心的移动和屈服面的增大,并改进背应力演化方程,使模型可以产生强烈的包辛格效应,从理论上具备模拟SMAs独特变形行为的能力;其次,构造全过程下的统一硬化函数显式表达式,代入本构方程后能得到符合要求的形函数;再次,利用选定的数据点构造统一光滑的上屈服函数,再利用上下屈服应力之间的一种线性关系,推导得到下屈服阶段的形函数;最后,只需要给定一个参数就可以得到单个循环结果,利用拉格朗日插值方法构建参数随循环次数变化的函数,就可以模拟任意循环荷载下的变形行为．通过模型结果和实验数据对比证明新方法的有效性．     

结构动力响应分析的三阶显隐式时程积分方法

基于泰勒级数展开式提出了一种用于结构动力响应分析的高精度时程积分方法,该方法假设t时刻的速度和加速度由t-Δt时刻、t时刻、t+Δt时刻的速度和加速度加权表示,并可根据求解需要调节权值,将积分算法构造成隐式格式或显式格式。通过理论分析和数值算例,计算讨论了该算法的稳定性和精度,确定了最佳的权值和允许的时间步长。结果表明:本文算法最高具有三阶精度,且具有振幅衰减率低、周期延长率极小等优点。最后结合一个铁道工程实例,表明本文算法适用于大型非线性动态响应的精确快速求解。     

李级数算法和显式辛算法的相位分析

以线性可分Hamilton动力学系统为例,研究了李级数算法和显式辛算法的相位精度,研究了李级数算法的保辛精度及其保辛精度的提高方法;指出了显式辛算法相位精度与算法阶次的不协调性,印辛算法的阶次高并不意味着其相位精度也高,李级数算法不存在这种问题,指出了一个算法的相位可能超前也可能滞后.分析结果表明三阶显式辛算法具有比较高的相位精度.     

一类带高频耗散的动力学显式算法

带数值耗散的算法因其能有效过滤虚假高频响应的影响而备受关注。基于离散控制理论Z变换,提出一类带数值耗散的结构动力学显式算法。该算法采用CR法的速度和位移递推式,满足零振幅衰减率,且对线性系统和非线性刚度软化系统为无条件稳定,对刚度硬化系统则是条件稳定的;该算法由单个参数ρ控制数值耗散能力。通过对振幅衰减和周期延长的理论分析表明,系数a可调节算法的精度和非线性稳定区间,给出精度最优时系数的取值。对特定的系数取值该算法可转变为CR法。通过算例对线性系统和非线性系统的分析验证了新算法具有良好的精度、稳定性和数值耗散,表明新算法是正确有效的。     

动、静摩擦对可接触型裂纹动态影响的一种算法

本文提出了一种分析考虑裂纹面摩擦影响的可接触型裂纹动力反应的时域显式计算方法，这一方法不但可以较好描述裂纹面之间碰撞过程，也能反映当裂纹面发生接触－滑移时，裂纹面间动、静摩擦力的影响．文中给出了使用这一方法的基本步骤，计算了无限大裂纹的动接触问题和直线型裂纹对入射Ｐ波的散射问题，初步验证了本文提出的方法．     

圆柱壳撞水时的弹塑性动力屈曲研究

本文研究圆柱壳在流固中冲击载荷下的弹塑性动力曲问题。建立了液体－气体－固体三相的数学模型。其中结构部分控制方程由弹塑性力学中关于加速度的最小原理获得，本构关系采用增量理论，液体假设为不可压缩，空气层认为最等熵压缩。分别讨论了不同冲击高度时压力变化规律，屈曲对载荷的影响。屈曲沿壳长的分布及发展规律等。     

一种基于SHTB的Ⅱ型动态断裂实验技术

冲击剪切载荷作用下动态断裂韧性的测定是材料力学性能和断裂行为研究中重要组成部分.为了测定材料的Ⅱ型动态断裂韧性,许多学者采用不同的试样与实验方法进行了实验,但限于实验条件,裂纹断裂模式往往是I+Ⅱ复合型,而不是纯Ⅱ型,因而不能准确测得材料的Ⅱ型动态断裂韧性.鉴于此,本文基于分离式霍普金森拉杆(split Hopkinson tension bar,SHTB)实验技术,提出一种改进的紧凑拉伸剪切(modified compact tension shear,MCTS)试样,通过夹具对MCTS试样施加约束,从而保证试样按照纯Ⅱ型模式断裂.采用实验-数值方法对MCTS试样动态加载过程进行分析,将实验测得的波形输入有限元软件ANSYS-LSDYNA,得到了裂纹尖端应力强度因子-时间曲线,并与紧凑拉伸剪切(compact tension shear,CTS)试样进行了对比.同时采用数字图像相关法进行了实验,验证了有限元分析结果.结果表明,MCTS试样在整个加载过程中K_I K_Ⅱ,裂纹没有张开;而CTS试样在同样的加载过程中K_IK_Ⅱ,出现裂纹张开现象.这说明MCTS试样能够准确地测定材料的Ⅱ型动态断裂韧性,为材料动态力学测试提供了一种有效的实验技术.     

由Gurtin变分原理求解一维动力响应的半解析法

基于Gurtin变分原理本文提出了求解动力响应问题的半解析法.该方法在空间域内进行有限元离散,在时间域上取级数。实例表明,该方法是求解动力响应问题的一种非常有效的方法。本文引入了一种节点位移函数,使所导出的位移模式可以求解各种初始条件和载荷的一维动力响应问题。     

随机有限元动力分析方法的研究进展

对随机结构动力分析方法,包括随机模拟法、摄动随机有限元法、动态随机有限元法和正交展开法等进行了评述,介绍了该领域内的研究进展,并讨论了研究中存在的问题以及发展方向．      

基于遗传算法的一种Tikhonov正则化改进方法

就Tikhonov正则化方法求解第一类算子方程进行了新的探讨,针对展平泛函的极值处理,采用了遗传算法来解决.结合遗传算法的优点,扩大了解决问题的范围.介绍了其基本原理及运算流程,作为此方法的应用,首先对一变截面悬臂梁模型进行截面积反求,另外解决了一个热传导反问题实例,证明了其方法的有效性.     

圆弧足被动行走器非光滑动力学仿真研究

传统的被动行走器动力学研究方法是将被动 行走器模型视为双足轮流铰接于地面(或纯滚动)的倒立摆模型, 并应用角动量守恒假设计 算支撑腿转换时足地碰撞后的系统运动状态. 但该方法不能分析足地之间的滑动现 象、足地之间的非完全塑性碰撞现象和两足同时着地的情况. 采用基于线性互补 的非光滑动力学方法, 给出了一种基于中点格式积分一步的时间步进算法, 可以有效地分析 足地间的摩擦和碰撞现象. 应用该方法对被动行走器进行数值仿真, 发现了两种新的周期性 步态, 一种是带摩擦滑动现象的步态, 另一种是支撑足转换时发生多次碰撞现象的步态. 数 值仿真结果还表明, 较大的摩擦系数和较小的法向恢复系数更有利于被动行走器实现稳定周 期步态行走.     

用“膜力因子法”分析简支刚塑性圆板的大挠度动力响应

在板的变形过程中,当其挠度与厚度同量级时,膜力的作用与弯矩的作用同样重要;当其挠度超过厚度甚至接近板的面内尺寸的量级时,相对于弯矩的作用,膜力的作用是更主要的。本文发展了一种基于能量平衡的简便的“膜力因子法”,有效地解决了当板的挠度达到厚度量级时弯矩和膜力联合作用问题;用这种方法并结合塑性膜阶段的分析成功地求解了简支刚塑性圆板在冲击载荷作用下的大挠度塑性响应问题。所得结果与Florence实验结果符合良好,比已有理论近似结果有所改进。     

径向基点插值法在旋转柔性梁动力学中的应用

将无网格径向基点插值法用于旋转柔性梁的动力学分析. 利用无网格方法对柔性梁的变形场进行离散,考虑梁的纵向拉伸变形和横向弯曲变形,并计入横向弯曲变形引起的纵向缩短,即非线性耦合项,运用第二类拉格朗日方程推导得到系统刚柔耦合动力学方程. 将无网格径向基点插值法的仿真结果有限元法和假设模态法进行比较分析,说明假设模态法的局限性,并表明其作为一种柔性体离散方法在刚柔耦合多体系统动力学的研究中具有可推广性,并讨论了径向基形状参数的影响. 同时运用3 种求解系统动力学方程的方法:纽马克方法、4阶龙格库塔法、亚当姆斯预报校正法,并比较各方法的计算效率, 结果表明纽马克方法最快.      

基于独立分量分析的多源动态载荷识别方法

提出了基于独立分量分析的多源动态载荷识别方法, 解决了在结构系统未知的情况下载荷波形的识别问题. 该方法基于结构在多源动态载荷作用下, 其响应是载荷与对应的结构脉冲响应卷积的原理, 并假设载荷源相互统计独立. 与既有的载荷识别方法相比,该识别方法特点表现在: 结构质量, 刚度及阻尼等信息可以完全未知, 但以实际载荷间的独立性为优化目标; 用互信息来度量识别载荷间的独立性, 通过梯度下降算法取消识别载荷间的各阶相关性, 使识别载荷间基本满足相互独立; 从波形的角度来进行载荷识别.通过数值仿真表明: 该方法对测点, 噪声, 不同载荷形式及不同结构有较好的鲁棒性; 识别载荷与实际载荷在归一化条件下, 识别载荷与实际载荷相关性系数约为1.