Dytran一般耦合流固耦合面失效实现方法与算例说明

1、Dytran一般耦合流固耦合面失效实现方法与算例说明一.建模用MSC.Patran建立舱室模型需要特殊处理的有两点1 .在孔的位置也要同样划分网格并将其单元属性定义为DUMMY即“哑元”这样可以将每个舱室构成一个闭合的空间，将来好用这个空间定义舱室内外的流体区域。2 .由于每个舱室都需要建成一个封闭的空间来定义流体区域，因此两个舱室之间的隔板需要重复划分单元，因为不同的流体区域的定义不能出现同一个单元。这样就形成了两层隔板，一层为结构一层为“哑元”。3 .建议将每个舱室作为一个组，这样划分网格会自动形成连续的单元编号，便于后面定义couple卡。划分流体网格流体，也就是舱室内外的空气，用eul

2、er单元划分。由于采用一般耦合(GENERALCOUPLE)定义舱室与空气冲击波之间的耦合，euler单元节点不用与lagrang单元节点重合，因此euler单元的划分比较随意，只要比所需要耦合的封闭空间大就行了。但是euler单元可能比较多，为了避免网格生成过程的麻烦,MESH,1,BOX,+,3.0,-0.5,2.0,4.5,6.0,3.5,+采用MESH命令自动生成euler网格。!MESH,1,box形网格划分，,+!+,box八个顶点中x坐标最小值，y坐标最小值，z,x!方向长度，y方向长度，z方向长度，-+y,z,，euler形单元，单+,20,25,15,EULER,201流体的

3、euler单元属性定义为!+,x方向划分单元数,!元属性号PEULER1PEULER1,201,2ndorder,流体区域属性定义101!PEULER1,属性号，二阶精度，流体初始条件号+A000003SPHERE8+A000004SPHERE5SPHERE85.25SPHERE55.25流体材料定义DMAT31.23流体状态方程定义EOSGAM31.4TICEUL10133流体初始条件定义2.52.5+A000003+A0000043.753.7511.015.0TICVALTICVALDENSITYDENSITY1.2120.0SIE0.1e6SIE5.0e6对于euler单元的说明可以参

4、考chapter_16_FeaturesofStandardEurer.ppt定义一般耦合dytran中有两种耦合方法，一种是ALE,另一种是一般耦合。ALE方法建模比较麻烦，结构和流体单元节点必须重合。因此这里用一般耦合COUPLE71OUTSIDEON+A000002+A00000211定义耦合面（SURFACE）SURFACE1ELEM4定义耦合面单元组SET141THRU716,定义COUP1FLCOUP1FL,1,1.2,0.1e6定义耦合之间的耦合COUP1INTON16+A000002+A000002COUP1INT,1,7,8把7号耦合（也就是“COUPLE7,”）和8号耦合，

5、耦合起来。这样舱室内爆炸冲击波就可以穿透孔（COUPOR卡定义）和失效单元传播到外面和其他舱室。下面两个命令定义两个耦合之间的孔PORFLCPL,COUPOR般耦合的说明可以参考chapter_17_GeneralCouple.ppt四.定义几个参数PARAM,FASTCOUP,FAILPARAM,INISTEP,1E-7PARAM,MINSTEP,1E-15PARAM,LIMITER,ROEPARAM,RKSCHEME,3PARAM,FAILOUT,NOPARAM,NZEROVEL,NOPARAM,STEPFCT,0.5!考虑耦合面失效的快速耦合法!初始步长!最小步长!欧拉流体单元采用ROE求解器!三阶龙个苦他积分!时间步长系数五.参考文献1. MSC.DytranExampleproblems中的4