急求一个用双椭球热源模型做焊接热模拟的例子

wei~wei

希望提供建模过程，尤其是step，interaction，load模块的设置过程的录像，子程序用Fortran语言编写，模型可以小一些，只要把具体操作讲清楚就好。为了使焊接数值模拟过程更接近于实际，热源的添加应该是体热源和表面热源同时进行，我只知到怎样添加一种热源，但不知道怎样同时添加，哪位大神做过这方面的模拟，请指教一下，在此先感谢了。我只知到怎样添加一种热源。

申博文1122

在Abaqus模拟焊接的方法(图文操作)_百度文库 http://wenku.baidu.com/link?url= ... AdkmfHI4iAaWoUr6jjW

Abaqus焊接模拟分析解析_百度文库 http://wenku.baidu.com/view/25857bc533687e21ae45a948.html
热源的话以前应该有人发过

wei~wei

申博文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

blswjs

wei~wei 发表于 2016-4-2 14:38
谢谢你的链接，挺有帮助的，问题已经解决了，但又遇到了新的问题，就是进行多层焊时，不知道怎么在子程序 ...

查看log文件把里面的错误解决了就行了

申博文1122

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啊，干嘛要这样呢

wei~wei

申博文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

问题已解决

wh1132566151

wei~wei 发表于 2016-4-6 19:13
问题已解决

请问方便留个联系方式么，我的QQ是1132566151

liu19930001

请问一下，用ABUSE加载双椭球热源和用ansys加载步骤相差大吗

liu19930001

请问一下，用ABUSE加载双椭球热源和用ansys加载步骤相差大吗

wei~wei

liu19930001 发表于 2016-4-19 19:58
请问一下，用ABUSE加载双椭球热源和用ansys加载步骤相差大吗

这个我不太清楚，我不会用ANSYS，你在论坛里再找一下吧，应该有。

liu19930001

好的，多谢！

b9638032

請問一下樓主   

後段
ELSE
          FLUX(1)=heat2*shape2
這邊的後橢球會作用嗎

ELSE  不是前面的 if  判斷有的話就不作用嗎

對程式語言這邊不熟   想請教請教

1099

楼主能不能解释一下两个heat的表达式的含义呢？