..:: ANSYS专区-斑竹,帮忙,耦合,热力,指点-热力耦合，请斑竹帮忙指点，{Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛-(手机验证注册)}

石头公子 发表于 2007-5-2 14:47:24

热力耦合，请斑竹帮忙指点

斑竹及各位高手，我是做焊接热应力耦合得，程序运行得时候出现以下问题，请各位帮我指点一下，我该怎么该。
谢谢！
Possible brick to tetrahedron interface detected at area 3 on volume 7.
The current element (type 2) does not support the required transition
shape (pyramids).No transition between bricks and tetrahedra will be
defined.
建模方面得命令是这样得
i=180
u=25
v=0.01
x1=-0.002
rb=0.002
length=0.12
wideth=0.06
height=0.006
et,1,solid5
et,2,solid98
block,0,length,0,wideth,0,height
block,0,length,0,3*rb,0,height
block,0,length,0,6*rb,0,height
block,0,length,0,12*rb,0,height
vovlap,all
numcmp,area
/view,1,1,1,1

lesize,4,rb
lesize,5,rb
lesize,9,rb
lesize,14,rb
lesize,50,rb
lesize,10,rb
lesize,52,rb
lesize,23,rb
lesize,51,rb
lesize,49,rb
lesize,24,rb
lesize,19,rb
lesize,31,2*rb
lesize,26,2*rb

lesize,43,4*rb
lesize,38,4*rb
lesize,2,4*rb
lesize,7,4*rb
lesize,53,4*rb
lesize,54,4*rb
lesize,55,4*rb
lesize,56,4*rb
lesize,11,1.5*rb
lesize,12,1.5*rb
lesize,47,1.5*rb
lesize,48,1.5*rb
save
type,1
mshkey,1
mshape,0
vmesh,5
vmesh,6
type,2
mshkey,0
mshape,1
vmesh,7
vmesh,8
save

icst_li 发表于 2007-5-2 21:23:41

好象说的是要有两种网格形式之间的过渡.有些单元对网格形式有要求.

[ 本帖最后由 icst_li 于 2007-5-2 21:52 编辑 ]

icst_li 发表于 2007-5-2 21:30:30

完整的命令流应该是这样的吧;
/prep7
mptemp,1,20,250,500,750,1000,1500,1700,2500
mpdate,kxx,1,1,50,47,40,27,30,3 5,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.4,0.45,0.48,0.5
mpdata,alpx,1,1,1.10E-05,1.22E-05,1.39E-05,1.48E-05,1.34E-05,1.33E-05,1.32E-05,1.31E-05
mpdata,ex,1,1,2.05E+11,1.87E+11,1.50E+11,7.00E+10,2.00E+10,1.90E+10,1.80E+10,1.20E+10

mptemp,,,,,,,,
mptemp,1,20
mptemp,2,1553
mptemp,3,1590
mptemp,4,1670
mpde,enth,1
mpdata,enth,1,,0
mpdata,enth,1,,7.5e+009
mpdata,enth,1,,9.6e+009
mpdata,enth,1,,1.1e+010

tb,biso,1,6,2,
tbtemp,20
TBDATA,,2.20E+08,2.O5E+10,,,,
TBTEMP,500
TBDATA,,1.30E+08,1.50E+10,,,,
TBTEMP,750
TBDATA,,4.00E+07,7.00E+09,,,,
TBTEMP,1000
TBDATA,,2.50E+07,2.00E+09,,,,
TBTEMP,1700
TBDATA,,1.00E+06,1.80E+09,,,,
TBTEMP,2500
TBDATA,,1.00E+05,1.20E+09,,,,
SAVE
FINI

/prep7
/units, si
rb=0.002
length=0.12
width=0.06
height=0.006
et,1,solid5
et,2,solid98
block,0,length,0,width,0,height
Block,0,length,0,3*Rb,0,height

Block,0,length,0,Rb*6,0,height
Block,0,length,0,Rb*12,0,height
VOVLAP, ALL
/view,1,1,1,1

Lesize,4,rb
Lesize,5,rb
Lesize,9,rb
Lesize,10,rb
Lesize,14,rb
Lesize,19,rb
Lesize,23,rb
Lesize,24,rb
Lesize,49,rb
Lesize,50,rb
Lesize,51,rb
Lesize,52,rb

Lesize,31,2*rb
Lesize,26,2*rb

Lesize,43,4*rb
Lesize,38,4*rb
Lesize,2,4*rb
Lesize,7,4*rb
Lesize,53,4*rb
Lesize,54,4*rb
Lesize,55,4*rb
Lesize,56,4*rb
Lesize,11,1.5*rb
Lesize,12,1.5*rb
Lesize,47,1.5*rb
Lesize,48,1.5*rb
save

type,1
mshkey,1
mshape,0
vmesh,5
vmesh,6
type,2
mshkey,0
mshape,1
vmesh,7
vmesh,8
Save

/SOLU
DA,3,SYMM
/GO
DL,53,,UX,
DL,61,,UX,
DL,57,,UX,
DL,49,,UX,
DL,4,,UZ,
DL,9,,UY,

/solu
antype,trans,new
nlgeom,on
timint,0,struct
Timint,1,thern
timint,0,mag
timint,0,elect
Tintp,0.005,,,-1,0.5,0.2
tref,25
nropt,auto,,
autos,on
outres,,all
v=0.008
tinc=rb/v

U=25
I=180
ef=0.75
p=25*180*ef

r=0.006
qmax=p/(3.14159*r**2)

max_time=length/rb+1
max_x=length/rb+1
max_y=width/rb+1
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*rb
*enddo

!加载移动热源
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/exp(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc !存储时间
nsel,s,loc,z,height
sf,all,hflux,%flux2%

!加载对流载荷
allsel,all
nsel,s,loc,x,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,x,length
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,width
sf,all,conv,10,20
allsel,all
nsel,s,loc,z,0
sf,all,conv,10,20
allsel,all
time,k*rb/v
nsubst,5
solve
*enddo

allsel,all
nsel,s,loc,z,height
sf,all,conv,10,20
allsel,all
save

石头公子 发表于 2007-5-3 22:17:49

谢谢楼上的，能和我联系一下吗？想请教几个问题
QQ：93446546
谢谢了

icst_li 发表于 2007-5-4 09:33:24

好啊!我也正打算做关于焊接过程的模拟. 最近看了一些资料.
QQ:396375268

石头公子 发表于 2007-5-4 11:18:23

运行了一下，出现很多错误
没有定义单元关键字
我也不知道怎么搞

lixiaodong163 发表于 2011-5-27 10:52:01

页: [1]

Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册)'s Archiver

热力耦合，请斑竹帮忙指点