Télécharger la plaquette Cast3M


fr en

Télécharger TirantLAB.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : tirantlab.dgibi
  2.  
  3. ************************************************************************
  4. ************************************************************************
  5. *
  6. * CAS TEST DU 13/11/14     PROVENANCE : TEST                                    
  7. SAUT PAGE;                                                                      
  8. *                                                                               
  9. *******************************************************        
  10. *                                                     *            
  11. *  TEST ELEMENT COAXIAL COS2 POUR MODELE DE LIAISON   * 
  12. *  ACIER-BETONEN REGIME LINEAIRE                      *            
  13. *                                                     *            
  14. *  TIRANT 3 m SECTION CARRE 0.1 * 0.1                 *            
  15. *                                                     *             
  16. *                                                     *            
  17. *      Le tirant est maille en 3D a l'aide d'elements *
  18. *  massif CUB8. L'acier est maille a l'aide d'elements*
  19. *  barre. Le tirant est soumise sous un chargemet     *                       
  20. *  monotone.                                          *           
  21. *******************************************************            
  22. *     
  23.  
  24.  
  25. OPTI ECHO 0; 
  26. ***parametres
  27. xp1 = 3.;
  28. xp2 = xp1 + 0.1;
  29. yp1 = 0.1;
  30. yp2 = yp1/2.;
  31. nelx = 100;
  32. nely = 1;
  33. nelz = 1;
  34. *
  35. opti dime 3 elem cub8;
  36. ***geometrie et maillage
  37. b1 = 0. 0. 0.; 
  38. b2 = xp1 0. 0.;
  39. lib1 = b1 droi nelx  b2 ;
  40. sb1 = lib1 tran nely (0. yp1 0.) ; 
  41. vol1 = sb1 volu tran nelz (0. 0. yp1) ; 
  42. a1 = 0. yp2 yp2; a2 = xp1 yp2 yp2;
  43. PA1 = (-0.1 yp2 yp2) ;
  44. PA2 = (xp2 yp2 yp2) ;
  45. *
  46. lia1 = droi nelx a1 a2 ;
  47. lia2 = lia1 plus (0. 0. 0.);
  48. *
  49. D1 = DROI 1 PA1 A1 ;
  50. D2 = DROI 1 A2 PA2 ;
  51. ACI = D1 et D2 et LIA1 ;
  52. *
  53. acibet = racc lia1 lia2 1.e-4;
  54. vtot = vol1 et acibet;
  55. ***loi d adherence
  56. pulo1 = prog 0. 1.;
  57. pulo2 = prog 0. 1.e11;
  58. *
  59. pulop = evol manu pulo1 pulo2;
  60. *   
  61. ***materiaux
  62. ***acier
  63. yga = 210.e9;
  64. nua = 0.3;
  65. secc1 = 7.854e-5;
  66. secc2 = 7.854e-5;
  67. ***beton
  68. ygb = 30.e9;
  69. nub = 0.2;
  70. ***model
  71. ***acier
  72. moa= model aci mecanique elastique isotrope barre;
  73. maa = mate moa young yga nu nua SECT secc1;
  74. ***beton
  75. mob = model vol1 mecanique elastique isotrope;
  76. mab =MATE mob youn ygb nu nub;
  77. ***beton acier
  78. moab = model acibet mecanique elastique isotrope plastique 
  79. liaison_acbe cos2;
  80. maab = mate moab 'PULO' pulop 'KN' 1e15 'KS' 1.e11 'SECT' secc2 ;
  81. *
  82. motot = moa et mob et moab;
  83. matot = maa et mab et maab;
  84. ***conditions aux limites
  85. cl1 = rela accro lia2 vol1; 
  86. cl2 = bloqu UZ sb1;
  87. cl3 = bloqu depl PA1;
  88. cl4 = bloqu ux PA2;
  89. cl5 = bloqu uz pa2;
  90. cl6 = bloqu uy pa2;
  91. pp1 = 0. 0. 0.;
  92. pp2 = 1. 0. 0.;
  93. pp3 = 1. 0. yp1;
  94. sy1 = vol1 poin plan pp1 pp2 pp3 1.e-4;
  95. cl7 = bloq uy sy1;
  96. *
  97. cltot = cl1 et cl2 et cl3 et cl4 et cl5 et cl6 et cl7;
  98. ***chargement
  99. evo1 = evol manu (prog 0. 1.) (prog 0. 2.e-3);     
  100. for1 = depi cl4 1.;
  101. char1 = chargement (for1 ) evo1 dimp;
  102. ***table
  103. tab1 = table;
  104. tab1.modele = motot;
  105. tab1.caracteristiques = matot;
  106. tab1.blocages_mecaniques = cltot;
  107. tab1.chargement= char1;        
  108. tab1.'PRECISION' = 1.e-5;
  109. tab1.temps_calcules =
  110. prog 0. pas 0.5 1.;
  111. *
  112. pasapas tab1;
  113. *
  114. *********************
  115. ***post traitement***
  116. *********************
  117. ***force 
  118. rr0 = reac cl4 tab1.deplacements.1;
  119. rr1 = resu rr0;
  120. fff = extr rr1 'FX' (poin ( extr rr0 mail ) initial );
  121. *
  122. errf = abs(fff - 33.75e3)/33.75e3;
  123. *
  124. si (> errf 1e-3);
  125. erreur (5);
  126. finsi;
  127.  
  128. fin;
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  

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