Simufact.Welding 6.0 热源二次开发初步教程

1、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 20133 基本步骤3.1 软件设置安装完之后需要配置path、Lib和include三个环境变量。以下列出环境变量值作为参考（需要根据VS和fortran的安装位置不同进行相应的修改）。-注：以下部分为参考网上内容，需要根据具体实际修改，如下面采用的是VS2012，X

2、E2015，根据需要修改Path：C:Program Files (x86)Microsoft Visual Studio 12.0Common7ToolsC:Program Files (x86)IntelComposer XE 2015redistintel64mklC:Program Files (x86)IntelComposer XE 2015binintel64C:Program Files (x86)Microsoft Visual Studio 12.0Common7IDEC:Program Files (x86)Windows Kits8.1Windows Performan

3、ce ToolkitC:Program Files (x86)Microsoft Visual Studio 12.0VCbinamd64C:Program Files (x86)Microsoft Visual Studio 12.0VCbinC:Program Files (x86)Microsoft SDKsWindowsv7.0ABinx64C:Program Files (x86)Microsoft SDKsWindowsv7.1ABinC:Program Files (x86)IntelMPI5.0.1.037intel64binC:Program Files (x86)Intel

4、Trace Analyzer and Collector9.0.1.035binLib：C:Program Files (x86)IntelComposer XE 2015compilerlibintel64C:Program Files (x86)IntelComposer XE 2015mkllibintel64C:Program Files (x86)Microsoft Visual Studio 12.0VClibamd64C:Program Files (x86)Microsoft Visual Studio 12.0VClibC:Program Files (x86)Microso

5、ft SDKsWindowsv7.1ALibx64C:Program Files (x86)Microsoft SDKsWindowsv7.1ALibINCLUDE：C:Program Files (x86)IntelComposer XE 2015compilerincludeC:Program Files (x86)IntelComposer XE 2015mklincludeC:Program Files (x86)Microsoft Visual Studio 12.0VCincludeC:Program Files (x86)Microsoft SDKsWindowsv7.1AInc

6、ludea3.2 常规建立并设置simufact.welding工程 建立如下常规的simufact.welding工程（常规建立过程省略） 图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 te

7、mperaturec (to be defined; optional, this might improve thec convergence behavior)cc temflu(1) estimated temperaturec temflu(2) previous volumetric fluxc temflu(3) temperature at beginning of incrementc temflu(4,5,6)integration point coordinatesc mibody(1) element numberc mibody(2) flux typec mibody

8、(3) integration point numberc mibody(4) flux index - not used if table inputc mibody(5) not usedc mibody(6) =1 : heat transferc =2 : joulec =3 : bearingc =4 : electrostaticc =5 : magnetostaticc =6 : acousticc mibody(7) internal element numberc mibody(8) layer number for heat transfer shells elements

9、c and volume fluxc mibody(9) Not usedc mibody(10) boundary condition number if table inputc time timec dimension mibody(*),temflu(*),welddim(*) real*8 f integer mibody real*8 temflu,time,welddimc* * * * * *c OPEN(UNIT=10,FILE='HELLO.TXT',POSITION='APPEND') a=0.003 b=0.003 c=0.003 v=0

10、.01T0=0.125temp=0.03xis=1 zq=ceiling(time/T0)if(abs(MOD(time-temp,T0).LT.1e-6) thenxis=0cwrite(10,*) time, 111else if(abs(MOD(time,T0).LT.1e-6) then xis=1cwrite(10,*) time,222 else if (time.LT.(zq*T0).and.(time.GT.(zq-1)*T0+temp) then xis=0c write(10,*) time, 333 else xis=1cwrite(10,*) time, 444end

11、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*xiscwrite(10,*) time, xis,fcclose(10) return end并存为文件名rey.f33.4 开始子函数运行1）点击常规项目的运行按钮图2 点击运行按钮点击后出现运行对话框，不用

12、管它，直接关闭。图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.b

13、at中的call后面为你的simufact.welding6.0安装后run_sfMarc.bat的位置 我的上述内容改为：call "E:softwarespecialsimufact6simufactwelding6.0solversimufactsfMarcsf_toolsrun_sfMarc.bat" -nthread_solver 2 -nthread_elem 2 -j F:studysoftwarestudysoftwaresimufactmarc_linkmarcProcess_Run_Process.dat -u F:studysoftwarestudysoftwaremarcudf_weldtoolrey.f -save yes 需要说明的是：-