Télécharger @ortho.procedur

Retour à la liste

Numérotation des lignes :

  1. * @ORTHO PROCEDUR FANDEUR 14/10/10 21:15:01 8178
  2. *-------------------------------------------------
  3. ****** PROCEDURE @ORTHO ******
  4. *-------------------------------------------------
  5. *
  6. * CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
  7. * A DISPOSOTION DE LA COMMUNAUTE CASTEM2000
  8. * PAR M. LIBEYRE ( CEA/DSM/DRFC )
  9. *
  10. * TEL : ( 33 1 ) 42 25 46 03
  11. *
  12. * PROCEDURE MISE A JOUR PAR OF LE 09/03/2009
  13. *-------------------------------------------------
  14. DEBPROC @ORTHO MAIL1*'MAILLAGE' LIG1*'MAILLAGE'
  15. CH1*'CHPOINT' CH2*'CHPOINT'
  16. ALPHA*'FLOTTANT' LL*'LISTREEL'
  17. TYPEMOD/'MOT' TYPEELT/'MOT' PTGENE/'POINT' ;
  18.  
  19. 'SI' ('<' ('DIMENSION' LL) 13) ;
  20. 'ERREUR' 'La liste des caracteristiques doit avoir 13 valeurs' ;
  21. 'FINSI' ;
  22.  
  23. DIMGEO = 'VALEUR' 'DIME' ; MODGEO = 'VALEUR' 'MODE' ;
  24. 'SI' ('NON' ('EXISTE' TYPEMOD)) ;
  25. TYPEMOD = MODGEO ;
  26. 'FINSI' ;
  27. 'SI' ('NON' ('EXISTE' TYPEELT)) ;
  28. TYPEELT = 'MOT' 'COQ3' ;
  29. 'FINSI' ;
  30. Mot_DPGE = 'TEXTE' ' ' ; Pt_DPGE = 'TEXTE' ' ' ;
  31.  
  32. 'SI' (('NEG' TYPEMOD 'COEP') 'ET' ('NEG' TYPEMOD 'COMI') 'ET'
  33. ('NEG' TYPEMOD MODGEO)) ;
  34. 'ERREUR' ('CHAINE' 'Mot cle ' TYPEMOD ' non reconnu') ;
  35. 'FINSI' ;
  36. 'SI' ('EGA' TYPEMOD 'COMI') ;
  37. 'SI' ('NON' ('EXISTE' ('MOTS' 'COQ3' 'DKT ') TYPEELT)) ;
  38. 'ERREUR' 'Element de COque MInce a support TRI3 inconnu' ;
  39. 'FINSI' ;
  40. 'FINSI' ;
  41. 'SI' ('EGA' TYPEMOD 'COEP') ;
  42. 'SI' ('NEG' DIMGEO 3) ;
  43. 'ERREUR' 'Modelisation COEP incompatible avec la DIMEnsion' ;
  44. 'FINSI' ;
  45. 'FINSI' ;
  46.  
  47. 'SI' ('EGA' DIMGEO 1) ;
  48. 'ERREUR' 'Cas DIME=1 non prevu actuellement' ;
  49. 'FINSI' ;
  50. 'SI' ('EGA' DIMGEO 3) ;
  51. 'SI' ('EGA' TYPEMOD 'TRID') ;
  52. lm_geo = 'MOTS' 'TET4' 'PRI6' 'PYR5' 'CUB8'
  53. 'TE10' 'PR15' 'PY14' 'CU20' ;
  54. lm_elt = lm_geo ;
  55. 'FINSI' ;
  56. 'SI' ('EGA' TYPEMOD 'COMI') ;
  57. lm_geo = 'MOTS' 'TRI3' ;
  58. lm_elt = 'MOTS' TYPEELT ;
  59. 'FINSI' ;
  60. 'SI' ('EGA' TYPEMOD 'COEP') ;
  61. lm_geo = 'MOTS' 'TRI3' 'QUA4' 'TRI6' 'QUA8' ;
  62. lm_elt = 'MOTS' 'DST ' 'COQ4' 'COQ6' 'COQ8' ;
  63. 'FINSI' ;
  64. 'FINSI' ;
  65. 'SI' ('EGA' DIMGEO 2) ;
  66. 'SI' ('EGA' MODGEO 'PLANGENE') ;
  67. 'SI' ('NON' ('EXISTE' PTGENE)) ;
  68. 'ERREUR' 'Le point support de la DPGE n est pas defini' ;
  69. 'FINSI' ;
  70. Mot_DPGE = 'MOT' 'DPGE' ; Pt_DPGE = PTGENE ;
  71. 'FINSI' ;
  72. lm_geo = 'MOTS' 'TRI3' 'QUA4' 'TRI6' 'QUA8' ;
  73. lm_elt = lm_geo ;
  74. 'SI' ('EGA' TYPEMOD 'COMI') ;
  75. lm_geo = 'MOTS' 'SEG2' ;
  76. lm_elt = 'MOTS' 'COQ2' ;
  77. 'FINSI' ;
  78. 'FINSI' ;
  79.  
  80. lm_mai = 'ELEM' MAIL1 'TYPE' ;
  81. 'REPETER' Boucle ('DIMENSION' lm_mai) ;
  82. str_z = 'EXTRAIRE' lm_mai &Boucle ;
  83. 'SI' ('NON' ('EXISTE' lm_geo str_z)) ;
  84. 'ERREUR' 'Element geometrique incompatible avec la modelisation' ;
  85. 'FINSI' ;
  86. 'FIN' Boucle ;
  87.  
  88. i_z = 0 ;
  89. 'REPETER' Boucle ('DIMENSION' lm_geo) ;
  90. str_z = 'EXTRAIRE' lm_geo &Boucle ;
  91. 'SI' ('EXISTE' lm_mai str_z) ;
  92. mail_z = MAIL1 'ELEM' str_z ;
  93. str_elt = 'EXTRAIRE' lm_elt &Boucle ;
  94. modl_z = 'MODE' mail_z 'MECANIQUE' 'ELASTIQUE' 'ORTHOTROPE'
  95. str_elt 'CONS' 'MECA' Mot_DPGE Pt_DPGE ;
  96. i_z = i_z + 1 ;
  97. 'SI' ('EGA' i_z 1) ;
  98. MODL1 = modl_z ;
  99. 'SINON' ;
  100. MODL1 = MODL1 'ET' modl_z ;
  101. 'FINSI' ;
  102. 'FINSI' ;
  103. 'FIN' Boucle ;
  104.  
  105. GEO = 'CHANGER' 'POI1' MAIL1 ;
  106. NBPT1 = 'NBNO' GEO ;
  107.  
  108. L1 = 'EXTRAIRE' CH1 'COMP' ;
  109. L2 = 'EXTRAIRE' CH2 'COMP' ;
  110.  
  111. 'SI' (DIMGEO 'EGA' 3) ;
  112. L11 = 'EXTRAIRE' 1 L1 ;
  113. L12 = 'EXTRAIRE' 2 L1 ;
  114. L13 = 'EXTRAIRE' 3 L1 ;
  115. L21 = 'EXTRAIRE' 1 L2 ;
  116. L22 = 'EXTRAIRE' 2 L2 ;
  117. L23 = 'EXTRAIRE' 3 L2 ;
  118. LCOMP_K = 'MOTS' 'KX ' 'KY ' 'KZ ' ;
  119. LCOMP_V1 = 'MOTS' 'V1X ' 'V1Y ' 'V1Z ' ;
  120. LCOMP_V2 = 'MOTS' 'V2X ' 'V2Y ' 'V2Z ' ;
  121. 'REPETER' Boucle NBPT1 ;
  122. POINMAIL = GEO 'POINT' &Boucle ;
  123. POINCOUR = LIG1 'POINT' 'PROCHE' POINMAIL ;
  124. CH1MAIL = 'MANU' 'CHPO' POINMAIL 3
  125. 'KX ' ('EXTRAIRE' CH1 L11 POINCOUR)
  126. 'KY ' ('EXTRAIRE' CH1 L12 POINCOUR)
  127. 'KZ ' ('EXTRAIRE' CH1 L13 POINCOUR)
  128. 'NATURE' 'DIFFUS' ;
  129. CH2MAIL = 'MANU' 'CHPO' POINMAIL 3
  130. 'KX ' ('EXTRAIRE' CH2 L21 POINCOUR)
  131. 'KY ' ('EXTRAIRE' CH2 L22 POINCOUR)
  132. 'KZ ' ('EXTRAIRE' CH2 L23 POINCOUR)
  133. 'NATURE' 'DIFFUS' ;
  134. 'SI' ('EGA' &Boucle 1) ;
  135. CH1GEO = CH1MAIL ;
  136. CH2GEO = CH2MAIL ;
  137. 'SINON' ;
  138. CH1GEO = CH1GEO 'ET' CH1MAIL ;
  139. CH2GEO = CH2GEO 'ET' CH2MAIL ;
  140. 'FINSI' ;
  141. 'FIN' Boucle ;
  142. 'FINSI' ;
  143. 'SI' (DIMGEO 'EGA' 2) ;
  144. L11 = 'EXTRAIRE' 1 L1 ;
  145. L12 = 'EXTRAIRE' 2 L1 ;
  146. L21 = 'EXTRAIRE' 1 L2 ;
  147. L22 = 'EXTRAIRE' 2 L2 ;
  148. LCOMP_K = 'MOTS' 'KX ' 'KY ' ;
  149. LCOMP_V1 = 'MOTS' 'V1X ' 'V1Y ' ;
  150. LCOMP_V2 = 'MOTS' 'V2X ' 'V2Y ' ;
  151. 'REPETER' Boucle NBPT1 ;
  152. POINMAIL = GEO 'POINT' &Boucle ;
  153. POINCOUR = LIG1 'POINT' 'PROCHE' POINMAIL ;
  154. CH1MAIL = 'MANU' 'CHPO' POINMAIL 2
  155. 'KX ' ('EXTRAIRE' CH1 L11 POINCOUR)
  156. 'KY ' ('EXTRAIRE' CH1 L12 POINCOUR)
  157. 'NATURE' 'DIFFUS' ;
  158. CH2MAIL = 'MANU' 'CHPO' POINMAIL 2
  159. 'KX ' ('EXTRAIRE' CH2 L21 POINCOUR)
  160. 'KY ' ('EXTRAIRE' CH2 L22 POINCOUR)
  161. 'NATURE' 'DIFFUS' ;
  162. 'SI' ('EGA' &Boucle 1) ;
  163. CH1GEO = CH1MAIL ;
  164. CH2GEO = CH2MAIL ;
  165. 'SINON' ;
  166. CH1GEO = CH1GEO 'ET' CH1MAIL ;
  167. CH2GEO = CH2GEO 'ET' CH2MAIL ;
  168. 'FINSI' ;
  169. 'FIN' Boucle ;
  170. 'FINSI' ;
  171.  
  172. CHP1 = (('COS' ALPHA)*CH1GEO) + (('SIN' ALPHA)*CH2GEO) ;
  173. CHP2 = (((-1.*('SIN' ALPHA))*CH1GEO)) + (('COS' ALPHA)*CH2GEO) ;
  174.  
  175. CHPV1 = 'NOMC' LCOMP_K LCOMP_V1 CHP1 ;
  176. CHMV1 = 'CHANGER' 'CHAM' CHPV1 MODL1 'RIGIDITE' 'CARACTERISTIQUES' ;
  177. CHPV2 = 'NOMC' LCOMP_K LCOMP_V2 CHP2 ;
  178. CHMV2 = 'CHANGER' 'CHAM' CHPV2 MODL1 'RIGIDITE' 'CARACTERISTIQUES' ;
  179.  
  180. YOUNG1S = 'EXTRAIRE' LL 1 ;
  181. YOUNG2S = 'EXTRAIRE' LL 2 ;
  182. YOUNG3S = 'EXTRAIRE' LL 3 ;
  183. NU1S = 'EXTRAIRE' LL 4 ;
  184. NU2S = 'EXTRAIRE' LL 5 ;
  185. NU3S = 'EXTRAIRE' LL 6 ;
  186. CIS1S = 'EXTRAIRE' LL 7 ;
  187. CIS2S = 'EXTRAIRE' LL 8 ;
  188. CIS3S = 'EXTRAIRE' LL 9 ;
  189. ALPH1S = 'EXTRAIRE' LL 10 ;
  190. ALPH2S = 'EXTRAIRE' LL 11 ;
  191. ALPH3S = 'EXTRAIRE' LL 12 ;
  192. RHOS = 'EXTRAIRE' LL 13 ;
  193.  
  194. YOUNG1 = 'MANU' 'CHML' MODL1 'YG1 ' YOUNG1S
  195. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  196. YOUNG2 = 'MANU' 'CHML' MODL1 'YG2 ' YOUNG2S
  197. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  198. YOUNG3 = 'MANU' 'CHML' MODL1 'YG3' YOUNG3S
  199. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  200. NU1 = 'MANU' 'CHML' MODL1 'NU12' NU1S
  201. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  202. NU2 = 'MANU' 'CHML' MODL1 'NU23' NU2S
  203. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  204. NU3 = 'MANU' 'CHML' MODL1 'NU13' NU3S
  205. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  206. CIS1 = 'MANU' 'CHML' MODL1 'G12 ' CIS1S
  207. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  208. CIS2 = 'MANU' 'CHML' MODL1 'G23 ' CIS2S
  209. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  210. CIS3 = 'MANU' 'CHML' MODL1 'G13 ' CIS3S
  211. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  212. ALPH1 = 'MANU' 'CHML' MODL1 'ALP1' ALPH1S
  213. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  214. ALPH2 = 'MANU' 'CHML' MODL1 'ALP2' ALPH2S
  215. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  216. ALPH3 = 'MANU' 'CHML' MODL1 'ALP3' ALPH3S
  217. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  218. RHO = 'MANU' 'CHML' MODL1 'RHO ' RHOS
  219. 'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
  220.  
  221. 'SI' ('EGA' TYPEMOD 'COMI') ;
  222. MAT11 = YOUNG1 'ET' YOUNG2 'ET' NU1 ;
  223. MAT22 = CIS1 ;
  224. 'SI' ('EGA' DIMGEO 3) ;
  225. MAT33 = ('EXCO' ('MOTS' 'V1X ' 'V1Y ') CHMV1) 'ET' ;
  226. MAT33 = ('EXCO' ('MOTS' 'V2X ' 'V2Y ') CHMV2) ;
  227. 'SINON' ;
  228. MAT33 = CHMV1 'ET' CHMV2 ;
  229. 'FINSI' ;
  230. MAT44 = ALPH1 'ET' ALPH2 'ET' RHO ;
  231. 'FINSI' ;
  232. 'SI' ('EGA' TYPEMOD 'COEP') ;
  233. MAT11 = YOUNG1 'ET' YOUNG2 'ET' NU1 ;
  234. MAT22 = CIS1 'ET' CIS2 'ET' CIS3 ;
  235. 'SI' ('EGA' DIMGEO 3) ;
  236. MAT33 = ('EXCO' ('MOTS' 'V1X ' 'V1Y ') CHMV1) 'ET' ;
  237. MAT33 = ('EXCO' ('MOTS' 'V2X ' 'V2Y ') CHMV2) ;
  238. 'SINON' ;
  239. MAT33 = CHMV1 'ET' CHMV2 ;
  240. 'FINSI' ;
  241. MAT44 = ALPH1 'ET' ALPH2 'ET' RHO ;
  242. 'FINSI' ;
  243. 'SI' ('EGA' TYPEMOD 'PLANCONT') ;
  244. MAT11 = YOUNG1 'ET' YOUNG2 'ET' YOUNG3 'ET' NU1 'ET' NU2 'ET' NU3 ;
  245. MAT22 = CIS1 ;
  246. MAT33 = CHMV1 'ET' CHMV2 ;
  247. MAT44 = ALPH1 'ET' ALPH2 'ET' RHO ;
  248. 'FINSI' ;
  249. 'SI' (('EGA' TYPEMOD 'PLANDEFO') 'OU' ('EGA' TYPEMOD 'PLANGENE') 'OU'
  250. ('EGA' TYPEMOD 'AXIS')) ;
  251. MAT11 = YOUNG1 'ET' YOUNG2 'ET' YOUNG3 'ET' NU1 'ET' NU2 'ET' NU3 ;
  252. MAT22 = CIS1 ;
  253. MAT33 = CHMV1 'ET' CHMV2 ;
  254. MAT44 = ALPH1 'ET' ALPH2 'ET' ALPH3 'ET' RHO ;
  255. 'FINSI' ;
  256. 'SI' (('EGA' TYPEMOD 'FOUR') 'OU' ('EGA' TYPEMOD 'TRID')) ;
  257. MAT11 = YOUNG1 'ET' YOUNG2 'ET' YOUNG3 'ET' NU1 'ET' NU2 'ET' NU3 ;
  258. MAT22 = CIS1 'ET' CIS2 'ET' CIS3 ;
  259. MAT33 = CHMV1 'ET' CHMV2 ;
  260. MAT44 = ALPH1 'ET' ALPH2 'ET' ALPH3 'ET' RHO ;
  261. 'FINSI' ;
  262.  
  263. MATGEO = MAT11 'ET' MAT22 'ET' MAT33 'ET' MAT44 ;
  264.  
  265. 'FINPROC' MATGEO MODL1 ;
  266.  
  267.  
  268.  
  269.  

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