Télécharger test_iwprd3D_sol.dgibi

Retour à la liste

Numérotation des lignes :

  1. *
  2. **********************************************************
  3. * Test pour la loi de comportement IWPR3D_SOL
  4. * On considere un essai numerique cyclique en imposant
  5. * un tensuer de deformation donne:
  6. *
  7. * eps = gam1[(ex^ez)+(ez^ex)]
  8. * avec gam1=1e-4
  9. *
  10. * L'etat de contraite initiale est:
  11. * sig0 = sig0[(ex^ex)+(ey^ey)+(ez^ez)]
  12. * avec sig0 = -3.0e4 Pa
  13. * les parametres de la loi IWPR3D_SOL sont:
  14. * n1 = 0.1 k1 = 0.75
  15. * et les seuil de chaque criteres sont:
  16. * tau1 = 200. m1 = (sqrt(3)*tau1)/sig0
  17. * tau2 = 400. m2 = (sqrt(3)*tau2)/sig0
  18. * tau3 = 600. m3 = (sqrt(3)*tau3)/sig0
  19. * tau4 = 800. m4 = (sqrt(3)*tau4)/sig0
  20. * tau5 = 900. m5 = (sqrt(3)*tau5)/sig0
  21. * tau6 = 1000. m6 = (sqrt(3)*tau6)/sig0
  22. * tau7 = 1200. m7 = (sqrt(3)*tau7)/sig0
  23. * tau8 = 1400. m8 = (sqrt(3)*tau8)/sig0
  24. * tau9 = 1600. m9 = (sqrt(3)*tau9)/sig0
  25. *
  26. * et les caracteristiques associées à chaque critere
  27. * H1 = 4000. D1=0.1 eta1=0.142 eta1=0.142
  28. * H2 = 3500. D2=0.1 eta2=0.142 eta2=0.142
  29. * H3 = 3000. D3=0.1 eta3=0.142 eta3=0.142
  30. * H4 = 2800. D4=0.1 eta4=0.142 eta4=0.142
  31. * H5 = 2600. D5=0.1 eta5=0.142 eta5=0.142
  32. * H6 = 2200. D6=0.1 eta6=0.142 eta6=0.142
  33. * H7 = 1800. D7=0.1 eta7=0.142 eta7=0.142
  34. * H8 = 1000. D8=0.1 eta8=0.142 eta8=0.142
  35. * H9 = 800. D9=0.1 eta9=0.142 eta9=0.142
  36. *
  37. * Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  38. *
  39. **********************************************************
  40. *
  41. OPTI DIME 3 MODE TRID ELEM CUB8;
  42. OPTI TRAC PSC;
  43. OPTI EPSI LINEAIRE;
  44. OPTI EPTR 12;
  45. *
  46. VSSOL1 = 200.;
  47. NU1 = 0.3;
  48. RHO1 = 2000.;
  49. G1 = (VSSOL1*VSSOL1)*(RHO1);
  50. E1 = ((G1)*((2.0)*(1.0 + NU1)));
  51. *
  52. DS1 = 1.0;
  53. *
  54. DT1 = 0.01;
  55. *
  56. D_IMP = 1.E-4;
  57. *
  58. S_XX_INI = -3.E4;
  59. S_YY_INI = -3.E4;
  60. S_ZZ_INI = -3.E4;
  61. *
  62. VAR_GRAP = 0;
  63. *
  64. *********************************************************************
  65. * Maillage
  66. *********************************************************************
  67. *
  68. P0 = 0. 0. 0.;
  69. P1 = P0 PLUS (DS1 0. 0.);
  70. P2 = P0 PLUS (DS1 DS1 0.);
  71. P3 = P0 PLUS (0. DS1 0.);
  72. *
  73. L1 = D 1 P0 P1;
  74. L2 = D 1 P1 P2;
  75. L3 = D 1 P2 P3;
  76. L4 = D 1 P3 P0;
  77. *
  78. SUR1 = DALL L1 L2 L3 L4;
  79. *
  80. VOL1 = SUR1 VOLU 1 TRANS (0. 0. DS1);
  81. *
  82. FLX_1 = L1 TRANS 1 (0. 0. DS1);
  83. FLX_2 = L3 TRANS 1 (0. 0. DS1);
  84. FLY_1 = L4 TRANS 1 (0. 0. DS1);
  85. FLY_2 = L2 TRANS 1 (0. 0. DS1);
  86. FB_IN = SUR1;
  87. FB_SU = SUR1 PLUS (0. 0. DS1);
  88. *
  89. VOL1 = VOL1 COUL VERT;
  90. FLX_1 = FLX_1 COUL ROSE;
  91. FLX_2 = FLX_2 COUL ROSE;
  92. FLY_1 = FLY_1 COUL BLEU;
  93. FLY_2 = FLY_2 COUL BLEU;
  94. FB_IN = FB_IN COUL ROUG;
  95. FB_SU = FB_SU COUL ROUG;
  96. *
  97. ELIM 0.001 (VOL1 ET FLX_1 ET FLX_2 ET FLY_1 ET FLY_2
  98. ET FB_IN ET FB_SU);
  99. *
  100. SI (VAR_GRAP EGA 1);
  101. TRAC (VOL1);
  102. TRAC (VOL1 ET FLX_1 ET FLX_2);
  103. TRAC (VOL1 ET FLX_1 ET FLX_2 ET FLY_1 ET FLY_2);
  104. TRAC (VOL1 ET FLX_1 ET FLX_2 ET FLY_1 ET FLY_2
  105. ET FB_IN ET FB_SU);
  106. FINSI;
  107. *
  108. *********************************************************************
  109. * Modele et Materiau
  110. *********************************************************************
  111. *
  112. *
  113. MOD1 = MODE VOL1 MECANIQUE ELASTIQUE PLASTIQUE IWPR3D_SOL CUB8;
  114. *
  115. XM1 = (((3.0)**(0.5))*(200.0))/(30000.0);
  116. XM2 = (((3.0)**(0.5))*(400.0))/(30000.0);
  117. XM3 = (((3.0)**(0.5))*(600.0))/(30000.0);
  118. XM4 = (((3.0)**(0.5))*(800.0))/(30000.0);
  119. XM5 = (((3.0)**(0.5))*(900.0))/(30000.0);
  120. XM6 = (((3.0)**(0.5))*(1000.0))/(30000.0);
  121. XM7 = (((3.0)**(0.5))*(1200.0))/(30000.0);
  122. XM8 = (((3.0)**(0.5))*(1400.0))/(30000.0);
  123. XM9 = (((3.0)**(0.5))*(1600.0))/(30000.0);
  124. *
  125. MAT1 = MATE MOD1 'YOUNG' E1 'NU' NU1 'RHO' RHO1;
  126. MAT1 = MAT1 ET
  127. (MATE MOD1 'E0' ((E1)/((30000)**(0.1))) 'ALP0' 1.0 'C' 0. 'N1' 0.1
  128. 'PREF' -1.0 'K1' 0.75);
  129. MAT1 = MAT1 ET
  130. (MATE MOD1 'M1' XM1 'M2' XM2 'M3' XM3 'M4' XM4 'M5' XM5
  131. 'M6' XM6 'M7' XM7 'M8' XM8 'M9' XM9 'M10' 19.0);
  132. MAT1 = MAT1 ET
  133. (MATE MOD1 'HC1' 4000. 'HC2' 3500. 'HC3' 3000. 'HC4' 2800.
  134. 'HC5' 2600. 'HC6' 2200. 'HC7' 1800. 'HC8' 1000.
  135. 'HC9' 800.);
  136. MAT1 = MAT1 ET
  137. (MATE MOD1 'HE1' 4000. 'HE2' 3500. 'HE3' 3000. 'HE4' 2800.
  138. 'HE5' 2600. 'HE6' 2200. 'HE7' 1800. 'HE8' 1000.
  139. 'HE9' 800.);
  140. MAT1 = MAT1 ET
  141. (MATE MOD1 'DA1' 0.1 'DA2' 0.1 'DA3' 0.1 'DA4' 0.1 'DA5' 0.1
  142. 'DA6' 0.1 'DA7' 0.1 'DA8' 0.1 'DA9' 0.1);
  143. MAT1 = MAT1 ET
  144. (MATE MOD1 'E_C1' 0.142 'E_C2' 0.142 'E_C3' 0.142 'E_C4' 0.142
  145. 'E_C5' 0.142 'E_C6' 0.142 'E_C7' 0.142 'E_C8' 0.142 'E_C9' 0.142);
  146. MAT1 = MAT1 ET
  147. (MATE MOD1 'E_E1' 0.142 'E_E2' 0.142 'E_E3' 0.142 'E_E4' 0.142
  148. 'E_E5' 0.142 'E_E6' 0.142 'E_E7' 0.142 'E_E8' 0.142 'E_E9' 0.142);
  149. *
  150. *********************************************************************
  151. * Condictions initiales
  152. *********************************************************************
  153. *
  154. SIG_INI = MANU 'CHML' MOD1 'SMXX' S_XX_INI
  155. 'SMYY' S_YY_INI
  156. 'SMZZ' S_ZZ_INI
  157. 'SMXY' 0.0 'SMXZ' 0.0 'SMYZ' 0.0;
  158. F_EX_INI = BSIGM MOD1 SIG_INI;
  159. *
  160. *********************************************************************
  161. * Chargement
  162. *********************************************************************
  163. *
  164. BL1 = BLOQUER 'DEPLA' FB_IN;
  165. BL2 = BLOQUER 'UY' FB_SU;
  166. BL3 = BLOQUER 'UZ' FB_SU;
  167. BL4 = BLOQUER 'UX' FB_SU;
  168. *
  169. DP1 = DEPI BL4 1.0;
  170. *
  171. L_T1 = PROG 0. PAS DT1 13.;
  172. L_E1 = IPOL L_T1 (PROG 0. 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13.)
  173. (PROG 0. 1. 0. -1. 0. 1. 0. -1. 0. 1. 0. -1. 0. 1.);
  174. L_E2 = PROG (DIME L_T1) * 1.0;
  175. *
  176. EV1 = EVOL MANU L_T1 ((L_E1)*(D_IMP));
  177. EV2 = EVOL MANU L_T1 ((L_E2));
  178. *
  179. CHA1 = CHAR 'DIMP' DP1 EV1;
  180. CHA2 = CHAR 'FORCE' F_EX_INI EV2;
  181. *
  182. *********************************************************************
  183. * Calcul PASAPAS
  184. *********************************************************************
  185. *
  186. CH_INI = MANU CHPO VOL1 3 'UX' 0. 'UY' 0. 'UZ' 0.;
  187. *
  188. LLT3 = PROG 0. PAS DT1 13.;
  189. LLT4 = PROG 0. PAS ((DT1)/(1.0)) 13.;
  190. *
  191. TAB1 = TABLE;
  192. TAB1.'MODELE' = MOD1;
  193. TAB1.'CARACTERISTIQUES' = MAT1;
  194. TAB1.'CHARGEMENT' = CHA1 ET CHA2;
  195. TAB1.'TEMPS_CALCULES' = LLT3;
  196. TAB1.'TEMPS_SAUVES' = LLT4;
  197. TAB1.'CONVERGENCE_FORCEE' = FAUX;
  198. TAB1.'BLOCAGES_MECANIQUES' = BL1 ET BL2 ET BL3 ET BL4;
  199. TAB1.'MES_SAUVEGARDES' = TABLE;
  200. TAB1.'MES_SAUVEGARDES'.'DEFIN' = VRAI;
  201. TAB1.'MOVA' = 'GAM1';
  202. TAB1.'DEPLACEMENTS' = TABLE;
  203. TAB1.'DEPLACEMENTS'. 0 = CH_INI;
  204. TAB1.'CONTRAINTES' = TABLE;
  205. TAB1.'CONTRAINTES'. 0 = SIG_INI;
  206. TAB1.'PRECISION' = 1E-4;
  207. *
  208. PASAPAS TAB1;
  209. *
  210. *********************************************************************
  211. * Posttraitement
  212. *********************************************************************
  213. *
  214. L_T1 = PROG;
  215. L_EXX = PROG;
  216. L_EYY = PROG;
  217. L_EZZ = PROG;
  218. L_EXY = PROG;
  219. L_EXZ = PROG;
  220. L_EYZ = PROG;
  221. L_SXX = PROG;
  222. L_SYY = PROG;
  223. L_SZZ = PROG;
  224. L_SXY = PROG;
  225. L_SXZ = PROG;
  226. L_SYZ = PROG;
  227. *
  228. I = 0;
  229. NN1 = DIME (TAB1.DEPLACEMENTS);
  230. *
  231. REPE IND1 (NN1);
  232. *
  233. L_T1 = L_T1 ET (PROG (TAB1.TEMPS. I));
  234. *
  235. DEP1 = TAB1.DEPLACEMENTS. I;
  236. EPS1 = EPSI DEP1 MOD1;
  237. SIG1 = TAB1.CONTRAINTES. I;
  238. *
  239. L_EXX = L_EXX ET (PROG (EXTR EPS1 'EPXX' 1 1 1));
  240. L_EYY = L_EYY ET (PROG (EXTR EPS1 'EPYY' 1 1 1));
  241. L_EZZ = L_EZZ ET (PROG (EXTR EPS1 'EPZZ' 1 1 1));
  242. L_EXY = L_EXY ET (PROG (EXTR EPS1 'GAXY' 1 1 1));
  243. L_EXZ = L_EXZ ET (PROG (EXTR EPS1 'GAXZ' 1 1 1));
  244. L_EYZ = L_EYZ ET (PROG (EXTR EPS1 'GAYZ' 1 1 1));
  245. *
  246. L_SXX = L_SXX ET (PROG (EXTR SIG1 'SMXX' 1 1 1));
  247. L_SYY = L_SYY ET (PROG (EXTR SIG1 'SMYY' 1 1 1));
  248. L_SZZ = L_SZZ ET (PROG (EXTR SIG1 'SMZZ' 1 1 1));
  249. L_SXY = L_SXY ET (PROG (EXTR SIG1 'SMXY' 1 1 1));
  250. L_SXZ = L_SXZ ET (PROG (EXTR SIG1 'SMXZ' 1 1 1));
  251. L_SYZ = L_SYZ ET (PROG (EXTR SIG1 'SMYZ' 1 1 1));
  252. *
  253. I = I + 1
  254. FIN IND1;
  255. *
  256. EV_EXX = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Exx []' L_EXX;
  257. EV_EYY = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Eyy []' L_EYY;
  258. EV_EZZ = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Ezz []' L_EZZ;
  259. EV_EXY = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Exy []' L_EXY;
  260. EV_EXZ = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Exz []' L_EXZ;
  261. EV_EYZ = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Eyz []' L_EYZ;
  262. *
  263. EV_SXX = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Sxx [Pa]' L_SXX;
  264. EV_SYY = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Syy [Pa]' L_SYY;
  265. EV_SZZ = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Szz [Pa]' L_SZZ;
  266. EV_SXY = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Sxy [Pa]' L_SXY;
  267. EV_SXZ = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Sxz [Pa]' L_SXZ;
  268. EV_SYZ = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Syz [Pa]' L_SYZ;
  269. *
  270. EV_1 = EVOL (ROSE) MANU 'Exz []' L_EXZ 'Sxz [Pa]' L_SXZ;
  271. EV_2 = EVOL (ROSE) MANU 'Exz []' L_EXZ 'p [Pa]'
  272. ((L_SXX + L_SYY + L_SZZ)/(3.0));
  273. *
  274. SI (VAR_GRAP EGA 1);
  275. DESS EV_EXX 'TITR' 'Exx';
  276. DESS EV_EYY 'TITR' 'Eyy';
  277. DESS EV_EZZ 'TITR' 'Ezz';
  278. DESS EV_EXY 'TITR' 'Exy';
  279. DESS EV_EXZ 'TITR' 'Exz';
  280. DESS EV_EYZ 'TITR' 'Eyz';
  281. *
  282. DESS EV_SXX 'TITR' 'Sxx';
  283. DESS EV_SYY 'TITR' 'Syy';
  284. DESS EV_SZZ 'TITR' 'Szz';
  285. DESS EV_SXY 'TITR' 'Sxy';
  286. DESS EV_SXZ 'TITR' 'Sxz';
  287. DESS EV_SYZ 'TITR' 'Syz';
  288. *
  289. DESS EV_1 'TITR' 'Exz - Sxz';
  290. DESS EV_2 'TITR' 'Exz - p';
  291. *
  292. FINSI;
  293. *
  294. *********************************************************************
  295. * Contruction de la solution de reference
  296. *********************************************************************
  297. *
  298. L_SXZT = (PROG
  299. 0.00000000e+00 8.00000000e+01 1.60000000e+02
  300. 2.12839070e+02 2.60203100e+02 3.05291070e+02
  301. 3.52371660e+02 3.96740470e+02 4.40781330e+02
  302. 4.84086590e+02 5.27096150e+02 5.68406540e+02
  303. 6.08937740e+02 6.49230570e+02 6.88908370e+02);
  304. L_SXZT = L_SXZT ET (PROG
  305. 7.28652200e+02 7.67707310e+02 8.05779040e+02
  306. 8.43343830e+02 8.80292300e+02 9.16612880e+02
  307. 9.50543070e+02 9.84382610e+02 1.01808480e+03
  308. 1.05141200e+03 1.08465590e+03 1.11536850e+03
  309. 1.14532130e+03 1.17519080e+03 1.20491190e+03);
  310. L_SXZT = L_SXZT ET (PROG
  311. 1.23433030e+03 1.26367880e+03 1.28733980e+03
  312. 1.30749550e+03 1.32755070e+03 1.34750650e+03
  313. 1.36735080e+03 1.38697860e+03 1.40652480e+03
  314. 1.42600100e+03 1.44539400e+03 1.46449510e+03
  315. 1.48114470e+03 1.49774810e+03 1.51431050e+03);
  316. L_SXZT = L_SXZT ET (PROG
  317. 1.53083650e+03 1.54733000e+03 1.56379460e+03
  318. 1.58023360e+03 1.59664980e+03 1.61304570e+03
  319. 1.62942370e+03 1.64578580e+03 1.66213370e+03
  320. 1.67846920e+03 1.69479360e+03 1.71110830e+03
  321. 1.72741440e+03 1.74371300e+03 1.76000510e+03);
  322. L_SXZT = L_SXZT ET (PROG
  323. 1.77629150e+03 1.79257290e+03 1.80885010e+03
  324. 1.82512380e+03 1.84139450e+03 1.85766270e+03
  325. 1.87392890e+03 1.89019360e+03 1.90645720e+03
  326. 1.92272000e+03 1.93898240e+03 1.95524480e+03
  327. 1.97150730e+03 1.98777040e+03 2.00403420e+03);
  328. L_SXZT = L_SXZT ET (PROG
  329. 2.02029900e+03 2.03656490e+03 2.05283240e+03
  330. 2.06910140e+03 2.08537220e+03 2.10164500e+03
  331. 2.11791990e+03 2.13419710e+03 2.15047680e+03
  332. 2.16675900e+03 2.18304400e+03 2.19933170e+03
  333. 2.21562250e+03 2.23191620e+03 2.24821320e+03);
  334. L_SXZT = L_SXZT ET (PROG
  335. 2.26451340e+03 2.28081710e+03 2.29712420e+03
  336. 2.31343480e+03 2.32974920e+03 2.34606720e+03
  337. 2.36238910e+03 2.37871500e+03 2.39504470e+03
  338. 2.41137860e+03 2.42771660e+03 2.34778880e+03
  339. 2.26786090e+03 2.18793310e+03 2.10800520e+03);
  340. L_SXZT = L_SXZT ET (PROG
  341. 2.04283480e+03 1.99707210e+03 1.95141180e+03
  342. 1.90572340e+03 1.85994860e+03 1.81407540e+03
  343. 1.76811290e+03 1.72334380e+03 1.67878060e+03
  344. 1.63569170e+03 1.59266830e+03 1.54966200e+03
  345. 1.50663980e+03 1.46357950e+03 1.42046400e+03);
  346. L_SXZT = L_SXZT ET (PROG
  347. 1.37722510e+03 1.33653190e+03 1.29600770e+03
  348. 1.25607350e+03 1.21616880e+03 1.17627200e+03
  349. 1.13636690e+03 1.09644030e+03 1.05647860e+03
  350. 1.01643240e+03 9.78694220e+02 9.39895240e+02
  351. 9.02095010e+02 8.64372590e+02 8.25892880e+02);
  352. L_SXZT = L_SXZT ET (PROG
  353. 7.90553610e+02 7.53565560e+02 7.17343480e+02
  354. 6.81188090e+02 6.44171820e+02 6.10099840e+02
  355. 5.75797600e+02 5.42109140e+02 5.08434210e+02
  356. 4.74767080e+02 4.41103010e+02 4.07437980e+02
  357. 3.73768270e+02 3.40089410e+02 3.06390860e+02);
  358. L_SXZT = L_SXZT ET (PROG
  359. 2.73862150e+02 2.41176290e+02 2.10389400e+02
  360. 1.80620380e+02 1.50863640e+02 1.21115570e+02
  361. 9.13731080e+01 6.16336650e+01 3.18949020e+01
  362. 2.15442210e+00 -2.75911010e+01 -5.73508650e+01
  363. -8.59856250e+01 -1.14697550e+02 -1.39112180e+02);
  364. L_SXZT = L_SXZT ET (PROG
  365. -1.58847760e+02 -1.78571410e+02 -1.98285320e+02
  366. -2.17991280e+02 -2.37690740e+02 -2.57384890e+02
  367. -2.77074670e+02 -2.96760870e+02 -3.16444120e+02
  368. -3.36125010e+02 -3.55804130e+02 -3.75482230e+02
  369. -3.95160660e+02 -4.14843900e+02 -4.33736680e+02);
  370. L_SXZT = L_SXZT ET (PROG
  371. -4.52617800e+02 -4.71540760e+02 -4.90675650e+02
  372. -5.06641090e+02 -5.22597350e+02 -5.38546040e+02
  373. -5.54488440e+02 -5.70425560e+02 -5.86358220e+02
  374. -6.02287050e+02 -6.18212540e+02 -6.34135100e+02
  375. -6.50055000e+02 -6.65972480e+02 -6.81887700e+02);
  376. L_SXZT = L_SXZT ET (PROG
  377. -6.97800790e+02 -7.13711830e+02 -7.29620890e+02
  378. -7.45528010e+02 -7.61433210e+02 -7.77336520e+02
  379. -7.93237930e+02 -8.09137470e+02 -8.25035130e+02
  380. -8.40930920e+02 -8.56824840e+02 -8.72716900e+02
  381. -8.88607110e+02 -9.04495490e+02 -9.20382050e+02);
  382. L_SXZT = L_SXZT ET (PROG
  383. -9.36266820e+02 -9.52149830e+02 -9.68031110e+02
  384. -9.83910680e+02 -9.99788610e+02 -1.01566490e+03
  385. -1.03153970e+03 -1.04741290e+03 -1.06328470e+03
  386. -1.07915510e+03 -1.09502410e+03 -1.11089190e+03
  387. -1.12675850e+03 -1.14262390e+03 -1.15848830e+03);
  388. L_SXZT = L_SXZT ET (PROG
  389. -1.17435170e+03 -1.19021430e+03 -1.20607590e+03
  390. -1.22193690e+03 -1.23779720e+03 -1.25365690e+03
  391. -1.26951610e+03 -1.28537500e+03 -1.30123350e+03
  392. -1.31709170e+03 -1.33294990e+03 -1.34880790e+03
  393. -1.36466600e+03 -1.38052430e+03 -1.39638270e+03);
  394. L_SXZT = L_SXZT ET (PROG
  395. -1.41224150e+03 -1.42810060e+03 -1.44396020e+03
  396. -1.45982040e+03 -1.47568130e+03 -1.49154290e+03
  397. -1.50740530e+03 -1.52326870e+03 -1.53913310e+03
  398. -1.55499850e+03 -1.57086520e+03 -1.58673320e+03
  399. -1.60260250e+03 -1.61847320e+03 -1.63434550e+03);
  400. L_SXZT = L_SXZT ET (PROG
  401. -1.65021940e+03 -1.66609500e+03 -1.68197240e+03
  402. -1.69785160e+03 -1.71373280e+03 -1.72961600e+03
  403. -1.74550130e+03 -1.76138880e+03 -1.77727860e+03
  404. -1.79317070e+03 -1.80906530e+03 -1.82496230e+03
  405. -1.84086190e+03 -1.85676410e+03 -1.87266910e+03);
  406. L_SXZT = L_SXZT ET (PROG
  407. -1.88857690e+03 -1.90448750e+03 -1.92040100e+03
  408. -1.93631760e+03 -1.95223720e+03 -1.96816000e+03
  409. -1.98408600e+03 -2.00001530e+03 -2.01594790e+03
  410. -2.03188390e+03 -2.04782340e+03 -2.06376640e+03
  411. -2.07971310e+03 -2.09566340e+03 -2.11161740e+03);
  412. L_SXZT = L_SXZT ET (PROG
  413. -2.12757530e+03 -2.14353690e+03 -2.15950250e+03
  414. -2.17547210e+03 -2.19144560e+03 -2.20742330e+03
  415. -2.22340510e+03 -2.23939110e+03 -2.25538130e+03
  416. -2.27137580e+03 -2.28737470e+03 -2.30337800e+03
  417. -2.31938580e+03 -2.33539800e+03 -2.35141480e+03);
  418. L_SXZT = L_SXZT ET (PROG
  419. -2.36743630e+03 -2.28770520e+03 -2.20797410e+03
  420. -2.12824300e+03 -2.04851200e+03 -1.98715270e+03
  421. -1.94203080e+03 -1.89699580e+03 -1.85191550e+03
  422. -1.80673460e+03 -1.76144450e+03 -1.71604710e+03
  423. -1.67222980e+03 -1.62870880e+03 -1.58623760e+03);
  424. L_SXZT = L_SXZT ET (PROG
  425. -1.54383720e+03 -1.50145170e+03 -1.45904280e+03
  426. -1.41658420e+03 -1.37404960e+03 -1.33258950e+03
  427. -1.29085750e+03 -1.25150480e+03 -1.21218240e+03
  428. -1.17286970e+03 -1.13355120e+03 -1.09421490e+03
  429. -1.05485040e+03 -1.01544110e+03 -9.77155120e+02);
  430. L_SXZT = L_SXZT ET (PROG
  431. -9.38238380e+02 -9.00744510e+02 -8.63320810e+02
  432. -8.25867240e+02 -7.89527150e+02 -7.52675570e+02
  433. -7.16763710e+02 -6.80914190e+02 -6.45012980e+02
  434. -6.10394570e+02 -5.75611400e+02 -5.42467990e+02
  435. -5.09337000e+02 -4.76212450e+02 -4.43089440e+02);
  436. L_SXZT = L_SXZT ET (PROG
  437. -4.09963820e+02 -3.76831860e+02 -3.43689320e+02
  438. -3.10528210e+02 -2.77286810e+02 -2.46141910e+02
  439. -2.15038800e+02 -1.85769100e+02 -1.56509860e+02
  440. -1.27257250e+02 -9.80080750e+01 -6.87596320e+01
  441. -3.95094940e+01 -1.02552060e+01 1.90065120e+01);
  442. L_SXZT = L_SXZT ET (PROG
  443. 4.82831510e+01 7.66897660e+01 1.05112840e+02
  444. 1.31754240e+02 1.51123650e+02 1.70483010e+02
  445. 1.89834660e+02 2.09180480e+02 2.28522030e+02
  446. 2.47860530e+02 2.67196970e+02 2.86532150e+02
  447. 3.05866750e+02 3.25201370e+02 3.44536640e+02);
  448. L_SXZT = L_SXZT ET (PROG
  449. 3.63873390e+02 3.83213100e+02 4.02559970e+02
  450. 4.21212020e+02 4.39844700e+02 4.58502200e+02
  451. 4.77271770e+02 4.93693170e+02 5.09384420e+02
  452. 5.25069860e+02 5.40750690e+02 5.56427860e+02
  453. 5.72102140e+02 5.87774110e+02 6.03444230e+02);
  454. L_SXZT = L_SXZT ET (PROG
  455. 6.19112860e+02 6.34780260e+02 6.50446620e+02
  456. 6.66112100e+02 6.81776800e+02 6.97440800e+02
  457. 7.13104130e+02 7.28766830e+02 7.44428910e+02
  458. 7.60090380e+02 7.75751250e+02 7.91411500e+02
  459. 8.07071140e+02 8.22730160e+02 8.38388560e+02);
  460. L_SXZT = L_SXZT ET (PROG
  461. 8.54046330e+02 8.69703490e+02 8.85360040e+02
  462. 9.01015980e+02 9.16671330e+02 9.32326120e+02
  463. 9.47980360e+02 9.63634080e+02 9.79287300e+02
  464. 9.94940070e+02 1.01059240e+03 1.02624440e+03
  465. 1.04189600e+03 1.05754740e+03 1.07319850e+03);
  466. L_SXZT = L_SXZT ET (PROG
  467. 1.08884950e+03 1.10450030e+03 1.12015110e+03
  468. 1.13580180e+03 1.15145270e+03 1.16710360e+03
  469. 1.18275470e+03 1.19840610e+03 1.21405780e+03
  470. 1.22970990e+03 1.24536250e+03 1.26101560e+03
  471. 1.27666940e+03 1.29232380e+03 1.30797890e+03);
  472. L_SXZT = L_SXZT ET (PROG
  473. 1.32363500e+03 1.33929190e+03 1.35494980e+03
  474. 1.37060870e+03 1.38626880e+03 1.40193020e+03
  475. 1.41759280e+03 1.43325670e+03 1.44892210e+03
  476. 1.46458900e+03 1.48025760e+03 1.49592770e+03
  477. 1.51159960e+03 1.52727340e+03 1.54294900e+03);
  478. L_SXZT = L_SXZT ET (PROG
  479. 1.55862660e+03 1.57430620e+03 1.58998790e+03
  480. 1.60567180e+03 1.62135800e+03 1.63704640e+03
  481. 1.65273730e+03 1.66843070e+03 1.68412650e+03
  482. 1.69982500e+03 1.71552620e+03 1.73123010e+03
  483. 1.74693680e+03 1.76264640e+03 1.77835890e+03);
  484. L_SXZT = L_SXZT ET (PROG
  485. 1.79407440e+03 1.80979300e+03 1.82551470e+03
  486. 1.84123970e+03 1.85696790e+03 1.87269940e+03
  487. 1.88843430e+03 1.90417260e+03 1.91991450e+03
  488. 1.93565990e+03 1.95140900e+03 1.96716170e+03
  489. 1.98291820e+03 1.99867850e+03 2.01444260e+03);
  490. L_SXZT = L_SXZT ET (PROG
  491. 2.03021070e+03 2.04598270e+03 2.06175870e+03
  492. 2.07753880e+03 2.09332310e+03 2.10911150e+03
  493. 2.12490410e+03 2.14070110e+03 2.15650230e+03
  494. 2.17230800e+03 2.18811810e+03 2.20393260e+03
  495. 2.21975170e+03 2.23557540e+03 2.25140370e+03);
  496. L_SXZT = L_SXZT ET (PROG
  497. 2.26723670e+03 2.28307440e+03 2.29891680e+03
  498. 2.31476410e+03 2.33061620e+03 2.34647320e+03
  499. 2.26686710e+03 2.18726110e+03 2.10765510e+03
  500. 2.02804910e+03 1.96912920e+03 1.92442230e+03
  501. 1.87979080e+03 1.83510320e+03 1.79030690e+03);
  502. L_SXZT = L_SXZT ET (PROG
  503. 1.74539470e+03 1.70035570e+03 1.65720720e+03
  504. 1.61436420e+03 1.57229270e+03 1.53029560e+03
  505. 1.48831170e+03 1.44629900e+03 1.40422740e+03
  506. 1.36205970e+03 1.32127700e+03 1.28050720e+03
  507. 1.24154620e+03 1.20261670e+03 1.16369580e+03);
  508. L_SXZT = L_SXZT ET (PROG
  509. 1.12476620e+03 1.08581420e+03 1.04682690e+03
  510. 1.00777620e+03 9.70239200e+02 9.31996960e+02
  511. 8.95022090e+02 8.58123340e+02 8.21076800e+02
  512. 7.86004300e+02 7.49958440e+02 7.14576920e+02
  513. 6.79260390e+02 6.43705940e+02 6.10360700e+02);
  514. L_SXZT = L_SXZT ET (PROG
  515. 5.76674250e+02 5.43878600e+02 5.11094190e+02
  516. 4.78314260e+02 4.45533230e+02 4.12746320e+02
  517. 3.79948990e+02 3.47135530e+02 3.14292570e+02
  518. 2.82532720e+02 2.50665450e+02 2.20351380e+02
  519. 1.91388230e+02 1.62433480e+02 1.33482850e+02);
  520. L_SXZT = L_SXZT ET (PROG
  521. 1.04532750e+02 7.55800860e+01 4.66219880e+01
  522. 1.76553310e+01 -1.13246870e+01 -4.03327340e+01
  523. -6.81639330e+01 -9.60781200e+01 -1.19788730e+02
  524. -1.38951110e+02 -1.58105570e+02 -1.77254600e+02
  525. -1.96400220e+02 -2.15544050e+02 -2.34687400e+02);
  526. L_SXZT = L_SXZT ET (PROG
  527. -2.53831310e+02 -2.72976650e+02 -2.92124160e+02
  528. -3.11274530e+02 -3.30428550e+02 -3.49587380e+02
  529. -3.68753360e+02 -3.87935860e+02 -4.06257040e+02
  530. -4.24571520e+02 -4.42940410e+02 -4.61432070e+02
  531. -4.77005150e+02 -4.92573330e+02 -5.08138020e+02);
  532. L_SXZT = L_SXZT ET (PROG
  533. -5.23700300e+02 -5.39261030e+02 -5.54820900e+02
  534. -5.70380420e+02 -5.85939980e+02 -6.01499890e+02
  535. -6.17060360e+02 -6.32621560e+02 -6.48183590e+02
  536. -6.63746530e+02 -6.79310410e+02 -6.94875260e+02
  537. -7.10441080e+02 -7.26007860e+02 -7.41575600e+02);
  538. L_SXZT = L_SXZT ET (PROG
  539. -7.57144260e+02 -7.72713840e+02 -7.88284300e+02
  540. -8.03855620e+02 -8.19427790e+02 -8.35000780e+02
  541. -8.50574580e+02 -8.66149180e+02 -8.81724580e+02
  542. -8.97300770e+02 -9.12877750e+02 -9.28455530e+02
  543. -9.44034110e+02 -9.59613510e+02 -9.75193760e+02);
  544. L_SXZT = L_SXZT ET (PROG
  545. -9.90774850e+02 -1.00635680e+03 -1.02193970e+03
  546. -1.03752360e+03 -1.05310840e+03 -1.06869420e+03
  547. -1.08428100e+03 -1.09986900e+03 -1.11545810e+03
  548. -1.13104840e+03 -1.14663980e+03 -1.16223260e+03
  549. -1.17782670e+03 -1.19342220e+03 -1.20901910e+03);
  550. L_SXZT = L_SXZT ET (PROG
  551. -1.22461740e+03 -1.24021730e+03 -1.25581880e+03
  552. -1.27142200e+03 -1.28702680e+03 -1.30263350e+03
  553. -1.31824190e+03 -1.33385220e+03 -1.34946450e+03
  554. -1.36507880e+03 -1.38069510e+03 -1.39631350e+03
  555. -1.41193420e+03 -1.42755700e+03 -1.44318220e+03);
  556. L_SXZT = L_SXZT ET (PROG
  557. -1.45880970e+03 -1.47443970e+03 -1.49007210e+03
  558. -1.50570710e+03 -1.52134460e+03 -1.53698480e+03
  559. -1.55262780e+03 -1.56827350e+03 -1.58392200e+03
  560. -1.59957350e+03 -1.61522780e+03 -1.63088520e+03
  561. -1.64654570e+03 -1.66220930e+03 -1.67787600e+03);
  562. L_SXZT = L_SXZT ET (PROG
  563. -1.69354600e+03 -1.70921930e+03 -1.72489590e+03
  564. -1.74057590e+03 -1.75625940e+03 -1.77194640e+03
  565. -1.78763700e+03 -1.80333110e+03 -1.81902900e+03
  566. -1.83473050e+03 -1.85043580e+03 -1.86614500e+03
  567. -1.88185800e+03 -1.89757490e+03 -1.91329580e+03);
  568. L_SXZT = L_SXZT ET (PROG
  569. -1.92902080e+03 -1.94474980e+03 -1.96048290e+03
  570. -1.97622020e+03 -1.99196170e+03 -2.00770750e+03
  571. -2.02345750e+03 -2.03921200e+03 -2.05497080e+03
  572. -2.07073410e+03 -2.08650180e+03 -2.10227410e+03
  573. -2.11805100e+03 -2.13383240e+03 -2.14961860e+03);
  574. L_SXZT = L_SXZT ET (PROG
  575. -2.16540940e+03 -2.18120500e+03 -2.19700540e+03
  576. -2.21281060e+03 -2.22862060e+03 -2.24443560e+03
  577. -2.26025550e+03 -2.27608040e+03 -2.29191030e+03
  578. -2.30774530e+03 -2.32358530e+03 -2.24402140e+03
  579. -2.16445750e+03 -2.08489350e+03 -2.00532960e+03);
  580. L_SXZT = L_SXZT ET (PROG
  581. -1.94724250e+03 -1.90268340e+03 -1.85819590e+03
  582. -1.81364970e+03 -1.76899320e+03 -1.72421920e+03
  583. -1.67930860e+03 -1.63641850e+03 -1.59387320e+03
  584. -1.55194360e+03 -1.51008990e+03 -1.46824920e+03
  585. -1.42637790e+03 -1.38444480e+03 -1.34240720e+03);
  586. L_SXZT = L_SXZT ET (PROG
  587. -1.30189260e+03 -1.26146370e+03 -1.22264030e+03
  588. -1.18384860e+03 -1.14506490e+03 -1.10627130e+03
  589. -1.06745340e+03 -1.02859700e+03 -9.89668210e+02
  590. -9.52430140e+02 -9.14374280e+02 -8.77582230e+02
  591. -8.40867130e+02 -8.03888640e+02 -7.69333970e+02);
  592. L_SXZT = L_SXZT ET (PROG
  593. -7.33442020e+02 -6.98245910e+02 -6.63113670e+02
  594. -6.27509130e+02 -5.94579260e+02 -5.61227250e+02
  595. -5.28539390e+02 -4.95861270e+02 -4.63185870e+02
  596. -4.30507380e+02 -3.97820770e+02 -3.65121180e+02
  597. -3.32402260e+02 -2.99647780e+02 -2.68105190e+02);
  598. L_SXZT = L_SXZT ET (PROG
  599. -2.36402690e+02 -2.06501080e+02 -1.77617570e+02
  600. -1.48740500e+02 -1.19865460e+02 -9.09887150e+01
  601. -6.21070360e+01 -3.32173970e+01 -4.31641630e+00
  602. 2.46013210e+01 5.35544060e+01 8.11996960e+01
  603. 1.08966280e+02 1.31829330e+02 1.50957680e+02);
  604. L_SXZT = L_SXZT ET (PROG
  605. 1.70080070e+02 1.89199010e+02 2.08316540e+02
  606. 2.27434320e+02 2.46553680e+02 2.65675690e+02
  607. 2.84801230e+02 3.03931080e+02 3.23065960e+02
  608. 3.42206720e+02 3.61354610e+02 3.80512310e+02
  609. 3.99692140e+02 4.17947610e+02 4.36204660e+02);
  610. L_SXZT = L_SXZT ET (PROG
  611. 4.54530790e+02 4.72728290e+02 4.88315540e+02
  612. 5.03900400e+02 5.19484110e+02 5.35067680e+02
  613. 5.50651890e+02 5.66237330e+02 5.81824450e+02
  614. 5.97413600e+02 6.13005030e+02 6.28598900e+02
  615. 6.44195350e+02 6.59794440e+02 6.75396210e+02);
  616. L_SXZT = L_SXZT ET (PROG
  617. 6.91000670e+02 7.06607820e+02 7.22217630e+02
  618. 7.37830080e+02 7.53445120e+02 7.69062710e+02
  619. 7.84682800e+02 8.00305360e+02 8.15930350e+02
  620. 8.31557710e+02 8.47187410e+02 8.62819430e+02
  621. 8.78453720e+02 8.94090260e+02 9.09729030e+02);
  622. L_SXZT = L_SXZT ET (PROG
  623. 9.25370020e+02 9.41013200e+02 9.56658580e+02
  624. 9.72306140e+02 9.87955880e+02 1.00360780e+03
  625. 1.01926190e+03 1.03491820e+03 1.05057680e+03
  626. 1.06623750e+03 1.08190050e+03 1.09756580e+03
  627. 1.11323340e+03 1.12890330e+03 1.14457550e+03);
  628. L_SXZT = L_SXZT ET (PROG
  629. 1.16025010e+03 1.17592710e+03 1.19160650e+03
  630. 1.20728840e+03 1.22297280e+03 1.23865970e+03
  631. 1.25434930e+03 1.27004140e+03 1.28573620e+03
  632. 1.30143370e+03 1.31713400e+03 1.33283700e+03
  633. 1.34854280e+03 1.36425160e+03 1.37996320e+03);
  634. L_SXZT = L_SXZT ET (PROG
  635. 1.39567780e+03 1.41139540e+03 1.42711600e+03
  636. 1.44283970e+03 1.45856660e+03 1.47429660e+03
  637. 1.49002990e+03 1.50576640e+03 1.52150620e+03
  638. 1.53724940e+03 1.55299600e+03 1.56874610e+03
  639. 1.58449970e+03 1.60025670e+03 1.61601740e+03);
  640. L_SXZT = L_SXZT ET (PROG
  641. 1.63178170e+03 1.64754970e+03 1.66332140e+03
  642. 1.67909690e+03 1.69487610e+03 1.71065930e+03
  643. 1.72644630e+03 1.74223720e+03 1.75803220e+03
  644. 1.77383110e+03 1.78963410e+03 1.80544130e+03
  645. 1.82125250e+03 1.83706800e+03 1.85288770e+03);
  646. L_SXZT = L_SXZT ET (PROG
  647. 1.86871170e+03 1.88453990e+03 1.90037260e+03
  648. 1.91620960e+03 1.93205110e+03 1.94789700e+03
  649. 1.96374740e+03 1.97960240e+03 1.99546200e+03
  650. 2.01132620e+03 2.02719510e+03 2.04306860e+03
  651. 2.05894690e+03 2.07483000e+03 2.09071790e+03);
  652. L_SXZT = L_SXZT ET (PROG
  653. 2.10661070e+03 2.12250830e+03 2.13841080e+03
  654. 2.15431830e+03 2.17023080e+03 2.18614830e+03
  655. 2.20207080e+03 2.21799850e+03 2.23393120e+03
    2.24986910e+03 2.26581220e+03 2.28176060e+03
  656. 2.29771410e+03 2.31367300e+03 2.32963720e+03);
  657. L_SXZT = L_SXZT ET (PROG
  658. 2.34560670e+03 2.26600390e+03 2.18640120e+03
  659. 2.10679840e+03 2.02719570e+03 1.96836960e+03
  660. 1.92369750e+03 1.87910140e+03 1.83445220e+03
  661. 1.78969740e+03 1.74482860e+03 1.69983100e+03
  662. 1.65678260e+03 1.61403870e+03 1.57199880e+03);
  663. L_SXZT = L_SXZT ET (PROG
  664. 1.53003390e+03 1.48808340e+03 1.44610490e+03
  665. 1.40406830e+03 1.36193570e+03 1.32124930e+03
  666. 1.28054770e+03 1.24161310e+03 1.20270940e+03
  667. 1.16381380e+03 1.12490870e+03 1.08598040e+03
  668. 1.04701570e+03 1.00798600e+03 9.70539260e+02);
  669. L_SXZT = L_SXZT ET (PROG
  670. 9.32441050e+02 8.95467630e+02 8.58569690e+02
  671. 8.21507750e+02 7.86565740e+02 7.50602190e+02
  672. 7.15206550e+02 6.79874470e+02 6.44277870e+02
  673. 6.11031600e+02 5.77417810e+02 5.44603120e+02
  674. 5.11796710e+02 4.78991810e+02 4.46182840e+02);
  675. L_SXZT = L_SXZT ET (PROG
  676. 4.13364990e+02 3.80533700e+02 3.47683300e+02
  677. 3.14800650e+02 2.83037460e+02 2.51147750e+02
  678. 2.20775160e+02 1.91752230e+02 1.62734090e+02
  679. 1.33716480e+02 1.04695840e+02 7.56690770e+01
  680. 4.66333630e+01 1.75856180e+01 -1.14787850e+01);
  681. L_SXZT = L_SXZT ET (PROG
  682. -4.05736840e+01 -6.84719410e+01 -9.64657050e+01
  683. -1.20267270e+02 -1.39526190e+02 -1.58780670e+02
  684. -1.78033180e+02 -1.97285680e+02 -2.16539800e+02
  685. -2.35796840e+02 -2.55057840e+02 -2.74323680e+02
  686. -2.93595080e+02 -3.12872730e+02 -3.32157410e+02);
  687. L_SXZT = L_SXZT ET (PROG
  688. -3.51450180e+02 -3.70753160e+02 -3.90074380e+02
  689. -4.08535530e+02 -4.26997830e+02 -4.45519570e+02
  690. -4.64186830e+02 -4.79917950e+02 -4.95648440e+02
  691. -5.11379550e+02 -5.27112260e+02 -5.42847320e+02
  692. -5.58585300e+02 -5.74326650e+02 -5.90071690e+02);
  693. L_SXZT = L_SXZT ET (PROG
  694. -6.05820640e+02 -6.21573650e+02 -6.37330830e+02
  695. -6.53092220e+02 -6.68857850e+02 -6.84627700e+02
  696. -7.00401750e+02 -7.16179950e+02 -7.31962240e+02
  697. -7.47748580e+02 -7.63538890e+02 -7.79333110e+02
  698. -7.95131170e+02 -8.10933020e+02 -8.26738580e+02);
  699. L_SXZT = L_SXZT ET (PROG
  700. -8.42547800e+02 -8.58360630e+02 -8.74177020e+02
  701. -8.89996920e+02 -9.05820280e+02 -9.21647080e+02
  702. -9.37477270e+02 -9.53310830e+02 -9.69147740e+02
  703. -9.84987980e+02 -1.00083150e+03 -1.01667840e+03
  704. -1.03252850e+03 -1.04838190e+03 -1.06423860e+03);
  705. L_SXZT = L_SXZT ET (PROG
  706. -1.08009860e+03 -1.09596190e+03 -1.11182840e+03
  707. -1.12769830e+03 -1.14357150e+03 -1.15944790e+03
  708. -1.17532780e+03 -1.19121100e+03 -1.20709760e+03
  709. -1.22298750e+03 -1.23888090e+03 -1.25477770e+03
  710. -1.27067800e+03 -1.28658180e+03 -1.30248910e+03);
  711. L_SXZT = L_SXZT ET (PROG
  712. -1.31840000e+03 -1.33431440e+03 -1.35023240e+03
  713. -1.36615400e+03 -1.38207940e+03 -1.39800840e+03
  714. -1.41394110e+03 -1.42987760e+03 -1.44581780e+03
  715. -1.46176190e+03 -1.47770980e+03 -1.49366160e+03
  716. -1.50961740e+03 -1.52557700e+03 -1.54154070e+03);
  717. L_SXZT = L_SXZT ET (PROG
  718. -1.55750840e+03 -1.57348010e+03 -1.58945590e+03
  719. -1.60543580e+03 -1.62141990e+03 -1.63740820e+03
  720. -1.65340060e+03 -1.66939740e+03 -1.68539840e+03
  721. -1.70140370e+03 -1.71741340e+03 -1.73342750e+03
  722. -1.74944600e+03 -1.76546890e+03 -1.78149630e+03);
  723. L_SXZT = L_SXZT ET (PROG
  724. -1.79752830e+03 -1.81356480e+03 -1.82960580e+03
  725. -1.84565150e+03 -1.86170190e+03 -1.87775690e+03
  726. -1.89381660e+03 -1.90988110e+03 -1.92595030e+03
  727. -1.94202430e+03 -1.95810320e+03 -1.97418700e+03
  728. -1.99027560e+03 -2.00636920e+03 -2.02246770e+03);
  729. L_SXZT = L_SXZT ET (PROG
  730. -2.03857120e+03 -2.05467980e+03 -2.07079330e+03
  731. -2.08691200e+03 -2.10303580e+03 -2.11916470e+03
  732. -2.13529870e+03 -2.15143800e+03 -2.16758240e+03
  733. -2.18373210e+03 -2.19988710e+03 -2.21604740e+03
  734. -2.23221300e+03 -2.24838400e+03 -2.26456040e+03);
  735. L_SXZT = L_SXZT ET (PROG
  736. -2.28074220e+03 -2.29692940e+03 -2.31312200e+03
  737. -2.32932020e+03 -2.34552390e+03 -2.36173310e+03
  738. -2.28202410e+03 -2.20231510e+03 -2.12260620e+03
  739. -2.04289720e+03 -1.98202240e+03 -1.93701810e+03
  740. -1.89209980e+03 -1.84713990e+03 -1.80208420e+03);
  741. L_SXZT = L_SXZT ET (PROG
  742. -1.75692230e+03 -1.71164920e+03 -1.66810010e+03
  743. -1.62473930e+03 -1.58237790e+03 -1.54008890e+03
  744. -1.49781650e+03 -1.45552160e+03 -1.41317730e+03
  745. -1.37075480e+03 -1.32954940e+03 -1.28811140e+03
  746. -1.24885820e+03 -1.20963430e+03 -1.17041880e+03);
  747. L_SXZT = L_SXZT ET (PROG
  748. -1.13119570e+03 -1.09195280e+03 -1.05267870e+03
  749. -1.01335530e+03 -9.75310730e+02 -9.36713540e+02
  750. -8.99287670e+02 -8.61931690e+02 -8.24528580e+02
  751. -7.88531980e+02 -7.51927490e+02 -7.16057670e+02
  752. -6.80248230e+02 -6.44362450e+02 -6.10066870e+02);
  753. L_SXZT = L_SXZT ET (PROG
  754. -5.75710520e+02 -5.42577660e+02 -5.09451500e+02
  755. -4.76325960e+02 -4.43195980e+02 -4.10057300e+02
  756. -3.76906030e+02 -3.43737790e+02 -3.10544340e+02
  757. -2.78450010e+02 -2.45976840e+02 -2.14661020e+02
  758. -1.85322180e+02 -1.55986210e+02 -1.26649390e+02);
  759. L_SXZT = L_SXZT ET (PROG
  760. -9.73085700e+01 -6.79610740e+01 -3.86045080e+01
  761. -9.23644540e+00 2.01463280e+01 4.95512180e+01
  762. 7.80079940e+01 1.06509910e+02 1.32712670e+02
  763. 1.52234580e+02 1.71753460e+02 1.91271570e+02
  764. 2.10790740e+02 2.30312490e+02 2.49838010e+02);
  765. L_SXZT = L_SXZT ET (PROG
  766. 2.69368290e+02 2.88904130e+02 3.08446190e+02
  767. 3.27995070e+02 3.47551360e+02 3.67115820e+02
  768. 3.86689760e+02 4.06276810e+02 4.25144500e+02
  769. 4.44008510e+02 4.62911230e+02 4.81940030e+02
  770. 4.98703860e+02 5.14683880e+02 5.30666070e+02);
  771. L_SXZT = L_SXZT ET (PROG
  772. 5.46651380e+02 5.62640560e+02 5.78634160e+02
  773. 5.94632610e+02 6.10636200e+02 6.26645160e+02
  774. 6.42659610e+02 6.58679640e+02 6.74705280e+02
  775. 6.90736530e+02 7.06773370e+02 7.22815740e+02
  776. 7.38863580e+02 7.54916830e+02 7.70975400e+02);
  777. L_SXZT = L_SXZT ET (PROG
  778. 7.87039220e+02 8.03108190e+02 8.19182230e+02
  779. 8.35261260e+02 8.51345210e+02 8.67433990e+02
  780. 8.83527530e+02 8.99625760e+02 9.15728630e+02
  781. 9.31836060e+02 9.47948020e+02 9.64064430e+02
  782. 9.80185270e+02 9.96310490e+02 1.01244010e+03);
  783. L_SXZT = L_SXZT ET (PROG
  784. 1.02857390e+03 1.04471210e+03 1.06085450e+03
  785. 1.07700120e+03 1.09315200e+03 1.10930710e+03
  786. 1.12546640e+03 1.14162990e+03 1.15779760e+03
  787. 1.17396940e+03 1.19014540e+03 1.20632560e+03
  788. 1.22251000e+03 1.23869860e+03 1.25489140e+03);
  789. L_SXZT = L_SXZT ET (PROG
  790. 1.27108840e+03 1.28728960e+03 1.30349500e+03
  791. 1.31970470e+03 1.33591860e+03 1.35213690e+03
  792. 1.36835940e+03 1.38458620e+03 1.40081730e+03
  793. 1.41705280e+03 1.43329260e+03 1.44953690e+03
  794. 1.46578550e+03 1.48203860e+03 1.49829610e+03);
  795. L_SXZT = L_SXZT ET (PROG
  796. 1.51455800e+03 1.53082450e+03 1.54709550e+03
  797. 1.56337100e+03 1.57965110e+03 1.59593580e+03
  798. 1.61222510e+03 1.62851900e+03 1.64481760e+03
  799. 1.66112090e+03 1.67742880e+03 1.69374150e+03
  800. 1.71005900e+03 1.72638120e+03 1.74270830e+03);
  801. L_SXZT = L_SXZT ET (PROG
  802. 1.75904020e+03 1.77537690e+03 1.79171850e+03
  803. 1.80806510e+03 1.82441650e+03 1.84077290e+03
  804. 1.85713430e+03 1.87350070e+03 1.88987210e+03
  805. 1.90624860e+03 1.92263020e+03 1.93901680e+03
  806. 1.95540860e+03 1.97180560e+03 1.98820770e+03);
  807. L_SXZT = L_SXZT ET (PROG
  808. 2.00461500e+03 2.02102750e+03 2.03744530e+03
  809. 2.05386830e+03 2.07029670e+03 2.08673030e+03
  810. 2.10316930e+03 2.11961370e+03 2.13606340e+03
  811. 2.15251860e+03 2.16897910e+03 2.18544520e+03
  812. 2.20191670e+03 2.21839370e+03 2.23487620e+03);
  813. L_SXZT = L_SXZT ET (PROG
  814. 2.25136430e+03 2.26785790e+03 2.28435710e+03
  815. 2.30086190e+03 2.31737240e+03 2.33388850e+03
  816. 2.35041020e+03 2.36693770e+03 2.38347080e+03
  817. 2.40000970e+03 2.41655430e+03);
  818. *
  819. *********************************************************************
  820. * Test
  821. *********************************************************************
  822. *
  823. EV_TEST = EVOL (VERT) MANU 'Exz []' L_EXZ 'Sxz [Pa]' L_SXZT;
  824. SI (VAR_GRAP EGA 1);
  825. DESS (EV_1 ET EV_TEST);
  826. FINSI;
  827. VAL_TEST = MAXI (ABS(L_SXZ - L_SXZT));
  828. VAL_TEST = VAL_TEST/(MAXI (ABS(L_SXZ)));
  829. *
  830. SI (VAL_TEST > 1.E-5);
  831. ERRE 5;
  832. FINSI;
  833. *
  834. FIN;
  835.  
  836.  
  837.  
  838.  
  839.  
  840.  
  841.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales