Strosailor 发表于 2007-5-31 16:04:10

我做的三维焊接温度场模拟例子，可以运行

感觉温度场分布不对，似乎不怎么传热，欢迎大家交流。

/prep7
/units,si
toffst,273
*set,rb,0.001
length=0.1
width=0.05
height=0.005
et,1,solid70
et,2,solid90
et,3,solid87
mptemp,1,20,250,500,750,1000,1500,1700,2500
mpdata,kxx,1,1,50,47,40,27,30,35,40,300
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,hf,1,1,24,46,73,140,125,224,400,450,500,500
block,0,0.1,0,0.05,-0.003,0
block,0,0.1,0,0.004,-0.003,0
block,0,0.1,0,0.01,-0.003,0
vovlap, all                     
numcmp,node
numcmp,area
numcmp,volume

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,37,rb
Lesize,38,rb
Lesize,39,rb
Lesize,40,rb

Lesize,45,rb！定义过渡区单元尺寸
Lesize,46,rb
Lesize,47,rb
Lesize,48,rb

Lesize,26,3*rb   ！定义远缝区单元尺寸
Lesize,35,3*rb
Lesize,31,3*rb
Lesize,36,3*rb
Lesize,41,3*rb
Lesize,42,3*rb
Lesize,43,3*rb
Lesize,44,3*rb
Lesize,2,3*rb
Lesize,7,3*rb
Lesize,11,3*rb
Lesize,12,3*rb
Save

type,1！划分网格
mat,1
mshkey,1
mshape,0
vmesh,1
vmesh,2

type,2
mat,1
mshape,1,3D
mopt,pyra,1
mshkey,0
mshape,1
vmesh,3
mopt,split,2
vimp,3,1,3
save

/solu

v=0.0045
tm=length/v
tinc=rb/v
U=25
I=100
ef=0.75
p=25*100*ef
r=0.003
qmax=3*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
a=(i-1)*rb-xcenter
b=(j-1)*rb-ycenter
distance=sqrt(a**2+b**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
*enddo

/solu

antype,trans
ic,all,temp,20
*do,t,1e-6,tm,tinc
time,t
autos,0
deltim,,,,1
kbc,1
nsel,s,loc,x,0,0.1   ！选择要加载热源的节点
nsel,r,loc,y,0,0.05
nsel,r,loc,z,0,0
sf,all,hflux,%flux2%

FLST,5,3,5,ORDE,3!选择要加载对流载荷的表面
FITEM,5,1   
FITEM,5,3   
FITEM,5,-4
ASEL,U, , ,P51X
ALLSEL,BELOW,AREA   
NSLA,S,1
sf,all,conv,-1,20 !在节点上施加对流载荷(与空气的对流换热系数）
allsel,below,volu
timint,on                                        ！打开时间积分
tintp,,,,1                                        ！将THETA设定为1
outres,all,all                                        ！输入每个子步的结果
solve
*enddo
finish

shen2mzy 发表于 2007-6-1 18:23:59

正好需要，十分感谢

eli-lu 发表于 2009-4-1 20:36:19

我也遇到了这样的问题，中心温度场过高

fangdian 发表于 2009-7-30 20:14:52

谢谢楼主。下了

acailaogu 发表于 2009-9-2 19:24:34

ansys的非线性分析能力够呛，热应力够呛，转向msc。marc了

tangxiushan 发表于 2009-9-13 21:59:49

参考下，谢谢分享

crruose 发表于 2009-9-21 11:07:38

问题解决没原因在呢

shezhang 发表于 2009-12-6 18:51:51

命令流有什么用啊

chutianyao 发表于 2010-3-9 16:16:35

是什么原因？

wxpand2004 发表于 2010-3-10 09:10:17

我做完也感觉不对，不知道楼主解决问题了没！

wyksuda 发表于 2010-4-1 10:29:33

:)学习中，谢谢楼主。不过可惜，暂时还帮不上忙

守护蓝天 发表于 2010-4-11 17:28:04

有没有做激光切削的啊？激光切削和焊接有什么区别呀？

niulian614 发表于 2010-4-12 20:00:18

学习中，谢谢Lz

kakaroteqiqi 发表于 2010-4-13 10:51:20

来看看哈！！

181809216 发表于 2010-7-10 15:46:22

楼主解决没啊

200802172 发表于 2010-7-14 18:03:28

可能是热源施加的问题

cumt410 发表于 2010-7-22 20:18:59

谢谢啦，正好需要瞬态传热的东西

守护蓝天 发表于 2010-9-3 00:11:49

17# cumt410 你好，能告诉下不同温度下材料的密度如何查吗？

mmmm123 发表于 2010-9-4 20:02:30

谢谢楼主。

永恒 发表于 2010-9-5 07:19:03

感谢楼主，我也在做这方面的。
页: [1] 2
查看完整版本: 我做的三维焊接温度场模拟例子，可以运行