超声流固耦合模拟求助！～～
小弟在用ansys做一个流固耦合的超声波无损检测的模拟。模型大概就是上面是液体，下面是固体，中间耦合面，超声波激励（用Fluid 中的 pressure加载）加在水中，希望得到从水中传到固体中的声波传播反射情况。程序是参考网上搜到的一个大坝的例子改的，现在的问题是水中的激励传到固液耦合面上就不对固体产生影响而反射了，因此无法观测到任何声波传播到固体中的现象。 我的程序如下：自己觉得应该是耦合面上边界条件的问题，若能使耦合面水和固体位移形变一样，激励就可以传下去了，不知道对不对？如何设定这样的边界条件。也有可能是其他方面的问题，请各位高手指点！/filname,frank
/units,si
/prep7
et,2,fluid29
mp,dens,2,1000
mp,sonc,2,1440
rectng,0,20e-3,10e-3,20e-3
allsel,all
aatt,2,,3
et,1,plane42,,,2
et,3,fluid29,,1
mp,dens,1,2700
mp,ex,1,70e9
mp,nuxy,1,0.34
rectng,0,20e-3,0,10e-3
asel,s,loc,y,0,10e-3
aatt,1,,1
esize,1e-4
mshape,0,2d
mshkey,0
allsel,all
amesh,all
esel,s,type,,2
nsle,s
nsel,r,loc,y,9.8e-3,10.2e-3
esln,r
emodif,all,type,3
allsel,all
/pnum,type,1
/number,1
nsel,s,loc,y,9.8e-3,10.2e-3
esln,s
esel,r,mat,,2
sf,all,fsi
esel,s,mat,,1
nsle,s
nsel,r,loc,x,19.8e-3,20.2e-3
d,all,uy,0
d,all,ux,0
esel,s,mat,,2
nsle,s
nsel,r,loc,y,19.8e-3,20.2e-3
d,all,pres,0
/pbc,all,,1
/psf,fsi,,2
allsel,all
gplot
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'myfunc'
*SET,_FNCCSYS,0
! /INPUT,H:\myfunc1.func,,,1
*DIM,%_FNCNAME%,TABLE,6,15,3,,,,%_FNCCSYS%
!
! Begin of equation: 0
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0
*SET,%_FNCNAME%(0,3,1), 0
*SET,%_FNCNAME%(0,4,1), 0
*SET,%_FNCNAME%(0,5,1), 0
*SET,%_FNCNAME%(0,6,1), 0
*SET,%_FNCNAME%(0,7,1), 0
*SET,%_FNCNAME%(0,8,1), 0
*SET,%_FNCNAME%(0,9,1), 0
*SET,%_FNCNAME%(0,10,1), 0
*SET,%_FNCNAME%(0,11,1), 0
*SET,%_FNCNAME%(0,12,1), 0
*SET,%_FNCNAME%(0,13,1), 0
*SET,%_FNCNAME%(0,14,1), 0
*SET,%_FNCNAME%(0,15,1), 0
! End of equation: 0
!
! Begin of equation: {TIME}/((4e-8)*(4e-8))*exp(-{TIME}/4e-8)*4e-6
*SET,%_FNCNAME%(0,0,2), 1e-8, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 4e-8, 0, 0, 0
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 4e-8, 0, 0, -1
*SET,%_FNCNAME%(0,3,2), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,2), 0.0, -1, 0, 1, 1, 4, -3
*SET,%_FNCNAME%(0,5,2), 0.0, -2, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,6,2), 0.0, -3, 0, 1, 0, 0, -2
*SET,%_FNCNAME%(0,7,2), 0.0, -4, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,8,2), 0.0, -2, 0, 1, -4, 3, 1
*SET,%_FNCNAME%(0,9,2), 0.0, -3, 0, 4e-8, 0, 0, -2
*SET,%_FNCNAME%(0,10,2), 0.0, -4, 0, 1, -2, 4, -3
*SET,%_FNCNAME%(0,11,2), 0.0, -2, 7, 1, -4, 0, 0
*SET,%_FNCNAME%(0,12,2), 0.0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,13,2), 0.0, -1, 0, 4e-6, 0, 0, -3
*SET,%_FNCNAME%(0,14,2), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,15,2), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: {TIME}/((4e-8)*(4e-8))*exp(-{TIME}/4e-8)*4e-6
!
! Begin of equation: 0
*SET,%_FNCNAME%(0,0,3), 3e-6, -999
*SET,%_FNCNAME%(2,0,3), 0.0
*SET,%_FNCNAME%(3,0,3), 0.0
*SET,%_FNCNAME%(4,0,3), 0.0
*SET,%_FNCNAME%(5,0,3), 0.0
*SET,%_FNCNAME%(6,0,3), 0.0
*SET,%_FNCNAME%(0,1,3), 1.0, 99, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,3), 0
*SET,%_FNCNAME%(0,3,3), 0
*SET,%_FNCNAME%(0,4,3), 0
*SET,%_FNCNAME%(0,5,3), 0
*SET,%_FNCNAME%(0,6,3), 0
*SET,%_FNCNAME%(0,7,3), 0
*SET,%_FNCNAME%(0,8,3), 0
*SET,%_FNCNAME%(0,9,3), 0
*SET,%_FNCNAME%(0,10,3), 0
*SET,%_FNCNAME%(0,11,3), 0
*SET,%_FNCNAME%(0,12,3), 0
*SET,%_FNCNAME%(0,13,3), 0
*SET,%_FNCNAME%(0,14,3), 0
*SET,%_FNCNAME%(0,15,3), 0
! End of equation: 0
!-->
FINISH
/SOL
FLST,5,1,1,ORDE,1
FITEM,5,4568
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
!*
/GO
!*
!*
D,_Y1,PRES, %MYFUNC%
!*
CMDELE,_Y1
!*
/solu
antype,trans
nsel,s,loc,y,0 !选择nodes子集
dsym,symm,y !指定节点上的约束是否为对称的
nsel,s,loc,x,20e-3
dsym,symm,x
allsel
time,8e-7
tintp,,,,1.0
autots,off
deltim,2e-9
outres,nsol,1
solve 自己顶，这问题已经纠结快半个月了，工作就卡在这里，老板都要发飙了……:'(:'(:'( 高手怎么还不出现啊 水和固体结合的以前单元要删除改为流固耦合单元
页:
[1]