simufactwelding熱源二次開發(fā)初步教程

1、Simufact.Welding 6.0 熱源二次開發(fā)初步教程1 主要目的：開始隨時間變化的焊接熱源，如脈沖熱源 主要目的：開始隨時間變化的焊接熱源，如脈沖熱源2 基本要求： 基本要求：1）脈沖熱源頻率：8Hz2）脈寬 30ms3）其它略3 軟件安裝 軟件安裝1）simufact.welding 6.0 2）VS20103）Intel Visual Fortran Composer XE 20133 基本步驟 基本步驟3.1 3.1 軟

2、件設(shè)置 軟件設(shè)置安裝完之后需要配置 path、Lib 和 include 三個環(huán)境變量。以下列出環(huán)境變量值作為參考（需要根據(jù) VS 和 fortran 的安裝位置不同進(jìn)行相應(yīng)的修改） 。---注：以下部分為參考網(wǎng)上內(nèi)容，需要根據(jù)具體實際修改，如下面采用的是 注：以下部分為參考網(wǎng)上內(nèi)容，需要根據(jù)具體實際修改，如下面采用的是 VS2012，XE2015，根據(jù)需要修改 ，根據(jù)需要修改Path：C:\Program Files (x86)\Mi

3、crosoft Visual Studio 12.0\Common7\ToolsC:\Program Files (x86)\Intel\Composer XE 2015\redist\intel64\mklC:\Program Files (x86)\Intel\Composer XE 2015\bin\intel64C:\Program Files (x86)\Microsoft Visual Studio 12.0\Common7

4、\IDEC:\Program Files (x86)\Windows Kits\8.1\Windows Performance Toolkit\C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\bin\amd64C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\binC:\Program Files (x86)\Mic

5、rosoft SDKs\Windows\v7.0A\Bin\x64C:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\BinC:\Program Files (x86)\Intel\MPI\5.0.1.037\intel64\binC:\Program Files (x86)\Intel\Trace Analyzer and Collector\9.0.1.035\binLib：C:\

6、Program Files (x86)\Intel\Composer XE 2015\compiler\lib\intel64C:\Program Files (x86)\Intel\Composer XE 2015\mkl\lib\intel64C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\lib\amd64C:\Program Files (x86)\Microsoft

7、 Visual Studio 12.0\VC\libC:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\Lib\x64C:\Program Files (x86)\Microsoft SDKs\Windows\v7.1A\Libc mibody(7) internal element numberc mibody(8) layer number for he

8、at transfer shells elementsc and volume fluxc mibody(9) Not usedc mibody(10) boundary condition number if table inputc time timecdimension mibody(*),temflu(*),welddim(*)real*8 fi

9、nteger mibodyreal*8 temflu,time,welddimc* * * * * *c OPEN(UNIT=10,FILE='HELLO.TXT',POSITION='APPEND')a=0.003b=0.003c=0.003v=0.01T0=0.125temp=0.03xis=1zq=ceiling(time/T0)if(abs(MOD(time-temp,T0)).LT.1

10、e-6) thenxis=0c write(10,*) time, 111else if(abs(MOD(time,T0)).LT.1e-6) thenxis=1c write(10,*) time,222else if ((time.LT.(zq*T0)).and.(time.GT.((zq-1)*T0+temp))) thenxis=0c write(10,*) time, 333elsexis=1c write(10,*)

11、 time, 444end ifd=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,*) t