- 积分
- 0
- 注册时间
- 2010-9-8
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2011-3-4 16:00:24
|
显示全部楼层
来自 大连理工大学
4# xccaii
这是其中的一段,热源的移动,但是发现焊了一半后发生了变化,前面一半比较好,后面就不知道扎回事,我觉得是不是移动有问题
max_time=b2/(v*ti)
max_y=b2/(v*ti)
max_x=t/(v*ti)
*dim,fl1,table,max_y,max_x,max_time,y,x,time
*do,k,1,max_time,1
*do,j,1,max_y,1
fl1(j,0,k)=(j-1)*t1
*enddo
*do,i,1,max_x,1
fl1(0,i,k)=(i-1)*t1
*enddo
*do,i,1,max_y,1
*do,j,1,max_x,1
y=0.5*b2-(i-1)*t1
x=0.5*t-(j-1)*t1
yce=(k-1)*t1!************************热源中心
xce=0
di=sqrt((y-yce)**2+(x-xce)**2)
*if,di,le,R,then
fl1(i,j,k)=qmax*EXP(-3*di**2/(R**2))
*else
fl1(i,j,k)=0
*endif
*enddo
*enddo
fl1(0,0,k)=(k-1)*ti
esel,s,type,,all
esel,u,type,,4
nsle,s
nsel,r,loc,x,-0.5*t,0.5*t
nsel,r,loc,z,0
sf,all,hflux,%fl1%
time,k*ti !载荷步结束的时间值
nsubst,1 !子步数和时间步长
allsel,all
solve
!sfedele,all,6,hflux
sfdele,all,hflux
allsel
*enddo |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|