Télécharger vari-youn-2-auto.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : vari-youn-2-auto.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **************************************************************
  5. *
  6. * C'est un k-test dont le but est de tester le fonctionnement
  7. * de PASAPAS avec la variation des données matérielles.
  8. *
  9. * Il s'agit d'une barre en traction simple. Son module de Young
  10. * dépend linéairement du paramètre d'évolution (à la fin du
  11. * chargement on descend à 10% de la valeur initiale).
  12. *
  13. * Pour la résolution, on utilise le mécanisme de traitement
  14. * des matériaux variables implanté dans PASAPAS. La variable
  15. * utilisée par cette mécanique porte le nom LDPL, ce chargement
  16. * n'est donc pas considéré comme un chargement thermique. Le
  17. * pilotage dans PASAPAS est automatique.
  18. *
  19. **************************************************************
  20.  
  21. *** Options ...
  22.  
  23. opti dime 2 mode plan defo elem seg2 echo 1 ;
  24. graph = faux ;
  25.  
  26. *** Paramètres ...
  27.  
  28. L = 10. ;
  29. H = 1. ;
  30. S = H * 1. ;
  31.  
  32. E_0 = 2.E+11 ;
  33. * E_1 = (2.0 * E_0) / (0.01 * L) ;
  34. E_1 = (-0.9 * E_0) / (1.0 * L) ;
  35.  
  36. * F = 1.E+10 ;
  37. F = 9.9E+9 ;
  38.  
  39. nbdiv = 1 ;
  40.  
  41. *** Points ...
  42.  
  43. dens (L / nbdiv) ;
  44.  
  45. p1 = 0 0 ;
  46. p2 = 0 H ;
  47.  
  48. *** Ligne ...
  49.  
  50. li1 = p1 d 1 p2 ;
  51.  
  52. *** Surface ...
  53.  
  54. opti elem qua4 ;
  55.  
  56. vectx = L 0 ;
  57. su1 = li1 tran nbdiv vectx ;
  58. pu = su1 poin proc vectx ;
  59.  
  60. titr 'Le maillage du modele' ;
  61. si(graph) ;
  62. trac su1 ;
  63. finsi ;
  64.  
  65. *** Modèle ...
  66.  
  67. listdepl = prog (-0.9 * L) (0.9 * L) ;
  68. listyoun = prog ;
  69. lon = dime listdepl ;
  70. repeter surdpl lon ;
  71. valdep = extr listdepl &surdpl ;
  72. valyou = E_0 + (E_1 * valdep) ;
  73. listyoun = listyoun et (prog valyou) ;
  74. fin surdpl ;
  75. titr 'Module de Young (deplacement)' ;
  76. evyou = evol manu LDPL listdepl YOUN listyoun ;
  77. si(graph) ;
  78. dess evyou ;
  79. finsi ;
  80.  
  81. mo1 = mode su1 mecanique elastique ;
  82. ma1 = mate mo1 NU 0.0 YOUN evyou ;
  83.  
  84. *** Préparation de la table pour PASAPAS ...
  85.  
  86. tabnl = table ;
  87. * tabnl . GRANDS_DEPLACEMENTS = vrai ;
  88. tabnl . HYPOTHESE_DEFORMATIONS = 'LINEAIRE';
  89. tabnl . MODELE = mo1 ;
  90. tabnl . CARACTERISTIQUES = ma1 ;
  91. cl1 = (bloq ux li1) et (bloq uy p1) ;
  92. tabnl . BLOCAGES_MECANIQUES = cl1 ;
  93. tabnl . PAS_AJUSTE = vrai ;
  94.  
  95. ordmax = 8.0 ;
  96. ldpl0 = manu CHML mo1 LDPL 1.0 STRESSES ;
  97. evconst = evol manu 't' (prog 0. 1.)
  98. 'f(t)' (prog 0. ordmax) ;
  99. chavar0 = char LDPL ldpl0 evconst ;
  100.  
  101. fobid = forc fx F (cote 3 su1) ;
  102. lr1 = prog 0. 1. ;
  103. lr2 = prog 0. 1. ;
  104. evbid = evol manu 't' lr1 'f(t)' lr2 ;
  105. chabid = char fobid evbid MECA ;
  106. tabnl . CHARGEMENT = chavar0 et chabid ;
  107.  
  108. nb_dt = 5 ;
  109. dt = 1.0 / nb_dt ;
  110. tabnl . TEMPS_SAUVES = prog 0. PAS dt 1. ;
  111.  
  112. *** Le calcul ...
  113.  
  114. pasapas tabnl ;
  115.  
  116. *** Petit post-traitement ...
  117.  
  118. listt = prog ;
  119. listux = prog ;
  120. listuxt = prog ;
  121. listfx = prog ;
  122. nbpas = dime (tabnl . TEMPS) ;
  123. repeter surpas nbpas ;
  124. lindice = &surpas - 1 ;
  125. listt = listt et (prog (tabnl . TEMPS . lindice)) ;
  126. laforc = F * (tabnl . TEMPS . lindice) ;
  127. listfx = listfx et (prog laforc) ;
  128. fauxdepl = ipol (tabnl . TEMPS . lindice) (extr evconst
  129. ABSC) (extr evconst ORDO) ;
  130. mess '-------------------------------------' ;
  131. mess 'fauxdepl = ' fauxdepl ;
  132. E_act = E_0 + (E_1 * fauxdepl) ;
  133. mess 'E_act = ' E_act ;
  134. uxbid = (laforc * L) / (E_act * S) ;
  135. listuxt = listuxt et (prog uxbid) ;
  136. depact = tabnl . DEPLACEMENTS . lindice ;
  137. ledepl = extr depact UX pu ;
  138. listux = listux et (prog ledepl) ;
  139. fin surpas ;
  140.  
  141. titr 'La deformee' ;
  142. defo0 = defo su1 depact 0.0 ;
  143. defo1 = defo su1 depact ROUG ;
  144. si(graph) ;
  145. trac (defo0 et defo1) ;
  146. finsi ;
  147.  
  148. titr 'Courbe de traction provenant de PASAPAS' ;
  149. courbetr = evol manu 'FX' listfx 'UX' listux ;
  150. si(graph) ;
  151. dess courbetr ;
  152. finsi ;
  153.  
  154. titr 'Courbes de traction PASAPAS + théorique' ;
  155. courbteo = evol ROUG manu 'FX' listfx 'UX' listuxt ;
  156. si(graph) ;
  157. dess (courbetr et courbteo) ;
  158. finsi ;
  159.  
  160. *** Vérification de la précision du calcul ...
  161.  
  162. chtit = chai 'Difference entre le depl. fourni par PASAPAS'
  163. ' et le depl. analytique' ;
  164. titr chtit ;
  165. ladiff = courbetr - courbteo ;
  166. ladiff = evol manu (extr ladiff ABSC) (abs (extr ladiff ORDO)) ;
  167. si(graph) ;
  168. dess ladiff ;
  169. finsi ;
  170. precis = maxi (extr ladiff ORDO) ;
  171. mess 'Écart maxi = ' precis ;
  172.  
  173. si(precis > 1.e-14) ;
  174. erre 5 ;
  175. finsi ;
  176.  
  177. *** Bye ...
  178.  
  179. fin ;
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  

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