资源描述：
《dytran一般耦合流固耦合面失效实现方法与算例说明》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、Dytran一般耦合流固耦合面失效实现方法与算例说明一．建模用MSC.Patran建立舱室模型需要特殊处理的有两点1.在孔的位置也要同样划分网格并将其单元属性定义为DUMMY即“哑元”这样可以将每个舱室构成一个闭合的空间，将来好用这个空间定义舱室内外的流体区域。2.由于每个舱室都需要建成一个封闭的空间来定义流体区域，因此两个舱室之间的隔板需要重复划分单元，因为不同的流体区域的定义不能出现同一个单元。这样就形成了两层隔板，一层为结构一层为“哑元”。3.建议将每个舱室作为一个组，这样划分网格会自动形成连续的单元编号，便
2、于后面定义couple卡。二．划分流体网格流体，也就是舱室内外的空气，用euler单元划分。由于采用一般耦合（GENERALCOUPLE）定义舱室与空气冲击波之间的耦合，euler单元节点不用与lagrang单元节点重合，因此euler单元的划分比较随意，只要比所需要耦合的封闭空间大就行了。但是euler单元可能比较多，为了避免网格生成过程的麻烦，采用MESH命令自动生成euler网格。MESH,1,BOX,,,,,,,+！MESH，1，box形网格划分，，，，，，，++,3.0,-0.5,2.0,4.5,6.0,
3、3.5,,,+！+，box八个顶点中x坐标最小值，y坐标最小值，z，x！方向长度，y方向长度，z方向长度，，，++,20,25,15,,,,EULER,201！+，x方向划分单元数，y…，z…，，，，，euler形单元，单！元属性号流体的euler单元属性定义为PEULER1PEULER1,201,,2ndorder,101！PEULER1，属性号，，二阶精度，流体初始条件号流体区域属性定义TICEUL101+A000003+A000003SPHERE8352+A000004+A000004SPHERE5341SP
4、HERE85.252.53.751.0SPHERE55.252.53.7515.0流体材料定义DMAT31.23流体状态方程定义EOSGAM31.4流体初始条件定义TICVAL4DENSITY1.2SIE0.1e6TICVAL5DENSITY120.0SIE5.0e6对于euler单元的说明可以参考chapter_16_FeaturesofStandardEurer.ppt一．定义一般耦合dytran中有两种耦合方法，一种是ALE，另一种是一般耦合。ALE方法建模比较麻烦，结构和流体单元节点必须重合。因此这里用一般
5、耦合COUPLE71页：2定义耦合的SURFACE号OUTSIDE页：2定义耦合的SURFACE外部EULER域不计算ON页：2SURFACE所包含的单元自动转换法向方向ON页：2检查单元法向是否正确16页：2定义孔（COUPOR）命令卡号+A000002+A000002+A000002+A0000021页：2COUP1FL号，COUP1FL是定义SURFACE破坏后周围流体性质1定义耦合面（SURFACE）页：2COUP1INT号，COUP1INT定义两个耦合之间的耦合SURFACE1ELEM页：2用单元定义SU
6、RFACE4页：2用于定义SURFACE的单元组SET1的号码定义耦合面单元组SET141THRU716页：2从1到716号单元…………定义COUP1FLCOUP1FL,1,1.2页：2密度,0.1e6页：2内能定义耦合之间的耦合COUP1INTCOUP1INT,1,7页：2耦合7,8把7号耦合（也就是“COUPLE7……”）和8号耦合，耦合起来。这样舱室内爆炸冲击波就可以穿透孔（COUPOR卡定义）和失效单元传播到外面和其他舱室。下面两个命令定义两个耦合之间的孔PORFLCPL，COUPOR一般耦合的说明可以参考
7、chapter_17_GeneralCouple.ppt二．定义几个参数PARAM,FASTCOUP,,FAIL！考虑耦合面失效的快速耦合法PARAM,INISTEP,1E-7！初始步长PARAM,MINSTEP,1E-15！最小步长PARAM,LIMITER,ROE！欧拉流体单元采用ROE求解器PARAM,RKSCHEME,3！三阶龙个苦他积分PARAM,FAILOUT,NO！PARAM,NZEROVEL,NO！PARAM,STEPFCT,0.5！时间步长系数三．参考文献1.MSC.DytranExamplepr
8、oblems中的4-7：bunker2.chapter_16_FeaturesofStandardEurer.ppt3.chapter_17_GeneralCouple.ppt页：2耦合8
