急求一个用双椭球热源模型做焊接热模拟的例子
本帖最后由 wei~wei 于 2016-4-9 09:39 编辑希望提供建模过程,尤其是step,interaction,load模块的设置过程的录像,子程序用Fortran语言编写,模型可以小一些,只要把具体操作讲清楚就好。为了使焊接数值模拟过程更接近于实际,热源的添加应该是体热源和表面热源同时进行,我只知到怎样添加一种热源,但不知道怎样同时添加,哪位大神做过这方面的模拟,请指教一下,在此先感谢了。我只知到怎样添加一种热源。
在Abaqus模拟焊接的方法(图文操作)_百度文库 http://wenku.baidu.com/link?url=CRiZKmsTkBglyXlIud0FEjs9Jo4SDxNgFKxIZ_iSX-K-ae3uuaWA3WkwIER4Fzlrs2Pwg3-xg5Z-F_elX18byvSwAdkmfHI4iAaWoUr6jjW
Abaqus焊接模拟分析解析_百度文库 http://wenku.baidu.com/view/25857bc533687e21ae45a948.html
热源的话以前应该有人发过 本帖最后由 wei~wei 于 2016-4-2 14:51 编辑
申博文1122 发表于 2016-3-29 09:21
在Abaqus模拟焊接的方法(图文操作)_百度文库 http://wenku.baidu.com/link?url=CRiZKmsTkBglyXlIud0FEjs9Jo ...
谢谢你的链接,挺有帮助的,问题已经解决了,但又遇到了新的问题,就是进行多层焊时,不知道怎么在子程序中设置循环,尝试着改了几个子程序,还是不行,总提醒程序编译错误,不知道你能不能解决一下,非常感谢,下面是焊单层时的子程序:
SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
1 TEMP,PRESS,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
DIMENSION COORDS(3),FLUX(2),TIME(2)
CHARACTER*80 SNAME
wu=380.0
wi=2.0
effi=0.8
v=0.006
q=wu*wi*effi
d=v*TIME(2)
x=COORDS(1)
y=COORDS(2)
z=COORDS(3)
x0=0.004
y0=0.0
z0=0.004
a=0.003
b=0.003
c=0.002
aa=0.006
f1=1.0
PI=3.1415926
beta=0.0
aacos=a/cos(beta/180.0*PI)
cccos=c*cos(beta/180.0*PI)
aaacos=aa/cos(beta/180.0*PI)
cccos=c*cos(beta/180.0*PI)
heat1=6.0*sqrt(3.0)*q/(a*b*c*PI*sqrt(PI))*f1
heat2=6.0*sqrt(3.0)*q/(aa*b*c*PI*sqrt(PI))*(2.0-f1)
shape1=exp(-3.0*(x-x0-d)**2/aacos**2-3.0*(y-y0)**2/b**2
$ -3.0*(z-z0)**2/cccos**2)
shape2=exp(-3.0*(x-x0-d)**2/aaacos**2-3.0*(y-y0)**2/b**2
$ -3.0*(z-z0)**2/cccos**2)
C JLTYP=1,表示为体热源
JLTYP=1
IF(x .GE.(x0+d)) THEN
FLUX(1)=heat1*shape1
ELSE
FLUX(1)=heat2*shape2
ENDIF
RETURN
END wei~wei 发表于 2016-4-2 14:38
谢谢你的链接,挺有帮助的,问题已经解决了,但又遇到了新的问题,就是进行多层焊时,不知道怎么在子程序 ...
查看log文件把里面的错误解决了就行了 wei~wei 发表于 2016-4-2 14:38
谢谢你的链接,挺有帮助的,问题已经解决了,但又遇到了新的问题,就是进行多层焊时,不知道怎么在子程序 ...
beta=0.0
aacos=a/cos(beta/180.0*PI)
cccos=c*cos(beta/180.0*PI)
aaacos=aa/cos(beta/180.0*PI)
cccos=c*cos(beta/180.0*PI)
一直不是很明白这几行是干嘛用的啊,beta=0.0,那aacos就是等于a啊,干嘛要这样呢 申博文1122 发表于 2016-4-6 09:59
beta=0.0
aacos=a/cos(beta/180.0*PI)
cccos=c*cos(beta/180.0*PI)
这个子程序我也是在网上找的,这几个赋值语句好像是控制焊丝与焊接方向的倾角的,beta=0.0,好像是说此时焊丝与焊接方向的倾角是90°,这是一个特殊情况,aa=a。如果倾角不是90°的话,aa就不等有a了,这是我通过网上的介绍,自己这样理解的,不知道对不对。比如要进行双面焊的话,当焊接背面时,加热方向肯定与正面不一样,这就需要有一个角度进行调节。 问题已解决 wei~wei 发表于 2016-4-6 19:13
问题已解决
请问方便留个联系方式么,我的QQ是1132566151 请问一下,用ABUSE加载双椭球热源和用ansys加载步骤相差大吗
请问一下,用ABUSE加载双椭球热源和用ansys加载步骤相差大吗 liu19930001 发表于 2016-4-19 19:58
请问一下,用ABUSE加载双椭球热源和用ansys加载步骤相差大吗
这个我不太清楚,我不会用ANSYS,你在论坛里再找一下吧,应该有。 好的,多谢! 請問一下樓主
後段
ELSE
FLUX(1)=heat2*shape2
這邊的後橢球會作用嗎
ELSE不是前面的 if判斷有的話就不作用嗎
對程式語言這邊不熟 想請教請教 :L 楼主能不能解释一下两个heat的表达式的含义呢?
页:
[1]