流固耦合错误，求救（附命令流）

shuizs

我现在做一个管道和水的耦合，水入口端加压，想得到当入口压力变化时，管壁的位移，可老是报错，都好几天了都没有一点进展，能帮忙看看么？我碰到一个错误，说：Fluid soild interface boundary incorrectly type,
*** WARNING ***                         CP =      90.562   TIME= 10:37:53
Flow solution requested for model containing only solid material        
elements.  This request will be ignored and flow solution will be      
turned off.                                                            
*** WARNING ***                         CP =      90.562   TIME= 10:37:53
Neither flow or thermal analyses are activated.  Analysis will be      
processed but the solution will not change.
命令流如下：
fini
/clear
/filname,xiexie,1
/prep7  
et,1,45
et,2,142,,,,1
mp,ex,1,2e11
mp,dens,1,7800   
mp,prxy,1,0.3   
mp,dens,2,1000   
mp,visc,2,0.001
CYL4, , ,0.3, ,0.32, ,5  ！创建外面的管道，固体
TYPE,   1   
MAT,  1
vsel,s,volu,,1
Vsweep,all

CYL4, , ,0.3, , , ,5    ！流体
TYPE,  2  
MAT,   2
vsel,s,volu,,2
vsweep,all

/solu
csys,1                  ！转换为柱坐标系
vsel,s,volu,,1          ！定义固体单元
allsel,below,volu
nslv,s,1
cm,solid,node

vsel,s,volu,,2          ！定义流体单元
allsel,below,volu
nslv,s,1
cm,fluid,node
alls                    ！外壳进行约束
nsel,s,loc,x,0.3
nsel,r,loc,z,5          ！出口
cmsel,r,solid
d,all,uz,0
alls
nsel,s,loc,x,0.3
nsel,r,loc,z,0          ！入口
cmsel,r,solid
d,all,uz,0

alls                    ！定义耦合面
cmsel,s,solid
nsel,r,loc,x,0.3
sf,all,fsin,1
alls                    ！流体入口压 0.2E6
cmsel,s,fluid
nsel,r,loc,x,0,0.3
nsel,r,loc,z,0
d,all,pres,0.2e6
alls                    ！出口压 0      
cmsel,s,fluid
nsel,r,loc,x,0,0.3
nsel,r,loc,z,5
d,all,pres,0
                        ！定义流体侧的耦合面
alls
cmsel,s,fluid
nsel,r,loc,x,0.3
d,all,vx,0
d,all,vy,0
d,all,vz,0
sf,all,fsin,1
alls
                           ！流体参数设置
FLDATA1,SOLU,TRAN,1        
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,0
FLDATA1,SOLU,TURB,0
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0  
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,1  
FLDATA1,SOLU,RDSF,1

FLDATA12,PROP,DENS,0
FLDATA13,VARY,DENS,0
FLDATA12,PROP,VISC,0
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0
FLDATA7,PROT,DENS,CONSTANT  
FLDATA8,NOMI,DENS,1000,
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,CONSTANT  
FLDATA8,NOMI,VISC,1e-3,
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS  
FLDATA7,PROT,COND,CONSTANT  
FLDATA8,NOMI,COND,-1,   
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT  
FLDATA8,NOMI,SPHT,-1,   
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
FLDATA15,PRES,REFE,101350,  
FLDATA16,BULK,BETA,36,  
FLDATA17,GAMM,COMP,1.4,
FLDATA14,TEMP,NOMI,293,
FLDATA14,TEMP,TTOT,293,
FLDATA14,TEMP,BULK,293,
TOFFSET,0,  
FLDATA4,TIME,STEP,0.001,
FLDATA4,TIME,ISTEP,0
FLDATA4,TIME,NUMB,10,   
FLDATA4,TIME,TEND,1.0e06
FLDATA4,TIME,GLOB,15,   
FLDATA4,TIME,VX,0.01,   
FLDATA4,TIME,VY,0.01,   
FLDATA4,TIME,VZ,0.01,   
FLDATA4,TIME,PRES,1e-004,   
FLDATA4,TIME,TEMP,1e-006,   
FLDATA4,TIME,ENKE,0.01,
FLDATA4,TIME,ENDS,0.01,
FLDATA4A,STEP,OVER,0,   
FLDATA4,TIME,OVER,0
FLDATA4A,STEP,APPE,0,   
FLDATA4,TIME,APPE,1.0e6
FLDATA4A,STEP,SUMF,10,  
FLDATA4,TIME,SUMF,1.0e6
FLDATA4,TIME,BC,0   
FLDATA4,TIME,TEND,1000000,  
FSAN,on  
FSOR,solid, FLUID
FSTR,FLUID,TRAN
FSTR,SOLID,TRAN
FSIN,1   
FSTI,0.025,0.001,   
FSDT,0.001,
FSIT,20,
FSCO,ALL,0.01,  
FSOU,0.1E-01
FSCO,ALL,0.01,  
FSOU,1  
FSRE,ALL,0.8,   
solve


版主大侠帮忙看看啊，不胜感激！

chenyinqiang

在流场分析中不能用固体单元,这就是你在求解的过程产生上述问题的原因

文枫

流体单元142必须是type,1