求助：双椭球热源表格加载中其后面数据怎么得来的？
浏览：29028 回答：2
这是我从论文《铝合金激光熔敷流场和温度场的数值模拟》中复制的程序，如下： finish /clear /BATCH WPSTYLE,,,,,,,,0 /FILNAME,3800-5-5,0 /TITLE,3800-5-5 /units,si !!! 参数定义 L=0.02 B=0.01 H1=0.002 H2=0.003+H1 Ni=0.07 Yi=0.1 D=0.005 ! 焊缝宽度 ah=D/2 ! 双椭球热源参数定义 chf=2*ah chb=4*ah bh=0.0005 rf=2*chf/(chf+chb) rb=2*chb/(chf+chb) fyl=chf/ah byl=chb/ah fbzw=bh/ah K=0.4 ! 功率有效系数 W0=3500 ! 功率 W=W0*K ! 有效功率 F1=0.6 F2=1-F1 /PREP7 et,1,fluid142 keyopt,1,1,3 ! 设定组份数 MSDATA,1,8314.29 mat,1 MSPROP,1,DENS,LIQUID,1024,933,0.3,0, ! 各组份的热物性 MSVARY,1,DENS,t MSPROP,1,COND,constant,220,0,0,0, MSVARY,1,COND,0 MSPROP,1,SPHT,CONSTANT,3000, MSVARY,1,COND,0 MSPROP,2,DENS,CONSTANT,8900 MSVARY,2,DENS,0 MSPROP,2,COND,CONSTANT,60 MSVARY,2,COND,0 MSVARY,2,SPHT,CONSTANT,427 MSVARY,2,SPHT,0 MSPROP,3,DENS,CONSTANT,4470 MSVARY,3,DENS,0 MSPROP,3,COND,CONSTANT,17.2 MSVARY,3,COND,0 MSPROP,3,SPHT,CONSTANT,300 MSVARY,3,SPHT,0 MPTEMP,1,273,933,934,950,975,1000 MPTEMP,7,1050,1075,1100,1150,1175,1200 MPDATA,VISC,1,1,1,1,1.35E-006,1.298E-006,1.235E-006,1.178E-006 MPDATA,VISC,1,7,1.079E-006,1.035E-006,9.96E-007,9.25E-007,8.94E-007,8.65E-007 MPTEMP,,,,, MPTEMP,1,273,933,934,2000 MPDATA,HF,1,1,20,20,40,40 type,1 ! 几何模型建立和网格划分 BLOCK,0,B,0,L,0,-H1 BLOCK,0,B,0,L,-H1,-H2 VGLUE,ALL lesize,all,0.0003 VMESH,ALL NUMCMP,ALL NUMMRG,ALL,,,,LOW ALLSEL VSEL,S,,,1 NSLV,S,1 IC,ALL,SP02,Ni IC,ALL,SP03,Yi IC,ALL,SP01,1-Ni-Yi ALLSEL VSEL,S,,,2 NSLV,S,1 IC,ALL,SP01,1 IC,ALL,SP02,0 IC,ALL,SP03,0 asel,all SAVE FINISH /SOLU ASEL,ALL ! 施加载荷和边界条件 ASEL,U,,,5 ASEL,U,,,8 NSLA,S,1 SF,ALL,CONV,40,293 ALLSEL,ALL ASEL,ALL ASEL,U,,,1 ASEL,U,,,5 ASEL,U,,,8 ASEL,U,,,2 NSLA,S,1 D,ALL,VX,0 D,ALL,VY,0 D,ALL,VZ,0 ASEL,S,,,2 NSLA,S,1 D,ALL,VZ,0 ASEL,S,,,5 ASEL,A,,,8 NSLA,S,1 D,ALL,VX,0 ALLSEL,ALL FINISH /SOLU ! 流体性质设置 FLDATA12,PROP,VISC,16 FLDATA13,VARY,VISC,1 FLDATA12,PROP,COND,CMIX FLDATA13,VARY,COND,1 FLDATA12,PROP,SPHT,CMIX FLDATA13,VARY,SPHT,1 FLDATA12,PROP,DENS,CMIX FLDATA13,VARY,DENS,1 FLDATA7,PROT,VISC,TABLE FLDATA8,NOMI,VISC,-1 FLDATA9,COF1,VISC,0 FLDATA10,COF2,VISC,0 FLDATA11,COF3,VISC,0 FLDATA12,PROP,IVIS FLDATA7,PROT,COND,CMIX FLDATA7,PROT,SPHT,CMIX FLDATA7,PROT,DENS,CMIX FLDATA7,PROT,SFTS,LIQUID FLDATA8,NOMI,SFTS,2400, FLDATA9,COF1,SFTS,933 FLDATA10,COF2,SFTS,0.274 FLDATA11,COF3,SFTS,0 FLDATA26,STAB,MOME,1, ! 收敛控制 FLDATA26,STAB,PRES,1, FLDATA26,STAB,TEMP,1, FLDATA26,STAB,TURB,1, FLDATA26,STAB,VISC,0, FLDATA25,RELX,VX,0.5, FLDATA25,RELX,VY,0.4, FLDATA25,RELX,VZ,0.5, FLDATA25,RELX,PRES,0.5, FLDATA25,RELX,TEMP,0.8, FLDATA25,RELX,ENKE,0.4, FLDATA25,RELX,ENDS,0.5, FLDATA1,SOLU,TRAN,1 FLDATA1,SOLU,FLOW,1 FLDATA1,SOLU,TEMP,1 FLDATA1,SOLU,TURB,0 FLDATA1,SOLU,COMP,0 FLDATA1,SOLU,VOF,0 FLDATA1,SOLU,SFTS,0 FLDATA1,SOLU,IVSH,0 FLDATA1,SOLU,SWRL,0 FLDATA1,SOLU,SPEC,0 FLDATA1,SOLU,ALE,0 FLDATA1,SOLU,RDSF,0 FLDATA1,SOLU,VOF,0 FLDATA15,PRES,REFE,101350, ! 求解参考条件 FLDATA16,BULK,BETA,1e+015, FLDATA17,GAMM,COMP,1.4, FLDATA14,TEMP,NOMI,293, FLDATA14,TEMP,TTOT,293, FLDATA14,TEMP,BULK,293, TOFFSET,0, ACEL,0,0,-9.81, STEP=0.001 ! 求解控制 N=L/STEP V=0.005 ! 速度 *DO,I,0,N YY=STEP*I TINC=STEP/V T=TINC*(I+1) FLDATA4,TIME,STEP,T FLDATA4,TIME,NUMB,3, FLDATA4,TIME,GLOB,10, FLDATA4,TIME,VX,1e-6, FLDATA4,TIME,VY,1e-6, FLDATA4,TIME,VZ,1e-6, FLDATA4,TIME,PRES,1e-006, FLDATA4,TIME,TEMP,1e-006, FLDATA4,TIME,ENKE,0.0001, FLDATA4,TIME,ENDS,0.0001, FLDATA4A,STEP,OVER,1, FLDATA4A,STEP,APPE,1, FLDATA4A,STEP,SUMF,10, FLDATA4,TIME,BC,0 FLDATA4,TIME,TEND,T FLDATA5,OUTP,PTOT,1 FLDATA5,OUTP,TTOT,1 FLDATA5,OUTP,HFLU,1 FLDATA5,OUTP,HFLM,1 FLDATA5,OUTP,STRM,1 FLDATA5,OUTP,PCOE,0 FLDATA5,OUTP,MACH,0 FLDATA5,OUTP,YPLU,0 FLDATA5,OUTP,TAUW,0 FLDATA5,OUTP,RDFL,0 /SOLU NROPT,FULL,,ON lnsrch,1 !!! TIMINT,1,ALL !!! 添加第一个载荷步（焊接阶段，有对流和热源） !!! 双椭球热源的加载 !!! 前半球节点的选取 LOCAL,11,2,0,YY,0.000, , , ,fyl,fbzw, NSEL,S,LOC,X,0,ah NSEL,R,LOC,Y,0,90 NSEL,R,LOC,Z,-90,0 CM,FRONT,NODE LOCAL,13,0,0,YY,0.000, , , ,1,1, !!! 前半球函数、表格定义 *DEL,TFRONT *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNC_C6 *DEL,_FNCCSYS *SET,_FNCNAME,'TFRONT' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *DIM,_FNC_C6,,1 *SET,_FNC_C1(1),AH *SET,_FNC_C2(1),CHF *SET,_FNC_C3(1),BH *SET,_FNC_C4(1),RF *SET,_FNC_C5(1),F1 *SET,_FNC_C6(1),W *SET,_FNCCSYS,13 *DIM,%_FNCNAME%,TABLE,8,41,1,,,,%_FNCCSYS% *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)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 4, 17 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, 3, 4, 18 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,9,1), 0.0, -5, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, 4, 4, 19 *SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 2, 0, 0, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -2, 17, -4 *SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -1, 1, -5 *SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 3, 0, 0, -2 *SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 1, -3, 3, -4 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 6, 0, 0, -1 *SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -1, 0, 1, 0, 0, 0 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 1, -1, 4, -2 *SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 3, 0, 0, -4 *SET,%_FNCNAME%(0,25,1), 0.0, -2, 0, 1, -1, 17, -4 *SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,27,1), 0.0, -2, 0, 1, -1, 3, 20 *SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -2, 3, 21 *SET,%_FNCNAME%(0,29,1), 0.0, -2, 0, 1, -1, 3, 22 *SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, 0 *SET,%_FNCNAME%(0,31,1), 0.0, -3, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,32,1), 0.0, -4, 0, 1, -1, 17, -3 *SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, 0, 0, 0 *SET,%_FNCNAME%(0,34,1), 0.0, -3, 0, 3, 0, 0, -1 *SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -3 *SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 1, -4, 17, -5 *SET,%_FNCNAME%(0,37,1), 0.0, -3, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -3, 3, 19 *SET,%_FNCNAME%(0,39,1), 0.0, -3, 0, 1, -1, 3, 18 *SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 1, -2, 4, -3 *SET,%_FNCNAME%(0,41,1), 0.0, 99, 0, 1, -1, 0, 0 BF,FRONT,HGEN, %TFRONT% CMDELE,FRONT ALLSEL,ALL NSEL,S,LOC,Z,0 NSEL,R,LOC,X,0,0.5*ah NSEL,R,LOC,Y,0,0.5*chf BF,all,FORC,1,1,0 ALLSEL,ALL !!! 后半球节点的选取 LOCAL,12,2,0,YY,0, , , ,byl,fbzw, NSEL,S,LOC,X,0,ah NSEL,R,LOC,Y,-90,0 NSEL,R,LOC,Z,-90,0 CM,BEHIND,NODE !!! 后半球加载 LOCAL,14,0,0,YY,0, , , ,1,1, !!! 后半球函数、表格定义 !!! 后半球函数、表格定义 *DEL,TBEHIND *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNC_C6 *DEL,_FNCCSYS *SET,_FNCNAME,'TBEHIND' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *DIM,_FNC_C6,,1 *SET,_FNC_C1(1),AH *SET,_FNC_C2(1),CHB *SET,_FNC_C3(1),BH *SET,_FNC_C4(1),RB *SET,_FNC_C5(1),F2 *SET,_FNC_C6(1),W *SET,_FNCCSYS,14 *DIM,%_FNCNAME%,TABLE,8,41,1,,,,%_FNCCSYS% *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)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 4, 17 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, 3, 4, 18 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,9,1), 0.0, -5, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, 4, 4, 19 *SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 2, 0, 0, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -2, 17, -4 *SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -1, 1, -5 *SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 3, 0, 0, -2 *SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 1, -3, 3, -4 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 6, 0, 0, -1 *SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -1, 0, 1, 0, 0, 0 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 1, -1, 4, -2 *SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 3, 0, 0, -4 *SET,%_FNCNAME%(0,25,1), 0.0, -2, 0, 1, -1, 17, -4 *SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,27,1), 0.0, -2, 0, 1, -1, 3, 20 *SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -2, 3, 21 *SET,%_FNCNAME%(0,29,1), 0.0, -2, 0, 1, -1, 3, 22 *SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, 0 *SET,%_FNCNAME%(0,31,1), 0.0, -3, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,32,1), 0.0, -4, 0, 1, -1, 17, -3 *SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, 0, 0, 0 *SET,%_FNCNAME%(0,34,1), 0.0, -3, 0, 3, 0, 0, -1 *SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -3 *SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 1, -4, 17, -5 *SET,%_FNCNAME%(0,37,1), 0.0, -3, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -3, 3, 19 *SET,%_FNCNAME%(0,39,1), 0.0, -3, 0, 1, -1, 3, 18 *SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 1, -2, 4, -3 *SET,%_FNCNAME%(0,41,1), 0.0, 99, 0, 1, -1, 0, 0 BF,BEHIND,HGEN, %TBEHIND% CMDELE,BEHIND ALLSEL,ALL NSEL,S,LOC,Z,0 NSEL,R,LOC,X,0,0.5*ah NSEL,R,LOC,Y,0,0.5*chb BF,all,FORC,1,1,0 ALLSEL,ALL SOLVE BFDELE,all,HGEN *ENDDO SAVE 其中热源加载表格*DIM,%_FNCNAME%,TABLE,8,41,1,,,,%_FNCCSYS%之后，它后面的数据怎么得来的，如*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1 求各位大神释疑？非常感谢！！ |
其实就是想问这表格怎么加载。你知道吗，求不吝赐教！