高斯热源3D打印APDL
多组元金属粉末选区激光烧结三维瞬态温度场模拟
多组元金属粉末选区激光烧结三维瞬态温度场模拟
沈以赴;顾冬冬;李守卫;李昊
【期刊名称】《南京航空航天大学学报》
【年(卷),期】2008(040)005
【摘要】利用Ansys有限元软件,建立了多组元金属粉末选区激光烧结过程中三维瞬态温度场计算模型.在考虑了相变潜热、辐射对流和随温度变化的热物性参数条件下,使用APDL参数化语言实现了高斯热源的施加和热源按一定速率的移动.通过该计算模型可掌握成形过程中烧结温度场随时间的变化规律,进而控制多组元金属粉末选区激光烧结成形机制,避免"球化"效应、翘曲变形等缺陷,为合理选取激光工艺参数提供理论依据.对Ni和Cu-10Sn多组元混合粉末进行了选区激光烧结实验,验证了模拟结果的正确性.
【总页数】6页(P611-616)
【作者】沈以赴;顾冬冬;李守卫;李昊
【作者单位】南京航空航天大学材料科学与技术学院,南京,210016;南京航空航天大学材料科学与技术学院,南京,210016;南京航空航天大学材料科学与技术学院,南京,210016;南京航空航天大学材料科学与技术学院,南京,210016
【正文语种】中文
【中图分类】TG14;TG665
【相关文献】
1.多组分铜基金属粉末选区激光烧结试验研究 [J], 顾冬冬;沈以赴;杨家林;王洋
2.多组分铜基金属粉末选区激光烧结致密化机理 [J], 顾冬冬;沈以赴;杨家林;王洋
3.钛粉选区激光烧结三维瞬态温度场模拟 [J], 段亮亮;党新安;杨立军;张勇
4.选择性激光烧结金属粉末瞬态温度场模拟 [J], 任继文;殷金菊
5.选区激光烧结瞬态温度场模拟与测试方法研究 [J], 邢键;孙晓刚;高益庆
因版权原因，仅展示原文概要，查看原文内容请购买。
逆向工程下的快速成型数据精度分析
逆向工程下的快速成型数据精度分析计三有;尹业腾;陈东【摘要】Put forward three-phase error analysis method for rapid prototyping based on reverse engineering technology.The accuracy factors of reverse engineering and rapid prototyping are analyzed.Based on error synthesis theory,respectively evaluate reverse engineering precision,rapid prototyping and whole precision and give evaluate calculation formula.%提出基于误差合成理论的逆向工程下快速成型技术的“三段式”精度评价方法.分析逆向工程及快速成型的精度影响因素.基于误差合成理论分别进行逆向工程精度评价、快速成型精度评价、整体精度评价,得出精度评价计算公式.【期刊名称】《自动化与仪表》【年(卷),期】2018(033)005【总页数】4页(P5-8)【关键词】快速成型;精度分析;逆向工程;增材制造【作者】计三有;尹业腾;陈东【作者单位】武汉理工大学物流工程学院,武汉430063;武汉理工大学物流工程学院,武汉430063;武汉理工大学物流工程学院,武汉430063【正文语种】中文【中图分类】TH16目前，针对逆向工程与快速成型集成应用的数据处理与精度分析研究相对较少，未对逆向工程和快速成型集成应用下的精度进行整体评价，对于两种技术的数据交换方式研究相对较少。
逆向工程中精度分析可在曲面拟合时进行，将测量点的相应位置和测量值对比得到点的偏差。
文献[1]研究出基于逆向工程的快速成型评价的方法，开发逆向工程系统对工件进行虚拟加工，将虚拟加工结果与原型对比，用表面积量化差异，从而对快速成型的加工工艺进行了评价。
ansys高斯热源施加
!平板对接开v型坡口多道焊成型(温度场) !采用guass表面热源finish/clear,nostart/filname,1.11thermal_analysis/title,3D_weld_thermal_analysis/units,si !定义国际制单位/prep7!**************************************************!定义焊接尺寸及焊接参数!**************************************************!******************************************构件尺寸len=100e-3 !长0.1mwid=50e-3 !宽0.05mthic=6e-3 !厚0.006m!b=0.002 !根部预留(m)!h=0 !钝边(m)wweld=15e-3 !焊接热影响区宽度lnum=10 !沿焊缝长度方向划分份数!******************************************焊接参数I=200 !焊接电流(A)U=30 !焊接电压(V)v=0.5/60 !焊接速度(m/s)arf=60 !坡口角度（°）pi=3.1415926 !pi为圆周率eta=0.7 !焊接热效率,手工电弧焊0.7,埋弧焊0.8K=1.2e4 !热流集中程度系数(/m^2) ?them0=25 !环境温度及母材初始温度,取室温25℃!**********************************************定义高斯表面移动热源参数!GAUSS热源模型qr=qmax*exp(-k*r*r) !(单位W/m2)!假设热源中心坐标为(a,b,c),a=0,b=每个焊缝的最高点的坐标y,c=v*dt!每个焊缝上表面任一点的坐标为(x,y,z) ,x>0, y>0, z>0!r^2=x^2+(y-b)^2+(z-c)^2，每个焊缝上表面任一点距热源中心距离qmax=0.24*k*eta*I*u/pi !加热斑点中心热流(W/m2)percent=0.99radius=abs(sqrt(1/k*log(1/(1-percent)))) !有效热半径，单位为m!**********************************************************************!定义单元类型及材料属性!**********************************************************************et,1,solid70 !计算单元类型!nothing !定义实常数mptemp,1,0,100,200,300,400,500 !定义材料属性温度范围mptemp,,1200,1250,1350,1400,1450,1500mptemp,,1505,1705,1905,2105,2305,2500mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定义材料密度(kg/m3)mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K))mpdata,kxx,1,,30,30,30,30,30,30mpdata,kxx,1,,30,30,30,30,30,30mpdata,c,1,1,450,469,481,508.5,536,569 !考虑相变潜热的比热容(J/kg*K)mpdata,c,1,,700,2172,5116,6589,8061,9533mpdata,c,1,,9533,7757,5982,4206,2431,700!****************************************************************!建立有限元计算模型!*****************************************************************afun,deg !如果涉及到角度，将使用“度”(系统默认是弧度)K,1,0,0,0K,2,wweld,0,0K,3,wid,0,0K,4,wid,thic,0K,5,thic*tan(arf/2)+wweld,thic,0K,6,thic*tan(arf/2),thic,0K,7,0,0,lenA,1,3,4,6A,1,2,5,6cyl4,0,0,thic/cos(arf/2),arf,0,90cyl4,0,0,thic/cos(arf/2)*sqrt(2)/2,arf,0,90L,1,7 !形成拖拉路径L14 Vdrag,all,,,,,,14 !沿线拖拉面Vptn,all !分割体Vglue,all !粘贴体？Nummrg,all!*************************************************************!划分网格!*************************************************************aadd,24,32 !面相加，以便生成映射网格lccat,40,49lccat,42,52!*************************************沿焊缝z方向划分，均为2mmlsel,s,line,,17,21,2lsel,a,line,,24,26,2lsel,a,line,,35,36lsel,a,line,,14,30,16lesize,all,2e-3!***************************************沿厚度y方向划分1-1.5mmlsel,s,line,,6,25,19lsel,a,line,,2,18,16lesize,all,,,6lsel,s,line,,41,43,2lsel,a,line,,50,53,3lesize,all,1.5e-3!***********************************************沿x方向划分1mmlsel,s,line,,45,47,2 !母材lsel,a,line,,44,46,2lesize,all,,,18,0.8lsel,s,line,,7,27,20 !母材lesize,all,,,12,0.7lsel,a,line,,5,23,18lesize,all,,,12,10/7lsel,s,line,,11,34,23 !焊缝lsel,a,line,,48,51,3lesize,all,,,4!****************************************划分体，生成六面体单元allsmshkey,1 !映射网格mshape,0,3d !六面体形状vsel,allvatt,,,1aslv,saatt,,,1type,1 !单元类型为1vmesh,all!*******************************************存档备份有限元模型allssave,'1.11thermal_analysis_meshed','db'finish!================================================================= !重新载入有限元模型进行加载求解!================================================================= /solu/view,,-3,2,-5 !调整窗口视角!*********************************************************************** ! 杀死焊缝区单元!***********************************************************************v1=5 !焊缝1所在的体v2=7 !焊缝2所在的体a1=18 !焊缝1的上表面a2=31 !焊缝2的上表面esel,s,type,,1 !选取1类单元*get,nemax,elem,,num,max !得到所有1类单元号码的最大值*get,nemin,elem,,num,min!********************************************分段杀死第2道焊缝vsel,s,,,v2eslvekill,allesel,s,liveeplot!***************************************************************!边界条件(暂时没有考虑对流)!***************************************************************nsel,s,loc,x,0 !在x=0处施加对称边界条件dsym,symm,x!*****************************************************************!设置非线性求解选项!*****************************************************************allsantype,trans !瞬态分析trnopt,full !瞬态分析选项，完全分析nropt,full,,on !定义完全牛顿--拉普森方法，激活自适应下降（默认）pred,on !打开预测校正timint,on !时间积分设置tintp,0.005,,,1,0.5,0.2 !定义瞬态综合参数tref,25 !参考温度25摄氏度!***************************************************************!设置载荷步参数!***************************************************************t=0 !求解时间初始值dt0=1e-6 !建立初始条件的小时间段dt1=0.5 !起始加热点对应的时间tinc=len/(lnum*v) !载荷步时间间隔,T =len/v,tinc =T/lnum!tsub=0.2 !求解温度场时子步时间!***************************************************************!稳态分析确定初始温度场!***************************************************************!初始条件time,dt0 !loadstep=1,time=dt0timint,off !时间积分选项，关闭（稳态分析）kbc,1 !阶越加载!nsubst,1 !子步为1步esel,s,type,,1 !选择类型1的单元节点nsle,sic,all,temp,25 !初始状态，温度25摄氏度allsoutres,all,allsolvesave,'1.11thermal_analysis_initial_temp','db'!查看初始温度场，应该全为25摄氏度!***************************************************************!按顺序模拟温度场，以第一道为例!***************************************************************!**************************************************第1道焊缝填充esel,s,type,,1 !在第一类单元中加热流载荷esel,r,liveeplot*do,im,0,lnum,1!****************************选择有效半径内节点为当前节点vsel,s,,,v1 !选择焊缝1上表面单元eslv,rnsleasel,s,,,a1nsla,r,1esln,r,0!热源加载, 将各段后点（0-lnum）为热源中心,先加载后删除热源!即先消除上段所加高斯热源,并将上段的温度值作为下段的初始值!*******************************************包括起始加热点c=V*tinc*im !热源位置b=thic/cos(arf/2)*sqrt(2)/2tm=dt1+im*tinctime,tmantype,4,resttimint,onautots,onkbc,1deltim,0.01,0.01,0.2 !时间步长!***********************************以下为施加热流密度载荷*do,i,nemin,nemax,1*if,esel(i),eq,1,then !如果该单元在上面的集合里，则xsy=centrx(i) !读取该单元的中心坐标ysy=centry(i)zsy=centrz(i)rr=abs(sqrt(xsy*xsy+(ysy-b)*(ysy-b)+(zsy-c)*(zsy-c))) !该单元中心距离热源中心*if,rr,le,radius,then !该单元中心在加热半径范围内qr=qmax*exp(-k*rr*rr) !该单元中心处的热流大小sfe,i,nmface(i),hflux,,qr !在每个单元所指定的面上施加热源载荷*endif*endif*enddoallsel !全选求解outres,all,allsolveesel,s,type,,1esel,r,livevsel,s,,,v1 !选择焊缝1上表面单元eslv,rnsleasel,s,,,a1nsla,r,1esln,r,0!!!!!!!!!!!!!!输入先暂时停止一下!插入"/post1$plnsol,temp$/psf,hflux,,1,1,on$/replot"的内容进行查看!问题：热流载荷的分布和温度场的分布不一致（该完全一致才对），有点偏移？!若要继续，读入"finish$/solu"*do,i,1,6sfedele,all,i,hflux !删除每个单元六个面上的热载荷*enddo*enddoallsesel,s,liveeplot !查看温度场save,'1.11thermal_analysis_alive1','db'!=================================================================!生成节点温度分布云图动画文件的程序段!=================================================================/post1/seg,dele !允许图形数据存储到局部终端存储器/dscale,1,1.0 !显示位移时，设置放大系数avprin,0,0 !规定主项量和总向量如何被计算avres,1 !规定当powergraphies可用时，结果数据如何被平均/seg,multi,1.11temp-10,0.2 !存动画,分节存储后来的显示,动画名字,每个画面间延迟系数esel,s,live!*************************************************读取初始稳态温度场结果set,,,1,,dt0 !从结果文件读,载荷步,子步,比例因子实部,比例因子虚部,时间plnsol,temp!*****************************************读取第1道焊缝过程中的温度场结果*do,im,0,lnum,1tm=dt1+im*tincset,,,1,,tmplnsol,temp*enddo/seg,off,1.11temp-10,0.2anim,1,1 !*********************************生成动画完毕!****************************************存储第1道焊缝过程中的温度场结果图片esel,s,live*do,im,0,lnum,1*if,im,eq,0,thentm=dt1set,,,1,,tmplnsol,temp/image,save,1.11temp0,jpeg*endiftm=dt1+im*tincset,,,1,,tmplnsol,temp*if,tm,eq,dt1+lnum/2*tinc,thenset,,,1,,tmplnsol,temp6.5,jpeg*endif*enddo/image,save,1.11temp12.5,jpeg。
金属激光3D打印过程数值模拟应用及研究现状
材料工程Journal of Materials Engineering第4 9卷 第4期2021年4月第52-62页Vol. 4 9 No. 4Apr. 2021 pp. 52―62金属激光3D 打印过程数值 模拟应用及研究现状Application and research status of numerical simulation of metallaser 3D printing process杨 鑫1，王 犇】，谷文萍2,张兆洋】，刘世锋3,武 涛1(1西安理工大学材料科学与工程学院，西安710048；2长安大学材料科学与工程学院，西安710061 ；3西安建筑科技大学冶金学院，西安710055) YANG Xin 1, WANG Ben 1 ,GU Wen-ping 2 , ZHANG Zhao-yang 1 , LIU Shi-feng 3 ,WU Tao 1(1 Department, of Materials Science and Engineering ,Xi ?an University ofTechnology, Xi an 71 0048, China ； 2 Department, of Materials Scienceand Engineering , Chang ? an University , Xi ? an 710061 , China ;3 School of Metallurgical and Engineering ,Xi ?an Universityof Architecture & Technology, Xi an 71 0055, China)摘要：数值模拟可以高效、有针对性地对金属激光选区熔化成型过程中的温度场、熔池形状、残余应力和变形、凝固过程 微观组织演变等过程建立相应的模型并对成形件的相关性能做出准确预测，为工艺优化提供科学的依据，显著降低工艺开发成本和缩短工艺开发周期，有力推动金属增材制造向工业级应用的转变。
间接选区激光多道烧结的三维温度场模拟
间接选区激光多道烧结的三维温度场模拟
贾蕾;王典;葛小兵
【期刊名称】《机械设计与制造工程》
【年(卷),期】2013(042)006
【摘要】为了得到间接选区激光烧结的温度场分布规律,基于ANSYS参数化设计语言APDL(ANSYS Parametric Design Language),建立了单层多道烧结的有限元模型.通过分析热传导方程、初始条件和边界条件,选取高斯热源模型来模拟激光加热.考虑到粉末的热物性参数随温度的非线性变化,在激光功率14 W、扫描速度2000mm/s、预热温度100℃、扫描间距0.1 mm、光斑直径0.2mm、铺粉厚度0.1 mm的条件下,进行了温度场的数值模拟.模拟结果显示:热源中心前方的温度梯度比后方的大；随着烧结的进行,热源中心的最高温度逐渐变大.
【总页数】3页(P17-19)
【作者】贾蕾;王典;葛小兵
【作者单位】江苏大学机械工程学院,江苏镇江212013;江苏大学机械工程学院,江苏镇江212013;江苏大学机械工程学院,江苏镇江212013
【正文语种】中文
【中图分类】TF124
【相关文献】
1.间接性选区激光烧结铁粉材料及成型工艺研究 [J], 刘越华;张坚;徐志锋;郑海忠
2.直接金属选区激光多道烧结温度场有限元模拟 [J], 沈显峰;王洋;姚进;杨家林
3.不锈钢316L间接选区激光烧结工艺参数优化研究 [J], 任乃飞;葛小兵;王权;李阳;赵新宏;杭雅慧;徐玮
4.聚合物选区激光烧结粉末加热模拟 [J], 蔡令令;丁浩亮;孟娟;温馨;严波
5.聚合物选区激光烧结粉末熔合介观模拟 [J], 蔡令令;丁浩亮;孟娟;温馨;严波
因版权原因，仅展示原文概要，查看原文内容请购买。
基于Sysweld的焊接接头热源模型二次开发
1．2 热源模型 Sysweld 软件内置了 3 种结构简单的热源模型［6］，
20 ·试验与研究·
分别为：
（1） 二维高斯面热源模型， 适用于表面热处理；
（2） 三维高斯锥形热源 模 型 ， 适 用 于 激 光 焊 、
电子束焊等高能束流焊接；
（3） 三维双椭球热源模型， 适用于 TIG， MIG 焊
1．3 热源计算
基于 Sysweld 平台开发自定义的异型焊 接 接 头 ，
建立三维有限元模型， 其单元类型、 单元组的数量
和种类必须和内置的热源模型一致， 具体包括：
（1） 建立包含焊缝和母材的 3D 单元组； （2） 从 3D
单元抽取表面网格， 生成 2D 单元组， 表征换热面；
（3） 从焊缝和母材的 3D 单元沿着焊接方向抽取 2 条
摘要： 焊接接头的结构形式和尺寸精度直接影响热源模型计算结果的准确性 ， 工程实际应用的焊接接头形式多样、 结构复杂， Sysweld
内置的焊接接头远不能满足实际需求。 基于 Hypermesh， Sysweld 软件平台开发了大坡口角焊缝和双侧多层角焊缝热源模型， 首先利用
Hypermesh 建立有限元模型生成可执行的内嵌文件， 利用 Sysweld 的 HSF 工具反复调整 双 椭 球 热 源 高 斯 参 数 ， 并 将 校 核 结 果 与 试 验 结
1 大坡口角焊缝的热源计算 1．1 几何模型
焊接结构的接头形式和形状多种多样， 只利用 专用商业焊接有限元软件所提供的自定义接头往往 不够， 这就要求必须根据工程实际自定义接头形式 和形状。 如图 1 所示的大坡口角焊缝， 其坡口位置、 坡口角度直接影响传热过程， 继而影响焊接结构件 最终的变形量和组织性能。
完整word版,ANSYS焊接高斯热源APDL参考.doc
不准/UNITS,SI/CONFIG ,NRES,10000/PREP7ET,1,SOLID70MP,DENS,1,7930MP,C,1,502mptemp,1,20,100,500mpdata,kxx,1,1,12.1,16.3,21.4/VIEW,1,1,1,1!*************定义基板尺寸***********lx=0.1lz=0.1ly=0.006!************* 定义焊接参数***********weld=0.01 ! 每道焊道宽度wheight=0.004 !焊道的高度LSIZE=0.005V=0.008!速度pi=3.1415926U=25I=180Q=U*I!电源有效功率R=0.006! 电源有效热半径effect=0.8Qmax=effect*Q/(pi*R*R) !中心处最大热流tinc=LSIZE/V! 每小段的时间间隔!*****************建立模型，生成网格************** block,0,lx,0,ly,0,lzwpoff,0,ly,lz/2block,0,lx,0,wheight,-weld/2,weld/2wpoff,0,0,-weld/2vsbw,all,,deletewpoff,0,0,weldvsbw,all,,deletevglue,alllsel,s,,,plsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35lsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35lesize,all,0.005,,,0.2,,,,1lsel,s,line,,9,10lsel,a,line,,12lsel,a,line,,38lesize,all,0.005,,,0.2,,,,1lsel,s,line,,2lsel,a,line,,4,5lsel,a,line,,7lsel,a,line,,14lsel,a,line,,16lsel,a,line,,19lsel,a,line,,39lesize,all,0.002,,,1,,,,1lsel,s,line,,23,24lsel,a,line,,41,42lsel,a,line,,45lsel,a,line,,47lsel,s,line,,23,24lsel,a,line,,41,42lsel,a,line,,45lsel,a,line,,47lesize,all,0.002,,,1,,,,1 vmesh,all/PNUM,DEFAEPLOTfini/soluantype,trans,new!nlgeom,ontimint,0,structtimint,1,thermtimint,0,magtimint,0,electtref,25nropt,autoautos,onkbc,0pred,onlnsrch,on!************杀死焊缝去单元***********nsel,s,loc,y,ly,ly+wheightesln,s,1ekill,allallsel,allesel,s,liveeplot!**********施加对流载荷*************esel,s,livensel,s,loc,z,0nsel,a,loc,z,lznsel,a,loc,x,0nsel,a,loc,x,lxsf,all,conv,10,25!************定义数组维数 **********MAX_X=1+lx/LSIZE!*************定义 table 数组 ****************************************** *do,i,1,MAX_X,1local,12,0,(i-1)*lsize,ly,lz/2,,,,,,*del,_FNCNAME*del,_FNCMTID*del,_FNC_C1*del,_FNC_C2*del,_FNC_C3*del,_FNC_C4*del,_FNCCSYS*del,'heatflux'!定义表格各行*set,_FNCNAME,'heatflux'*dim,_FNC_C1,,1*dim,_FNC_C2,,1*dim,_FNC_C3,,1*dim,_FNC_C4,,1!表格各行赋值*set,_FNC_C1(1),effect*set,_FNC_C2(1),U*set,_FNC_C3(1),I*set,_FNC_C4(1),R*set,_FNCCSYS,12*DIM,%_FNCNAME%,TABLE,6,26,1,,,,%_FNCCSYS%!!Begin of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) %_FNCNAME%(0,0,1)= 0.0, -999 %_FNCNAME%(2,0,1)= 0.0%_FNCNAME%(3,0,1)= %_FNC_C1(1)% %_FNCNAME%(4,0,1)= %_FNC_C2(1)% %_FNCNAME%(5,0,1)= %_FNC_C3(1)% %_FNCNAME%(6,0,1)= %_FNC_C4(1)% %_FNCNAME%(0,1,1)= 1.0, -1, 0, 3, 0, 0, 17%_FNCNAME%(0,2,1)= 0.0, -2, 0, 1, -1, 3, 17%_FNCNAME%(0,3,1)=0, -1, 0, 1, -2, 3, 18%_FNCNAME%(0,4,1)= 0.0, -2, 0, 1, -1, 3, 19%_FNCNAME%(0,5,1)= 0.0, -1, 0, 0, 0, 0, 0%_FNCNAME%(0,6,1)= 0.0, -3, 0, 1, 0, 0, -1%_FNCNAME%(0,7,1)= 0.0, -4, 0, 1, -1, 2, -3%_FNCNAME%(0,8,1)= 0.0, -1, 0, 3, 0, 0, -4%_FNCNAME%(0,9,1)= 0.0, -3, 0, 1, -4, 3, -1%_FNCNAME%(0,10,1)= 0.0, -1, 0, 2, 0, 0, 2%_FNCNAME%(0,11,1)= 0.0, -4, 0, 1, 2, 17, -1%_FNCNAME%(0,12,1)= 0.0, -1, 0, 2, 0, 0, 4%_FNCNAME%(0,13,1)= 0.0, -5, 0, 1, 4, 17, -1%_FNCNAME%(0,14,1)= 0.0, -1, 0, 1, -4, 1, -5%_FNCNAME%(0,15,1)= 0.0, -4, 0, 1, -3, 3, -1%_FNCNAME%(0,16,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,17,1)= 0.0, -3, 0, 1, 20, 17, -1%_FNCNAME%(0,18,1)= 0.0, -1, 0, 1, -4, 4, -3%_FNCNAME%(0,19,1)= 0.0, -1, 7, 1, -1, 0, 0%_FNCNAME%(0,20,1)= 0.0, -3, 0, 1, -2, 3, -1%_FNCNAME%(0,21,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,22,1)= 0.0, -2, 0, 1, 20, 17, -1%_FNCNAME%(0,23,1)= 0.0, -1, 0, 3.14159265358979310, 0, 0, -2 %_FNCNAME%(0,24,1)= 0.0, -4, 0, 1, -1, 3, -2%_FNCNAME%(0,25,1)= 0.0, -1, 0, 1, -3, 4, -4%_FNCNAME%(0,26,1)= 0.0, 99, 0, 1, -1, 0, 0!End of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) !-->!**********激活单元 *********esel,s,livensel,s,loc,x,(i-1)*lsize-0.002,(i-1)*lsize+0.002nsel,r,loc,z,-weld/2,weld/2nsel,r,loc,y,0,wheightesln,s,0ealive,allallsel,all!**********施加热流载荷********** esel,s,livensel,s,loc,Y,wheightsf,all,hflux,%heatflux%allsel,alltime,i*tincnsubst,2SOLVE!*************删除热流载荷*********** nsel,s,loc,y,wheightsfdele,all,hfluxallsel,allesel,s,liveeplotOUTRES,ALL,ALL,/PSF,HFLUX,,2/REPLOTOUTRES,ALL,ALL,。
高斯移动热源作用下的选区激光熔化近场动力学模型
高斯移动热源作用下的选区激光熔化近场动力学模型高斯移动热源作用下的选区激光熔化近场动力学模型1. 引言选区激光熔化（Selective Laser Melting，SLM）是一种具有广泛应用前景的先进制造技术，被广泛应用于快速原型制造、复杂结构制造以及定制产品制造等领域。
在SLM过程中，高功率密度的激光束将粉末材料选区加热到熔点，从而实现材料的局部熔化和熔固交替，最终形成一个完全密实的三维构件。
然而，SLM过程中复杂的相变和热传导现象往往导致部件内部的应力和形变问题。
为了更好地理解和优化SLM过程，科学家们提出了不同的模型和方法，其中包括高斯移动热源作用下的选区激光熔化近场动力学模型。
2. 高斯移动热源作用下的选区激光熔化动力学模型在SLM过程中，激光热源在焊接区域内移动，通过激光辐射与材料的相互作用，使材料局部加热并熔化。
高斯移动热源是描述激光热输入分布的数学模型之一，它通过高斯函数来表示热源的温度分布。
基于高斯移动热源模型，研究者们建立了选区激光熔化的动力学模型，用于解释和预测SLM过程中的熔化和凝固行为。
3. 动力学模型的基本原理高斯移动热源作用下的选区激光熔化近场动力学模型是基于质量守恒、能量守恒和动量守恒等基本原理建立的。
其中，质量守恒方程描述了选区激光熔化过程中粉末材料的熔化和凝固行为，能量守恒方程考虑了激光辐射和传导的热效应以及相变潜热的释放，动量守恒方程描述了热传导引起的流体流动效应。
综合考虑这些守恒方程，可以建立一个全面的选区激光熔化动力学模型。
4. 模型数学表达根据高斯移动热源模型，可以得到选区激光熔化动力学模型的数学表达式。
在该模型中，涉及到材料的熔化速率、温度分布、热应力分布以及热传导等关键参数。
通过对这些参数进行数值模拟和分析，可以更好地理解和预测SLM过程中的热、流和固态物理现象，为优化SLM过程提供理论支持。
5. 对模型的应用和优化高斯移动热源作用下的选区激光熔化近场动力学模型被广泛应用于SLM过程的设计和预测中。
金属激光立体成形热应力场参数化有限元模型_马良
第36卷 第12期中 国 激 光V ol.36,No.122009年12月CH INESE JOU RNAL OF LASERSDecember,2009文章编号:0258-7025(2009)12-3226-07金属激光立体成形热应力场参数化有限元模型马 良 黄卫东 于 君 王 波 杨海欧(西北工业大学凝固技术国家重点实验室,陕西西安710072)摘要 针对金属材料的激光立体成形(M LSF)工艺,利用大型有限元分析工具A nsy s 的二次开发语言A PDL ,开发了成形过程温度场、应力场的参数化有限元模型。
该模型综合考虑了随温度变化的材料非线性、高斯激光能量分布、对流/辐射换热边界条件、相变以及自由变形约束等一系列问题。
通过使用过渡网格划分技术,在提高计算精度的基础上,大幅减少了单元数目,从而实现了金属激光立体成形过程的整体建模。
采用移动热源和单元生死技术,对激光成形过程热应力场进行了有效仿真,在准确计算温度场演化规律的基础上,揭示了塑性压缩区、塑性拉伸区、卸载区等热应力场产生的原因。
使用该参数化模型便于研究金属激光立体成形工艺条件、不同材料等对成形过程热应力场的影响。
关键词 激光技术;金属激光立体成形;热应力场;有限元仿真中图分类号 T N 249;T G665 文献标识码 A doi:10.3788/CJL 20093612.3226Parametric Finite Eleme nt Model of Temperature/Stre ss Field Evolutionby Metal Laser Solid FormingMa Liang Huang Weidong Yu Jun Wang Bo Yang Haio u(St at e Key La bor a tor y of Solidif ica t ion Pr ocessing ,Nor t hwest er n Poly techn ical Un iver sity ,Xi c an,Shaanxi 710072,China )Abstract Based on Ansys secondary development language APDL,a parametric temperature/stress finite element m odel for metal laser solid forming (MLSF)proc ess was developed.The m odel a llows for non -linear behavior of m ateria l properties,Gaussian distribution of laser energy,thermal convection and radiation boundary c ondition,phase change,free deformation constraints c ondition et a l.A full finite element model of ML SF process is established by transition mesh tec hnology.Accuracy of simula tion is improved and the element number is reduced under the m odel.The temperature/stress evolut ion of ML SF proc ess is simulated based on moving laser beam and element birth and death technology.On the basis of the temperature distribution is evolution,the reason of inducing the compress plastic stress,tensile plastic stress and unloading zones is analyzed.The influence of ML SF proc ess c onditions on temperature/stress can be studied by utilizing the parametric finite elem ent model in order to optimize the MLSF process.Key wo rds la ser technique;metal laser solid forming;temperature/stress field;finite elem ent simulation收稿日期:2009-01-08;收到修改稿日期:2009-04-12 基金项目:国家自然科学基金(50331010)资助项目。
三维移动高斯分布表面热源焊接例子
/prep7
/UNITS,si
et,1,solid70
mp,c,1,544.28
mp,kxx,1,15.06
mp,dens,1,4500
block,0,30,0,15,0,3
esize,1
vmesh,all
FI,trans
tunif,20
kbc,1
*SET,delt,1
timint,on
autots,off
*do,tm,1e-6,30,0.5
/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST
FLST,2,1,5,ORDE,1
*DO,K,1,7,1
*DO,I,1,NUM_TIME,1
FLUX2(I,0,K)=(I-1)
*ENDDO
*DO,J,1,NUM_X,1
FLUX2(0,J,K)=(J*v-v)
*ENDDO
Q=1000
KAPA=1/4
qmax=((Q*KAPA)/3.14159)
WIDTH=W ! width of support bar
WHALF=WIDTH/2
LCYL=2.5 ! length of cylinder
LHALF=LCYL/2
RID=.50 ! inside radius of cylinder
ROD=.75 ! outside radius of cylinder
/DIST,1,.6
/FOCUS,1,.47,.59
WPCSYS,-1,13
CYLIND,,TWELD,-WHALF,WHALF,35,90+(ANG3+ANG4)/2
3d高斯溅射原理
3d高斯溅射原理哎呀，说起3D高斯溅射，这可真是个让人头大的话题。
不过呢，我得承认，这玩意儿确实挺有趣的。
你知道吗，我上次去朋友家玩，他那台3D打印机可真是让我大开眼界。
那天，我一进门，就看到那台3D打印机在角落里默默地工作着。
我心想，这玩意儿能打印出什么花来？结果，我朋友得意洋洋地给我展示了他刚打印出来的一个迷你小房子。
哇塞，那细节，简直了！窗户、门、屋顶，样样俱全，连屋顶上的瓦片都清晰可见。
我当时就惊呆了，这不就是传说中的3D高斯溅射技术吗？我朋友看我一脸懵逼，就给我解释了一下。
他说，3D高斯溅射啊，其实就是一种打印技术，通过控制打印头的移动和材料的喷射，一层一层地堆叠起来，最后就形成了立体的物体。
他一边说，一边还给我演示了一下。
只见那打印头在平台上来回移动，一丝丝的材料就像蜘蛛丝一样，慢慢地堆积起来，逐渐形成了一个立体的形状。
我当时就好奇了，这材料是怎么来的呢？我朋友告诉我，这些材料啊，有的是塑料，有的是金属，甚至还有生物材料。
他说，这些材料都是通过特殊的工艺制成的，可以精确控制材料的喷射量和喷射速度，这样才能保证打印出来的物体精确无误。
我看着那台3D打印机，心里想，这玩意儿真是太神奇了。
我还记得，小时候我们玩的是橡皮泥，那时候要是有这玩意儿，那得多好玩啊！我朋友说，现在3D打印技术已经越来越成熟了，不仅可以打印出各种小玩意儿，还能打印出复杂的机械零件，甚至是人体器官呢！我看着那台3D打印机，心里不禁感慨，科技真是日新月异啊。
想想看，以前我们只能想象的东西，现在都能通过3D打印技术变成现实了。
这3D高斯溅射技术，虽然听起来复杂，但其实就像我们小时候玩泥巴一样，一点一点地塑造出我们想要的形状。
所以啊，下次你再听到3D高斯溅射，可别觉得它遥不可及。
其实，它就在我们身边，说不定哪天，你也能拥有一台属于自己的3D打印机，想打印什么就打印什么。
到时候，你可得请我去看看你的作品哦！。
(整理)ANSYS焊接高斯热源APDL参考.
不准/UNITS,SI/CONFIG,NRES,10000/PREP7ET,1,SOLID70MP,DENS,1,7930MP,C,1,502mptemp,1,20,100,500mpdata,kxx,1,1,12.1,16.3,21.4/VIEW,1,1,1,1!*************定义基板尺寸***********lx=0.1lz=0.1ly=0.006!*************定义焊接参数***********weld=0.01 !每道焊道宽度wheight=0.004 !焊道的高度LSIZE=0.005V=0.008 !速度pi=3.1415926U=25I=180Q=U*I !电源有效功率R=0.006 !电源有效热半径effect=0.8Qmax=effect*Q/(pi*R*R) !中心处最大热流tinc=LSIZE/V !每小段的时间间隔!*****************建立模型，生成网格************** block,0,lx,0,ly,0,lzwpoff,0,ly,lz/2block,0,lx,0,wheight,-weld/2,weld/2wpoff,0,0,-weld/2vsbw,all,,deletewpoff,0,0,weldvsbw,all,,deletevglue,alllsel,s,,,plsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35lsel,s,line,,26lsel,a,line,,30lsel,a,line,,33lsel,a,line,,35 lesize,all,0.005,,,0.2,,,,1lsel,s,line,,9,10lsel,a,line,,12lsel,a,line,,38 lesize,all,0.005,,,0.2,,,,1lsel,s,line,,2lsel,a,line,,4,5lsel,a,line,,7lsel,a,line,,14lsel,a,line,,16lsel,a,line,,19lsel,a,line,,39 lesize,all,0.002,,,1,,,,1lsel,s,line,,23,24lsel,a,line,,41,42 lsel,a,line,,45lsel,a,line,,47lsel,s,line,,23,24lsel,a,line,,41,42 lsel,a,line,,45lsel,a,line,,47 lesize,all,0.002,,,1,,,,1vmesh,all/PNUM,DEFA EPLOTfini/soluantype,trans,new!nlgeom,ontimint,0,struct timint,1,therm timint,0,magtimint,0,electtref,25nropt,autoautos,onkbc,0pred,onlnsrch,on!************杀死焊缝去单元***********nsel,s,loc,y,ly,ly+wheightesln,s,1ekill,allallsel,allesel,s,liveeplot!**********施加对流载荷*************esel,s,livensel,s,loc,z,0nsel,a,loc,z,lznsel,a,loc,x,0nsel,a,loc,x,lxsf,all,conv,10,25!************定义数组维数**********MAX_X=1+lx/LSIZE!*************定义table数组****************************************** *do,i,1,MAX_X,1local,12,0,(i-1)*lsize,ly,lz/2,,,,,,*del,_FNCNAME*del,_FNCMTID*del,_FNC_C1*del,_FNC_C2*del,_FNC_C3*del,_FNC_C4*del,_FNCCSYS*del,'heatflux'!定义表格各行*set,_FNCNAME,'heatflux'*dim,_FNC_C1,,1*dim,_FNC_C2,,1*dim,_FNC_C3,,1*dim,_FNC_C4,,1!表格各行赋值*set,_FNC_C1(1),effect*set,_FNC_C2(1),U*set,_FNC_C3(1),I*set,_FNC_C4(1),R*set,_FNCCSYS,12*DIM,%_FNCNAME%,TABLE,6,26,1,,,,%_FNCCSYS%!! Begin of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) %_FNCNAME%(0,0,1)= 0.0, -999%_FNCNAME%(2,0,1)= 0.0%_FNCNAME%(3,0,1)= %_FNC_C1(1)%%_FNCNAME%(4,0,1)= %_FNC_C2(1)%%_FNCNAME%(5,0,1)= %_FNC_C3(1)%%_FNCNAME%(6,0,1)= %_FNC_C4(1)%%_FNCNAME%(0,1,1)= 1.0, -1, 0, 3, 0, 0, 17%_FNCNAME%(0,2,1)= 0.0, -2, 0, 1, -1, 3, 17%_FNCNAME%(0,3,1)= 0, -1, 0, 1, -2, 3, 18%_FNCNAME%(0,4,1)= 0.0, -2, 0, 1, -1, 3, 19%_FNCNAME%(0,5,1)= 0.0, -1, 0, 0, 0, 0, 0%_FNCNAME%(0,6,1)= 0.0, -3, 0, 1, 0, 0, -1%_FNCNAME%(0,7,1)= 0.0, -4, 0, 1, -1, 2, -3%_FNCNAME%(0,8,1)= 0.0, -1, 0, 3, 0, 0, -4%_FNCNAME%(0,9,1)= 0.0, -3, 0, 1, -4, 3, -1%_FNCNAME%(0,10,1)= 0.0, -1, 0, 2, 0, 0, 2%_FNCNAME%(0,11,1)= 0.0, -4, 0, 1, 2, 17, -1%_FNCNAME%(0,12,1)= 0.0, -1, 0, 2, 0, 0, 4%_FNCNAME%(0,13,1)= 0.0, -5, 0, 1, 4, 17, -1%_FNCNAME%(0,14,1)= 0.0, -1, 0, 1, -4, 1, -5%_FNCNAME%(0,15,1)= 0.0, -4, 0, 1, -3, 3, -1%_FNCNAME%(0,16,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,17,1)= 0.0, -3, 0, 1, 20, 17, -1%_FNCNAME%(0,18,1)= 0.0, -1, 0, 1, -4, 4, -3%_FNCNAME%(0,19,1)= 0.0, -1, 7, 1, -1, 0, 0%_FNCNAME%(0,20,1)= 0.0, -3, 0, 1, -2, 3, -1%_FNCNAME%(0,21,1)= 0.0, -1, 0, 2, 0, 0, 20%_FNCNAME%(0,22,1)= 0.0, -2, 0, 1, 20, 17, -1%_FNCNAME%(0,23,1)= 0.0, -1, 0, 3.14159265358979310, 0, 0, -2%_FNCNAME%(0,24,1)= 0.0, -4, 0, 1, -1, 3, -2%_FNCNAME%(0,25,1)= 0.0, -1, 0, 1, -3, 4, -4%_FNCNAME%(0,26,1)= 0.0, 99, 0, 1, -1, 0, 0! End of equation: 3*effect*U*I*exp(-3*({X}^2+{Z}^2)/R^2)/({PI}*R^2) !-->!**********激活单元*********esel,s,livensel,s,loc,x,(i-1)*lsize-0.002,(i-1)*lsize+0.002nsel,r,loc,z,-weld/2,weld/2nsel,r,loc,y,0,wheightesln,s,0ealive,allallsel,all!**********施加热流载荷********** esel,s,livensel,s,loc,Y,wheightsf,all,hflux,%heatflux%allsel,alltime,i*tincnsubst,2SOLVE!*************删除热流载荷*********** nsel,s,loc,y,wheightsfdele,all,hfluxallsel,allesel,s,liveeplotOUTRES,ALL,ALL,/PSF,HFLUX,,2/REPLOTOUTRES,ALL,ALL,*ENDDO。
热压条件下激光深熔焊接温度场的数值模拟
热压条件下激光深熔焊接温度场的数值模拟
熊智军;李永强;赵熹华;李民
【期刊名称】《焊接学报》
【年(卷),期】2007(028)008
【摘要】对铝合金LF3Y2搭接板进行了热压条件下激光深熔焊接.建立了该条件下的焊接热源模型,热源模型由作用在工件表面的高斯面热源和沿激光入射方向的旋转高斯体热源构成.使用该热源模型和ANSYS有限元分析软件对前述的试验进行了数值模拟.结果表明,热压使激光在高速焊接时依然有足够的熔深和熔宽;试件表面的等温线呈椭圆形,在移动热源前方温度梯度大,后方温度梯度小;计算所得的熔池截面成形与工艺试验结果吻合良好,验证了该热源模型的合理性.
【总页数】4页(P41-44)
【作者】熊智军;李永强;赵熹华;李民
【作者单位】吉林大学,材料科学与工程学院,长春,130022;吉林大学,材料科学与工程学院,长春,130022;吉林大学,材料科学与工程学院,长春,130022;吉林大学,材料科学与工程学院,长春,130022
【正文语种】中文
【中图分类】TG456.7
【相关文献】
1.深熔激光焊接熔池温度场的数值模拟 [J], 张瑞华;樊丁;片山聖二
2.基于铝合金激光深熔焊接温度场研究 [J], 钟荣华;金湘中
3.激光深熔温度场数值模拟热源模型分析 [J], 顾兰;薛忠明;张彦华
4.深熔激光焊接熔池温度场的数值模拟 [J], 张瑞华;陈磊;樊丁;片山聖二
5.奥氏体不锈钢激光焊接温度场数值模拟 [J], 郑旺旺
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
/clear
jbc=2.5e-3/2 Leabharlann jbk=1.6e-3/2
jbh=0.4e-3
fmc=0.9e-3
fmk=0.6e-3
fmh=0.12e-3
v=0.12
p=195
r=0.05e-3
sh=0.04e-3
sc=0.025e-3
sk=0.025e-3
!************************参数计算莵E
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 21, 17, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -1, 4, -3
TUNIF,298,
NSUBST,3,0,0
ESEL,S,MAT,,2
EKILL,all
alls
pred,off
alls
solve
QMX=2*0.72*P/3.1415926/R/R*1800
ZDZB=R
/SOLU
*DEL,_FNCNAME
*DEL,_FNCMTID
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,sh, , , , , , ,0
FLST,5,4,4,ORDE,4
FITEM,5,2
FITEM,5,4
FITEM,5,-5
FITEM,5,7
! 2/r^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
VSBW, all
wpro,,30.000000,
wpro,,30.000000,
wpro,,30.000000,
VSBW, all
KWPAVE, 3
VSBW, all
!************************材料属性的设置
TYPE,2
ESIZE,JBC/20,0,
VSEL,U,MAT,,2
VMESH,ALL
!************************************边界条件加载
/SOL
ANTYPE,4
time,0.00001
TRNOPT,FULL
ASEL,S,EXT
SFA,all,1,CONV,5,298
DT=SC/V
NNN=FMC/SC
ppp=fmk/r/2
/prep7
!************************单元选詠E
ET,1,SOLID70
ET,2,SOLID70
!************************模型建立
BLOCK,,fmc,,fmk,,fmh,
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 1, 4, 2, 20
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 2, 0, 0, -3
*SET,_FNCCSYS,11
! /INPUT,ZH3.func,,,1
*DIM,%_FNCNAME%,TABLE,8,25,1,,,,%_FNCCSYS%
!
! Begin of equation: qmx*exp(-2*((ZDZB-{y})^2+({TIME}*v-(fmc-{X}))^2)+({Z}-sh)^
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,sc, , , , , , ,0
!*
FLST,5,4,4,ORDE,4
FITEM,5,1
FITEM,5,3
FITEM,5,6
FITEM,5,8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,sk, , , , , , ,0
mat,2
TYPE,1
MSHAPE,0,3D
MSHKEY,1
VMESH,1
MAT,1
BLOCK,-(jbc-fmc)/2,fmc+(jbc-fmc)/2,-(jbk-fmk)/2,fmk+(jbk-fmk)/2,,-jbh,
VGLUE,all
wpro,,,30.000000
wpro,,,30.000000
wpro,,,30.000000
VSBW,all
KWPAVE, 4
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, 18, 2, 3
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 17, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, 1, 3, 19
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNC_C6
*DEL,_FNCCSYS
*SET,_FNCNAME,'ZHF'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
! End of equation: qmx*exp(-2*((ZDZB-{y})^2+({TIME}*v-{X})^2+({Z}-sh)^2)/r^2)
!-->
*DEL,ZHF
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -2, 0, 0
MPDATA,C,2,,412
MPDATA,C,2,,460
MPDATA,C,2,,544
MPDATA,C,2,,583
!*******************************划分网竵E
FLST,5,4,4,ORDE,2
FITEM,5,9
FITEM,5,-12
MPTEMP,1,0
MPDATA,DENS,2,,8240
MPTEMP,,,,,,,,
MPTEMP,1,273
MPTEMP,2,373
MPTEMP,3,473
MPTEMP,4,673
MPTEMP,5,873
MPTEMP,6,1073
MPTEMP,7,1573
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
! /INPUT,ZH2.func,,,1
*DIM,%_FNCNAME%,TABLE,7,24,1,,,,%_FNCCSYS%
!
! Begin of equation: qmx*exp(-2*((ZDZB-{y})^2+({TIME}*v-{X})^2+({Z}-sh)^2)/r^2)
MPDATA,KXX,2,,10.31
MPDATA,KXX,2,,11.88
MPDATA,KXX,2,,13.6
MPDATA,KXX,2,,16.6
MPDATA,KXX,2,,20.1
MPDATA,KXX,2,,26.3
MPDATA,KXX,2,,30.75
MPTEMP,,,,,,,,
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),qmx
*SET,_FNC_C2(1),ZDZB