Simufact.Welding 6.0 熱源二次開發(fā)初步教程

./Simufact.Welding6.0熱源二次開發(fā)初步教程1主要目的：開始隨時間變化的焊接熱源,如脈沖熱源2基本要求：1脈沖熱源頻率：8Hz2脈寬30ms3其它略3軟件安裝1simufact.welding6.02VS20103IntelVisualFortranComposerXE20133基本步驟3.1軟件設(shè)置安裝完之后需要配置path、Lib和include三個環(huán)境變量。以下列出環(huán)境變量值作為參考〔需要根據(jù)VS和fortran的安裝位置不同進行相應(yīng)的修改。注：以下部分為參考網(wǎng)上容,需要根據(jù)具體實際修改,如下面采用的是VS2012,XE2015,根據(jù)需要修改Path：C:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\Common7\ToolsC:\ProgramFiles<x86>\Intel\ComposerXE2015\redist\intel64\mklC:\ProgramFiles<x86>\Intel\ComposerXE2015\bin\intel64C:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\Common7\IDEC:\ProgramFiles<x86>\WindowsKits\8.1\WindowsPerformanceToolkit\C:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\VC\bin\amd64C:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\VC\binC:\ProgramFiles<x86>\MicrosoftSDKs\Windows\v7.0A\Bin\x64C:\ProgramFiles<x86>\MicrosoftSDKs\Windows\v7.1A\BinC:\ProgramFiles<x86>\Intel\MPI\5.0.1.037\intel64\bin.035\binLib：C:\ProgramFiles<x86>\Intel\ComposerXE2015\compiler\lib\intel64C:\ProgramFiles<x86>\Intel\ComposerXE2015\mkl\lib\intel64C:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\VC\lib\amd64C:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\VC\libC:\ProgramFiles<x86>\MicrosoftSDKs\Windows\v7.1A\Lib\x64C:\ProgramFiles<x86>\MicrosoftSDKs\Windows\v7.1A\LibINCLUDE：C:\ProgramFiles<x86>\Intel\ComposerXE2015\compiler\includeC:\ProgramFiles<x86>\Intel\ComposerXE2015\mkl\includeC:\ProgramFiles<x86>\MicrosoftVisualStudio12.0\VC\includeC:\ProgramFiles<x86>\MicrosoftSDKs\Windows\v7.1A\Includea3.2常規(guī)建立并設(shè)置simufact.welding工程建立如下常規(guī)的simufact.welding工程〔常規(guī)建立過程省略圖1常規(guī)工程主要坐標(biāo)系為：沿X方向?qū)?0mm,沿y方向板厚10mm,沿z方向長度50mm,焊接沿著z負(fù)方向,焊接起始點坐標(biāo)〔50,10,-50[單位mm],建立對稱的一半模型。3.3建立.f子函數(shù)subroutineuweldflux<f,temflu,mibody,welddim,time>cf<1>fluxvalue<tobedefined>cf<2>derivativeoffluxwithrespecttotemperaturec<tobedefined;optional,thismightimprovethecconvergencebehavior>cctemflu<1>estimatedtemperaturectemflu<2>previousvolumetricfluxctemflu<3>temperatureatbeginningofincrementctemflu<4,5,6>integrationpointcoordinatescmibody<1>elementnumbercmibody<2>fluxtypecmibody<3>integrationpointnumbercmibody<4>fluxindex-notusediftableinputcmibody<5>notusedcmibody<6>=1:heattransferc=2:joulec=3:bearingc=4:electrostaticc=5:magnetostaticc=6:acousticcmibody<7>internalelementnumbercmibody<8>layernumberforheattransfershellselementscandvolumefluxcmibody<9>Notusedcmibody<10>boundaryconditionnumberiftableinputctimetimecdimensionmibody<*>,temflu<*>,welddim<*>real*8fintegermibodyreal*8temflu,time,welddimc******cOPEN<UNIT=10,FILE='HELLO.TXT',POSITION='APPEND'>a=0.003b=0.003c=0.003v=0.01 T0=0.125 temp=0.03 xis=1zq=ceiling<time/T0> if<abs<MOD<time-temp,T0>>.LT.1e-6>then xis=0c write<10,*>time,111 elseif<abs<MOD<time,T0>>.LT.1e-6>thenxis=1c write<10,*>time,222elseif<<time.LT.<zq*T0>>.and.<time.GT.<<zq-1>*T0+temp>>>thenxis=0cwrite<10,*>time,333elsexis=1c write<10,*>time,444 endifd=0.05-v*timeq=7000*0.7pi=3.14x=temflu<4>y=temflu<5>z=temflu<6>mibody<6>=1y0=0.009x0=0.05heat=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*xisc write<10,*>time,xis,fc close<10>returnend并存為文件名rey.f33.4開始子函數(shù)運行1點擊常規(guī)項目的運行按鈕圖2點擊運行按鈕點擊后出現(xiàn)運行對話框,不用管它,直接關(guān)閉。圖3運行對話框〔直接關(guān)閉2此時找到你常規(guī)工程下的_Run_文件夾位置并打開該文件夾圖4該工程的_Run_文件夾會發(fā)現(xiàn)該文件夾下面有Process.dat文件和run.bat文件。3用記事本方式打開并修改Process.dat文件<此步不可少！>找到"WELDFLUX-DefinitionsofMotionandFluxParametersforWeldHeatSource"這一行,并修改這一行下面的第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_solver2-nthread_elem2-jF:\study\softwarestudy\software\simufact\marc_link\marc\Process\_Run_\Process.dat-uF:\study\softwarestudy\software\marc\udf_weldtool\rey.f-saveye