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

1、DytranDytran般耦合流固耦合面失效实现方法与算例说明一.建模用MSC.Patran建立舱室模型而令1.2.3.特殊处理的有两点1在孔的位置也要同样每个舱室构成一个闭合的空敗将来好好用这个空由于每个舱室都需要建成一个封要重复划分左隔板，一.,建议将每个舱室作为一个组，这样划分网格会自动couple卡。要哪划分网格并将其单元属性定为定来好丿，-.闭的空间来定义本区域的定义不能”。，因为不同的流丿层为结构划分流体网格层为“哑元而令1.2.3.特殊处理的有两点1在孔的位置也要同样每个舱室构成一个闭合的空敗将来好好用这个空由于每个舱室都需要建成一个封要重复划分左隔板，一.,建议将每个舱室作为一

2、个组，这样划分网格会自动couple卡。要哪划分网格并将其单元属性定为定来好丿，-.闭的空间来定义本区域的定义不能”。，因为不同的流丿层为结构划分流体网格层为“哑元DytranDytran般耦合流固耦合面失效实现方法与算例说明流体，也就是舱室内外的空气，用euler单元划分。由于采用一般耦合(GENERALCOUPLE)定义舱室与空气冲击波之间的耦合，euler单元节点不用与lagrang单元节点重合，因此euler单元的划分比较随意，只要比所需要耦合的封闭空间大就行了。但是euler单元可能比较多，为了避免网格生成过程的麻烦，采用MESH命令自动生成euler网格。MESH,1,BOX,+,

3、3.0,-0.5,2.0,4.5,6.0,3.5,+,20,25,15,EULER,201!MESH，1，box形网格划分,MESH,1,BOX,+,3.0,-0.5,2.0,4.5,6.0,3.5,+,20,25,15,EULER,201!+，box八个顶点中x坐标最小值，y坐标最小值，z,x!方向长度，y方向长度，z方向长度，+！+，x方向划分单元数，y，z,euler形单元，单！元属性号流体的euler单元属性定义为PEULER1PEULERl,201,2ndorder,101!PEULER1,属性号，二阶精度，流体初始条件号流体区域属性定义TICEUL101+A000003+A0000

4、03SPHERE8352+A000004+A000004SPHERE5341SPHERE85.252.53.751.0SPHERE55.252.53.7515.0流体材料定义DMAT31.23流体状态方程定义EOSGAM31.4流体初始条件定义TICVALTICVALDENSITY1.2TICVALTICVALDENSITY1.2DENSITY120.0SIE0.1e6SIE5.0e6对于euler单兀的说明可以参考Chapter_16_FeaturesofStandardEurer.ppt三定义一般耦合dytran中有两种耦合方法，一种是ALE,另一种是一般耦合。ALE方法建模比较麻烦，结构

5、和流体单兀节点必须重合。因此这里用一般耦合COUPLE71OUTSIDECOUPLE71OUTSIDE+A000002+A0000021定义耦合面（SURFACE）SURFACE1ELEM定义耦合面单元组SET141THRU定义COUP1FLCOUP1FL,1,1.2,0.1e6定义耦合之间的耦合COUP1INTCOUP1INT,1,7,8ONON16+A000002+A00000214716DytranDytran般耦合流固耦合面失效实现方法与算例说明！考虑耦合面失效的快速耦合法！考虑耦合面失效的快速耦合法！初始步长！最小步长!欧拉流体单元采用ROE求解器！三阶龙个苦他积分!时间步长系数把7

6、号耦合（也就是“COUPLE7”）和8号耦合，耦合起来。这样舱室内爆炸冲击波就可以穿透孔（COUPOR卡定义）和失效单元传播到外面和其他舱室。下面两个命令定义两个耦合之间的孔PORFLCPL，COUPOR一般耦合的说明可以参考hapter_17_GeneralCouple.ppt四.定义几个参数PARAM,FASTCOUP,FAILPARAM,INISTEP,1E-7PARAM,MINSTEP,1E-15PARAM,LIMITER,ROEPARAM,RKSCHEME,3PARAM,FAILOUT,NOPARAM,NZEROVEL,NOPARAM,STEPFCT,0.5五.参考文献MSC.Dytra