一、多步法的小时间步问题 二、回顾局部切断误差最小化 三、一个不收敛的实例 四、稳定+连续导出收敛 五、考察大时间步问题 六、两个时间刻度实例的绝对稳定性 七、振荡器

CFD the systematic application of computing systems and computational solution techniques to mathematical models formulated to describe and simulate fluid dynamic phenomena. Simulation is used by engineers and physicists to forecast or reconstruct the behaviour of an engineering product or physical situation under assumed or measured

结合北京地铁十号线一期光华路地铁站的工程和地质条件,运用FLAC3D软件对车站主体浅埋暗挖的五种不同施工顺序的方案进行模拟研究.在相同的地表荷载、固定约束和土体参数的情况下,模拟研究得出车站周围土体应力、塑性区和位移的基本变化规律,并确定群洞效应下开挖顺序的推荐方案.研究结果表明:(1)车站开挖过程中和完成后,破坏区主要集中在中洞与两侧洞之间的区域,即群洞之间部位,并从中洞两侧垂线以45°角向外扩散;(2)中洞的跨度对地表的沉降量有着明显的影响,对于群洞而言,先对左右洞室进行开挖可有效地减小开挖中洞的位移变化量;(3)中洞的拱形支撑能有效抑制地表的破坏区,考虑群洞效应时应先保证中洞洞顶部分的强度

为降低P110级石油套管淬火冷却过程中的内应力,提出\水淬—空冷—水淬\的优化冷却方式,并利用有限元方法对冷却过程中温度、应力场的变化规律和分布状态进行了模拟.模拟结果表明:冷却至7.5s出水时,横截面上最大温差为104℃,空冷结束时断面温度均匀;再次水冷的最大温差为80℃,与7.5s时相比,温差降低了24℃.对于应力,在最初的水冷阶段,从开始到2.5s,切向应力增大,2.5～5.5s,切向应力降低,冷却至5.5s时发生组织转变,此后热应力和组织应力共存,切向应力随冷却进行迅速升高,并在7.5s时达到最大,为563MPa;出水空冷阶段,热应力减小,组织应力消失,13s空冷结束时切向应力分布较均匀,为-11～27MPa;再次入水冷却至13.6s,切向应力再次达到最大,为451MPa,比7.5s时的563MPa降低了112MPa,达到了优化冷却工艺的目的

木建筑由于结构冗余性以及钉接节点超强的吸能耗能能力在抗震中表现良好.交叉层积木是一种新型的建筑材料.本文以交叉层积木3种柔性连接试验为基础,采用Open Sees中Pinching4自定义模型模拟连接滞回曲线的高度非线性、强度退化、刚度退化和捏拢现象.基于主次半循环累积能量的损伤模型,对交叉层积木连接进行损伤分析,并提出该连接的5种性能水平的损伤指数.Pinching4模型与连接试验结果吻合较好,进一步证明该模型模拟木节点连接性能的可行性和有效性.损伤因子对应的损伤程度基本符合试验规律,其平均值在合理范围内,计算结果离散性较低

在GTN模型基础上,考虑到微孔洞剪切变形对材料劣化的影响,建立适用于压应力状态的剪切修正模型.通过用户子程序接口VUMAT将与损伤耦合的弹塑性本构模型嵌入具有ALE法的有限元软件Abaqus/Explicit中.利用模拟拉伸与纯剪切试验拟合载荷-位移曲线以确定模型参数.将修正模型应用到辊冲工艺有限元模拟中预测断面质量,并进行试验验证.结果表明:前刃口为小间隙冲裁,塌角较小,光亮带较大但带有一定的倾角;后刃口为大间隙冲裁,塌角与断裂带较大,光亮带较小;裂纹会同时在前刃口凸模与凹模侧面萌发,而对于后刃口,会首先在凹模侧面产生

以直径20 mm, 包覆比50%的银包铝细棒为研究对象, 通过有限元数值模拟以及相应的实验验证, 得出了银包铝复合材料立式连铸复合成形工艺的边界条件.采用ProCAST软件模拟了立式连铸成形过程, 得出各工艺参数对连铸结果的影响规律, 给出了可行的连铸工艺参数范围及工艺调控策略, 以模拟结果为指导, 制备出表面质量高、复合界面效果良好的银包铝复合棒材.实验结果表明, 芯管长度、连铸速度对结果的影响最大, 芯管长度影响了芯管出口处双金属的接触温度、接触时间, 并直接改变了铝芯固液界面的相对位置.当芯管长度过短时, 银铝界面反应较强烈, 当芯管长度过长时, 芯棒冷却强度大, 芯部铝产生明显的冷隔.随着连铸速度的增大, 银的固液界面到芯管出口距离逐渐减小, 铝的固液界面距出口距离逐渐增大; 铝液铸造温度升高, 冷却水减少也会带来相似的作用.结果显示, 芯管长度30 mm, 速度37~67 mm·min-1, 银的铸造温度1225~1325℃, 铝的铸造温度800℃, 冷却水流量约300 L·h-1是可行的银包铝连铸工艺

利用Pro CAST软件对2400 mm×400 mm宽厚板坯结晶器建立三维动态模型,采用移动边界法实现结晶器内流场、温度场及应力场的耦合模拟.结果表明:考虑凝固坯壳的影响,下回流区位置向铸坯中心靠拢,真实反映了钢液在连铸结晶器内的流动情况.自由液面的钢液从窄面流向水口,速度先增大后减小,距水口约0.7 m处,出现最大表面流速,约为0.21 m·s-1.结晶器出口坯壳窄面中心厚度最小且由中心向两侧逐渐增大,最小厚度约为10.4 mm;受流股冲击影响较弱的宽面坯壳与窄面相比生长更均匀,宽面偏角部和中心的坯壳厚度分别为18.9 mm和27.6 mm.铸坯坯壳应力变化趋势与温度基本保持一致,表明初凝坯壳应力主要是热应力.结晶器内铸坯宽窄面上的等效应力均沿着结晶器高度下降方向呈增大趋势,铸坯角部、宽面中心及窄面中心位置的最大应力各约为200、100和25 MPa

对氢致钢内部疲劳裂纹的萌生和扩展进行了数值模拟.首先用有限元法分析了氢在疲劳载荷作用下向钢中缺陷处扩散富集的过程,然后计算得到氢含量分布结果.根据夹杂理论将氢富集区视为在缺陷附近分布的弹性夹杂,用有限元法计算得到的氢含量场求出夹杂处的应力强度因子,进而建立疲劳裂纹萌生和扩展的判据.比较了在不同加载条件下氢致疲劳裂纹萌生和扩展的规律.用梯形法修正了Sofronis和McMeeking的瞬态扩散有限元公式,发现用梯形法可以缓解加载初期较高的浓度梯度和应力梯度引起的计算结果震荡的情况,这对于计算开裂判据是十分重要的.最后讨论了提高模拟精度和改进模型的方法

通过理论分析和数值模拟的方法,研究了1873 K钢液中球状夹杂物颗粒运动过程中的受力情况.静止钢液和湍流场中,夹杂物颗粒所受的压力梯度力、虚拟质量力、Saffman力及Magnus力较小,可以忽略不计.静止钢液和均匀湍流场中,随着夹杂物颗粒尺寸的变大,颗粒所受Brown力逐渐衰减,当夹杂物颗粒直径大于20μm时颗粒所受Brown力基本不会对其运动造成影响,可以忽略不计.小于10μm的夹杂物颗粒布朗运动较为剧烈,单位质量夹杂物颗粒所受Basset力和Brown力较大,成为主导夹杂物颗粒运动的主要作用力;尺寸较大的夹杂物颗粒(直径D>50μm),Basset力、质量力和Stokes力共同作用影响夹杂物颗粒的运动、碰撞和去除