高斯激光热源下的n层材料瞬态热传导方程

xy112

高斯激光热源下的n层材料瞬态热传导方程
' W  o, H) Q5 }2 ~) p高斯, 激光, 热源, 方程, 热传导高斯, 激光, 热源, 方程, 热传导
# ]: k) [5 D6 e2 `' N2 e9 O出于需要， 写了个 n层材料的瞬态热传导方程， 由表面激光加热。每层的热源由其光学性质N= n + ik决定. 通过解Maxwell方程获得激光在材料中的强度分布， 由多项式 a0+a1*y+a2*y^2+a3*y^3+a4*y^4拟和获得。
; K5 \2 E2 C& P# eP=Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*( A0+A1*(-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4).
  
0 E( m; H; V, g+ e& ^' d& o, l高斯分布体热源在程序中由HG.mac加载。不同层的a0,a1,a2,a3,a4决定了其内部的热源分布。
1 L2 B# `! Q% B9 g; K9 t( v( h  
本程序可已根据问题选择所需的材料层数，所有参数都可以改变。  
9 w+ D6 s+ S7 P2 b6 N2 s  
以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！
  
FINISH  ! Make sure we are at BEGIN level   
/ X! s7 A2 X2 d0 D' w/Clear
/FILNAME,Media,0  
  
/PREP7
- q& p  K2 Z( m+ A' x; p) R/PNUM,VOLU,1
  
! In simulation, the X-Z is the media plane, while Y is the depth direction
  
2 Y  _+ s' ]# J6 e8 C% e, y4 wLayer=5                       !Total layer number
nm= 1e-9                    !length unit
( h& ~1 ^9 U3 N! p. X4 @Power=0.005               ! Laser incident power   W
r=5e-7                        ! Laser spot radius      m
$ U8 w! q5 n* e% O7 oV=0                           ! Laser moving veolcity, along z direction,  m/s
  
pw=Power*1000                 !For display Title
rh=r/1e-6                           !For display Title
  
! l' V9 ?! N' }4 A  e" R/TITLE,%Layer% layer media, Laser Power=%pw%mW, Heating radius=%rh%um, Velocity=%V%m/s  
  
  g# P+ X# d3 b3 Y; s1 L+ _*DIM,Thickness,ARRAY,Layer    ! Array input for each layer thickness
*DIM,k_L,ARRAY,Layer          ! Array input for each layer lateral thermal conductivity  
*DIM,k_N,ARRAY,Layer          ! Array input for each layer normal thermal conductivity
: K( N# ?0 g8 q6 E) n) u*DIM,Cp,ARRAY,Layer           ! Array input for each layer heat capacity
6 ~/ p) Z! Y5 u8 t9 p*DIM,Den,ARRAY,Layer          ! Array input for each layer density
*DIM,div,ARRAY,Layer          ! Array input for each layer thickness division
) `1 S* |# j1 N/ {*DIM,dscale,ARRAY,Layer       ! Array input for each layer thickness division scale
+ j1 b3 [: Y: q*DIM,A0,ARRAY,Layer           ! Array input for each layer a0 coefficient for heat generation curve
*DIM,A1,ARRAY,Layer           ! Array input for each layer a1 coefficient for heat generation curve
*DIM,A2,ARRAY,Layer           ! Array input for each layer a2 coefficient for heat generation curve
' `, h$ H" _" z' ~1 p*DIM,A3,ARRAY,Layer           ! Array input for each layer a3 coefficient for heat generation curve
0 T2 t' O4 v$ _" n0 M! I! m*DIM,A4,ARRAY,Layer           ! Array input for each layer a4 coefficient for heat generation curve
  
scale=10000                   !for vovlap operation
* @3 L  N2 A$ v3 a" L4 ^8 cPulse_Duration=1e-8           !Laser pulse duration
3 h+ C, \' X8 jEnd_Time=2e-8                 !Final time for simulation
& f; q, t# ]/ T& O2 VStep=20                       !Time interval  
1 `7 L7 R9 I' m4 u6 r- Q. R0 gT0=50                         !Environment temperature   
  
2 H. x. ^. _9 o  ?9 E: OLength=4000*nm                !Length of the block in z direction
Width=2000*nm                 !Width of the block in x direction
. D# ]' j- I$ i1 U5 z9 P1 Z! o  
!Layer properties, from top to bottom  
  
) I2 T, F% t' u1 J6 U4 J3 y) S!First layer thermal property
- J6 Z/ ?) G1 D4 |# s% q4 Z& ]( CNum=1                        
3 y' `2 G0 e: u( r8 J- Q2 Hk_L(Num)=1.4                  
0 {+ F% K- V; b- h" ak_N(Num)=1.4
Den(Num)=4000
Cp(Num)=500
5 @* n* ]% Z( ZThickness(Num)=50*nm
  v& S' x# t$ M  y/ Z9 C7 m# rdiv(Num)=5
) F" n" |: [: \+ c% V  Ddscale(Num)=1                 !last/first in -Y direction
A0(Num)=0
A1(Num)=0
+ Z+ t6 C( I. ~- p: {- }* o( vA2(Num)=0
) j/ ~2 B) N! J  y$ p6 w& uA3(Num)=0
6 i" {! K& b6 q' JA4(Num)=0
* ^& u; q4 i, B9 v4 S6 I  
!Second layer thermal property
Num=2
- H$ q- d& a) X- H1 ^k_L(Num)=16
# L9 w+ W' u" @: _k_N(Num)=16
' L5 Y% [0 a* e0 K( @& HDen(Num)=3000
Cp(Num)=900
$ h. |' J1 w8 }Thickness(Num)=20*nm
div(Num)=2
' L7 r9 r1 k" s2 q: Y$ s! s3 H. Gdscale(Num)=1
A0(Num)=807.929E+6
( w, r/ ?9 F& y# C4 `7 e4 M9 xA1(Num)=-43.896E+6
A2(Num)=934.529E+3
! ]( ^9 d: x  e* x5 M1 ~0 pA3(Num)=-8.864E+3
A4(Num)=31.703
  
!Third layer thermal property
Num=3
k_L(Num)=1.5
5 H, Q; O- E9 Y* ^! ik_N(Num)=1.5
( r8 B& \$ }6 R: n) p! X8 n) [; RDen(Num)=4000
Cp(Num)=500
+ W1 T2 E) H) n. Z1 q( uThickness(Num)=100*nm
3 ]" Q, F% _" R5 ]0 U6 \: z5 e% Mdiv(Num)=5
9 l4 u: ~, W# M+ q7 v0 d# v. idscale(Num)=1
A0(Num)=0
A1(Num)=0
+ [3 p( e6 j$ \' O- e! h5 U7 u) _A2(Num)=0
. C; K+ M6 u9 L! bA3(Num)=0
  S, x2 G% Y0 y$ ]A4(Num)=0
+ M" f& {# M* l4 }) M( |  
( u# D# f$ G* X!Fourth layer thermal property
Num=4
k_L(Num)=100
1 j$ @& y: C# ik_N(Num)=100
4 p( j6 h7 V) w/ G! H: X$ E2 mDen(Num)=3900
Cp(Num)=700
+ h* ]5 X  h! k: HThickness(Num)=50*nm
2 |/ }+ ^2 [2 M1 F" [$ O- h' tdiv(Num)=2
6 T/ b9 K. v) C6 V5 U! mdscale(Num)=2
% u* e+ X0 e  F3 t9 b- t" n; MA0(Num)=1.627E+9
, \3 R) w6 B3 f3 d  RA1(Num)=-32.146E+6
A2(Num)=238.195E+3
A3(Num)=-784.173
A4(Num)=967.716E-3
- o. y* y# x0 o4 p5 i  k  
!Fifth layer thermal property
Num=5
- g7 x! v7 b! g9 \k_L(Num)=0.2
; ?; q9 h; n$ o" a, ^k_N(Num)=0.2
Den(Num)=3900
" `0 O+ l7 A  O0 Z2 L# U5 R' zCp(Num)=450
" x$ X, x- s' JThickness(Num)=100*nm
3 G) c: U" e9 J# I! w; ]1 s& a6 Tdiv(Num)=5
dscale(Num)=2
+ t7 O1 V4 I1 h) ^# EA0(Num)=0
A1(Num)=0
3 W6 g5 B/ V( N2 b- D0 _A2(Num)=0
A3(Num)=0
A4(Num)=0
  
! Define the thermal element
# O6 h' C5 [: A7 u% S5 B6 VET,1,SOLID70  
KEYOPT,1,2,0  
KEYOPT,1,4,0  
) B5 M9 g1 B: p) @* HKEYOPT,1,7,0  
) a' w* ^1 e8 s6 R3 b& hKEYOPT,1,8,0  
  
5 F0 O( K. N) m; I- v!For post-precessing
nPoint=200
( o) d0 M% i2 G0 adisplay_r=1e-6
, h/ z" l# S1 K0 D  
4 j) r# ]3 n6 |6 H7 k# MTSum1=0
4 i, X: W0 f; W3 i! m0 [TSum2=0
9 m- k4 z  g* t: i! j  \*do,Num,1,Layer,1
  UIMP,Num,DENS, , ,Den(Num),                                        !Define density of material Num
) k: V' v0 C1 R) f- C  UIMP,Num,KXX,KYY,KZZ,k_L(Num),k_N(Num),k_L(Num),                   !Define thermal conductivity of material Num
  UIMP,Num,C, , ,Cp(Num),                                            !Define heat capacity of material Num   
  TSum1=TSum1+Thickness(Num)  
' b3 P: L8 b/ [+ P+ K4 P2 A  BLOCK,-Width/2,Width/2,-TSum1,-TSum2,-Length/2,Length/2,           !Draw the layer of Num  
. q# Y) i8 r" u4 T. F+ v  TSum2=TSum2+Thickness(Num)
*Enddo
  
d_total=TSum1                 !Total thickness of the multilayer                                                  
  
VLSCALE,ALL, , , scale, scale, scale,  ,1,1,
1 Y8 s: F; S6 k& n/ E6 V  
0 T0 }8 S) D2 M- \" v2 K*IF,Layer,GT,1,THEN
6 I5 v( }7 g# i" Z0 e# U   VGLUE,ALL
*ENDIF
% E$ Q" p9 s# c/ M" Y  
VLSCALE,ALL, , , 1/scale, 1/scale, 1/scale, ,1,1,
/ u+ x/ _7 W7 S6 ~* z  
TSum1=0
6 w9 f; X: S1 l; [$ R- qTSum2=0
$ P. f6 o: a$ [$ j  
Size=100*nm                             !Area mesh size, in X-Z plane
aesize,all,Size,  
8 g! ]7 O' X' _/ O' `: q6 A  
# b8 t: ]3 C: C" O; p% d3 p4 O8 d*do,Num,1,Layer,1
   TSum1=TSum1+THICKNESS(Num)
   lsel,all
  
& j1 P6 @9 e- x5 D8 t; O   !Determine the mesh scale and divison for each layer at Y direction
   *do,i,-1,1,2
" t, |5 `$ R( w+ C9 r1 L8 t      *do, j,-1,1,2
: l% K# H" V2 V2 L% X( J         lsel,s,loc,y,-TSum1,-TSum2
  O% [" v* y1 E         lsel,r,loc,x,i*Width/2,i*Width/2
         lsel,r,loc,z,j*Length/2,j*Length/2
         *GET,Line_Num,LINE, 0, NUM,MIN,
$ W# w# E& ?5 ?         *GET,KP1,LINE,Line_Num,KP,1,
         *GET,KP2,LINE,Line_Num,KP,2,
         *GET,KP1_Y,KP,KP1, LOC, Y,  
# s1 B: ]. \+ ]" Z         *GET,KP2_Y,KP,KP2, LOC, Y,
         *IF,KP1_Y,GT,KP2_Y,THEN
3 A. A. ]6 `9 d             lesize,all,THICKNESS(Num)/div(Num), , ,dscale(Num),1, , ,1,
         *ELSE
             lesize,all,THICKNESS(Num)/div(Num), , ,1/dscale(Num),1, , ,1,
8 |0 x) M5 \9 T1 t( Z3 i8 K1 t4 E         *ENDIF
      *ENDDO
9 _. p" f$ Y9 P   *ENDDO   
  
* z; L9 x. @% w; F2 x   asel,all
- N( |) b' A3 V4 Q+ Z  
   esize,Size,0,
  
   VSEL,S,LOC,Y,-TSum1,-TSum2
  
1 n; `* |% A1 H" ]# l) \. ^   TYPE,1
   MAT,Num
   mshkey,1  
7 W3 H, I2 r* n8 P* x   mshape,0,3d  
( W3 [2 o) J4 F1 o5 D( n   vmesh,ALL             !Mesh the layer Num
  
6 Z4 L: R& f# P/ [+ n8 L( C4 y   allsel,all
  
   TSum2=TSum2+Thickness(Num)
' r3 ^- ~9 Q9 v& z( ~  
1 w" ?1 Q! K( P$ l*enddo
  
/REP,FAST  
/psf,hflu,,2,0  
/PNUM,LINE,1  
$ }% [5 E- |/ [) v5 l% Jfinish  
/ c, n' u' o: o$ Z+ H3 a  X  q- v) b/solu  
, b" s6 \& w) @. v- O  
* u$ M# }1 X. A' Xantype,trans             !Transient simulation
3 y) x( S3 S9 j1 {) r8 M! m# \$ Ftunif,T0                 !Ambient temperature
- d3 Q) [9 v. `; p: d6 I% V3 ETREF,T0
! B$ O3 e" }( w; j2 [timint,on !time integration = on  
1 x) B) N' }' u2 Hkbc,1  
  e% T: N* P& V- F3 @OUTRES,ALL,ALL,          !Output all the result
  
/view, 1 ,1,1,1  
/ang, 1  
/rep,fast  
8 o3 k/ [" y4 U9 D  
*do,Num,1,Layer,1
   *USE,HG.mac,'Layer%Num%',Power,r,V,A0(Num),A1(Num),A2(Num),A3(Num),A4(Num)     !Load heat generation for each layer
. N% f' n- L# z* P) K7 ~; @: G*enddo
  
5 ^0 |( f* g1 ]' r# G; [0 s3 A!Start simulation
$ i. |6 g* P% W* v*do,tm,1e-15,End_Time,End_Time/Step  
# J* o0 t: L  z$ P7 _$ h: K" r  S  
9 q  L3 z! i0 i* @   time,tm  
+ A% N4 B+ e: W; @) }2 X- |3 C   autots,0  
* L/ I. z9 ?: x1 M   deltim, , , ,1  
   kbc,0  
( W( f9 ~) I. b! T' ^1 K  
; C) ^' R% f/ t( p5 J5 n   allsel,all  
3 K+ B( A% a1 z% S& B- Z( [7 w' J! L   TSum1=0
5 E! [2 k7 Q- m" n$ Q  L   TSum2=0
   Num=0
7 G# ?1 M9 S( p6 f$ Q) s  G  
# W4 [0 g$ k# D/ v8 I% E   *if,tm,LE,Pulse_Duration+1e-12,then  
        
      !Apply the heat generation in each layer  
6 t% k6 C/ `+ E8 K6 H      Num=1
* F5 P% b6 z) S      TSum1=TSum1+THICKNESS(Num)
, h2 U8 T5 B( n1 W# M      nsel,s,LOC,y,-TSum1,-TSum2   
) V. \& Q  z6 m9 Q) Y( y( B" p      BF,all,HGEN,%Layer1%
% m! F: V0 c. {0 p/ m      TSum2=TSum2+Thickness(Num)
      
8 u. t0 e' g  z( G( p. r- `      Num=2
- j+ h# \/ e' T8 A4 B- o% }      TSum1=TSum1+THICKNESS(Num)
! @) c  j. ?( s4 B+ L  r# h* X      nsel,s,LOC,y,-TSum1,-TSum2   
  U  o+ F# E  \0 ?      BF,all,HGEN,%Layer2%
      TSum2=TSum2+Thickness(Num)
$ e9 W4 U, ?  j' P1 E# L2 w   
      Num=3   
      TSum1=TSum1+THICKNESS(Num)
* I+ W  g/ u) z; q7 [* k% n      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer3%
      TSum2=TSum2+Thickness(Num)
   
      Num=4
      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer4%
/ b3 W  r8 G* V0 R      TSum2=TSum2+Thickness(Num)
   
      Num=5
      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
6 m/ N( X2 K0 R  k2 w2 r      BF,all,HGEN,%Layer5%
      TSum2=TSum2+Thickness(Num)
; F" j9 |. o) L% Z4 K( B1 W  
4 [! D7 v& [& N* u9 X. |9 [   *else
6 Q3 H/ e3 |9 Y5 S1 q  [      allsel,all  
      BFDELE, ALL,HGEN
   *endif  
  
2 q, ]+ P3 b; n3 M   allsel,all
   asel,s,LOC,y,-1e-9,0
!   DA,ALL,TEMP,T0  
# }# M( r8 }0 F5 O5 m5 S/ _; a   allsel,all
  
0 R) w- q- @, M  o/ B2 k   /status,solu  
* f) F" \$ }) d   solve  
$ N6 ?/ w2 j0 P* m9 a# j  Q*enddo  
+ J4 e3 ?6 \6 r9 Q% p  
2 P3 L- x. ~! [% h' Lfinish  
/POST1
2 x- u- O" s6 X) T9 @" |* f5 t- J; O  
!*ULIB, Post,mac
!*USE,TTV    !Temperature decay at each layer thickness
# W: w( ?2 D4 V, e2 E!*USE,TRV,1  !Normal Temperature profile at each layer at given time  
!*USE,TRL,1  !Plane Temperature profile at each layer at given time  
% \( j/ Q9 }: c" M8 o) }  
!Laser power absorption in each layer
  P  K4 @4 s6 R+ q+ G!Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*( A0+A1*(-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)      
' i1 n5 C$ z; X4 s; \4 G* }' f  
- J4 }; x, Z6 z//UG.mac
- Z0 _% ~  }# ^( w9 i2 x  
. u1 Q. Z9 f( \* _% }+ u1 |!FileName=ARG1
. @8 F. T3 y- [( B& `3 n!Power=ARG2
4 O$ e% \+ |" q9 d3 f!r=ARG3
5 r5 F/ Z, Z0 J!V=ARG4
2 G; U2 F& c+ e. F7 |; A3 S$ C/ T!A0=ARG5
!A1=ARG6
!A2=ARG7
' X# n. }  ]: q, B3 n!A3=ARG8
& F( q. L" k. h7 W- X!A4=ARG9
  o7 W7 O5 T2 n7 j! J: w1 ^  
2 \9 B( n$ A: f, w. S/ @! g*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
4 @8 i, H* h# a& _. S* x3 E*DEL,_FNC_C2
*DEL,_FNC_C3
. C' S7 ?$ P$ J4 \- i$ ^9 L*DEL,_FNC_C4
*DEL,_FNC_C5
0 ]4 K3 A( m6 Z! u*DEL,_FNC_C6
*DEL,_FNC_C7
*DEL,_FNC_C8
& Z  M- m4 V' D* m*SET,_FNCNAME,ARG1   
*DIM,_FNC_C1,,1  
*DIM,_FNC_C2,,1  
) J- r# T: Z8 y) h! h*DIM,_FNC_C3,,1  
*DIM,_FNC_C4,,1  
! T& B% T- c' l( Q) S2 [) m*DIM,_FNC_C5,,1  
*DIM,_FNC_C6,,1  
, o& z  x; N: A: [- C" {2 f*DIM,_FNC_C7,,1  
8 x+ U# d( P4 D+ P* B3 F7 t+ N$ |! [*DIM,_FNC_C8,,1  
7 j4 [- [' h8 L# P5 h/ b) D*SET,_FNC_C1(1),ARG2   
*SET,_FNC_C2(1),ARG3
*SET,_FNC_C3(1),ARG4
*SET,_FNC_C4(1),ARG5
, K& |8 V8 b% `*SET,_FNC_C5(1),ARG6
*SET,_FNC_C6(1),ARG7
*SET,_FNC_C7(1),ARG8
*SET,_FNC_C8(1),ARG9
! /INPUT,HG.func
4 R2 p/ E" e" c% E: P& @4 M. C*DIM,%_FNCNAME%,TABLE,10,61,1   
7 m8 c0 k& A7 R1 t  {; K1 f1 {3 i!   
! Begin of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*   
! (-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)   
$ M/ ^1 k- L3 K% }' `) Q*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0  
4 r  p' u% X* o; W*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
2 F0 C! c- D4 P) V; H8 C*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
& i( h" V7 K7 }  d1 d; _+ ?*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
  A9 M8 K+ f. ?! i% B, O0 N$ K*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
4 J4 u' D% X; N9 t8 O( S*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 18  
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 18, 17, -1   
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, 17, 4, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1   
6 s, s- C" ]# E8 g) [8 `9 L*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 4, -2
, M" S8 Q# ~+ u2 [*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0   
*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, 0, 0, -1  
( r. t! t+ [8 `! {*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -2
4 ?- p2 G2 R) \$ @) R. W7 Y3 ^*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2   
+ Z7 ]. P/ l4 ^" C4 U: k+ w*SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, 2, 17, -1   
9 U  }4 w8 s+ A1 A8 s*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, 19, 3, 1
5 [9 [8 Q' t( P4 P*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 4, 2, -1
7 O  k7 d- _- ?) R, Y: ~*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 2, 0, 0, -5
*SET,%_FNCNAME%(0,14,1), 0.0, -6, 0, 1, -5, 17, -1   
2 G# U  f" n9 a7 d3 L*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -2, 1, -6   
*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 1, -4, 3, -1   
! Q/ V2 A3 O. o*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 2, 0, 0, 18
2 B' I0 Y& Z$ m/ c*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, 18, 17, -1   
9 N; x; _5 u1 u1 ^3 G' ~/ O*SET,%_FNCNAME%(0,19,1), 0.0, -1, 0, 1, -2, 4, -4   
% E* w; q/ t  T, m! s# @: M*SET,%_FNCNAME%(0,20,1), 0.0, -1, 7, 1, -1, 0, 0
. ?# v$ g3 u" K9 D# H*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -3, 3, -1   
& ~  v' n8 S6 B; e*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 0, 0, 0, 0  
- p5 O8 P! O& l, Y  l*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, 0, 0, -1
4 c+ u' ~% l! R6 U( Z* j& j! [! z*SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, -1, 2, -3   
. H1 x- i: |! Z: x- z*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 1, -4, 3, 3
*SET,%_FNCNAME%(0,26,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
*SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, -1, 4, -3   
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, 21, 3, -4   
2 L* X) a& ], a*SET,%_FNCNAME%(0,29,1), 0.0, -3, 0, 1, 20, 1, -1   
5 i, u1 Z7 D1 F# L# Y7 e% ?0 R+ Y*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 0, 0, 0, 0  
3 u* S- X0 Z4 M$ w# S*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, 0, 0, -1
% P7 C, S% S1 y$ B/ `4 V  ~5 ~*SET,%_FNCNAME%(0,32,1), 0.0, -5, 0, 1, -1, 2, -4   
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, -5, 3, 3
- m# E  B- \  V; F*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
! R: }! V* a! M*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -4   
- k" N4 H8 G- L' v3 }$ W0 c+ H*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, -5
+ t& I1 {0 M3 ~: [0 S6 G5 K*SET,%_FNCNAME%(0,37,1), 0.0, -4, 0, 1, -5, 17, -1   
7 V( E7 V3 A) o' v" c8 }& \* y  @*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, 22, 3, -4   
*SET,%_FNCNAME%(0,39,1), 0.0, -4, 0, 1, -3, 1, -1   
" B8 m9 h* o1 S  P: G*SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, 0, 0, -1
; z' S2 [) w) S% s*SET,%_FNCNAME%(0,42,1), 0.0, -5, 0, 1, -1, 2, -3   
8 b$ Z  [' ]+ @0 x( j! F*SET,%_FNCNAME%(0,43,1), 0.0, -1, 0, 1, -5, 3, 3
*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
" l2 T( m: }9 R/ j/ Q0 ?*SET,%_FNCNAME%(0,45,1), 0.0, -5, 0, 1, -1, 4, -3   
2 o5 h5 E" n% Z0 W) W*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 3, 0, 0, -5
*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -5, 17, -1   
*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 1, 23, 3, -3   
*SET,%_FNCNAME%(0,49,1), 0.0, -3, 0, 1, -4, 1, -1   
*SET,%_FNCNAME%(0,50,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,51,1), 0.0, -4, 0, 1, 0, 0, -1
. B5 Y! P# f9 \( F+ \1 P0 g*SET,%_FNCNAME%(0,52,1), 0.0, -5, 0, 1, -1, 2, -4   
/ b/ u! B/ C6 `" p. B' B& g*SET,%_FNCNAME%(0,53,1), 0.0, -1, 0, 1, -5, 3, 3
6 [+ x3 _8 S' D$ N: T- s*SET,%_FNCNAME%(0,54,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
, p$ n  [: l3 X( g*SET,%_FNCNAME%(0,55,1), 0.0, -5, 0, 1, -1, 4, -4   
0 J4 y7 w0 _2 s' M) @*SET,%_FNCNAME%(0,56,1), 0.0, -1, 0, 4, 0, 0, -5
*SET,%_FNCNAME%(0,57,1), 0.0, -4, 0, 1, -5, 17, -1   
" M* t5 J* \( _*SET,%_FNCNAME%(0,58,1), 0.0, -1, 0, 1, 24, 3, -4   
*SET,%_FNCNAME%(0,59,1), 0.0, -4, 0, 1, -3, 1, -1   
& [9 x% d0 ~1 M' A*SET,%_FNCNAME%(0,60,1), 0.0, -1, 0, 1, -2, 3, -4   
*SET,%_FNCNAME%(0,61,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*(-   
! {Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)
( K7 R+ a+ q* R: k!-->