下面为定义的命令流。
如果是时间的函数就好些了,关键是没法转换成时间的函数。
/prep7
et,1,70
!定义导热系数表格及比热容表格随温度变换情况
*dim,K,table,75,1,1,temp
K(1,1)=1.96,1.94,1.92,1.90,1.88,1.86,1.85,1.83,1.81,1.79,1.77,1.75,1.74,1.72,1.70,1.68,1.67,1.65,1.63,1.62,1.60,1.58,1.57,1.55,1.54,1.52,1.51,1.49,1.48,1.46,1.45,1.43,1.42,1.40,1.39,1.37,1.36,1.35,1.33,1.32,1.31,1.29,1.28,1.27,1.26,1.24,1.23,1.22,1.21,1.20,1.19,1.17,1.16,1.15,1.14,1.13,1.12,1.11,1.10,1.09,1.08,1.07,1.06,1.05,1.04,1.03,1.03,1.02,1.01,0.97,0.93,0.90,0.88,0.85,0.83
K(1,0)=20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,280,290,300,310,320,330,340,350,360,370,380,390,400,410,420,430,440,450,460,470,480,490,500,510,520,530,540,550,560,570,580,590,600,610,620,630,640,650,660,670,680,690,700,750,800,850,900,950,1000
K(0,1)=1.0
!比热容
*dim,C,table,75,1,1,temp
C(1,1)=913.22,919.75,926.22,932.64,939.00,945.31,951.56,957.75,963.89,969.97,976.00,981.97,987.89,993.75,999.56,1005.31,1011.00,1016.64,1022.22,1027.75,1033.22,1038.64,1044.00,1049.31,1054.56,1059.75,1064.89,1069.97,1075.00,1079.97,1084.89,1089.75,1094.56,1099.31,1104.00,1108.64,1113.22,1117.75,1122.22,1126.64,1131.00,1135.31,1139.56,1143.75,1147.89,1151.97,1156.00,1159.97,1163.89,1167.75,1171.56,1175.31,1179.00,1182.64,1186.22,1189.75,1193.22,1196.64,1200.00,1203.31,1206.56,1209.75,1212.89,1215.97,1219.00,1221.97,1224.89,1227.75,1230.56,1243.75,1255.56,1265.97,1275.00,1282.64,1288.89
C(1,0)=20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,280,290,300,310,320,330,340,350,360,370,380,390,400,410,420,430,440,450,460,470,480,490,500,510,520,530,540,550,560,570,580,590,600,610,620,630,640,650,660,670,680,690,700,750,800,850,900,950,1000
C(0,1)=1.0
mp,kxx,1,%K%
mp,C,1,%C%
mp,dens,1,2450 |