应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » 其他 » Simufact.Welding 6.0 热源二次开发的小例子
141/1返回列表
查看: 2603  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 fangdian 的 3 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会

fangdian

木虫 (正式写手)

[交流] Simufact.Welding 6.0 热源二次开发的小例子

最近做了一个Simufact.Welding 6.0 热源二次开发的小例子，与坛友共勉：
1 主要目的：开始随时间变化的焊接热源，如脉冲热源
2 基本要求：
  1）脉冲热源频率：8Hz
  2）脉宽30ms
  3）其它略
3 软件安装
   1）simufact.welding 6.0
   2）VS2010
   3）Intel Visual Fortran Composer XE 2013
3 基本步骤
3.1 软件设置
安装完之后需要配置path、Lib和include三个环境变量。
以下列出环境变量值作为参考（需要根据VS和fortran的安装位置不同进行相应的修改）。---注：以下部分为参考网上内容，需要根据具体实际修改，如下面采用的是VS2012，XE2015，根据需要修改
Path：
C:\Program Files (x86)\Microsoft Visual Studio 12.0\Common7\Tools
C:\Program Files (x86)\Intel\Composer XE 2015\redist\intel64\mkl
C:\Program Files (x86)\Intel\Composer XE 2015\bin\intel64
C:\Program Files (x86)\Microsoft Visual Studio 12.0\Common7\IDE
C:\Program Files (x86)\Windows Kits\8.1\Windows Performance Toolkit\
C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\bin\amd64
C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\bin
C:\Program Files (x86)\Microsoft SDKs\Windows\v7.0A\Bin\x64
C:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\Bin
C:\Program Files (x86)\Intel\MPI\5.0.1.037\intel64\bin
C:\Program Files (x86)\Intel\Trace Analyzer and Collector\9.0.1.035\bin

Lib：
C:\Program Files (x86)\Intel\Composer XE 2015\compiler\lib\intel64
C:\Program Files (x86)\Intel\Composer XE 2015\mkl\lib\intel64
C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\lib\amd64
C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\lib
C:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\Lib\x64
C:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\Lib

INCLUDE：
C:\Program Files (x86)\Intel\Composer XE 2015\compiler\include
C:\Program Files (x86)\Intel\Composer XE 2015\mkl\include
C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\include
C:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\Includea

3.2 常规建立并设置simufact.welding工程
  建立如下常规的simufact.welding工程
（常规建立过程省略）

Simufact.Welding 6.0 热源二次开发的小例子
  图1 常规工程
主要坐标系为：沿X方向宽50mm，沿y方向板厚10mm，沿z方向长度50mm，焊接沿着z负方向，焊接起始点坐标（50,10，-50）【单位mm】,建立对称的一半模型。

3.3 建立.f子函数
      subroutine uweldflux(f,temflu,mibody,welddim,time)
c     f(1)         flux value (to be defined)
c     f(2)         derivative of flux with respect to temperature
c                  (to be defined; optional, this might improve the
c                   convergence behavior)
c
c     temflu(1)    estimated temperature
c     temflu(2)    previous volumetric flux
c     temflu(3)    temperature at beginning of increment
c     temflu(4,5,6)integration point coordinates
c     mibody(1)    element number
c     mibody(2)    flux type
c     mibody(3)    integration point number
c     mibody(4)    flux index - not used if table input
c     mibody(5)    not used
c     mibody(6)    =1 : heat transfer
c                  =2 : joule
c                  =3 : bearing
c                  =4 : electrostatic
c                  =5 : magnetostatic
c                  =6 : acoustic
c     mibody(7)    internal element number
c     mibody(8)    layer number for heat transfer shells elements
c                  and volume flux
c     mibody(9)    Not used
c     mibody(10)   boundary condition number if table input
c     time         time
c
      dimension mibody(*),temflu(*),welddim(*)
      real*8 f
      integer mibody
      real*8 temflu,time,welddim
c* * * * * *
c      OPEN(UNIT=10,FILE='HELLO.TXT',POSITION='APPEND')
      a=0.003
      b=0.003
      c=0.003
      v=0.01
        T0=0.125
        temp=0.03
        xis=1
      zq=ceiling(time/T0)
        if(abs(MOD(time-temp,T0)).LT.1e-6) then
                xis=0
c        write(10,*) time, 111
        else if(abs(MOD(time,T0)).LT.1e-6) then
          xis=1
c        write(10,*) time,222
      else if ((time.LT.(zq*T0)).and.(time.GT.((zq-1)*T0+temp))) then
          xis=0
c     write(10,*) time, 333
      else
          xis=1
c        write(10,*) time, 444
        end if
       

      d=0.05-v*time
      q=7000*0.7
      pi=3.14
      x=temflu(4)
      y=temflu(5)
      z=temflu(6)
      mibody(6)=1
      y0=0.009
      x0=0.05
      heat=6*sqrt(3.0)*q/(a*b*c*pi)
      ex=exp(-3*(z-d)**2/c**2-3*(y-y0)**2/b**2-3*(x-x0)**2/a**2)
      f=heat*ex*xis
c        write(10,*) time, xis,f
c        close(10)
      return
      end

并存为文件名rey.f

3.4 开始子函数运行

1）点击常规项目的运行按钮
略
图2 点击运行按钮
点击后出现运行对话框，不用管它，直接关闭。
略

图3 运行对话框（直接关闭）
2）此时找到你常规工程下的_Run_文件夹位置并打开该文件夹


图4 该工程的_Run_文件夹

会发现该文件夹下面有Process.dat文件和run.bat文件。

3）用记事本方式打开并修改Process.dat文件(此步不可少！)
找到“WELD FLUX - Definitions of Motion and Flux Parameters for Weld Heat Source”这一行，并修改这一行下面的第5行第三列的值改为3，如下图中的6改为3

图5  修改前

图6 修改后

修改完成后保存。

4）用记事本打开并修改run.bat文件

  
图7 修改前
run.bat中的call后面为你的simufact.welding6.0安装后run_sfMarc.bat的位置
我的上述内容改为：
call "E:\software\special\simufact6\simufact\welding\6.0\solver\simufact\sfMarc\sf_tools\run_sfMarc.bat" -nthread_solver 2 -nthread_elem 2 -j F:\study\softwarestudy\software\simufact\marc_link\marc\Process\_Run_\Process.dat -u F:\study\softwarestudy\software\marc\udf_weldtool\rey.f -save yes


需要说明的是：
-j 后面增加刚才修改的Process.dat的位置
-u 后面增加刚才rey.f文件的位置
增加
-save yes

修改完成后并保存退出

5）命令运行run.bat

                 图8 点击红色框中的黑色小三角并选择“Open a shell”


图9 点击红色框中的黑色小三角并选择“Open a shell”

在出现的对话框中输入“run.bat”并回车





运行一段时间后可以看到Results出现，即可得到想要的结果。
Simufact.Welding 6.0 热源二次开发的小例子-1

热源的脉冲效果如下：
Simufact.Welding 6.0 热源二次开发的小例子-2

  具体内容详见附件的文件“Simufact.Welding 6.0 热源二次开发初步教程.doc”
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
  • 附件 1 : Simufact.Welding_6.0_热源二次开发初步教程.doc
    • 2019-01-22 10:37:49, 3.06 M

» 收录本帖的淘贴专辑推荐

simufact仿真分析

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2019-01-22 10:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangdaqing

新虫 (正式写手)
★
fangdian(金币+1): 谢谢参与
连专门做焊接仿真的软件都要自定义热源？

发自小木虫IOS客户端
赞一下 回复此楼
5楼2019-01-22 11:55:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿七不笨

金虫 (正式写手)
★
fangdian(金币+1): 谢谢参与
楼主有用abaqus做么

发自小木虫Android客户端
赞一下 回复此楼
7楼2019-01-22 19:21:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

美丽人生1940

铜虫 (初入文坛)
★
fangdian(金币+1): 谢谢参与
附件怎么下载呢，觉得资料比较有价值
赞一下 回复此楼
14楼2021-06-10 21:53:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
tzynew2楼
2019-01-22 10:49   回复  
fangdian(金币+1): 谢谢参与
miaojiabing3楼
2019-01-22 11:14   回复  
fangdian(金币+1): 谢谢参与
youngen4楼
2019-01-22 11:26   回复  
fangdian(金币+1): 谢谢参与
发自小木虫Android客户端
nono20096楼
2019-01-22 17:30   回复  
fangdian(金币+1): 谢谢参与
·
miaojiabing8楼
2019-01-23 08:04   回复  
Biotage9楼
2019-01-23 08:46   回复  
fangdian(金币+1): 谢谢参与
sz_kebaoyuan10楼
2019-01-23 08:59   回复  
fangdian(金币+1): 谢谢参与
miaojiabing11楼
2019-01-25 08:17   回复  
youngen12楼
2019-01-25 17:14   回复  
youngen13楼
2019-01-28 13:10   回复  
相关版块跳转 我要订阅楼主 fangdian 的主题更新
141/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请基金代表性成果 +14 lancet0903 2024-05-17 19/950 2024-05-20 22:12 by 星火12
[找工作] 绍兴文理学院怎么样？有没有坑啊 +9 zhaojiang427 2024-05-16 23/1150 2024-05-20 21:44 by 高敖曹
[考博] 海南大学国家级人才田新龙教授团队氢燃料电池与器件方向招收科研助理 (重点支持方向) +6 Dreammzp 2024-05-17 9/450 2024-05-20 18:40 by Dreammzp
[基金申请] 基金委也是打工的 +24 wsjing 2024-05-16 32/1600 2024-05-20 11:29 by dd921ww
[硕博家园] 博三一直没文章怎么办 +24 133456 2024-05-17 35/1750 2024-05-20 11:29 by zlx040628
[硕博家园] 被导师抓着一个褥羊毛该咋办 +12 未来可期简直不� 2024-05-15 24/1200 2024-05-20 11:26 by zlx040628
[教师之家] 今天才知道硕士考进大专比博士进本科工资高得多 +9 ZHONGWU_U 2024-05-14 11/550 2024-05-20 11:24 by 少年太美
[职场人生] 硕士毕业四年工作不顺心 +5 noora2018 2024-05-14 5/250 2024-05-20 10:29 by umvyn
[找工作] 浙江的高校现在门槛都这么高吗 +10 dadqweq_qw 2024-05-16 10/500 2024-05-19 17:02 by wenwen0825
[考博] 想被211以上高校课题组接收 +8 风起沧澜 2024-05-16 10/500 2024-05-19 06:52 by Guo369
[基金申请] 这个教授的科研做得真好 +6 zju2000 2024-05-14 7/350 2024-05-18 15:18 by shl2112501
[基金申请] 河北省基金 +11 3101251516 2024-05-16 18/900 2024-05-18 11:41 by 晓晓爱翠翠
[论文投稿] 中文投稿 +7 liujinyu80 2024-05-13 7/350 2024-05-18 11:09 by dgfj564
[论文投稿] 投稿成功后又想撤回 +5 otani 2024-05-16 5/250 2024-05-17 16:02 by topedit
[电化学] 锂离子电池石墨负极用 1M LiPF6 in DEC:EC=1:1 Vol% 可以吗？ 50+3 fffhhhhh 2024-05-15 8/400 2024-05-17 14:57 by 多点关心多点i
[硕博家园] 五氯化铌怎么溶解啊 +3 南南枝枝 2024-05-17 5/250 2024-05-17 11:37 by ad_fish
[硕博家园] 中国科学院大学李海艳老师课题组诚招硕/博研究生和博士后，诚聘客座本/硕/博学生 +4 ucaszmh 2024-05-14 5/250 2024-05-16 19:51 by 我是小班
[基金申请] 系统导出申报书中参与人那一页的格式，跨入两页，要通过调整行间距调整到一页吗？ +3 danisha 2024-05-15 5/250 2024-05-16 18:21 by danisha
[考博] 深度学习申博 +4 酸菜鱼． 2024-05-15 4/200 2024-05-16 11:12 by 安塔瓦拉多
[论文投稿] 七个月了，还在selecting for review +3 g9522 2024-05-14 6/300 2024-05-14 19:00 by lizhengke06
信息提示
关闭
请填处理意见
关闭