Télécharger dependance.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dependance.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************
  5. * TEST DE CONDENSATION
  6. *
  7. * poutre des section rectangulaire en appui simples
  8. * face inferieure renforcee par une plaque en acier
  9. * un renfort interne de type poutre
  10. * un renfort interne de type barr
  11. *
  12. * les renforts sont disjoints du maillage massif
  13. *
  14. * comparaison des solutions obtenues par :
  15.  
  16. * rela accro / resou
  17. * dependance / resou
  18. * cmct / resou + mrem
  19. *
  20. * test de fonctionnement dans pasapas
  21. *
  22. *********************************************************************
  23. opti echo 1;
  24. graph='N' ;
  25. opti dime 3 elem cub8 ;
  26. *
  27. * GEOMETRIE
  28. *
  29. long = 10.0;
  30. *---BETON
  31. h = 1.;
  32. h1= .9 ;
  33. *h1 = 1. - (1./6.1);
  34. pa=0.0 0. 0.0;
  35. pb=long 0. 0.0;
  36. pc=long 0. h ;
  37. pd=0.0 0. h;pn = 0. 0. (h - h1) ;
  38.  
  39. *---ARMATURE
  40. prof = 1. ;
  41. pm = long (prof/2) (h - h1) ;
  42. pn = 0. (prof/2) (h - h1) ;
  43. nb1 = 20;
  44. nb2 = 4 ;
  45. nb3=17 ;
  46. c1=pa droi nb1 pb;
  47. c2=pb droi nb2 pc;
  48. c3=pc droi nb1 pd;
  49. c4=pd droi nb2 pa;
  50. sur1=(dall c1 c2 c3 c4) ;
  51.  
  52.  
  53. sur1 = sur1 volu 2 tran (0 prof 0) ;
  54. d1 = d nb3 pm pn coul roug ;
  55. d2 = d1 plus ( 0 0 h1 ) ;
  56. finf = sur1 points plan pa pb (pa plus (0 1 0)) .01 ;
  57. finf = (( sur1 envel) elem appu stric finf ) ;
  58. fsup = finf plus (0 0 h) ;
  59. elim .1 sur1 fsup ;
  60. *----- on decolle la plaque pour pouvoir l ACCROCHER
  61. plaq = finf plus ( 0 0 0 ) coul bleu ;
  62.  
  63. 'SI' ( 'NEG' graph 'N' ) ;
  64. trac (1000 1000 1000 ) ( sur1 et d1) ;
  65. 'FINSI' ;
  66.  
  67. *
  68. * MODELE
  69. *------------------------- beton ---------------------------------
  70. Eb = 0.4e11 ;
  71. modb= 'MODE' sur1 mecanique elastique ;
  72. matb= 'MATE' modb youn Eb nu 0.3;
  73. ribe = rigi modb matb;
  74.  
  75. *----------------------- barre ----------------------------
  76. YOUCAB = 1495.D6/0.718D-2 ; ;
  77. leps1 = 1.D-2 * (PROG 0. .718 .804 .890 1.00 1.278 2.957) ;
  78. lsig1 = 1.D6 * (PROG 0. 1495. 1553. 1623. 1669. 1727. 1855.) ;
  79. TRACAB = EVOL MANU 'EPS' leps1 'SIG' lsig1 ;
  80. ym1 = (extr lsig1 2) / (extr leps1 2) ;
  81. lsig1 = lsig1 enle 1 ;
  82. lepsp1 = (leps1 enle 1) - (lsig1 / ym1) ;
  83. trecro = evol vert manu 'EPS' lepsp1 'SIG' lsig1 ;
  84. 'SI' ( 'NEG' graph 'N' ) ;
  85. dess (TRACAB et trecro) titr ' Courbes de traction et d ecrouissage (vert)' ;
  86. fins ;
  87. mod1 = MODE d1 MECANIQUE ELASTIQUE ISOTROPE PLASTIQUE ISOTROPE BARR ;
  88. mat1= MATE mod1 YOUN youcab NU 0.3 SECT 0.04 'ECRO' trecro ALPH 12.5E-6 ;
  89. riba = rigi mod1 mat1 ;
  90.  
  91. *------------------- - poutre -----------------------
  92.  
  93. mod2= 'MODE' d2 mecanique elastique pout ;
  94. mat2=( 'MATE' mod2 youn Eb nu 0.0 INRY 10. INRZ 10. TORS 3.)
  95. et ( cara mod2 sect .02) ;
  96. ripo = rigi mod2 mat2 ;
  97.  
  98. *--------------------- coque -------------------------
  99.  
  100. *---- Courbe d'écrouissage de l'acier A42
  101. EPS1 = PROG 0. 0.002 0.018 0.02 0.036 0.05 0.068 0.084 0.102 ;
  102. SIG1 = PROG 0. 312. 325. 335. 400. 425. 450. 465. 475. ;
  103. SIG1 = SIG1*1.E6 ;
  104. EVSIG = EVOL MANU 'EPS' EPS1 'SIG' SIG1 ;
  105. ym1 = (extr sig1 2) / (extr eps1 2) ;
  106. sig1 = sig1 enle 1 ;
  107. eps1 = (eps1 enle 1) - (sig1 / ym1) ;
  108. EVECR = EVOL vert MANU 'EPS' EPS1 'SIG' SIG1 ;
  109. 'SI' ( 'NEG' graph 'N' ) ;
  110. dess (EVSIG et EVECR) titr ' Courbes de traction et d ecrouissage (vert)' ;
  111. fins ;
  112. *---- Module de Young de la Peau d'étanchéité
  113. YOUNT = 312.E6 / 0.002 ;
  114. mod3 = mode plaq MECANIQUE ELASTIQUE ISOTROPE
  115. PLASTIQUE ISOTROPE coq4 ;
  116. Mat3 = MATE mod3 YOUN YOUNT NU 0.3 EPAI 0.03 ecro EVecr ALPH 12.5E-6 ;
  117. rico = rigi mod3 mat3 ;
  118. *----------------------------------------------------------
  119. *
  120. * CONDITIONS AUX LIMITES
  121. *
  122. aa = sur1 points droite pa ( pa plus (0 1 0)) .01 ;
  123. bb = sur1 points droite pb ( pb plus (0 1 0)) .01 ;
  124. cl1 = bloq aa ux ;
  125. cl2 = bloq aa uz ;
  126. cl3 = bloq bb uz ;
  127. cl4 = bloquer uy ( pa et pb) ;
  128. clb = cl1 et cl2 et cl3 et cl4 ;
  129. *
  130. * on veut bloquer la rotation sur elle meme de la poutre
  131. *
  132. dp2 = d2 elem ( lect 1 pas 1 ( nb3 - 1)) ;
  133. pf = d2 point final ;
  134. elim ( sur1 point proche pf ) pf ;
  135. *
  136. clt = clb et ( bloquer RX pf ) ;
  137. *
  138.  
  139. racc = rela accro (plaq et d1 et dp2 ) sur1 ;
  140. radd = chan DEPE racc ;
  141.  
  142. * ------ concatenation des raideurs ------------
  143. rign = ribe et riba et ripo et rico ;
  144. riga = rign et racc et clt ;
  145. rigd = rign et radd et clt ;
  146. *
  147. * --- pression sur face superieure + force sur la plaque -----
  148. MOP = 'MODE' fsup 'CHARGEMENT' 'PRESSION' ;
  149. PRZ = 'PRES' MOP pres 1.e8 ;
  150. f1 = bsig mop PRZ ;
  151. f2 = force 1.E8 FZ plaq ;
  152. ff = f1 + f2 ;
  153. *-------------- solution par depend-------------------------
  154. dep2 = resou rigd ff ;
  155. *-------------- solution par accro -------------------------
  156. dep1 = resou riga ff ;
  157. *--------------------------------------------------------
  158. *------------- solution par cmct puis mrem ------------
  159. * --- condensation des matrices liees aux noeuds esclaves ---
  160. rige = cmct (rign et clt ) radd ;
  161. * --- condensation des la force liees aux noeuds esclaves -----
  162. racd radu = chan COND racc ;
  163. fad = f1 et ( radu * f2 ) ;
  164. * ----- solution reduite ---------------------------------------
  165. * dep30 = resou rige fad ;
  166. *-- remeontee a la solution complete --------------------------
  167. * on n'utilise plus tout cela, c'est intégré dans résou
  168. * dep3 = mrem dep30 (rign et clt et radd ) ff ;
  169. dep3 = resou rige fad ;
  170. *--------------------------------------------------------------
  171. mde1 = mini dep1 avec ( mots UZ) ;
  172. mde2 = mini dep2 avec ( mots UZ) ;
  173. mde3 = mini dep3 avec ( mots UZ) ;
  174. surd = (sur1 envel ) et d1 et d2 ;
  175.  
  176. * ----------- reactions aux appuis et liaisons ---------------------
  177. rea1 = reac (clt et racc ) dep1 ;
  178. rea2 = reac (clt et radd ) dep2 ;
  179. *----------examen des differences entre les deux methodes ----------
  180. aaa = rea1 - rea2 ;
  181. reamo1 = (psca rea1 rea1 ( mots FX FY FZ) ( mots FX FY FZ)) ** .5 ;
  182. reamo2 = (psca rea2 rea2 ( mots FX FY FZ) ( mots FX FY FZ)) ** .5 ;
  183. depmod = (psca dep1 dep1 ( mots UX UY UZ) ( mots UX UY UZ)) ** .5 ;
  184. si ( neg GRAPH N) ;
  185. av = maxi ( abs rea1 ) avec (mots FZ) ;
  186. amv = 4./av ;
  187. vr1 = vecteur rea1 amv FX FY FZ roug ;
  188. vr2 = vecteur rea2 amv FX FY FZ vert ;
  189. amp = .3/ ( maxi depmod) ;
  190. amd = 4./(maxi reamo1) ;
  191. vd1 = vecteur rea1 amd FX FY FZ roug ;
  192. vd2 = vecteur rea2 amd FX FY FZ vert ;
  193. titre 'UZ mini accro ' mde1 ' dependence ' mde2 'cmct' mde3 ;
  194. dde0 = defo surd dep1 0. blan ;
  195. dde1 = defo surd dep1 amp vr1 roug ;
  196. dde2 = defo surd dep2 amp vr2 vert ;
  197. trac ( dde1 et dde2 ) ;
  198. dur1 = defo surd dep1 amp vd1 roug ;
  199. dur2 = defo surd dep2 amp vd2 vert ;
  200. trac ( dur1 et dur2 ) ;
  201. vdif = vecteur aaa FX FY FZ amd vert ;
  202. titre ' difference sur reactions accro depen ' ;
  203. trac surd vdif ;
  204. finsi ;
  205. rdif = reamo1 - reamo2 ;
  206. a3= rdif point maxi ;
  207. av1 = (abs (maxi rdif)) /( extr reamo1 SCAL (a3 point 1 )) ;
  208. ruz = abs (((mde1 - mde2) * 100 ) / mde2 ) ;
  209.  
  210. mess ' differences relatives max sur la fleche maximum ' ruz '%' ;
  211. mess ' differences relatives max sur les reactions ' av1 '%' ;
  212. 'SI' (( ruz > 1.e-7) OU ( av1 > 1.e-7));
  213. erre 5 ;
  214. 'FINSI' ;
  215.  
  216. * suite uniquement pour tester que tout se passe correctement ds pasapas
  217.  
  218. evf = EVOL MANU 'Sec.' (prog 0. 1000.) 'Pa' (prog 1. 1. ) ;
  219. chforce = char meca f2 evf ;
  220. *
  221. chpress = char pres prz evf ;
  222.  
  223. TABMC = TABLE ;
  224. TABMC.MODELE = modb et mod1 et mod2 et mod3 et mop ;
  225. TABMC.BLOCAGES_MECANIQUES = radd et clt ;
  226. TABMC.CARACTERISTIQUES = matb et mat1 et mat2 et mat3 ;
  227. TABMC.CHARGEMENT = chpress et chforce ;
  228. evtt = prog 0. pas 1. 5.;
  229. TABMC.TEMPS_CALCULES = evtt ;
  230. TABMC.TEMPS_SAUVES = evtt ;
  231. TABMC.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  232. PASAPAS TABMC ;
  233.  
  234. ttt = tabmc.'DEPLACEMENTS' ;
  235. i1 = index ttt ;nt = dime i1 ;
  236. dep = ttt . (i1.nt ) ;uref = -8.94153E-02 ;
  237. uymi = mini dep avec (mots UZ) ;
  238. * fleche attendue 5 iterations =-8.94153E-02 ;
  239. ruz = abs (((uref - uymi) * 100 ) / (abs uymi )) ;
  240. mess ' difference relative ' ruz '%' ;
  241.  
  242. 'SI' ( ruz > 5.e-1) ;
  243. erre 5 ;
  244. 'FINSI' ;
  245. 'FIN' ;
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  
  261.  
  262.  
  263.  
  264.  
  265.  
  266.  
  267.  
  268.  
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  
  275.  

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