Télécharger mrcshell.procedur

Retour à la liste

Numérotation des lignes :

  1. * MRCSHELL PROCEDUR AF221230 16/11/03 21:15:12 9139
  2. DEBPROC MRCSHELL TYP_CAL1*'MOT '
  3. TAB_IN*'TABLE' TOL_1*'FLOTTANT'
  4. L_ELE1/'LISTENTI';
  5. *
  6. * --------------------------------------------------------
  7. * Procedure MRCSHELL
  8. * Computation of security coef Lambda
  9. *
  10. * Develeppers:
  11. * Alberto FRAU
  12. * DEN/DANS/DM2S/SEMT/EMSI
  13. * Nicolas ILE
  14. * DEN/DANS/DM2S/SEMT/EMSI
  15. *
  16. * INPUT:
  17. * TYP_CAL1: Calculation option
  18. * 'RECTANGLE'
  19. * 'ELLIPSOIDE'
  20. * TAB_IN.'MAILLAGE': Maillage du voile [MAILLAGE]
  21. * TAB_IN.'EFFORT_SEISME': MCHAML de la matrice
  22. * representant l'enveloppe des
  23. * efforts sismiques (voir SISSIB)
  24. * TAB_IN.'EFFORT_STATIQUE': MCHAML des efforts statiques qui
  25. * agissent sur le voile
  26. * (issu EFFMARTI)
  27. * TAB_IN.'CARACTERISTIQUES_EXTERNE':MCHAML contenant les
  28. * caracteristiques de la couche
  29. * externe selon le modele de MARTI
  30. * (voir EFFMARTI):
  31. * 'RHO1' taux d'acier direction 1
  32. * 'RHO2' taux d'acier direction 2
  33. * 'ENRO' Enrobage
  34. * TAB_IN.'CARACTERISTIQUES_INTERNE':MCHAML contenant les
  35. * caracteristiques de la couche
  36. * interne selon le modele de MARTI
  37. * (voir EFFMARTI):
  38. * 'RHO1' taux d'acier direction 1
  39. * 'RHO2' taux d'acier direction 2
  40. * 'ENRO' Enrobage
  41. * TAB_IN.'CARACTERISTIQUES_CORE': MCHAML contenant les
  42. * caracteristiques et le coeficient
  43. * pour la couche intermediaire
  44. * selon le modele de MARTI
  45. * 'H' Epaisseur Totale
  46. * 'RHOT' taux d'acier transversale
  47. * 'COTH' Facteur de cisaillement
  48. * TAB_IN.'FC_BET': Resistance caracteristique beton
  49. * [FLOTTANT]
  50. * TAB_IN.'GAM_C': Coef gammac 1.5/1.2 EC2/EC8
  51. * [FLOTTANT]
  52. * TAB_IN.'ALP_C': Coef alpa 1.0 EC2
  53. * [FLOTTANT]
  54. * TAB_IN.'NU_C': Coef nu 0.6*(1-250/fck) EC2
  55. * [FLOTTANT]
  56. * TAB_IN.'FS_ACI': resistance caracteristique de l'acier
  57. * [FLOTTANT]
  58. * TAB_IN.'GAM_S': Coef gammas 1.15 EC2
  59. * [FLOTTANT]
  60. * TOL_1: Tollerance [FLOTTANT]
  61. * L_ELE1: Liste des elements sur lesquels on veut sortir les
  62. * graphiques des surfaces limites et des enveloppes
  63. * (sans ou avec covariance - RECTANGLE ou ELLISPOIDE)
  64. * pour les couches externes, interne et
  65. * intermediaire (pas necessaire) [LISTREEL]
  66. * OUTPUT:
  67. * TT2: Table Contenant:
  68. * .'CH_LAMBDA_E': MCHAML les valeurs des marges pour la couche
  69. * externe (composent LAME)
  70. * .'CH_LAMBDA_I': MCHAML des valeurs des marges pour la couche
  71. * interne (composent LAMI)
  72. * .'CH_LAMBDA_C': MCHAML des valeurs des marges pour la couche
  73. * intermediaire (composent LAMC)
  74. * .'CARTE_E': Maillage avec deux colorations pour indiquer les
  75. * element avec une marge superieure à 1 ou
  76. * inferieure (couche externe):
  77. * ROUGE elements -> Lambda < 1.0
  78. * VERT elements -> Lambda > 1.0
  79. * .'CARTE_I': Maillage avec deux colorations pour indiquer les
  80. * element avec une marge superieure à 1 ou
  81. * inferieure (couche interne):
  82. * ROUGE elements -> Lambda < 1.0
  83. * VERT elements -> Lambda > 1.0
  84. * .'CARTE_C': Maillage avec deux colorations pour indiquer les
  85. * element avec une marge superieure à 1 ou
  86. * inferieure (couche intermediaire):
  87. * ROUGE elements -> Lambda < 1.0
  88. * VERT elements -> Lambda > 1.0
  89. * .'GRAPHIQUES': Sous table pour les outils de graphique:
  90. * .I. Ieme element
  91. * .'LIMITE_E': Surface limite de la couche externe
  92. * [MAILLAGE]
  93. * .'LIMITE_E': Surface limite de la couche interne
  94. * [MAILLAGE]
  95. * .'LIMITE_C': Surface limite de la couche
  96. * intermediaire [MAILLAGE]
  97. * .'RECTANGLE_E': enveloppe couche externe sans prise
  98. * en compte des covariances
  99. * (methode RECTANGLE) [MAILLAGE]
  100. * .'RECTANGLE_I': enveloppe couche interne sans prise
  101. * en compte des covariances
  102. * (methode RECTANGLE) [MAILLAGE]
  103. * .'RECTANGLE_I': enveloppe couche intermediaire
  104. * sans prise en compte des covariances
  105. * (methode RECTANGLE) [MAILLAGE]
  106. * .'ELLIPSOIDE_E': enveloppe couche externe avec prise
  107. * en compte des covariances
  108. * (methode ELLIPSOIDE) [MAILLAGE]
  109. * .'ELLIPSOIDE_I': enveloppe couche interne avec prise
  110. * en compte des covariances
  111. * (methode ELLIPSOIDE) [MAILLAGE]
  112. * .'ELLIPSOIDE_I': enveloppe couche intermediaire
  113. * avec prise en compte des covariances
  114. * (methode ELLIPSOIDE) [MAILLAGE]
  115. * .'RECTANGLE_E_AUG': enveloppe couche externe sans prise
  116. * en compte des covariances amplifie par
  117. * la marge LAME
  118. * (methode RECTANGLE) [MAILLAGE]
  119. * .'RECTANGLE_I_AUG': enveloppe couche interne sans prise
  120. * en compte des covariances amplifie par
  121. * la marge LAMI
  122. * (methode RECTANGLE) [MAILLAGE]
  123. * .'RECTANGLE_C_AUG': enveloppe couche intermediaire sans prise
  124. * en compte des covariances amplifie par
  125. * la marge LAMC
  126. * (methode RECTANGLE) [MAILLAGE]
  127. * .'ELLIPSOIDE_E_AUG': enveloppe couche externe avec prise
  128. * en compte des covariances amplifie par
  129. * la marge LAME
  130. * (methode ELLIPSOIDE) [MAILLAGE]
  131. * .'ELLIPSOIDE_I_AUG': enveloppe couche interne avec prise
  132. * en compte des covariances amplifie par
  133. * la marge LAMI
  134. * (methode ELLIPSOIDE) [MAILLAGE]
  135. * .'ELLIPSOIDE_C_AUG': enveloppe couche intermediaire avec
  136. * prise en compte des covariances
  137. * amplifie par la marge LAMC
  138. * (methode ELLIPSOIDE) [MAILLAGE]
  139. *
  140. * --------------------------------------------------------
  141. *
  142. MAI_V = TAB_IN.'MAILLAGE';
  143. CH_SEI1 = TAB_IN.'EFFORT_SEISME';
  144. CH_STA1 = TAB_IN.'EFFORT_STATIQUE';
  145. CH_CARE = TAB_IN.'CARACTERISTIQUES_EXTERNE';
  146. CH_CARI = TAB_IN.'CARACTERISTIQUES_INTERNE';
  147. CH_CARC = TAB_IN.'CARACTERISTIQUES_CORE';
  148. *D_VOI1 = TAB_IN.'D';
  149. FCBET1 = TAB_IN.'FC_BET';
  150. GAMC = TAB_IN.'GAM_C';
  151. ALPC1 = TAB_IN.'ALP_C';
  152. NUC1 = TAB_IN.'NU_C';
  153. FSACI1 = TAB_IN.'FS_ACI';
  154. GAMS1 = TAB_IN.'GAM_S';
  155. *
  156. 'MESS' ;
  157. 'MESS' '**************************************************************';
  158. 'MESS' ' Calculs des coef de securite ';
  159. 'MESS' '**************************************************************';
  160. 'MESS';
  161. *
  162. * controle of input data
  163. 'SI' (NEG (TYP_CAL1) 'ELLIPSOIDE');
  164. 'SI' (NEG (TYP_CAL1) 'RECTANGLE');
  165. 'MESS' ' ERROR....Type de calcul erronee!!!';
  166. 'QUIT' LAM_RCM;
  167. 'FINSI';
  168. 'FINSI';
  169. 'SI' (('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 1) 'C11E') 'ET'
  170. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 2) 'C22E') 'ET'
  171. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 3) 'C33E') 'ET'
  172. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 4) 'C12E') 'ET'
  173. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 5) 'C13E') 'ET'
  174. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 6) 'C23E'));
  175. 'SINON';
  176. 'SI' (('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 7) 'C11I') 'ET'
  177. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 8) 'C22I') 'ET'
  178. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 9) 'C33I') 'ET'
  179. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 10) 'C12I') 'ET'
  180. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 11) 'C13I') 'ET'
  181. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 12) 'C23I'));
  182. 'SINON';
  183. 'SI' (('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 13) 'CV11') 'ET'
  184. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 14) 'CV22') 'ET'
  185. ('NEG' ('EXTR' ('EXTR' CH_SEI1 'COMP') 15) 'CV12'));
  186. 'MESS' ' ERROR....MCHAML seisme';
  187. 'MESS' ' composantes erronee!!!';
  188. 'QUIT' LAM_RCM;
  189. 'FINSI';
  190. 'FINSI';
  191. 'FINSI';
  192. 'SI' (('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 1) 'N11E') 'ET'
  193. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 2) 'N22E') 'ET'
  194. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 3) 'N12E') 'ET'
  195. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 4) 'N11I') 'ET'
  196. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 5) 'N22I') 'ET'
  197. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 6) 'N12I'));
  198. 'MESS' ' ERROR....MCHAML statique';
  199. 'MESS' ' composantes erronee!!!';
  200. 'FINSI';
  201. 'SI' (('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 1) 'N11E') 'ET'
  202. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 2) 'N22E') 'ET'
  203. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 3) 'N12E') 'ET'
  204. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 4) 'N11I') 'ET'
  205. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 5) 'N22I') 'ET'
  206. ('NEG' ('EXTR' ('EXTR' CH_STA1 'COMP') 6) 'N12I'));
  207. 'MESS' ' ERROR....MCHAML seisme';
  208. 'MESS' ' composantes erronee!!!';
  209. 'FINSI';
  210. 'SI' (('NEG' ('EXTR' ('EXTR' CH_CARE 'COMP') 1) 'RHO1') 'ET'
  211. ('NEG' ('EXTR' ('EXTR' CH_CARE 'COMP') 2) 'RHO2') 'ET'
  212. ('NEG' ('EXTR' ('EXTR' CH_CARE 'COMP') 3) 'ENRO'));
  213. 'MESS' ' ERROR....MCHAML caracteristiques couche externe';
  214. 'MESS' ' composantes erronee!!!';
  215. 'FINSI';
  216. 'SI' (('NEG' ('EXTR' ('EXTR' CH_CARI 'COMP') 1) 'RHO1') 'ET'
  217. ('NEG' ('EXTR' ('EXTR' CH_CARI 'COMP') 2) 'RHO2') 'ET'
  218. ('NEG' ('EXTR' ('EXTR' CH_CARI 'COMP') 3) 'ENRO'));
  219. 'MESS' ' ERROR....MCHAML caracteristiques couche interne';
  220. 'MESS' ' composantes erronee!!!';
  221. 'FINSI';
  222. *
  223. * creation of the table output
  224. TT2 = 'TABLE';
  225. *
  226. * creation of the MCHAML field for the output
  227. CH_LAME = 'VIDE' MCHAML;
  228. CH_LAMI = 'VIDE' MCHAML;
  229. CH_LAMC = 'VIDE' MCHAML;
  230. MAI_REE = 'VIDE' MAILLAGE;
  231. MAI_REI = 'VIDE' MAILLAGE;
  232. MAI_REC = 'VIDE' MAILLAGE;
  233. MODTOT = 'VIDE' MMODEL;
  234. *
  235. * number of element shear wall
  236. N_ELE1 = 'NBEL' MAI_V;
  237. *
  238. * Ellipsoide methode
  239. 'SI' ('EGA' TYP_CAL1 'ELLIPSOIDE');
  240. *------------------ Calcul par l'ellipsoide ---------------------------*
  241. 'MESS' ' ------------------------------------------------------------';
  242. 'MESS' ' Methode Ellipsoide';
  243. 'MESS' ' ------------------------------------------------------------';
  244. 'MESS';
  245. *
  246. * loop for the layer (external or intenal)
  247. I_IE1 = 1;
  248. 'REPE' IND_IE1 (3);
  249. *
  250. * I_IE1 = 1 External layer
  251. 'SI' (I_IE1 'EGA' 1);
  252. 'MESS';
  253. 'MESS' ' Membrane Exterieure';
  254. 'FINSI';
  255. *
  256. * I_IE1 = 2 Internal layer
  257. 'SI' (I_IE1 'EGA' 2);
  258. 'MESS';
  259. 'MESS' ' Membrane Interieure';
  260. 'FINSI';
  261. *
  262. 'SI' (I_IE1 'EGA' 3);
  263. 'MESS';
  264. 'MESS' ' Core';
  265. 'FINSI';
  266. *
  267. * loop for the element
  268. I_ELE1 = 1;
  269. 'REPE' IND_ELE1 (N_ELE1);
  270. * element I
  271. ELE_VOI1 = MAI_V ELEM I_ELE1;
  272. * extracting the forces and characteristics for the external layer
  273. 'SI' (I_IE1 'EGA' 1);
  274. RHOXS = 'EXTR' ('REDU' CH_CARE ELE_VOI1) 'RHO1' 1 1 1;
  275. RHOYS = 'EXTR' ('REDU' CH_CARE ELE_VOI1) 'RHO2' 1 1 1;
  276. ENR1 = 'EXTR' ('REDU' CH_CARE ELE_VOI1) 'ENRO' 1 1 1;
  277. *
  278. N11 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C11E' 1 1 1)**(0.5);
  279. N11N22 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C12E' 1 1 1);
  280. N11N12 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C13E' 1 1 1);
  281. N22 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C22E' 1 1 1)**(0.5);
  282. N22N12 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C23E' 1 1 1);
  283. N12 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C33E' 1 1 1)**(0.5);
  284. N110 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'N11E' 1 1 1);
  285. N220 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'N22E' 1 1 1);
  286. N120 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'N12E' 1 1 1);
  287. *
  288. * seismic stress
  289. N11 = ((N11)/((2.0)*(ENR1)));
  290. N22 = ((N22)/((2.0)*(ENR1)));
  291. N12 = ((N12)/((2.0)*(ENR1)));
  292. N11N22 = ((N11N22)/(((2.0)*(ENR1))**(2.0)));
  293. N11N12 = ((N11N12)/(((2.0)*(ENR1))**(2.0)));
  294. N22N12 = ((N22N12)/(((2.0)*(ENR1))**(2.0)));
  295. *
  296. * static stress
  297. N110 = ((N110)/((2.0)*(ENR1)));
  298. N220 = ((N220)/((2.0)*(ENR1)));
  299. N120 = ((N120)/((2.0)*(ENR1)));
  300. *
  301. 'FINSI';
  302. * extracting the force and characteristics for the internal layer
  303. 'SI' (I_IE1 'EGA' 2);
  304. RHOXS = 'EXTR' ('REDU' CH_CARI ELE_VOI1) 'RHO1' 1 1 1;
  305. RHOYS = 'EXTR' ('REDU' CH_CARI ELE_VOI1) 'RHO2' 1 1 1;
  306. ENR1 = 'EXTR' ('REDU' CH_CARI ELE_VOI1) 'ENRO' 1 1 1;
  307. *
  308. N11 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C11I' 1 1 1)**(0.5);
  309. N11N22 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C12I' 1 1 1);
  310. N11N12 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C13I' 1 1 1);
  311. N22 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C22I' 1 1 1)**(0.5);
  312. N22N12 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C23I' 1 1 1);
  313. N12 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'C33I' 1 1 1)**(0.5);
  314. N110 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'N11I' 1 1 1);
  315. N220 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'N22I' 1 1 1);
  316. N120 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'N12I' 1 1 1);
  317. *
  318. * seismic stress
  319. N11 = ((N11)/((2.0)*(ENR1)));
  320. N22 = ((N22)/((2.0)*(ENR1)));
  321. N12 = ((N12)/((2.0)*(ENR1)));
  322. N11N22 = ((N11N22)/(((2.0)*(ENR1))**(2.0)));
  323. N11N12 = ((N11N12)/(((2.0)*(ENR1))**(2.0)));
  324. N22N12 = ((N22N12)/(((2.0)*(ENR1))**(2.0)));
  325. *
  326. * static stress
  327. N110 = ((N110)/((2.0)*(ENR1)));
  328. N220 = ((N220)/((2.0)*(ENR1)));
  329. N120 = ((N120)/((2.0)*(ENR1)));
  330. *
  331. 'FINSI';
  332. *
  333. * extracting the force and characteristics for the internal layer
  334. 'SI' (I_IE1 'EGA' 3);
  335. RHOT = 'EXTR' ('REDU' CH_CARC ELE_VOI1) 'RHOT' 1 1 1;
  336. ENR1 = 'EXTR' ('REDU' CH_CARE ELE_VOI1) 'ENRO' 1 1 1;
  337. ENR2 = 'EXTR' ('REDU' CH_CARI ELE_VOI1) 'ENRO' 1 1 1;
  338. H_11 = 'EXTR' ('REDU' CH_CARC ELE_VOI1) 'H' 1 1 1;
  339. D_11 = H_11 - ((2.0)*(ENR1)) - ((2.0)*(ENR2));
  340. CO_TH1 = 'EXTR' ('REDU' CH_CARC ELE_VOI1) 'COTH' 1 1 1;
  341. *
  342. V1 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'CV11' 1 1 1)**(0.5);
  343. V1V2 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'CV12' 1 1 1);
  344. V2 = ('EXTR' ('REDU' CH_SEI1 ELE_VOI1) 'CV22' 1 1 1)**(0.5);
  345. *
  346. V10 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'V1T' 1 1 1);
  347. V20 = ('EXTR' ('REDU' CH_STA1 ELE_VOI1) 'V2T' 1 1 1);
  348. *
  349. * seismic stress
  350. T1 = ((V1)/((D_11)));
  351. T2 = ((V2)/((D_11)));
  352. T1T2 = ((V1V2)/(((D_11))**(2.0)));
  353. *
  354. * static stress
  355. T10 = ((V10)/((D_11)));
  356. T20 = ((V20)/((D_11)));
  357. *
  358. 'FINSI';
  359. *
  360. * Verification Core -----------------------------------------
  361. 'SI' (I_IE1 'EGA' 3);
  362. * --------------------------------------------------------------
  363. TTL1 = TABLE;
  364. TTL1.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  365. TTL1.'T1' = T10;
  366. TTL1.'T2' = T20;
  367. TTL1.'FCK' = FCBET1;
  368. TTL1.'FSK_T' = FSACI1;
  369. TTL1.'RHO_T' = RHOT;
  370. TTL1.'GAMMA_C' = GAMC;
  371. TTL1.'ALPA_C' = ALPC1;
  372. TTL1.'GAMMA_S' = GAMS1;
  373. TTL1.'NU_C' = NUC1;
  374. TTL1.'COT_THETA' = CO_TH1;
  375. T_ST1 = G_ULTI2D TTL1;
  376. *
  377. 'SI' (T_ST1 '&lt;EG' 0.0);
  378. LAM_FIN = 0.0;
  379. TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda='
  380. LAM_FIN '.....Sig0 externe';
  381. 'MESS' TIT1;
  382. MOD1 = 'MODE' ELE_VOI1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';
  383. MODTOT = MODTOT 'ET' MOD1;
  384. CH_LAMC = CH_LAMC 'ET' ('MANU' 'CHML' ELE_VOI1 'LAMC' LAM_FIN);
  385. MAI_REC = MAI_REC 'ET' (ELE_VOI1 'COUL' 'ROUG');
  386. OUBL LAM_FIN;
  387. 'SINON';
  388.  
  389. *
  390. * Start of the iteration
  391. P1 = 0. 0. 0.;
  392. P2 = 1. 0. 0.;
  393. P3 = 2. 0. 0.;
  394. *
  395. MAI1 = P1;
  396. * Matrix X
  397. MATX1 = VIDE 'RIGIDITE'/'RIGIDITE';
  398. MATX1 = MATX1 'ET' ('MANU' 'RIGIDITE' MAI1 (MOTS 'UX' 'UY')
  399. ('PROG' (T1*T1) T1T2
  400. T1T2 (T2*T2)));
  401. *
  402. * Static state stress
  403. EF_ST0 = 'MANU' 'CHPO' MAI1 2 'FX' T10 'FY' T20;
  404. EF_ST0 = 'CHAN' 'ATTRIBUT' EF_ST0 'NATURE' 'DISCRET';
  405. *
  406. * Initial Vecteur
  407. ALP_0 = 'MANU' 'CHPO' MAI1 2 'UX' 2.0 'UY' 2.0;
  408. XALP_0 = (MATX1*ALP_0)/(('XTMX' MATX1 ALP_0)**(0.5));
  409. XALP_0 = ('CHAN' ATTRIBUT XALP_0 NATURE DISCRET) + EF_ST0;
  410. SB0 = ('XTX' (XALP_0 - EF_ST0))**(0.5);
  411. RB0 = SB0;
  412. LAMB0 = RB0/SB0;
  413. *
  414. * incremental step
  415. H_VAL1 = (SB0)/(100000.0);
  416. H_VAL2 = (SB0)/(100000.0);
  417. *
  418. * direction used for the ellispoide methode
  419. TB_DIR = 'TABLE';
  420. TB_DIR. 1 = 'MANU' 'CHPO' MAI1 2 'UX' 1.0 'UY' 0.0 ;
  421. TB_DIR. 2 = 'MANU' 'CHPO' MAI1 2 'UX' 0.0 'UY' 1.0 ;
  422. TB_DIR. 3 = 'MANU' 'CHPO' MAI1 2 'UX' -1.0 'UY' 0.0 ;
  423. TB_DIR. 4 = 'MANU' 'CHPO' MAI1 2 'UX' 0.0 'UY' -1.0;
  424. *
  425. * loop for each direction
  426. I_K1 = 1;
  427. 'REPE' IND3 (4);
  428. *
  429. TAB_ITE = 'TABLE';
  430. TAB_ITE. 0 = 'TABLE';
  431. TAB_ITE. 0 . 'ALPA' = ALP_0;
  432. TAB_ITE. 0 . 'LAMBDA' = LAMB0;
  433. TAB_ITE. 0 . 'RB' = RB0;
  434. *
  435. ALP_K = (TB_DIR. I_K1);
  436. I_M1 = 1;
  437. * first loop in order to have n_elllips=n_limite
  438. 'REPE' IND2 (200);
  439. * sauving Lambda for j iteraction
  440. TAB_ITE. I_M1 = 'TABLE';
  441. TAB_ITE. I_M1 . 'ALPA' = ALP_K;
  442. * determination of the stress state for the chosen direction
  443. XALP_J = (MATX1*ALP_K)/(('XTMX' MATX1 ALP_K)**(0.5));
  444. XALP_J = ('CHAN' ATTRIBUT XALP_J NATURE DISCRET) + EF_ST0;
  445. * nomalisez vector of hte rayon of the ellipsoide for the
  446. * chosen direction
  447. SBJ = ('XTX' (XALP_J - EF_ST0))**(0.5);
  448. BJ = (XALP_J - EF_ST0)/(SBJ);
  449. * rayon of the ellipsoide for the chosen direction
  450. RBVAL1 = (TAB_ITE.(I_M1 - 1). 'RB');
  451. * second loop in order to have f=0
  452. I_I1 = 1;
  453. REPE IND1 (200);
  454. * strees state x0, x0+h, x0-h
  455. VV1 = (RBVAL1 + H_VAL1)*(BJ);
  456. VV2 = (RBVAL1 - H_VAL1)*(BJ);
  457. VV3 = (RBVAL1)*(BJ);
  458. VV1 = ('CHAN' 'ATTRIBUT' VV1 'NATURE' 'DISCRET') + EF_ST0;
  459. VV2 = ('CHAN' 'ATTRIBUT' VV2 'NATURE' 'DISCRET') + EF_ST0;
  460. VV3 = ('CHAN' 'ATTRIBUT' VV3 'NATURE' 'DISCRET') + EF_ST0;
  461. *
  462. T1V1 = 'EXTR' VV1 'FX' MAI1;
  463. T2V1 = 'EXTR' VV1 'FY' MAI1;
  464.  
  465. T1V2 = 'EXTR' VV2 'FX' MAI1;
  466. T2V2 = 'EXTR' VV2 'FY' MAI1;
  467.  
  468. T1V3 = 'EXTR' VV3 'FX' MAI1;
  469. T2V3 = 'EXTR' VV3 'FY' MAI1;
  470.  
  471. 'OUBL' VV1;
  472. 'OUBL' VV2;
  473. 'OUBL' VV3;
  474. * calculation of f and D terms
  475. TTL1 = TABLE;
  476. TTL1.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  477. TTL1.'T1' = T1V3;
  478. TTL1.'T2' = T2V3;
  479. TTL1.'FCK' = FCBET1;
  480. TTL1.'FSK_T' = FSACI1;
  481. TTL1.'RHO_T' = RHOT;
  482. TTL1.'GAMMA_C' = GAMC;
  483. TTL1.'ALPA_C' = ALPC1;
  484. TTL1.'GAMMA_S' = GAMS1;
  485. TTL1.'NU_C' = NUC1;
  486. TTL1.'COT_THETA' = CO_TH1;
  487. *
  488. TTL2 = TABLE;
  489. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  490. TTL2.'T1' = T1V1;
  491. TTL2.'T2' = T2V1;
  492. TTL2.'FCK' = FCBET1;
  493. TTL2.'FSK_T' = FSACI1;
  494. TTL2.'RHO_T' = RHOT;
  495. TTL2.'GAMMA_C' = GAMC;
  496. TTL2.'ALPA_C' = ALPC1;
  497. TTL2.'GAMMA_S' = GAMS1;
  498. TTL2.'NU_C' = NUC1;
  499. TTL2.'COT_THETA' = CO_TH1;
  500. *
  501. TTL3 = TABLE;
  502. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  503. TTL3.'T1' = T1V2;
  504. TTL3.'T2' = T2V2;
  505. TTL3.'FCK' = FCBET1;
  506. TTL3.'FSK_T' = FSACI1;
  507. TTL3.'RHO_T' = RHOT;
  508. TTL3.'GAMMA_C' = GAMC;
  509. TTL3.'ALPA_C' = ALPC1;
  510. TTL3.'GAMMA_S' = GAMS1;
  511. TTL3.'NU_C' = NUC1;
  512. TTL3.'COT_THETA' = CO_TH1;
  513. *
  514. F_VAL1 = G_ULTI2D TTL1;
  515. D_VAL1 = ((G_ULTI2D TTL2) -
  516. (G_ULTI2D TTL3))/((2.0)*(H_VAL1));
  517. * controle
  518. 'SI' (('ABS'(F_VAL1)) '<' TOL_1);
  519. * convergence!!!
  520. RBJ = RBVAL1;
  521. 'QUIT' IND1;
  522. 'SINON';
  523. 'SI' (I_I1 'EGA' 200);
  524. * end of loop, we take the last value!!!
  525. RBJ = RBVAL1;
  526. 'QUIT' IND1;
  527. 'FINSI';
  528. * Check beetwen 2 values of R (radius)!!!
  529. 'SI' (I_I1 > 2);
  530. 'SI' ((ABS((RBVAL1 - RBVAL2)/(RBVAL1))) < ((TOL_1)*(1.e-10)));
  531. RBJ = RBVAL1;
  532. 'QUIT' IND1;
  533. 'FINSI';
  534. 'FINSI';
  535. * case when D is equal to zero
  536. 'SI' (('ABS'(D_VAL1)) < 1.E-16);
  537. RBVAL2 = RBVAL1;
  538. RBVAL1 = ((1.1)*(RBVAL1));
  539. 'SINON';
  540. * otherwise
  541. RBVAL2 = RBVAL1;
  542. RBVAL1 = RBVAL1 - ((F_VAL1)/(D_VAL1));
  543. 'FINSI';
  544. 'FINSI';
  545. *
  546. I_I1 = I_I1 + 1;
  547. 'FIN' IND1;
  548. * computation of lambda
  549. LAMBJ = RBJ/SBJ;
  550. TAB_ITE. I_M1 . 'LAMBDA' = LAMBJ;
  551. TAB_ITE. I_M1 . 'RB' = RBJ;
  552. XCOUB1 = (((TAB_ITE. I_M1 .'LAMBDA')*(TAB_ITE. I_M1 .'LAMBDA'))*
  553. (MATX1));
  554. * amplified stress state using the lambda coef
  555. XALP_L = (XCOUB1*ALP_K)/(('XTMX' XCOUB1 ALP_K)**(0.5));
  556. XALP_L = ('CHAN' ATTRIBUT XALP_L NATURE DISCRET) + EF_ST0;
  557. * Dh*ex ; Dh*ey ; Dh*ez
  558. DE_HX = 'MANU' 'CHPO' MAI1 2 'FX' H_VAL2 'FY' 0.0 ;
  559. DE_HY = 'MANU' 'CHPO' MAI1 2 'FX' 0.0 'FY' H_VAL2 ;
  560. * x0*lam + Dh*ex ; x0*lam + Dh*ey ; x0*lam + Dh*ez
  561. VV_HXP = XALP_L + DE_HX;
  562. VV_HYP = XALP_L + DE_HY;
  563. * x0*lam + Dh*ex ; x0*lam + Dh*ey ; x0*lam + Dh*ez
  564. VV_HXN = XALP_L - DE_HX;
  565. VV_HYN = XALP_L - DE_HY;
  566. * x0*lam
  567. VV_0 = XALP_L;
  568. *
  569. VV_HXP1 = 'EXTR' VV_HXP 'FX' MAI1;
  570. VV_HXP2 = 'EXTR' VV_HXP 'FY' MAI1;
  571. VV_HYP1 = 'EXTR' VV_HYP 'FX' MAI1;
  572. VV_HYP2 = 'EXTR' VV_HYP 'FY' MAI1;
  573. *
  574. VV_HXN1 = 'EXTR' VV_HXN 'FX' MAI1;
  575. VV_HXN2 = 'EXTR' VV_HXN 'FY' MAI1;
  576. VV_HYN1 = 'EXTR' VV_HYN 'FX' MAI1;
  577. VV_HYN2 = 'EXTR' VV_HYN 'FY' MAI1;
  578. *
  579. VV_01 = 'EXTR' VV_HXP 'FX' MAI1;
  580. VV_02 = 'EXTR' VV_HXP 'FY' MAI1;
  581. *
  582. TTL2 = TABLE;
  583. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  584. TTL2.'T1' = VV_HXP1;
  585. TTL2.'T2' = VV_HXP2;
  586. TTL2.'FCK' = FCBET1;
  587. TTL2.'FSK_T' = FSACI1;
  588. TTL2.'RHO_T' = RHOT;
  589. TTL2.'GAMMA_C' = GAMC;
  590. TTL2.'ALPA_C' = ALPC1;
  591. TTL2.'GAMMA_S' = GAMS1;
  592. TTL2.'NU_C' = NUC1;
  593. TTL2.'COT_THETA' = CO_TH1;
  594. TTL3 = TABLE;
  595. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  596. TTL3.'T1' = VV_HXN1;
  597. TTL3.'T2' = VV_HXN2;
  598. TTL3.'FCK' = FCBET1;
  599. TTL3.'FSK_T' = FSACI1;
  600. TTL3.'RHO_T' = RHOT;
  601. TTL3.'GAMMA_C' = GAMC;
  602. TTL3.'ALPA_C' = ALPC1;
  603. TTL3.'GAMMA_S' = GAMS1;
  604. TTL3.'NU_C' = NUC1;
  605. TTL3.'COT_THETA' = CO_TH1;
  606. VAL_X1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  607. *
  608. TTL2 = TABLE;
  609. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  610. TTL2.'T1' = VV_HYP1;
  611. TTL2.'T2' = VV_HYP2;
  612. TTL2.'FCK' = FCBET1;
  613. TTL2.'FSK_T' = FSACI1;
  614. TTL2.'RHO_T' = RHOT;
  615. TTL2.'GAMMA_C' = GAMC;
  616. TTL2.'ALPA_C' = ALPC1;
  617. TTL2.'GAMMA_S' = GAMS1;
  618. TTL2.'NU_C' = NUC1;
  619. TTL2.'COT_THETA' = CO_TH1;
  620. TTL3 = TABLE;
  621. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_COR';
  622. TTL3.'T1' = VV_HYN1;
  623. TTL3.'T2' = VV_HYN2;
  624. TTL3.'FCK' = FCBET1;
  625. TTL3.'FSK_T' = FSACI1;
  626. TTL3.'RHO_T' = RHOT;
  627. TTL3.'GAMMA_C' = GAMC;
  628. TTL3.'ALPA_C' = ALPC1;
  629. TTL3.'GAMMA_S' = GAMS1;
  630. TTL3.'NU_C' = NUC1;
  631. TTL3.'COT_THETA' = CO_TH1;
  632. VAL_Y1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  633. *
  634. * costruction n_limit
  635. DELTG1 = 'MANU' 'CHPO' MAI1 'UX' ((-1.0)*(VAL_X1))
  636. 'UY' ((-1.0)*(VAL_Y1));
  637. *
  638. * variable for the convergence test
  639. ID1_CON = 0;
  640. * check on the n_limit
  641. 'SI' ((('XTX' DELTG1)**(0.5)) < 1.E-7);
  642. * if |n_limit|=0 rotate the chosen direction by 10 degrees
  643. MATR1 = 'VIDE' 'RIGIDITE'/'RIGIDITE';
  644. MATR1 = MATR1 'ET'
  645. ('MANU' 'RIGIDITE' MAI1 ('MOTS' 'UX' 'UY')
  646. ('PROG' (COS(10.0)) ((-1.0)*(SIN(10.0)))
  647. (SIN(10.0)) (COS(10.0))
  648. 0.0 0.0 ));
  649. AAA1 = MATR1*ALP_K;
  650. DELTG1 = MANU CHPO MAI1 'UX' (EXTR AAA1 MAI1 'FX')
  651. 'UY' (EXTR AAA1 MAI1 'FY');
  652. 'FINSI';
  653. * convergence test
  654. VET_TEST = (TAB_ITE. I_M1 . 'ALPA')
  655. - (TAB_ITE. (I_M1 - 1) . 'ALPA');
  656. *
  657. 'SI' ((('XTX' VET_TEST)**(0.5)) '<' 0.001);
  658. * deduction of the final value of Lambda
  659. 'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0);
  660. 'SI' ('NEG' ('TYPE' LAM_FIN) 'FLOTTANT');
  661. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  662. 'SINON';
  663. 'SI' (LAM_FIN '>' (TAB_ITE. I_M1 . 'LAMBDA'));
  664. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  665. 'FINSI';
  666. 'FINSI';
  667. 'FINSI';
  668. 'QUIT' IND2;
  669. 'SINON';
  670. 'SI' (I_M1 'EGA' 200);
  671. ID1_CON = 1;
  672. 'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0);
  673. 'SI' ('NEG' ('TYPE' LAM_FIN) 'FLOTTANT');
  674. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  675. 'SINON';
  676. 'SI' (LAM_FIN '>' (TAB_ITE. I_M1 . 'LAMBDA'));
  677. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  678. 'FINSI';
  679. 'FINSI';
  680. 'FINSI';
  681. 'SINON';
  682. * otherwise update the n_limit
  683. ALP_K = DELTG1/((XTX DELTG1)**(0.5));
  684. 'FINSI';
  685. 'FINSI';
  686. I_M1 = I_M1 + 1;
  687. 'FIN' IND2;
  688. I_K1 = I_K1 + 1;
  689. 'FIN' IND3;
  690.  
  691. *
  692. * screen message
  693. 'SI' (ID1_CON 'EGA' 0);
  694. TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda='
  695. LAM_FIN;
  696. 'MESS' TIT1;
  697. 'FINSI';
  698. 'SI' (ID1_CON 'EGA' 1);
  699. TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda='
  700. LAM_FIN '.....Attention non convergence';
  701.  
  702. 'SI' ('EGA' ('TYPE' LL_ERRE1) 'LISTENTI');
  703. LL_ERRE1 = LL_ERRE1 ET ('LECT' I_ELE1);
  704. 'SINON';
  705. LL_ERRE1 = ('LECT' I_ELE1);
  706. 'FINSI';
  707. 'FINSI';
  708. *
  709. 'SI' (I_IE1 'EGA' 3);
  710. MOD1 = 'MODE' ELE_VOI1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';
  711. MODTOT = MODTOT 'ET' MOD1;
  712. CH_LAMC = CH_LAMC 'ET' ('MANU' 'CHML' ELE_VOI1 'LAMC' LAM_FIN);
  713. 'SI' (LAM_FIN '<' 1.0);
  714. MAI_REC = MAI_REC 'ET' (ELE_VOI1 'COUL' 'ROUG');
  715. 'SINON';
  716. MAI_REC = MAI_REC 'ET' (ELE_VOI1 'COUL' 'VERT');
  717. 'FINSI';
  718. 'FINSI';
  719. 'OUBL' LAM_FIN;
  720. 'OUBL' LAM_FIN2;
  721. 'FINSI';
  722. * Verification Externe et interne-------------------------------
  723. 'SINON';
  724. * --------------------------------------------------------------
  725. *
  726. * Test about the static state stress
  727. TTL1 = TABLE;
  728. TTL1.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  729. TTL1.'S11' = N110;
  730. TTL1.'S22' = N220;
  731. TTL1.'S12' = N120;
  732. TTL1.'FCK' = FCBET1;
  733. TTL1.'FSK_1' = FSACI1;
  734. TTL1.'FSK_2' = FSACI1;
  735. TTL1.'RHO_1' = RHOXS;
  736. TTL1.'RHO_2' = RHOYS;
  737. TTL1.'GAMMA_C' = GAMC;
  738. TTL1.'ALPA_C' = ALPC1;
  739. TTL1.'NU_C' = NUC1;
  740. TTL1.'GAMMA_S' = GAMS1;
  741. *
  742. T_ST1 = G_ULTI2D TTL1;
  743. 'SI' (T_ST1 '&lt;EG' 0.0);
  744. LAM_FIN = 0.0;
  745. TIT1 = 'CHAINE' ' Element ' I_IE1 ' Lambda='
  746. LAM_FIN '.....Sig0 externe';
  747. 'MESS' TIT1;
  748. *
  749. * results stored in the MCHAML field
  750. 'SI' (I_IE1 'EGA' 1);
  751. MOD1 = 'MODE' ELE_VOI1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';
  752. MODTOT = MODTOT 'ET' MOD1;
  753. * CH_LAME = CH_LAME 'ET' ('MANU' 'CHML' MOD1 'LAME' LAM_FIN);
  754. CH_LAME = CH_LAME 'ET' ('MANU' 'CHML' ELE_VOI1 'LAME' LAM_FIN);
  755. 'SI' (LAM_FIN '<' 1.0);
  756. MAI_REE = MAI_REE 'ET' (ELE_VOI1 'COUL' 'ROUG');
  757. 'SINON';
  758. MAI_REE = MAI_REE 'ET' (ELE_VOI1 'COUL' 'VERT');
  759. 'FINSI';
  760. 'FINSI';
  761. 'SI' (I_IE1 'EGA' 2);
  762. MOD1 = 'MODE' ELE_VOI1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';
  763. MODTOT = MODTOT 'ET' MOD1;
  764. * CH_LAMI = CH_LAMI 'ET' ('MANU' 'CHML' MOD1 'LAMI' LAM_FIN);
  765. CH_LAMI = CH_LAMI 'ET' ('MANU' 'CHML' ELE_VOI1 'LAMI' LAM_FIN);
  766. 'SI' (LAM_FIN '<' 1.0);
  767. MAI_REI = MAI_REI 'ET' (ELE_VOI1 'COUL' 'ROUG');
  768. 'SINON';
  769. MAI_REI = MAI_REI 'ET' (ELE_VOI1 'COUL' 'VERT');
  770. 'FINSI';
  771. 'FINSI';
  772. *
  773. 'OUBL' LAM_FIN;
  774. *
  775. 'SINON';
  776. *
  777. * Start of the iteration
  778. P1 = 0. 0. 0.;
  779. P2 = 1. 0. 0.;
  780. P3 = 2. 0. 0.;
  781. *
  782. MAI1 = P1;
  783. * Matrix X
  784. MATX1 = VIDE 'RIGIDITE'/'RIGIDITE';
  785. MATX1 = MATX1 'ET' ('MANU' 'RIGIDITE' MAI1 (MOTS 'UX' 'UY' 'UZ')
  786. ('PROG' (N11*N11) N11N22 N11N12
  787. N11N22 (N22*N22) N22N12
  788. N11N12 N22N12 (N12*N12)));
  789. *
  790. * Static state stress
  791. EF_ST0 = 'MANU' 'CHPO' MAI1 3 'FX' N110 'FY' N220 'FZ' N120;
  792. EF_ST0 = 'CHAN' 'ATTRIBUT' EF_ST0 'NATURE' 'DISCRET';
  793. *
  794. * Initial Vecteur
  795. ALP_0 = 'MANU' 'CHPO' MAI1 3 'UX' 2.0 'UY' 2.0 'UZ' 2.0;
  796. XALP_0 = (MATX1*ALP_0)/(('XTMX' MATX1 ALP_0)**(0.5));
  797. XALP_0 = ('CHAN' ATTRIBUT XALP_0 NATURE DISCRET) + EF_ST0;
  798. SB0 = ('XTX' (XALP_0 - EF_ST0))**(0.5);
  799. RB0 = SB0;
  800. LAMB0 = RB0/SB0;
  801. *
  802. * incremental step
  803. H_VAL1 = (SB0)/(100000.0);
  804. H_VAL2 = (SB0)/(100000.0);
  805. *
  806. * direction used for the ellispoide methode
  807. TB_DIR = 'TABLE';
  808. TB_DIR. 1 = 'MANU' 'CHPO' MAI1 3 'UX' 1.0 'UY' 0.0 'UZ' 0.0;
  809. TB_DIR. 2 = 'MANU' 'CHPO' MAI1 3 'UX' 0.0 'UY' 1.0 'UZ' 0.0;
  810. TB_DIR. 3 = 'MANU' 'CHPO' MAI1 3 'UX' 0.0 'UY' 0.0 'UZ' 1.0;
  811. TB_DIR. 4 = 'MANU' 'CHPO' MAI1 3 'UX' -1.0 'UY' 0.0 'UZ' 0.0;
  812. TB_DIR. 5 = 'MANU' 'CHPO' MAI1 3 'UX' 0.0 'UY' -1.0 'UZ' 0.0;
  813. TB_DIR. 6 = 'MANU' 'CHPO' MAI1 3 'UX' 0.0 'UY' 0.0 'UZ' -1.0;
  814. *
  815. * loop for each direction
  816. I_K1 = 1;
  817. 'REPE' IND3 (6);
  818. *
  819. TAB_ITE = 'TABLE';
  820. TAB_ITE. 0 = 'TABLE';
  821. TAB_ITE. 0 . 'ALPA' = ALP_0;
  822. TAB_ITE. 0 . 'LAMBDA' = LAMB0;
  823. TAB_ITE. 0 . 'RB' = RB0;
  824. *
  825. ALP_K = (TB_DIR. I_K1);
  826. I_M1 = 1;
  827. * first loop in order to have n_elllips=n_limite
  828. 'REPE' IND2 (200);
  829. * sauving Lambda for j iteraction
  830. TAB_ITE. I_M1 = 'TABLE';
  831. TAB_ITE. I_M1 . 'ALPA' = ALP_K;
  832. * determination of the stress state for the chosen direction
  833. XALP_J = (MATX1*ALP_K)/(('XTMX' MATX1 ALP_K)**(0.5));
  834. XALP_J = ('CHAN' ATTRIBUT XALP_J NATURE DISCRET) + EF_ST0;
  835. * nomalisez vector of hte rayon of the ellipsoide for the
  836. * chosen direction
  837. SBJ = ('XTX' (XALP_J - EF_ST0))**(0.5);
  838. BJ = (XALP_J - EF_ST0)/(SBJ);
  839. * rayon of the ellipsoide for the chosen direction
  840. RBVAL1 = (TAB_ITE.(I_M1 - 1). 'RB');
  841. * second loop in order to have f=0
  842. I_I1 = 1;
  843. REPE IND1 (200);
  844. * strees state x0, x0+h, x0-h
  845. VV1 = (RBVAL1 + H_VAL1)*(BJ);
  846. VV2 = (RBVAL1 - H_VAL1)*(BJ);
  847. VV3 = (RBVAL1)*(BJ);
  848. VV1 = ('CHAN' 'ATTRIBUT' VV1 'NATURE' 'DISCRET') + EF_ST0;
  849. VV2 = ('CHAN' 'ATTRIBUT' VV2 'NATURE' 'DISCRET') + EF_ST0;
  850. VV3 = ('CHAN' 'ATTRIBUT' VV3 'NATURE' 'DISCRET') + EF_ST0;
  851. *
  852. N11V1 = 'EXTR' VV1 'FX' MAI1;
  853. N22V1 = 'EXTR' VV1 'FY' MAI1;
  854. N12V1 = 'EXTR' VV1 'FZ' MAI1;
  855. N11V2 = 'EXTR' VV2 'FX' MAI1;
  856. N22V2 = 'EXTR' VV2 'FY' MAI1;
  857. N12V2 = 'EXTR' VV2 'FZ' MAI1;
  858. N11V3 = 'EXTR' VV3 'FX' MAI1;
  859. N22V3 = 'EXTR' VV3 'FY' MAI1;
  860. N12V3 = 'EXTR' VV3 'FZ' MAI1;
  861. 'OUBL' VV1;
  862. 'OUBL' VV2;
  863. 'OUBL' VV3;
  864. * calculation of f and D terms
  865. TTL1 = TABLE;
  866. TTL1.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  867. TTL1.'S11' = N11V3;
  868. TTL1.'S22' = N22V3;
  869. TTL1.'S12' = N12V3;
  870. TTL1.'FCK' = FCBET1;
  871. TTL1.'FSK_1' = FSACI1;
  872. TTL1.'FSK_2' = FSACI1;
  873. TTL1.'RHO_1' = RHOXS;
  874. TTL1.'RHO_2' = RHOYS;
  875. TTL1.'GAMMA_C' = GAMC;
  876. TTL1.'ALPA_C' = ALPC1;
  877. TTL1.'NU_C' = NUC1;
  878. TTL1.'GAMMA_S' = GAMS1;
  879. *
  880. TTL2 = TABLE;
  881. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  882. TTL2.'S11' = N11V1;
  883. TTL2.'S22' = N22V1;
  884. TTL2.'S12' = N12V1;
  885. TTL2.'FCK' = FCBET1;
  886. TTL2.'FSK_1' = FSACI1;
  887. TTL2.'FSK_2' = FSACI1;
  888. TTL2.'RHO_1' = RHOXS;
  889. TTL2.'RHO_2' = RHOYS;
  890. TTL2.'GAMMA_C' = GAMC;
  891. TTL2.'ALPA_C' = ALPC1;
  892. TTL2.'NU_C' = NUC1;
  893. TTL2.'GAMMA_S' = GAMS1;
  894. *
  895. TTL3 = TABLE;
  896. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  897. TTL3.'S11' = N11V2;
  898. TTL3.'S22' = N22V2;
  899. TTL3.'S12' = N12V2;
  900. TTL3.'FCK' = FCBET1;
  901. TTL3.'FSK_1' = FSACI1;
  902. TTL3.'FSK_2' = FSACI1;
  903. TTL3.'RHO_1' = RHOXS;
  904. TTL3.'RHO_2' = RHOYS;
  905. TTL3.'GAMMA_C' = GAMC;
  906. TTL3.'ALPA_C' = ALPC1;
  907. TTL3.'NU_C' = NUC1;
  908. TTL3.'GAMMA_S' = GAMS1;
  909. *
  910. F_VAL1 = G_ULTI2D TTL1;
  911. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  912. ((2.0)*(H_VAL1));
  913. * controle
  914. 'SI' (('ABS'(F_VAL1)) '<' TOL_1);
  915. * convergence!!!
  916. RBJ = RBVAL1;
  917. 'QUIT' IND1;
  918. 'SINON';
  919. 'SI' (I_I1 'EGA' 200);
  920. * fin of loop, we take the last value!!!
  921. RBJ = RBVAL1;
  922. 'QUIT' IND1;
  923. 'FINSI';
  924. * Check beetwen 2 values of R (radius)!!!
  925. 'SI' (I_I1 > 2);
  926. 'SI' ((ABS((RBVAL1 - RBVAL2)/(RBVAL1))) < ((TOL_1)*(1.e-10)));
  927. RBJ = RBVAL1;
  928. 'QUIT' IND1;
  929. 'FINSI';
  930. 'FINSI';
  931. * Modification in case of the x0+h stress case is out of the
  932. * external boudary of the limit surface
  933. * [(-fc-rx*fs),(-fc-ry*fs)]x[(rx*fs),(ry*fs)]
  934. 'SI' (N11V3 > ((FSACI1*RHOXS)/(GAMS1)));
  935. RBVAL2 = RBVAL1;
  936. RBVAL1 = (((FSACI1*RHOXS)/(GAMS1)) -
  937. (EXTR EF_ST0 'FX' MAI1) -
  938. TOL_1)/
  939. (EXTR BJ 'FX' MAI1);
  940. 'SINON';
  941. 'SI' (N22V3 '>' ((FSACI1*RHOYS)/(GAMS1)));
  942. RBVAL2 = RBVAL1;
  943. RBVAL1 = (((FSACI1*RHOYS)/(GAMS1))
  944. - ('EXTR' EF_ST0 'FY' MAI1)
  945. - TOL_1)/
  946. ('EXTR' BJ 'FY' MAI1);
  947. 'SINON';
  948. 'SI' (N11V3 '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  949. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))));
  950. RBVAL2 = RBVAL1;
  951. RBVAL1 = ((((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  952. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1))))
  953. - ('EXTR' EF_ST0 'FX' MAI1)
  954. + TOL_1)/
  955. ('EXTR' BJ 'FX' MAI1);
  956. 'SINON';
  957. 'SI' (N22V3 < (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  958. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1)))));
  959. RBVAL2 = RBVAL1;
  960. RBVAL1 = ((((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  961. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1))))
  962. - ('EXTR' EF_ST0 'FY' MAI1)
  963. + TOL_1)/
  964. ('EXTR' BJ 'FY' MAI1);
  965. 'SINON';
  966. * case when D is equal to sero
  967. 'SI' (('ABS'(D_VAL1)) < 1.E-16);
  968. RBVAL2 = RBVAL1;
  969. RBVAL1 = ((1.1)*(RBVAL1));
  970. 'SINON';
  971. * otherwise
  972. RBVAL2 = RBVAL1;
  973. RBVAL1 = RBVAL1 - ((F_VAL1)/(D_VAL1));
  974. 'FINSI';
  975. 'FINSI';
  976. 'FINSI';
  977. 'FINSI';
  978. 'FINSI';
  979. 'FINSI';
  980. *
  981. I_I1 = I_I1 + 1;
  982. 'FIN' IND1;
  983. * computation of lambda
  984. LAMBJ = RBJ/SBJ;
  985. TAB_ITE. I_M1 . 'LAMBDA' = LAMBJ;
  986. TAB_ITE. I_M1 . 'RB' = RBJ;
  987. XCOUB1 = (((TAB_ITE. I_M1 .'LAMBDA')*(TAB_ITE. I_M1 .'LAMBDA'))*
  988. (MATX1));
  989. * amplified stress state using the lambda coef
  990. XALP_L = (XCOUB1*ALP_K)/(('XTMX' XCOUB1 ALP_K)**(0.5));
  991. XALP_L = ('CHAN' ATTRIBUT XALP_L NATURE DISCRET) + EF_ST0;
  992. * Dh*ex ; Dh*ey ; Dh*ez
  993. DE_HX = 'MANU' 'CHPO' MAI1 3 'FX' H_VAL2 'FY' 0.0 'FZ' 0.0;
  994. DE_HY = 'MANU' 'CHPO' MAI1 3 'FX' 0.0 'FY' H_VAL2 'FZ' 0.0;
  995. DE_HZ = 'MANU' 'CHPO' MAI1 3 'FX' 0.0 'FY' 0.0 'FZ' H_VAL2;
  996. * x0*lam + Dh*ex ; x0*lam + Dh*ey ; x0*lam + Dh*ez
  997. VV_HXP = XALP_L + DE_HX;
  998. VV_HYP = XALP_L + DE_HY;
  999. VV_HZP = XALP_L + DE_HZ;
  1000. * x0*lam + Dh*ex ; x0*lam + Dh*ey ; x0*lam + Dh*ez
  1001. VV_HXN = XALP_L - DE_HX;
  1002. VV_HYN = XALP_L - DE_HY;
  1003. VV_HZN = XALP_L - DE_HZ;
  1004. * x0*lam
  1005. VV_0 = XALP_L;
  1006. *
  1007. VV_HXP1 = 'EXTR' VV_HXP 'FX' MAI1;
  1008. VV_HXP2 = 'EXTR' VV_HXP 'FY' MAI1;
  1009. VV_HXP3 = 'EXTR' VV_HXP 'FZ' MAI1;
  1010. VV_HYP1 = 'EXTR' VV_HYP 'FX' MAI1;
  1011. VV_HYP2 = 'EXTR' VV_HYP 'FY' MAI1;
  1012. VV_HYP3 = 'EXTR' VV_HYP 'FZ' MAI1;
  1013. VV_HZP1 = 'EXTR' VV_HZP 'FX' MAI1;
  1014. VV_HZP2 = 'EXTR' VV_HZP 'FY' MAI1;
  1015. VV_HZP3 = 'EXTR' VV_HZP 'FZ' MAI1;
  1016. VV_HXN1 = 'EXTR' VV_HXN 'FX' MAI1;
  1017. VV_HXN2 = 'EXTR' VV_HXN 'FY' MAI1;
  1018. VV_HXN3 = 'EXTR' VV_HXN 'FZ' MAI1;
  1019. VV_HYN1 = 'EXTR' VV_HYN 'FX' MAI1;
  1020. VV_HYN2 = 'EXTR' VV_HYN 'FY' MAI1;
  1021. VV_HYN3 = 'EXTR' VV_HYN 'FZ' MAI1;
  1022. VV_HZN1 = 'EXTR' VV_HZN 'FX' MAI1;
  1023. VV_HZN2 = 'EXTR' VV_HZN 'FY' MAI1;
  1024. VV_HZN3 = 'EXTR' VV_HZN 'FZ' MAI1;
  1025. *
  1026. VV_01 = 'EXTR' VV_HXP 'FX' MAI1;
  1027. VV_02 = 'EXTR' VV_HXP 'FY' MAI1;
  1028. VV_03 = 'EXTR' VV_HXP 'FZ' MAI1;
  1029. * Calculation of n_limite if x0*lam + Dh*ex is out of the external
  1030. * boudary of limite surface
  1031. 'SI' (VV_HXP1 '>' ((FSACI1*RHOXS)/(GAMS1)));
  1032. *
  1033. TTL2 = TABLE;
  1034. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1035. TTL2.'S11' = VV_01;
  1036. TTL2.'S22' = VV_02;
  1037. TTL2.'S12' = VV_03;
  1038. TTL2.'FCK' = FCBET1;
  1039. TTL2.'FSK_1' = FSACI1;
  1040. TTL2.'FSK_2' = FSACI1;
  1041. TTL2.'RHO_1' = RHOXS;
  1042. TTL2.'RHO_2' = RHOYS;
  1043. TTL2.'GAMMA_C' = GAMC;
  1044. TTL2.'ALPA_C' = ALPC1;
  1045. TTL2.'NU_C' = NUC1;
  1046. TTL2.'GAMMA_S' = GAMS1;
  1047. TTL3 = TABLE;
  1048. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1049. TTL3.'S11' = VV_HXN1;
  1050. TTL3.'S22' = VV_HXN2;
  1051. TTL3.'S12' = VV_HXN3;
  1052. TTL3.'FCK' = FCBET1;
  1053. TTL3.'FSK_1' = FSACI1;
  1054. TTL3.'FSK_2' = FSACI1;
  1055. TTL3.'RHO_1' = RHOXS;
  1056. TTL3.'RHO_2' = RHOYS;
  1057. TTL3.'GAMMA_C' = GAMC;
  1058. TTL3.'ALPA_C' = ALPC1;
  1059. TTL3.'NU_C' = NUC1;
  1060. TTL3.'GAMMA_S' = GAMS1;
  1061. VAL_X1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL2));
  1062. *
  1063. TTL2 = TABLE;
  1064. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1065. TTL2.'S11' = VV_HYP1;
  1066. TTL2.'S22' = VV_HYP2;
  1067. TTL2.'S12' = VV_HYP3;
  1068. TTL2.'FCK' = FCBET1;
  1069. TTL2.'FSK_1' = FSACI1;
  1070. TTL2.'FSK_2' = FSACI1;
  1071. TTL2.'RHO_1' = RHOXS;
  1072. TTL2.'RHO_2' = RHOYS;
  1073. TTL2.'GAMMA_C' = GAMC;
  1074. TTL2.'ALPA_C' = ALPC1;
  1075. TTL2.'NU_C' = NUC1;
  1076. TTL2.'GAMMA_S' = GAMS1;
  1077. TTL3 = TABLE;
  1078. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1079. TTL3.'S11' = VV_HYN1;
  1080. TTL3.'S22' = VV_HYN2;
  1081. TTL3.'S12' = VV_HYN3;
  1082. TTL3.'FCK' = FCBET1;
  1083. TTL3.'FSK_1' = FSACI1;
  1084. TTL3.'FSK_2' = FSACI1;
  1085. TTL3.'RHO_1' = RHOXS;
  1086. TTL3.'RHO_2' = RHOYS;
  1087. TTL3.'GAMMA_C' = GAMC;
  1088. TTL3.'ALPA_C' = ALPC1;
  1089. TTL3.'NU_C' = NUC1;
  1090. TTL3.'GAMMA_S' = GAMS1;
  1091. VAL_Y1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  1092. *
  1093. TTL2 = TABLE;
  1094. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1095. TTL2.'S11' = VV_HZP1;
  1096. TTL2.'S22' = VV_HZP2;
  1097. TTL2.'S12' = VV_HZP3;
  1098. TTL2.'FCK' = FCBET1;
  1099. TTL2.'FSK_1' = FSACI1;
  1100. TTL2.'FSK_2' = FSACI1;
  1101. TTL2.'RHO_1' = RHOXS;
  1102. TTL2.'RHO_2' = RHOYS;
  1103. TTL2.'GAMMA_C' = GAMC;
  1104. TTL2.'ALPA_C' = ALPC1;
  1105. TTL2.'NU_C' = NUC1;
  1106. TTL2.'GAMMA_S' = GAMS1;
  1107. TTL3 = TABLE;
  1108. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1109. TTL3.'S11' = VV_HZN1;
  1110. TTL3.'S22' = VV_HZN2;
  1111. TTL3.'S12' = VV_HZN3;
  1112. TTL3.'FCK' = FCBET1;
  1113. TTL3.'FSK_1' = FSACI1;
  1114. TTL3.'FSK_2' = FSACI1;
  1115. TTL3.'RHO_1' = RHOXS;
  1116. TTL3.'RHO_2' = RHOYS;
  1117. TTL3.'GAMMA_C' = GAMC;
  1118. TTL3.'ALPA_C' = ALPC1;
  1119. TTL3.'NU_C' = NUC1;
  1120. TTL3.'GAMMA_S' = GAMS1;
  1121. VAL_Z1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  1122. 'SINON';
  1123. * Calculation of n_limit if x0*lam - Dh*ex is out of the external
  1124. * boudary of the limite surface
  1125. 'SI' (VV_HXN1 '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1126. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))));
  1127. *
  1128. TTL2 = TABLE;
  1129. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1130. TTL2.'S11' = VV_HXP1;
  1131. TTL2.'S22' = VV_HXP2;
  1132. TTL2.'S12' = VV_HXP3;
  1133. TTL2.'FCK' = FCBET1;
  1134. TTL2.'FSK_1' = FSACI1;
  1135. TTL2.'FSK_2' = FSACI1;
  1136. TTL2.'RHO_1' = RHOXS;
  1137. TTL2.'RHO_2' = RHOYS;
  1138. TTL2.'GAMMA_C' = GAMC;
  1139. TTL2.'ALPA_C' = ALPC1;
  1140. TTL2.'NU_C' = NUC1;
  1141. TTL2.'GAMMA_S' = GAMS1;
  1142. TTL3 = TABLE;
  1143. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1144. TTL3.'S11' = VV_01;
  1145. TTL3.'S22' = VV_02;
  1146. TTL3.'S12' = VV_03;
  1147. TTL3.'FCK' = FCBET1;
  1148. TTL3.'FSK_1' = FSACI1;
  1149. TTL3.'FSK_2' = FSACI1;
  1150. TTL3.'RHO_1' = RHOXS;
  1151. TTL3.'RHO_2' = RHOYS;
  1152. TTL3.'GAMMA_C' = GAMC;
  1153. TTL3.'ALPA_C' = ALPC1;
  1154. TTL3.'NU_C' = NUC1;
  1155. TTL3.'GAMMA_S' = GAMS1;
  1156. VAL_X1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL2));
  1157. *
  1158. TTL2 = TABLE;
  1159. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1160. TTL2.'S11' = VV_HYP1;
  1161. TTL2.'S22' = VV_HYP2;
  1162. TTL2.'S12' = VV_HYP3;
  1163. TTL2.'FCK' = FCBET1;
  1164. TTL2.'FSK_1' = FSACI1;
  1165. TTL2.'FSK_2' = FSACI1;
  1166. TTL2.'RHO_1' = RHOXS;
  1167. TTL2.'RHO_2' = RHOYS;
  1168. TTL2.'GAMMA_C' = GAMC;
  1169. TTL2.'ALPA_C' = ALPC1;
  1170. TTL2.'NU_C' = NUC1;
  1171. TTL2.'GAMMA_S' = GAMS1;
  1172. TTL3 = TABLE;
  1173. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1174. TTL3.'S11' = VV_HYN1;
  1175. TTL3.'S22' = VV_HYN2;
  1176. TTL3.'S12' = VV_HYN3;
  1177. TTL3.'FCK' = FCBET1;
  1178. TTL3.'FSK_1' = FSACI1;
  1179. TTL3.'FSK_2' = FSACI1;
  1180. TTL3.'RHO_1' = RHOXS;
  1181. TTL3.'RHO_2' = RHOYS;
  1182. TTL3.'GAMMA_C' = GAMC;
  1183. TTL3.'ALPA_C' = ALPC1;
  1184. TTL3.'NU_C' = NUC1;
  1185. TTL3.'GAMMA_S' = GAMS1;
  1186. VAL_Y1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  1187. *
  1188. TTL2 = TABLE;
  1189. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1190. TTL2.'S11' = VV_HZP1;
  1191. TTL2.'S22' = VV_HZP2;
  1192. TTL2.'S12' = VV_HZP3;
  1193. TTL2.'FCK' = FCBET1;
  1194. TTL2.'FSK_1' = FSACI1;
  1195. TTL2.'FSK_2' = FSACI1;
  1196. TTL2.'RHO_1' = RHOXS;
  1197. TTL2.'RHO_2' = RHOYS;
  1198. TTL2.'GAMMA_C' = GAMC;
  1199. TTL2.'ALPA_C' = ALPC1;
  1200. TTL2.'NU_C' = NUC1;
  1201. TTL2.'GAMMA_S' = GAMS1;
  1202. TTL3 = TABLE;
  1203. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1204. TTL3.'S11' = VV_HZN1;
  1205. TTL3.'S22' = VV_HZN2;
  1206. TTL3.'S12' = VV_HZN3;
  1207. TTL3.'FCK' = FCBET1;
  1208. TTL3.'FSK_1' = FSACI1;
  1209. TTL3.'FSK_2' = FSACI1;
  1210. TTL3.'RHO_1' = RHOXS;
  1211. TTL3.'RHO_2' = RHOYS;
  1212. TTL3.'GAMMA_C' = GAMC;
  1213. TTL3.'ALPA_C' = ALPC1;
  1214. TTL3.'NU_C' = NUC1;
  1215. TTL3.'GAMMA_S' = GAMS1;
  1216. VAL_Z1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  1217. 'SINON';
  1218. * Calculation of n_limit if x0*lam + Dh*ey is out of the external
  1219. * boudary of the limite surface
  1220. 'SI' (VV_HYP1 '>' ((FSACI1*RHOYS)/(GAMS1)));
  1221. *
  1222. TTL2 = TABLE;
  1223. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1224. TTL2.'S11' = VV_HXP1;
  1225. TTL2.'S22' = VV_HXP2;
  1226. TTL2.'S12' = VV_HXP3;
  1227. TTL2.'FCK' = FCBET1;
  1228. TTL2.'FSK_1' = FSACI1;
  1229. TTL2.'FSK_2' = FSACI1;
  1230. TTL2.'RHO_1' = RHOXS;
  1231. TTL2.'RHO_2' = RHOYS;
  1232. TTL2.'GAMMA_C' = GAMC;
  1233. TTL2.'ALPA_C' = ALPC1;
  1234. TTL2.'NU_C' = NUC1;
  1235. TTL2.'GAMMA_S' = GAMS1;
  1236. TTL3 = TABLE;
  1237. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1238. TTL3.'S11' = VV_HXN1;
  1239. TTL3.'S22' = VV_HXN2;
  1240. TTL3.'S12' = VV_HXN3;
  1241. TTL3.'FCK' = FCBET1;
  1242. TTL3.'FSK_1' = FSACI1;
  1243. TTL3.'FSK_2' = FSACI1;
  1244. TTL3.'RHO_1' = RHOXS;
  1245. TTL3.'RHO_2' = RHOYS;
  1246. TTL3.'GAMMA_C' = GAMC;
  1247. TTL3.'ALPA_C' = ALPC1;
  1248. TTL3.'NU_C' = NUC1;
  1249. TTL3.'GAMMA_S' = GAMS1;
  1250. VAL_X1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  1251. *
  1252. TTL2 = TABLE;
  1253. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1254. TTL2.'S11' = VV_01;
  1255. TTL2.'S22' = VV_02;
  1256. TTL2.'S12' = VV_03;
  1257. TTL2.'FCK' = FCBET1;
  1258. TTL2.'FSK_1' = FSACI1;
  1259. TTL2.'FSK_2' = FSACI1;
  1260. TTL2.'RHO_1' = RHOXS;
  1261. TTL2.'RHO_2' = RHOYS;
  1262. TTL2.'GAMMA_C' = GAMC;
  1263. TTL2.'ALPA_C' = ALPC1;
  1264. TTL2.'NU_C' = NUC1;
  1265. TTL2.'GAMMA_S' = GAMS1;
  1266. TTL3 = TABLE;
  1267. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1268. TTL3.'S11' = VV_HYN1;
  1269. TTL3.'S22' = VV_HYN2;
  1270. TTL3.'S12' = VV_HYN3;
  1271. TTL3.'FCK' = FCBET1;
  1272. TTL3.'FSK_1' = FSACI1;
  1273. TTL3.'FSK_2' = FSACI1;
  1274. TTL3.'RHO_1' = RHOXS;
  1275. TTL3.'RHO_2' = RHOYS;
  1276. TTL3.'GAMMA_C' = GAMC;
  1277. TTL3.'ALPA_C' = ALPC1;
  1278. TTL3.'NU_C' = NUC1;
  1279. TTL3.'GAMMA_S' = GAMS1;
  1280. VAL_Y1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL2));
  1281. *
  1282. TTL2 = TABLE;
  1283. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1284. TTL2.'S11' = VV_HZP1;
  1285. TTL2.'S22' = VV_HZP2;
  1286. TTL2.'S12' = VV_HZP3;
  1287. TTL2.'FCK' = FCBET1;
  1288. TTL2.'FSK_1' = FSACI1;
  1289. TTL2.'FSK_2' = FSACI1;
  1290. TTL2.'RHO_1' = RHOXS;
  1291. TTL2.'RHO_2' = RHOYS;
  1292. TTL2.'GAMMA_C' = GAMC;
  1293. TTL2.'ALPA_C' = ALPC1;
  1294. TTL2.'NU_C' = NUC1;
  1295. TTL2.'GAMMA_S' = GAMS1;
  1296. TTL3 = TABLE;
  1297. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1298. TTL3.'S11' = VV_HZN1;
  1299. TTL3.'S22' = VV_HZN2;
  1300. TTL3.'S12' = VV_HZN3;
  1301. TTL3.'FCK' = FCBET1;
  1302. TTL3.'FSK_1' = FSACI1;
  1303. TTL3.'FSK_2' = FSACI1;
  1304. TTL3.'RHO_1' = RHOXS;
  1305. TTL3.'RHO_2' = RHOYS;
  1306. TTL3.'GAMMA_C' = GAMC;
  1307. TTL3.'ALPA_C' = ALPC1;
  1308. TTL3.'NU_C' = NUC1;
  1309. TTL3.'GAMMA_S' = GAMS1;
  1310. VAL_Z1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL2));
  1311. 'SINON';
  1312. * Calculation of n_limit if x0*lam + Dh*ey is out of the external
  1313. * boudary of the limite surface
  1314. 'SI' (VV_HYN1 '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1315. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1)))));
  1316. *
  1317. TTL2 = TABLE;
  1318. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1319. TTL2.'S11' = VV_HXP1;
  1320. TTL2.'S22' = VV_HXP2;
  1321. TTL2.'S12' = VV_HXP3;
  1322. TTL2.'FCK' = FCBET1;
  1323. TTL2.'FSK_1' = FSACI1;
  1324. TTL2.'FSK_2' = FSACI1;
  1325. TTL2.'RHO_1' = RHOXS;
  1326. TTL2.'RHO_2' = RHOYS;
  1327. TTL2.'GAMMA_C' = GAMC;
  1328. TTL2.'ALPA_C' = ALPC1;
  1329. TTL2.'NU_C' = NUC1;
  1330. TTL2.'GAMMA_S' = GAMS1;
  1331. TTL3 = TABLE;
  1332. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1333. TTL3.'S11' = VV_HXN1;
  1334. TTL3.'S22' = VV_HXN2;
  1335. TTL3.'S12' = VV_HXN3;
  1336. TTL3.'FCK' = FCBET1;
  1337. TTL3.'FSK_1' = FSACI1;
  1338. TTL3.'FSK_2' = FSACI1;
  1339. TTL3.'RHO_1' = RHOXS;
  1340. TTL3.'RHO_2' = RHOYS;
  1341. TTL3.'GAMMA_C' = GAMC;
  1342. TTL3.'ALPA_C' = ALPC1;
  1343. TTL3.'NU_C' = NUC1;
  1344. TTL3.'GAMMA_S' = GAMS1;
  1345. VAL_X1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  1346. ((2.0)*(H_VAL2));
  1347. *
  1348. TTL2 = TABLE;
  1349. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1350. TTL2.'S11' = VV_HYP1;
  1351. TTL2.'S22' = VV_HYP2;
  1352. TTL2.'S12' = VV_HYP3;
  1353. TTL2.'FCK' = FCBET1;
  1354. TTL2.'FSK_1' = FSACI1;
  1355. TTL2.'FSK_2' = FSACI1;
  1356. TTL2.'RHO_1' = RHOXS;
  1357. TTL2.'RHO_2' = RHOYS;
  1358. TTL2.'GAMMA_C' = GAMC;
  1359. TTL2.'ALPA_C' = ALPC1;
  1360. TTL2.'NU_C' = NUC1;
  1361. TTL2.'GAMMA_S' = GAMS1;
  1362. TTL3 = TABLE;
  1363. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1364. TTL3.'S11' = VV_01;
  1365. TTL3.'S22' = VV_02;
  1366. TTL3.'S12' = VV_03;
  1367. TTL3.'FCK' = FCBET1;
  1368. TTL3.'FSK_1' = FSACI1;
  1369. TTL3.'FSK_2' = FSACI1;
  1370. TTL3.'RHO_1' = RHOXS;
  1371. TTL3.'RHO_2' = RHOYS;
  1372. TTL3.'GAMMA_C' = GAMC;
  1373. TTL3.'ALPA_C' = ALPC1;
  1374. TTL3.'NU_C' = NUC1;
  1375. TTL3.'GAMMA_S' = GAMS1;
  1376. VAL_Y1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  1377. ((1.0)*(H_VAL2));
  1378. *
  1379. TTL2 = TABLE;
  1380. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1381. TTL2.'S11' = VV_HZP1;
  1382. TTL2.'S22' = VV_HZP2;
  1383. TTL2.'S12' = VV_HZP3;
  1384. TTL2.'FCK' = FCBET1;
  1385. TTL2.'FSK_1' = FSACI1;
  1386. TTL2.'FSK_2' = FSACI1;
  1387. TTL2.'RHO_1' = RHOXS;
  1388. TTL2.'RHO_2' = RHOYS;
  1389. TTL2.'GAMMA_C' = GAMC;
  1390. TTL2.'ALPA_C' = ALPC1;
  1391. TTL2.'NU_C' = NUC1;
  1392. TTL2.'GAMMA_S' = GAMS1;
  1393. TTL3 = TABLE;
  1394. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1395. TTL3.'S11' = VV_HZN1;
  1396. TTL3.'S22' = VV_HZN2;
  1397. TTL3.'S12' = VV_HZN3;
  1398. TTL3.'FCK' = FCBET1;
  1399. TTL3.'FSK_1' = FSACI1;
  1400. TTL3.'FSK_2' = FSACI1;
  1401. TTL3.'RHO_1' = RHOXS;
  1402. TTL3.'RHO_2' = RHOYS;
  1403. TTL3.'GAMMA_C' = GAMC;
  1404. TTL3.'ALPA_C' = ALPC1;
  1405. TTL3.'NU_C' = NUC1;
  1406. TTL3.'GAMMA_S' = GAMS1;
  1407. VAL_Z1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  1408. ((2.0)*(H_VAL2));
  1409. 'SINON';
  1410. * Ordinary method
  1411. TTL2 = TABLE;
  1412. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1413. TTL2.'S11' = VV_HXP1;
  1414. TTL2.'S22' = VV_HXP2;
  1415. TTL2.'S12' = VV_HXP3;
  1416. TTL2.'FCK' = FCBET1;
  1417. TTL2.'FSK_1' = FSACI1;
  1418. TTL2.'FSK_2' = FSACI1;
  1419. TTL2.'RHO_1' = RHOXS;
  1420. TTL2.'RHO_2' = RHOYS;
  1421. TTL2.'GAMMA_C' = GAMC;
  1422. TTL2.'ALPA_C' = ALPC1;
  1423. TTL2.'NU_C' = NUC1;
  1424. TTL2.'GAMMA_S' = GAMS1;
  1425. TTL3 = TABLE;
  1426. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1427. TTL3.'S11' = VV_HXN1;
  1428. TTL3.'S22' = VV_HXN2;
  1429. TTL3.'S12' = VV_HXN3;
  1430. TTL3.'FCK' = FCBET1;
  1431. TTL3.'FSK_1' = FSACI1;
  1432. TTL3.'FSK_2' = FSACI1;
  1433. TTL3.'RHO_1' = RHOXS;
  1434. TTL3.'RHO_2' = RHOYS;
  1435. TTL3.'GAMMA_C' = GAMC;
  1436. TTL3.'ALPA_C' = ALPC1;
  1437. TTL3.'NU_C' = NUC1;
  1438. TTL3.'GAMMA_S' = GAMS1;
  1439. VAL_X1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  1440. ((2.0)*(H_VAL2));
  1441. *
  1442. TTL2 = TABLE;
  1443. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1444. TTL2.'S11' = VV_HYP1;
  1445. TTL2.'S22' = VV_HYP2;
  1446. TTL2.'S12' = VV_HYP3;
  1447. TTL2.'FCK' = FCBET1;
  1448. TTL2.'FSK_1' = FSACI1;
  1449. TTL2.'FSK_2' = FSACI1;
  1450. TTL2.'RHO_1' = RHOXS;
  1451. TTL2.'RHO_2' = RHOYS;
  1452. TTL2.'GAMMA_C' = GAMC;
  1453. TTL2.'ALPA_C' = ALPC1;
  1454. TTL2.'NU_C' = NUC1;
  1455. TTL2.'GAMMA_S' = GAMS1;
  1456. TTL3 = TABLE;
  1457. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1458. TTL3.'S11' = VV_HYN1;
  1459. TTL3.'S22' = VV_HYN2;
  1460. TTL3.'S12' = VV_HYN3;
  1461. TTL3.'FCK' = FCBET1;
  1462. TTL3.'FSK_1' = FSACI1;
  1463. TTL3.'FSK_2' = FSACI1;
  1464. TTL3.'RHO_1' = RHOXS;
  1465. TTL3.'RHO_2' = RHOYS;
  1466. TTL3.'GAMMA_C' = GAMC;
  1467. TTL3.'ALPA_C' = ALPC1;
  1468. TTL3.'NU_C' = NUC1;
  1469. TTL3.'GAMMA_S' = GAMS1;
  1470. VAL_Y1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  1471. ((2.0)*(H_VAL2));
  1472. *
  1473. TTL2 = TABLE;
  1474. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1475. TTL2.'S11' = VV_HZP1;
  1476. TTL2.'S22' = VV_HZP2;
  1477. TTL2.'S12' = VV_HZP3;
  1478. TTL2.'FCK' = FCBET1;
  1479. TTL2.'FSK_1' = FSACI1;
  1480. TTL2.'FSK_2' = FSACI1;
  1481. TTL2.'RHO_1' = RHOXS;
  1482. TTL2.'RHO_2' = RHOYS;
  1483. TTL2.'GAMMA_C' = GAMC;
  1484. TTL2.'ALPA_C' = ALPC1;
  1485. TTL2.'NU_C' = NUC1;
  1486. TTL2.'GAMMA_S' = GAMS1;
  1487. TTL3 = TABLE;
  1488. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1489. TTL3.'S11' = VV_HZN1;
  1490. TTL3.'S22' = VV_HZN2;
  1491. TTL3.'S12' = VV_HZN3;
  1492. TTL3.'FCK' = FCBET1;
  1493. TTL3.'FSK_1' = FSACI1;
  1494. TTL3.'FSK_2' = FSACI1;
  1495. TTL3.'RHO_1' = RHOXS;
  1496. TTL3.'RHO_2' = RHOYS;
  1497. TTL3.'GAMMA_C' = GAMC;
  1498. TTL3.'ALPA_C' = ALPC1;
  1499. TTL3.'NU_C' = NUC1;
  1500. TTL3.'GAMMA_S' = GAMS1;
  1501. VAL_Z1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/
  1502. ((2.0)*(H_VAL2));
  1503. 'FINSI';
  1504. 'FINSI';
  1505. 'FINSI';
  1506. 'FINSI';
  1507. *
  1508. * costruction n_limit
  1509. DELTG1 = 'MANU' 'CHPO' MAI1 'UX' ((-1.0)*(VAL_X1))
  1510. 'UY' ((-1.0)*(VAL_Y1))
  1511. 'UZ' ((-1.0)*(VAL_Z1));
  1512. *
  1513. * variable for the convergence test
  1514. ID1_CON = 0;
  1515. * check on the n_limit
  1516. 'SI' ((('XTX' DELTG1)**(0.5)) < 1.E-7);
  1517. * if |n_limit|=0 rotate the chosen direction by 10 degrees
  1518. MATR1 = 'VIDE' 'RIGIDITE'/'RIGIDITE';
  1519. MATR1 = MATR1 'ET'
  1520. ('MANU' 'RIGIDITE' MAI1 ('MOTS' 'UX' 'UY' 'UZ')
  1521. ('PROG' (COS(10.0)) ((-1.0)*(SIN(10.0))) 0.0
  1522. (SIN(10.0)) (COS(10.0)) 0.0
  1523. 0.0 0.0 1.0));
  1524. AAA1 = MATR1*ALP_K;
  1525. DELTG1 = MANU CHPO MAI1 'UX' (EXTR AAA1 MAI1 'FX')
  1526. 'UY' (EXTR AAA1 MAI1 'FY')
  1527. 'UZ' (EXTR AAA1 MAI1 'FZ');
  1528. 'FINSI';
  1529. * convergence test
  1530. VET_TEST = (TAB_ITE. I_M1 . 'ALPA')
  1531. - (TAB_ITE. (I_M1 - 1) . 'ALPA');
  1532. *
  1533. 'SI' ((('XTX' VET_TEST)**(0.5)) '<' 0.001);
  1534. * deduction of the final value of Lambda
  1535. 'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0);
  1536. 'SI' ('NEG' ('TYPE' LAM_FIN) 'FLOTTANT');
  1537. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  1538. 'SINON';
  1539. 'SI' (LAM_FIN '>' (TAB_ITE. I_M1 . 'LAMBDA'));
  1540. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  1541. 'FINSI';
  1542. 'FINSI';
  1543. 'FINSI';
  1544. 'QUIT' IND2;
  1545. 'SINON';
  1546. 'SI' (I_M1 'EGA' 200);
  1547. ID1_CON = 1;
  1548. 'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0);
  1549. 'SI' ('NEG' ('TYPE' LAM_FIN) 'FLOTTANT');
  1550. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  1551. 'SINON';
  1552. 'SI' (LAM_FIN '>' (TAB_ITE. I_M1 . 'LAMBDA'));
  1553. LAM_FIN = (TAB_ITE. I_M1 . 'LAMBDA');
  1554. 'FINSI';
  1555. 'FINSI';
  1556. 'FINSI';
  1557. 'SINON';
  1558. * otherwise updating the n_limit
  1559. ALP_K = DELTG1/((XTX DELTG1)**(0.5));
  1560. 'FINSI';
  1561. 'FINSI';
  1562. I_M1 = I_M1 + 1;
  1563. 'FIN' IND2;
  1564. I_K1 = I_K1 + 1;
  1565. 'FIN' IND3;
  1566. *
  1567. * screen message
  1568. 'SI' (ID1_CON 'EGA' 0);
  1569. TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda='
  1570. LAM_FIN;
  1571. 'MESS' TIT1;
  1572. 'FINSI';
  1573. 'SI' (ID1_CON 'EGA' 1);
  1574.  
  1575. TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda='
  1576. LAM_FIN '.....Attention non convergence';
  1577.  
  1578. * New computation of Lambda....rectangle method
  1579. * H_VAL3 = FCBET1/100000;
  1580. LAMB0 = 1.0;
  1581. H_VAL3 = LAMB0/10000;
  1582. TOL_2 = FCBET1/1000;
  1583. *
  1584. * reference list for the 8 point for rectangle solid
  1585. LL_I = PROG 1. 2. 3. 4. 5. 6. 7. 8.;
  1586. *
  1587. * creation of the resuming table
  1588. TAB_ITE = 'TABLE';
  1589. TAB_ITE. 0 = 'TABLE';
  1590. TAB_ITE. 0 .'LAMBDA' = LAMB0 ;
  1591. *
  1592. I_MM1 = 0;
  1593. ID1_CON = 0;
  1594. *
  1595. * loop for in order to find f=0
  1596. 'REPE' IND4 (200);
  1597. *
  1598. * stress state for the 8 points
  1599. X1_0 = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  1600. Y1_0 = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1601. Z1_0 = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1602. X2_0 = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1603. Y2_0 = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1604. Z2_0 = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1605. X3_0 = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1606. Y3_0 = ((N22)*(((1.0)*(LAMB0)))) + N220;
  1607. Z3_0 = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1608. X4_0 = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  1609. Y4_0 = ((N22)*(((1.0)*(LAMB0)))) + N220;
  1610. Z4_0 = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1611. X5_0 = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  1612. Y5_0 = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1613. Z5_0 = ((N12)*(((1.0)*(LAMB0)))) + N120;
  1614. X6_0 = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1615. Y6_0 = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1616. Z6_0 = ((N12)*(((1.0)*(LAMB0)))) + N120;
  1617. X7_0 = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1618. Y7_0 = ((N22)*(((1.0)*(LAMB0)))) + N220;
  1619. Z7_0 = ((N12)*(((1.0)*(LAMB0)))) + N120;
  1620. X8_0 = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  1621. Y8_0 = ((N22)*(((1.0)*(LAMB0)))) + N220;
  1622. Z8_0 = ((N12)*(((1.0)*(LAMB0)))) + N120;
  1623. *
  1624. * put the N11, N22 et N12 stress state in the lists
  1625. LL_X0 = 'PROG' (X1_0) (X2_0) (X3_0) (X4_0)
  1626. (X5_0) (X6_0) (X7_0) (X8_0);
  1627. LL_Y0 = 'PROG' (Y1_0) (Y2_0) (Y3_0) (Y4_0)
  1628. (Y5_0) (Y6_0) (Y7_0) (Y8_0);
  1629. LL_Z0 = 'PROG' (Z1_0) (Z2_0) (Z3_0) (Z4_0)
  1630. (Z5_0) (Z6_0) (Z7_0) (Z8_0);
  1631. LL_F0 = 'PROG';
  1632. *
  1633. * deduction of f for each point
  1634. III1 = 1;
  1635. 'REPE' IND2 (8);
  1636. TTL1 = TABLE;
  1637. TTL1.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1638. TTL1.'S11' = ('EXTR' LL_X0 III1);
  1639. TTL1.'S22' = ('EXTR' LL_Y0 III1);
  1640. TTL1.'S12' = ('EXTR' LL_Z0 III1);
  1641. TTL1.'FCK' = FCBET1;
  1642. TTL1.'FSK_1' = FSACI1;
  1643. TTL1.'FSK_2' = FSACI1;
  1644. TTL1.'RHO_1' = RHOXS;
  1645. TTL1.'RHO_2' = RHOYS;
  1646. TTL1.'GAMMA_C' = GAMC;
  1647. TTL1.'ALPA_C' = ALPC1;
  1648. TTL1.'NU_C' = NUC1;
  1649. TTL1.'GAMMA_S' = GAMS1;
  1650. LL_F0 = LL_F0 ET
  1651. ('PROG' (G_ULTI2D TTL1));
  1652. III1 = III1 + 1;
  1653. 'FIN' IND2;
  1654. * saving
  1655. TAB_ITE. I_MM1 .'L_X' = LL_X0;
  1656. TAB_ITE. I_MM1 .'L_Y' = LL_Y0;
  1657. TAB_ITE. I_MM1 .'L_Z' = LL_Z0;
  1658. TAB_ITE. I_MM1 .'L_F' = LL_F0;
  1659. * finding the minimum value of f and the corresponding point
  1660. EV1 = 'EVOL' 'MANU' LL_I LL_F0;
  1661. II1 II2 VAL_F = 'MINI' EV1;
  1662. * computation of D (gradiant) in case of Point_1=Argmin(fmin)
  1663. 'SI' (II1 'EGA' 1);
  1664. X1_DP = ((N11)*(((-1.0)*(LAMB0 + H_VAL3)))) + N110;
  1665. Y1_DP = ((N22)*(((-1.0)*(LAMB0 + H_VAL3)))) + N220;
  1666. Z1_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  1667. X1_DN = ((N11)*(((-1.0)*(LAMB0 - H_VAL3)))) + N110;
  1668. Y1_DN = ((N22)*(((-1.0)*(LAMB0 - H_VAL3)))) + N220;
  1669. Z1_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  1670. 'SI' (((X1_DP) '>' ((FSACI1*RHOXS)/(GAMS1)))
  1671. 'OU' ((Y1_DP) '>' ((FSACI1*RHOYS)/(GAMS1))));
  1672. X1_DP = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  1673. Y1_DP = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1674. Z1_DP = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1675. X1_DN = ((N11)*(((-1.0)*(LAMB0 - H_VAL3)))) + N110;
  1676. Y1_DN = ((N22)*(((-1.0)*(LAMB0 - H_VAL3)))) + N220;
  1677. Z1_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  1678. *
  1679. TTL2 = TABLE;
  1680. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1681. TTL2.'S11' = X1_DP;
  1682. TTL2.'S22' = Y1_DP;
  1683. TTL2.'S12' = Z1_DP;
  1684. TTL2.'FCK' = FCBET1;
  1685. TTL2.'FSK_1' = FSACI1;
  1686. TTL2.'FSK_2' = FSACI1;
  1687. TTL2.'RHO_1' = RHOXS;
  1688. TTL2.'RHO_2' = RHOYS;
  1689. TTL2.'GAMMA_C' = GAMC;
  1690. TTL2.'ALPA_C' = ALPC1;
  1691. TTL2.'NU_C' = NUC1;
  1692. TTL2.'GAMMA_S' = GAMS1;
  1693. TTL3 = TABLE;
  1694. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1695. TTL3.'S11' = X1_DN;
  1696. TTL3.'S22' = Y1_DN;
  1697. TTL3.'S12' = Z1_DN;
  1698. TTL3.'FCK' = FCBET1;
  1699. TTL3.'FSK_1' = FSACI1;
  1700. TTL3.'FSK_2' = FSACI1;
  1701. TTL3.'RHO_1' = RHOXS;
  1702. TTL3.'RHO_2' = RHOYS;
  1703. TTL3.'GAMMA_C' = GAMC;
  1704. TTL3.'ALPA_C' = ALPC1;
  1705. TTL3.'NU_C' = NUC1;
  1706. TTL3.'GAMMA_S' = GAMS1;
  1707. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  1708. 'SINON';
  1709. 'SI' (((X1_DN) '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1710. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))))
  1711. 'OU' ((Y1_DN) '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1712. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1))))));
  1713. X1_DP = ((N11)*(((-1.0)*(LAMB0 + H_VAL3)))) + N110;
  1714. Y1_DP = ((N22)*(((-1.0)*(LAMB0 + H_VAL3)))) + N220;
  1715. Z1_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  1716. X1_DN = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  1717. Y1_DN = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1718. Z1_DN = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1719. TTL2 = TABLE;
  1720. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1721. TTL2.'S11' = X1_DP;
  1722. TTL2.'S22' = Y1_DP;
  1723. TTL2.'S12' = Z1_DP;
  1724. TTL2.'FCK' = FCBET1;
  1725. TTL2.'FSK_1' = FSACI1;
  1726. TTL2.'FSK_2' = FSACI1;
  1727. TTL2.'RHO_1' = RHOXS;
  1728. TTL2.'RHO_2' = RHOYS;
  1729. TTL2.'GAMMA_C' = GAMC;
  1730. TTL2.'ALPA_C' = ALPC1;
  1731. TTL2.'NU_C' = NUC1;
  1732. TTL2.'GAMMA_S' = GAMS1;
  1733. TTL3 = TABLE;
  1734. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1735. TTL3.'S11' = X1_DN;
  1736. TTL3.'S22' = Y1_DN;
  1737. TTL3.'S12' = Z1_DN;
  1738. TTL3.'FCK' = FCBET1;
  1739. TTL3.'FSK_1' = FSACI1;
  1740. TTL3.'FSK_2' = FSACI1;
  1741. TTL3.'RHO_1' = RHOXS;
  1742. TTL3.'RHO_2' = RHOYS;
  1743. TTL3.'GAMMA_C' = GAMC;
  1744. TTL3.'ALPA_C' = ALPC1;
  1745. TTL3.'NU_C' = NUC1;
  1746. TTL3.'GAMMA_S' = GAMS1;
  1747. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  1748. 'SINON';
  1749. TTL2 = TABLE;
  1750. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1751. TTL2.'S11' = X1_DP;
  1752. TTL2.'S22' = Y1_DP;
  1753. TTL2.'S12' = Z1_DP;
  1754. TTL2.'FCK' = FCBET1;
  1755. TTL2.'FSK_1' = FSACI1;
  1756. TTL2.'FSK_2' = FSACI1;
  1757. TTL2.'RHO_1' = RHOXS;
  1758. TTL2.'RHO_2' = RHOYS;
  1759. TTL2.'GAMMA_C' = GAMC;
  1760. TTL2.'ALPA_C' = ALPC1;
  1761. TTL2.'NU_C' = NUC1;
  1762. TTL2.'GAMMA_S' = GAMS1;
  1763. TTL3 = TABLE;
  1764. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1765. TTL3.'S11' = X1_DN;
  1766. TTL3.'S22' = Y1_DN;
  1767. TTL3.'S12' = Z1_DN;
  1768. TTL3.'FCK' = FCBET1;
  1769. TTL3.'FSK_1' = FSACI1;
  1770. TTL3.'FSK_2' = FSACI1;
  1771. TTL3.'RHO_1' = RHOXS;
  1772. TTL3.'RHO_2' = RHOYS;
  1773. TTL3.'GAMMA_C' = GAMC;
  1774. TTL3.'ALPA_C' = ALPC1;
  1775. TTL3.'NU_C' = NUC1;
  1776. TTL3.'GAMMA_S' = GAMS1;
  1777. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL3));
  1778. 'FINSI';
  1779. 'FINSI';
  1780. 'FINSI';
  1781. * computation of D (gradiant) in case of Point_2=Argmin(fmin)
  1782. 'SI' (II1 'EGA' 2);
  1783. X2_DP = ((N11)*(((1.0)*(LAMB0 + H_VAL3)))) + N110;
  1784. Y2_DP = ((N22)*(((-1.0)*(LAMB0 + H_VAL3)))) + N220;
  1785. Z2_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  1786. X2_DN = ((N11)*(((1.0)*(LAMB0 - H_VAL3)))) + N110;
  1787. Y2_DN = ((N22)*(((-1.0)*(LAMB0 - H_VAL3)))) + N220;
  1788. Z2_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  1789. 'SI' (((X2_DP) '>' ((FSACI1*RHOXS)/(GAMS1)))
  1790. 'OU' ((Y2_DP) '>' ((FSACI1*RHOYS)/(GAMS1))));
  1791. X2_DP = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1792. Y2_DP = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1793. Z2_DP = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1794. X2_DN = ((N11)*(((1.0)*(LAMB0 - H_VAL3)))) + N110;
  1795. Y2_DN = ((N22)*(((-1.0)*(LAMB0 - H_VAL3)))) + N220;
  1796. Z2_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  1797. *
  1798. TTL2 = TABLE;
  1799. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1800. TTL2.'S11' = X2_DP;
  1801. TTL2.'S22' = Y2_DP;
  1802. TTL2.'S12' = Z2_DP;
  1803. TTL2.'FCK' = FCBET1;
  1804. TTL2.'FSK_1' = FSACI1;
  1805. TTL2.'FSK_2' = FSACI1;
  1806. TTL2.'RHO_1' = RHOXS;
  1807. TTL2.'RHO_2' = RHOYS;
  1808. TTL2.'GAMMA_C' = GAMC;
  1809. TTL2.'ALPA_C' = ALPC1;
  1810. TTL2.'NU_C' = NUC1;
  1811. TTL2.'GAMMA_S' = GAMS1;
  1812. TTL3 = TABLE;
  1813. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1814. TTL3.'S11' = X2_DN;
  1815. TTL3.'S22' = Y2_DN;
  1816. TTL3.'S12' = Z2_DN;
  1817. TTL3.'FCK' = FCBET1;
  1818. TTL3.'FSK_1' = FSACI1;
  1819. TTL3.'FSK_2' = FSACI1;
  1820. TTL3.'RHO_1' = RHOXS;
  1821. TTL3.'RHO_2' = RHOYS;
  1822. TTL3.'GAMMA_C' = GAMC;
  1823. TTL3.'ALPA_C' = ALPC1;
  1824. TTL3.'NU_C' = NUC1;
  1825. TTL3.'GAMMA_S' = GAMS1;
  1826. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  1827. 'SINON';
  1828. 'SI' (((X2_DN) '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1829. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))))
  1830. 'OU' ((Y2_DN) '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1831. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1))))));
  1832. X2_DP = ((N11)*(((1.0)*(LAMB0 + H_VAL3)))) + N110;
  1833. Y2_DP = ((N22)*(((-1.0)*(LAMB0 + H_VAL3)))) + N220;
  1834. Z2_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  1835. X2_DN = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1836. Y2_DN = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  1837. Z2_DN = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1838. *
  1839. TTL2 = TABLE;
  1840. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1841. TTL2.'S11' = X2_DP;
  1842. TTL2.'S22' = Y2_DP;
  1843. TTL2.'S12' = Z2_DP;
  1844. TTL2.'FCK' = FCBET1;
  1845. TTL2.'FSK_1' = FSACI1;
  1846. TTL2.'FSK_2' = FSACI1;
  1847. TTL2.'RHO_1' = RHOXS;
  1848. TTL2.'RHO_2' = RHOYS;
  1849. TTL2.'GAMMA_C' = GAMC;
  1850. TTL2.'ALPA_C' = ALPC1;
  1851. TTL2.'NU_C' = NUC1;
  1852. TTL2.'GAMMA_S' = GAMS1;
  1853. TTL3 = TABLE;
  1854. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1855. TTL3.'S11' = X2_DN;
  1856. TTL3.'S22' = Y2_DN;
  1857. TTL3.'S12' = Z2_DN;
  1858. TTL3.'FCK' = FCBET1;
  1859. TTL3.'FSK_1' = FSACI1;
  1860. TTL3.'FSK_2' = FSACI1;
  1861. TTL3.'RHO_1' = RHOXS;
  1862. TTL3.'RHO_2' = RHOYS;
  1863. TTL3.'GAMMA_C' = GAMC;
  1864. TTL3.'ALPA_C' = ALPC1;
  1865. TTL3.'NU_C' = NUC1;
  1866. TTL3.'GAMMA_S' = GAMS1;
  1867. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  1868. 'SINON';
  1869. TTL2 = TABLE;
  1870. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1871. TTL2.'S11' = X2_DP;
  1872. TTL2.'S22' = Y2_DP;
  1873. TTL2.'S12' = Z2_DP;
  1874. TTL2.'FCK' = FCBET1;
  1875. TTL2.'FSK_1' = FSACI1;
  1876. TTL2.'FSK_2' = FSACI1;
  1877. TTL2.'RHO_1' = RHOXS;
  1878. TTL2.'RHO_2' = RHOYS;
  1879. TTL2.'GAMMA_C' = GAMC;
  1880. TTL2.'ALPA_C' = ALPC1;
  1881. TTL2.'NU_C' = NUC1;
  1882. TTL2.'GAMMA_S' = GAMS1;
  1883. TTL3 = TABLE;
  1884. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1885. TTL3.'S11' = X2_DN;
  1886. TTL3.'S22' = Y2_DN;
  1887. TTL3.'S12' = Z2_DN;
  1888. TTL3.'FCK' = FCBET1;
  1889. TTL3.'FSK_1' = FSACI1;
  1890. TTL3.'FSK_2' = FSACI1;
  1891. TTL3.'RHO_1' = RHOXS;
  1892. TTL3.'RHO_2' = RHOYS;
  1893. TTL3.'GAMMA_C' = GAMC;
  1894. TTL3.'ALPA_C' = ALPC1;
  1895. TTL3.'NU_C' = NUC1;
  1896. TTL3.'GAMMA_S' = GAMS1;
  1897. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL3));
  1898. 'FINSI';
  1899. 'FINSI';
  1900. 'FINSI';
  1901. * computation of D (gradiant) in case of Point_3=Argmin(fmin)
  1902. 'SI' (II1 'EGA' 3);
  1903. X3_DP = ((N11)*(((1.0)*(LAMB0 + H_VAL3)))) + N110;
  1904. Y3_DP = ((N22)*(((1.0)*(LAMB0 + H_VAL3)))) + N220;
  1905. Z3_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  1906. X3_DN = ((N11)*(((1.0)*(LAMB0 - H_VAL3)))) + N110;
  1907. Y3_DN = ((N22)*(((1.0)*(LAMB0 - H_VAL3)))) + N220;
  1908. Z3_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  1909. 'SI' (((X3_DP) '>' ((FSACI1*RHOXS)/(GAMS1)))
  1910. 'OU' ((Y3_DP) '>' ((FSACI1*RHOYS)/(GAMS1))));
  1911. X3_DP = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1912. Y3_DP = ((N22)*(((1.0)*(LAMB0)))) + N220;
  1913. Z3_DP = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1914. X3_DN = ((N11)*(((1.0)*(LAMB0 - H_VAL3)))) + N110;
  1915. Y3_DN = ((N22)*(((1.0)*(LAMB0 - H_VAL3)))) + N220;
  1916. Z3_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  1917. *
  1918. TTL2 = TABLE;
  1919. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1920. TTL2.'S11' = X3_DP;
  1921. TTL2.'S22' = Y3_DP;
  1922. TTL2.'S12' = Z3_DP;
  1923. TTL2.'FCK' = FCBET1;
  1924. TTL2.'FSK_1' = FSACI1;
  1925. TTL2.'FSK_2' = FSACI1;
  1926. TTL2.'RHO_1' = RHOXS;
  1927. TTL2.'RHO_2' = RHOYS;
  1928. TTL2.'GAMMA_C' = GAMC;
  1929. TTL2.'ALPA_C' = ALPC1;
  1930. TTL2.'NU_C' = NUC1;
  1931. TTL2.'GAMMA_S' = GAMS1;
  1932. TTL3 = TABLE;
  1933. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1934. TTL3.'S11' = X3_DN;
  1935. TTL3.'S22' = Y3_DN;
  1936. TTL3.'S12' = Z3_DN;
  1937. TTL3.'FCK' = FCBET1;
  1938. TTL3.'FSK_1' = FSACI1;
  1939. TTL3.'FSK_2' = FSACI1;
  1940. TTL3.'RHO_1' = RHOXS;
  1941. TTL3.'RHO_2' = RHOYS;
  1942. TTL3.'GAMMA_C' = GAMC;
  1943. TTL3.'ALPA_C' = ALPC1;
  1944. TTL3.'NU_C' = NUC1;
  1945. TTL3.'GAMMA_S' = GAMS1;
  1946. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  1947. 'SINON';
  1948. 'SI' (((X3_DN) '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1949. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))))
  1950. 'OU' ((Y3_DN) '<' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  1951. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1))))));
  1952. X3_DP = ((N11)*(((1.0)*(LAMB0 + H_VAL3)))) + N110;
  1953. Y3_DP = ((N22)*(((1.0)*(LAMB0 + H_VAL3)))) + N220;
  1954. Z3_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  1955. X3_DN = ((N11)*(((1.0)*(LAMB0)))) + N110;
  1956. Y3_DN = ((N22)*(((1.0)*(LAMB0)))) + N220;
  1957. Z3_DN = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  1958. *
  1959. TTL2 = TABLE;
  1960. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1961. TTL2.'S11' = X3_DP;
  1962. TTL2.'S22' = Y3_DP;
  1963. TTL2.'S12' = Z3_DP;
  1964. TTL2.'FCK' = FCBET1;
  1965. TTL2.'FSK_1' = FSACI1;
  1966. TTL2.'FSK_2' = FSACI1;
  1967. TTL2.'RHO_1' = RHOXS;
  1968. TTL2.'RHO_2' = RHOYS;
  1969. TTL2.'GAMMA_C' = GAMC;
  1970. TTL2.'ALPA_C' = ALPC1;
  1971. TTL2.'NU_C' = NUC1;
  1972. TTL2.'GAMMA_S' = GAMS1;
  1973. TTL3 = TABLE;
  1974. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1975. TTL3.'S11' = X3_DN;
  1976. TTL3.'S22' = Y3_DN;
  1977. TTL3.'S12' = Z3_DN;
  1978. TTL3.'FCK' = FCBET1;
  1979. TTL3.'FSK_1' = FSACI1;
  1980. TTL3.'FSK_2' = FSACI1;
  1981. TTL3.'RHO_1' = RHOXS;
  1982. TTL3.'RHO_2' = RHOYS;
  1983. TTL3.'GAMMA_C' = GAMC;
  1984. TTL3.'ALPA_C' = ALPC1;
  1985. TTL3.'NU_C' = NUC1;
  1986. TTL3.'GAMMA_S' = GAMS1;
  1987. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  1988. 'SINON';
  1989. TTL2 = TABLE;
  1990. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  1991. TTL2.'S11' = X3_DP;
  1992. TTL2.'S22' = Y3_DP;
  1993. TTL2.'S12' = Z3_DP;
  1994. TTL2.'FCK' = FCBET1;
  1995. TTL2.'FSK_1' = FSACI1;
  1996. TTL2.'FSK_2' = FSACI1;
  1997. TTL2.'RHO_1' = RHOXS;
  1998. TTL2.'RHO_2' = RHOYS;
  1999. TTL2.'GAMMA_C' = GAMC;
  2000. TTL2.'ALPA_C' = ALPC1;
  2001. TTL2.'NU_C' = NUC1;
  2002. TTL2.'GAMMA_S' = GAMS1;
  2003. TTL3 = TABLE;
  2004. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2005. TTL3.'S11' = X3_DN;
  2006. TTL3.'S22' = Y3_DN;
  2007. TTL3.'S12' = Z3_DN;
  2008. TTL3.'FCK' = FCBET1;
  2009. TTL3.'FSK_1' = FSACI1;
  2010. TTL3.'FSK_2' = FSACI1;
  2011. TTL3.'RHO_1' = RHOXS;
  2012. TTL3.'RHO_2' = RHOYS;
  2013. TTL3.'GAMMA_C' = GAMC;
  2014. TTL3.'ALPA_C' = ALPC1;
  2015. TTL3.'NU_C' = NUC1;
  2016. TTL3.'GAMMA_S' = GAMS1;
  2017. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL3));
  2018. 'FINSI';
  2019. 'FINSI';
  2020. 'FINSI';
  2021. * computation of D (gradiant) in case of Point_4=Argmin(fmin)
  2022. 'SI' (II1 'EGA' 4);
  2023. X4_DP = ((N11)*(((-1.0)*(LAMB0 + H_VAL3)))) + N110;
  2024. Y4_DP = ((N22)*(((1.0)*(LAMB0 + H_VAL3)))) + N220;
  2025. Z4_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  2026. X4_DN = ((N11)*(((-1.0)*(LAMB0 - H_VAL3)))) + N110;
  2027. Y4_DN = ((N22)*(((1.0)*(LAMB0 - H_VAL3)))) + N220;
  2028. Z4_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  2029. 'SI' (((X4_DP) '>' ((FSACI1*RHOXS)/(GAMS1)))
  2030. 'OU' ((Y4_DP) '>' ((FSACI1*RHOYS)/(GAMS1))));
  2031. X4_DP = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  2032. Y4_DP = ((N22)*(((1.0)*(LAMB0)))) + N220;
  2033. Z4_DP = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  2034. X4_DN = ((N11)*(((-1.0)*(LAMB0 - H_VAL3)))) + N110;
  2035. Y4_DN = ((N22)*(((1.0)*(LAMB0 - H_VAL3)))) + N220;
  2036. Z4_DN = ((N12)*(((-1.0)*(LAMB0 - H_VAL3)))) + N120;
  2037. *
  2038. TTL2 = TABLE;
  2039. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2040. TTL2.'S11' = X4_DP;
  2041. TTL2.'S22' = Y4_DP;
  2042. TTL2.'S12' = Z4_DP;
  2043. TTL2.'FCK' = FCBET1;
  2044. TTL2.'FSK_1' = FSACI1;
  2045. TTL2.'FSK_2' = FSACI1;
  2046. TTL2.'RHO_1' = RHOXS;
  2047. TTL2.'RHO_2' = RHOYS;
  2048. TTL2.'GAMMA_C' = GAMC;
  2049. TTL2.'ALPA_C' = ALPC1;
  2050. TTL2.'NU_C' = NUC1;
  2051. TTL2.'GAMMA_S' = GAMS1;
  2052. TTL3 = TABLE;
  2053. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2054. TTL3.'S11' = X4_DN;
  2055. TTL3.'S22' = Y4_DN;
  2056. TTL3.'S12' = Z4_DN;
  2057. TTL3.'FCK' = FCBET1;
  2058. TTL3.'FSK_1' = FSACI1;
  2059. TTL3.'FSK_2' = FSACI1;
  2060. TTL3.'RHO_1' = RHOXS;
  2061. TTL3.'RHO_2' = RHOYS;
  2062. TTL3.'GAMMA_C' = GAMC;
  2063. TTL3.'ALPA_C' = ALPC1;
  2064. TTL3.'NU_C' = NUC1;
  2065. TTL3.'GAMMA_S' = GAMS1;
  2066. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  2067. 'SINON';
  2068. 'SI' (((X4_DN) '>' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  2069. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))))
  2070. 'OU' ((Y4_DN) '>' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  2071. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1))))));
  2072. X4_DP = ((N11)*(((-1.0)*(LAMB0 + H_VAL3)))) + N110;
  2073. Y4_DP = ((N22)*(((1.0)*(LAMB0 + H_VAL3)))) + N220;
  2074. Z4_DP = ((N12)*(((-1.0)*(LAMB0 + H_VAL3)))) + N120;
  2075. X4_DN = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  2076. Y4_DN = ((N22)*(((1.0)*(LAMB0)))) + N220;
  2077. Z4_DN = ((N12)*(((-1.0)*(LAMB0)))) + N120;
  2078. *
  2079. TTL2 = TABLE;
  2080. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2081. TTL2.'S11' = X4_DP;
  2082. TTL2.'S22' = Y4_DP;
  2083. TTL2.'S12' = Z4_DP;
  2084. TTL2.'FCK' = FCBET1;
  2085. TTL2.'FSK_1' = FSACI1;
  2086. TTL2.'FSK_2' = FSACI1;
  2087. TTL2.'RHO_1' = RHOXS;
  2088. TTL2.'RHO_2' = RHOYS;
  2089. TTL2.'GAMMA_C' = GAMC;
  2090. TTL2.'ALPA_C' = ALPC1;
  2091. TTL2.'NU_C' = NUC1;
  2092. TTL2.'GAMMA_S' = GAMS1;
  2093. TTL3 = TABLE;
  2094. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2095. TTL3.'S11' = X4_DN;
  2096. TTL3.'S22' = Y4_DN;
  2097. TTL3.'S12' = Z4_DN;
  2098. TTL3.'FCK' = FCBET1;
  2099. TTL3.'FSK_1' = FSACI1;
  2100. TTL3.'FSK_2' = FSACI1;
  2101. TTL3.'RHO_1' = RHOXS;
  2102. TTL3.'RHO_2' = RHOYS;
  2103. TTL3.'GAMMA_C' = GAMC;
  2104. TTL3.'ALPA_C' = ALPC1;
  2105. TTL3.'NU_C' = NUC1;
  2106. TTL3.'GAMMA_S' = GAMS1;
  2107. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  2108. 'SINON';
  2109. TTL2 = TABLE;
  2110. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2111. TTL2.'S11' = X4_DP;
  2112. TTL2.'S22' = Y4_DP;
  2113. TTL2.'S12' = Z4_DP;
  2114. TTL2.'FCK' = FCBET1;
  2115. TTL2.'FSK_1' = FSACI1;
  2116. TTL2.'FSK_2' = FSACI1;
  2117. TTL2.'RHO_1' = RHOXS;
  2118. TTL2.'RHO_2' = RHOYS;
  2119. TTL2.'GAMMA_C' = GAMC;
  2120. TTL2.'ALPA_C' = ALPC1;
  2121. TTL2.'NU_C' = NUC1;
  2122. TTL2.'GAMMA_S' = GAMS1;
  2123. TTL3 = TABLE;
  2124. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2125. TTL3.'S11' = X4_DN;
  2126. TTL3.'S22' = Y4_DN;
  2127. TTL3.'S12' = Z4_DN;
  2128. TTL3.'FCK' = FCBET1;
  2129. TTL3.'FSK_1' = FSACI1;
  2130. TTL3.'FSK_2' = FSACI1;
  2131. TTL3.'RHO_1' = RHOXS;
  2132. TTL3.'RHO_2' = RHOYS;
  2133. TTL3.'GAMMA_C' = GAMC;
  2134. TTL3.'ALPA_C' = ALPC1;
  2135. TTL3.'NU_C' = NUC1;
  2136. TTL3.'GAMMA_S' = GAMS1;
  2137. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL3));
  2138. 'FINSI';
  2139. 'FINSI';
  2140. 'FINSI';
  2141. * computation of D (gradiant) in case of Point_5=Argmin(fmin)
  2142. 'SI' (II1 'EGA' 5);
  2143. X5_DP = ((N11)*(((-1.0)*(LAMB0) + H_VAL3))) + N110;
  2144. Y5_DP = ((N22)*(((-1.0)*(LAMB0) + H_VAL3))) + N220;
  2145. Z5_DP = ((N12)*(((1.0)*(LAMB0) + H_VAL3))) + N120;
  2146. X5_DN = ((N11)*(((-1.0)*(LAMB0) - H_VAL3))) + N110;
  2147. Y5_DN = ((N22)*(((-1.0)*(LAMB0) - H_VAL3))) + N220;
  2148. Z5_DN = ((N12)*(((1.0)*(LAMB0) - H_VAL3))) + N120;
  2149. 'SI' (((X5_DP) '>' ((FSACI1*RHOXS)/(GAMS1)))
  2150. 'OU' ((Y5_DP) '>' ((FSACI1*RHOYS)/(GAMS1))));
  2151. X5_DP = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  2152. Y5_DP = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  2153. Z5_DP = ((N12)*(((1.0)*(LAMB0)))) + N120;
  2154. X5_DN = ((N11)*(((-1.0)*(LAMB0 - H_VAL3)))) + N110;
  2155. Y5_DN = ((N22)*(((-1.0)*(LAMB0 - H_VAL3)))) + N220;
  2156. Z5_DN = ((N12)*(((1.0)*(LAMB0 - H_VAL3)))) + N120;
  2157. *
  2158. TTL2 = TABLE;
  2159. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2160. TTL2.'S11' = X5_DP;
  2161. TTL2.'S22' = Y5_DP;
  2162. TTL2.'S12' = Z5_DP;
  2163. TTL2.'FCK' = FCBET1;
  2164. TTL2.'FSK_1' = FSACI1;
  2165. TTL2.'FSK_2' = FSACI1;
  2166. TTL2.'RHO_1' = RHOXS;
  2167. TTL2.'RHO_2' = RHOYS;
  2168. TTL2.'GAMMA_C' = GAMC;
  2169. TTL2.'ALPA_C' = ALPC1;
  2170. TTL2.'NU_C' = NUC1;
  2171. TTL2.'GAMMA_S' = GAMS1;
  2172. TTL3 = TABLE;
  2173. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2174. TTL3.'S11' = X5_DN;
  2175. TTL3.'S22' = Y5_DN;
  2176. TTL3.'S12' = Z5_DN;
  2177. TTL3.'FCK' = FCBET1;
  2178. TTL3.'FSK_1' = FSACI1;
  2179. TTL3.'FSK_2' = FSACI1;
  2180. TTL3.'RHO_1' = RHOXS;
  2181. TTL3.'RHO_2' = RHOYS;
  2182. TTL3.'GAMMA_C' = GAMC;
  2183. TTL3.'ALPA_C' = ALPC1;
  2184. TTL3.'NU_C' = NUC1;
  2185. TTL3.'GAMMA_S' = GAMS1;
  2186. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  2187. 'SINON';
  2188. 'SI' (((X5_DN) '>' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  2189. + ((-1.0)*((FSACI1*RHOXS)/(GAMS1)))))
  2190. 'OU' ((Y5_DN) '>' (((-1.0)*(((FCBET1*ALPC1)*(NUC1))/(ALPC1)))
  2191. + ((-1.0)*((FSACI1*RHOYS)/(GAMS1))))));
  2192. X5_DP = ((N11)*(((-1.0)*(LAMB0 + H_VAL3)))) + N110;
  2193. Y5_DP = ((N22)*(((-1.0)*(LAMB0 + H_VAL3)))) + N220;
  2194. Z5_DP = ((N12)*(((1.0)*(LAMB0 + H_VAL3)))) + N120;
  2195. X5_DN = ((N11)*(((-1.0)*(LAMB0)))) + N110;
  2196. Y5_DN = ((N22)*(((-1.0)*(LAMB0)))) + N220;
  2197. Z5_DN = ((N12)*(((1.0)*(LAMB0)))) + N120;
  2198. *
  2199. TTL2 = TABLE;
  2200. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2201. TTL2.'S11' = X5_DP;
  2202. TTL2.'S22' = Y5_DP;
  2203. TTL2.'S12' = Z5_DP;
  2204. TTL2.'FCK' = FCBET1;
  2205. TTL2.'FSK_1' = FSACI1;
  2206. TTL2.'FSK_2' = FSACI1;
  2207. TTL2.'RHO_1' = RHOXS;
  2208. TTL2.'RHO_2' = RHOYS;
  2209. TTL2.'GAMMA_C' = GAMC;
  2210. TTL2.'ALPA_C' = ALPC1;
  2211. TTL2.'NU_C' = NUC1;
  2212. TTL2.'GAMMA_S' = GAMS1;
  2213. TTL3 = TABLE;
  2214. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2215. TTL3.'S11' = X5_DN;
  2216. TTL3.'S22' = Y5_DN;
  2217. TTL3.'S12' = Z5_DN;
  2218. TTL3.'FCK' = FCBET1;
  2219. TTL3.'FSK_1' = FSACI1;
  2220. TTL3.'FSK_2' = FSACI1;
  2221. TTL3.'RHO_1' = RHOXS;
  2222. TTL3.'RHO_2' = RHOYS;
  2223. TTL3.'GAMMA_C' = GAMC;
  2224. TTL3.'ALPA_C' = ALPC1;
  2225. TTL3.'NU_C' = NUC1;
  2226. TTL3.'GAMMA_S' = GAMS1;
  2227. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((1.0)*(H_VAL3));
  2228. 'SINON';
  2229. TTL2 = TABLE;
  2230. TTL2.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2231. TTL2.'S11' = X5_DP;
  2232. TTL2.'S22' = Y5_DP;
  2233. TTL2.'S12' = Z5_DP;
  2234. TTL2.'FCK' = FCBET1;
  2235. TTL2.'FSK_1' = FSACI1;
  2236. TTL2.'FSK_2' = FSACI1;
  2237. TTL2.'RHO_1' = RHOXS;
  2238. TTL2.'RHO_2' = RHOYS;
  2239. TTL2.'GAMMA_C' = GAMC;
  2240. TTL2.'ALPA_C' = ALPC1;
  2241. TTL2.'NU_C' = NUC1;
  2242. TTL2.'GAMMA_S' = GAMS1;
  2243. TTL3 = TABLE;
  2244. TTL3.'TYPE_COUCHE' = 'CHAINE' 'COUCHE_EXT';
  2245. TTL3.'S11' = X5_DN;
  2246. TTL3.'S22' = Y5_DN;
  2247. TTL3.'S12' = Z5_DN;
  2248. TTL3.'FCK' = FCBET1;
  2249. TTL3.'FSK_1' = FSACI1;
  2250. TTL3.'FSK_2' = FSACI1;
  2251. TTL3.'RHO_1' = RHOXS;
  2252. TTL3.'RHO_2' = RHOYS;
  2253. TTL3.'GAMMA_C' = GAMC;
  2254. TTL3.'ALPA_C' = ALPC1;
  2255. TTL3.'NU_C' = NUC1;
  2256. TTL3.'GAMMA_S' = GAMS1;
  2257. D_VAL1 = ((G_ULTI2D TTL2) - (G_ULTI2D TTL3))/((2.0)*(H_VAL3));
  2258. 'FINSI';
  2259. 'FINSI';
  2260. 'FINSI';
  2261. * computation of D (gradiant) in case of Point_6=Argmin(fmin)
  2262. 'SI' (II1 'EGA' 6);
  2263. X6_DP = ((N11)*(((1.0)*