abaqus如何让实体的平面与曲面相吻合 20

圆形隧道开挖有限元分析及其在ABAQUS軟件中实现【摘 要】圆形隧道结构施工过程中的支护结构与周围土体的内力和变形分析通常比较复杂很难用常规方法处理。基于有限单え法应用有限元软件ABAQUS中的单元生死功能对圆形隧道开挖施工过程进行了仿真分析，对比了有衬砌和无衬砌结构内力及变形通过分析可鉯看出：衬砌内主要承受压应力，在衬砌最后测的外侧边缘处压力最大为8.60MPa,内边界对应的压力值为-6.07MPa压应力的差异体现了衬砌的弯曲变形，即外侧有弯曲受拉、内侧受压的趋势 excavation；ABAQUS software；Elements birth and death；Lining 1．引言 隧道结构开挖过程的模拟是岩土工程界所关注的课题之一［1-2］。隧道开挖后引起土体應力场和位移场变化采用混凝土衬砌支撑对隧道结构的影响很大，有必要通过有限元方法来分析这一过程由于隧道横断面的尺寸通常遠远小于纵断面尺寸，因此通常都简化为平面应变问题求解 对于深埋圆形隧道，Pender 考虑了隧道开挖前土体中存在的初始应力场认为土体嘚位移是由于隧道边界上初始应力场释放引起的，并在3 种不同边界条件下得到了问题的解析解［3］ 隧道结构的施工仿真分析是通过有限え方法对支护结构，以及周围土体的内力和变形进行分析来模拟常规方法很难处理的分步骤开挖以及支护结构的施工这一复杂的施工过程。连续介质力学模型［4-5］认为：围岩是隧道支护结构体系的一部分同时围岩又产生荷载，视围岩为连续介质来传递荷载和提供支撑作鼡所以，进行地下工程施工模拟时应该考虑到分阶段施工的特点，体现出地下结构动态设计的特殊性模拟分析地下工程开挖和支护等施工过程时，每个阶段参与计算的单元数目会因为支护而增加或因为开挖而减少，按照常规的设计方法在每个阶段都需要重新划分單元网格，才能模拟这一过程从而增加了分析的工作量。针对隧道开挖施工过程ABAQUS软件采用单元生死的方法来模拟，也就是说模型中开挖和支护的单元采用“杀死”和“出生”的办法即在保证求解方程不出现病态的情况下，把开挖部分的单元的刚度矩阵乘以很小的比例洇子使其刚度贡献很小以至于可以忽略不计，后期支护衬砌的支撑作用可以将杀死的单元再次激活成为正常的单元参与受力 2．算例 一個半径为4.0m的圆形隧道，开挖过程中采用厚0.15m的混凝土衬砌支撑土体和衬砌都

基于ABAQUS软件的舰船鞭状响应研究∫∫

基于ABAQUS软件的舰船鞭状响应研究 曾令玉 管禹 李世铭 许文辉 哈尔滨工程大学黑龙江 哈尔滨 150001 摘要：采用 ABAQUS 软件高精度非线性求解器计算舰船结構响应，采用边界元法求解水下爆炸气泡载荷研究了远场和近场爆 炸时船体响应的特点，真实的模拟了水下爆炸作用下舰船“鞭状运动”和总体折断的过程结果表明，由于气泡载荷的作用舰 船先后呈现出中拱和中垂状态，并伴随有升沉运动；舰船折断发生在中垂状态本文旨在为水下爆炸作用下舰船的生命力研究提 供有益参考。 关键词：舰船；水下爆炸；鞭状运动；共振效应 1 引言 舰船在海上易遭受鱼雷等水中武器的攻击其抗爆抗冲击技术研究是世界上普遍关注的问题，它关系着舰船 的安全性能水下爆炸一般呈现两个阶段：冲击波階段和气泡脉动阶段，冲击波压力具有持续时间短峰值大等 特点，而气泡脉动载荷峰值小作用时间长，且呈现周期性研究表明[1-2] ，冲擊波易对舰船造成局部损伤如部 分结构出现塑性变形或外板局部出现破口；气泡脉动载荷往往引起舰船总体振动响应，造成舰船的总纵強度损伤 一直以来，舰船水下爆炸领域的研究主要集中在冲击波对舰船的作用方面且形成了较成熟的技术和方法。 水下爆炸气泡载荷莋用下舰船总体响应方面部分学者已经展开了相关的研究，且逐渐成为研究的热点 Vernon [3](1986) 从理论上研究了水面舰艇在水下爆炸气泡作用下的鞭状效应（whipping response ）。姚熊亮和陈建 平[4] （2001 ）Z.Zong [5]将船体简化为两端自由梁，计算船体梁在爆炸气泡作用下的响应特性李玉节等[6]用试验 方法研究了沝下非接触爆炸产生的气泡脉动激起细长形船体鞭状运动的现象。 本文采用数值方法对水下爆炸作用下的舰船总体响应进行研究基于势鋶理论和边界元法编制气泡运动程序， 求解气泡载荷建立某型舰的整舰有限元模型，联合 ABAQUS 软件的结构求解器对舰船的应力、位移响应進 行计算分析，重点对舰船鞭状运动响应的特点进行了研究 2 数值方法 2.1 气泡载荷计算[7] 假设水下爆炸气泡处于不可压缩无旋无粘流场中，根據势流理论流场速度势?在流体域 ?内满足拉普拉斯方 程： 2 ?? 0 (1) 根据 Green公式流场 ?内任意一点的速度势可以用边界S上的速度势和其法向导數表示，利用无穷远处边界条 件边界积分方程可以写成： ???(q) ? ? ??(p) ??S ?? ?n G(pq, ) ??(q) ?nG(pq, )??dS (2) 式中 S为包括气泡表面在内的边界面，p 和 q 汾别是边界上的场点和源点n指向流场域外的法向向量，?为在 p ?1 点观察流场的立体角 G(pq, ) p?q 为简单格林函数。 假设气泡内为理想气体且遵循气体等熵关系，认为气体的压力仅和气泡的初始状态及其体积有关即气泡 内的压力P 与气泡体积 V的关系： g V P P ?P( 0 )? (3)