zxw19851015 发表于 2008-10-19 20:34:09

帮我看看焊接移动热源程序

以下是加载程序，最后一行中sf,all,hflux,%flux2% 中的%flux2% 是不是一次把移动热源全加上去还是什么意思，望知道的同行帮我解答
*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)*lsize
*enddo
   
*do,j,1,max_y,1
   flux2(0,j,k)=(j-1)*lsize
*enddo
   
*do,i,1,max_x,1
    *do,j,1,max_y,1
   xcenter=v*(k-1)*tinc
   ycenter=0
   distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
    *if,distance,le,r,then
      flux2(i,j,k)=Qmax*exp(-3*kk*distance**2/(r**2)) !高斯移动热源公式
    *else
      flux2(i,j,k)=0
    *endif
    *enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
/solu
antype,trans
tunif,15
kbc,1
timint,on
autots,on
esel,s,type,,1 !选择表面效应单元
nsle,s
sf,all,conv,-370 !施加对流边界条件
allsel
*do,t,1e-6,tm,tinc
time,t
autots,0
deltim,,,,2
kbc,0
outres,all,all
tsres,erase
sf,all,hflux,%flux2%
solve
*enddo

y1949b 发表于 2008-10-22 16:18:38

帮我看看焊接移动热源程序

也做这方面，但是是新手，帮你顶一下:)

半个苹果 发表于 2008-10-30 20:39:35

顶下！

呵呵能互相交流下么 qq：80976658

半个苹果 发表于 2008-10-30 22:23:34

！

我觉得应当是把热源按照时间的顺序加上去的！ 前面定义的三维表格数组就是起到这个作用的！

acailaogu 发表于 2008-11-1 20:58:47

正在探索中。。。。。。。。。。。。。。。。。。

eli-lu 发表于 2009-4-1 19:32:21

我也碰到过这样的热源，只是自己也不太明白

zxw19851015 发表于 2009-4-3 16:35:07

这些问题已经解决了

lengxuef 发表于 2009-4-3 17:08:31

能不能发全部代码？

ljian25 发表于 2009-4-30 17:11:56

呵呵能互相交流

ljian25 发表于 2009-5-7 14:48:50

谢谢
好像来过

tangxiushan 发表于 2009-9-13 22:45:09

学习了，不过好像还不是很清楚

qaz54809 发表于 2009-9-26 18:37:01

这个好像武汉理工某个硕士论文也用过这个，他这个是生死单元做的，如果我没记错的话，楼主没有全发。

wxjjuezhan13 发表于 2012-10-4 11:55:53

武汉理工几年前有个硕士论文是移动热源是、，至今只找到论文，至于怎么做的，命令流等等，无法求之，急急急
页: [1]
查看完整版本: 帮我看看焊接移动热源程序