Télécharger plas_coufdp.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas_coufdp.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * ESSAIS COUDE EN FLEXION DANS LE PLAN
  6. * COUDE MINCE ( e = 2.15 mm dext = 179.0 mm )
  7. *
  8. OPTI ECHO 0 ;
  9. OPTI DIME 3 ELEM SEG2 ;
  10. OEIL = 100. 0. 0. ;
  11. NELCOU = 10 ;
  12. LH1 = 1.450 ;
  13. LA1 = .179 ;
  14. RCO1 = .275 ;
  15. LARC1 = LA1 + RCO1 ;
  16. LB1 = LH1 - LARC1 ;
  17. *
  18. P1 = 0. 0. LH1 ;
  19. P2 = 0. 0. LARC1 ;
  20. P3 = 0. 0. RCO1 ;
  21. P4 = 0. RCO1 0. ;
  22. P5 = 0. LARC1 0. ;
  23. P6 = 0. ( LARC1 + 50.E-3 ) 0. ;
  24. *
  25. LREN1 = P1 D 5 P2 ;
  26. LDRO1 = P2 D 3 P3 ;
  27. COUD1 = P3 C NELCOU ( 0. RCO1 RCO1 ) P4 ;
  28. LDRO2 = P4 D 3 P5 ;
  29. LREN2 = P5 D 1 P6 ;
  30. *
  31. LIGN1 = LREN1 ET LDRO1 ET COUD1 ET LDRO2 ;
  32. LIGN1 = LIGN1 ET LREN2 ;
  33. ELIM LIGN1 1.E-4 ;
  34. TASS LIGN1 ;
  35. *
  36. * MODELE PLASTIQUE ECROUISSAGE ISOTROPE
  37. *
  38. OBLR1 = MODE ( LREN1 ET LREN2 ) 'MECANIQUE' 'ELASTIQUE'
  39. 'PLASTIQUE' 'TUYA' ;
  40. OBLD1 = MODE LDRO1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'TUYA' ;
  41. OBCO1 = MODE COUD1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'TUYA' ;
  42. OBLD2 = MODE LDRO2 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'TUYA' ;
  43. OBDRO1 = OBLD1 ET OBLD2 ;
  44. modlign = OBLR1 ET OBLD1 ET OBCO1 ET OBLD2 ;
  45. *
  46. REXT1 = 179.0E-3 / 2. ;
  47. ECOU1 = 2.15E-3 ;
  48. EDRO1 = 2.15E-3 ;
  49. ERIG1 = 10.0E-3 ;
  50. CAREN1 = 'CARA' OBLR1 'EPAI' ERIG1 'RAYO' REXT1 'VECT' ( 1. 0. 0. ) ;
  51. CARDRO1 = 'CARA' OBDRO1 'EPAI' EDRO1 'RAYO' REXT1 'VECT' ( 1. 0. 0. ) ;
  52. CARCOU1 = 'CARA' OBCO1 'EPAI' ECOU1 'RAYO' REXT1
  53. 'RACO' RCO1 'VECT' ( 1. 0. 0. ) ;
  54. CARLIGN = CAREN1 ET CARDRO1 ET CARCOU1 ;
  55. *
  56. * COURBE DE TRACTION
  57. *
  58. LSIG1 = 'PROG' 0. 220. 260.0 300.0 330. 360.0 390.0
  59. 420.0 450.0 480. 520.0 620.0 ;
  60. LSIG1 = LSIG1 * 1.E6 ;
  61. LEPS1 = 'PROG' 0. .118 .152 .203 .273 .398 1.044
  62. 4.331 8.55 12.78 20.5 67. ;
  63. LEPS1 = LEPS1 * 1.E-2 ;
  64. *
  65. SSTBG = 'EVOL' 'MANU' 'DEFO' LEPS1 'CONT' LSIG1 ;
  66. *
  67. SIGM1 = 'EXTR' SSTBG 'ORDO' ;
  68. EPSI1 = 'EXTR' SSTBG 'ABSC' ;
  69. YOUN1 = ( 'EXTR' SIGM1 2 ) / ( 'EXTR' EPSI1 2 ) ;
  70. 'MESS' 'MODULE D YOUNG = ' YOUN1 ;
  71. *
  72. MATLIGN = 'MATE' modlign 'YOUN' YOUN1 'NU' .3 'TRAC' SSTBG ;
  73. *
  74. *BLO1 = 'BLOQ' 'DEPL' 'ROTA' P6 ;
  75. *BLO2 = 'BLOQ' 'UX' 'RY' P1 ;
  76. BLO1 = 'BLOQ' 'UY' 'UZ' 'RX' 'RZ' P6 ;
  77. BLO2 = 'BLOQ' 'UX' 'RY' lign1 ;
  78. rigblo1 = BLO1 ET BLO2 ;
  79. *
  80. * PAS DE CHARGE
  81. *
  82. DFE1 = 1.E3 ;
  83. FELF = 2.E3 ;
  84. *
  85. LISCEF = 'PROG' 0. PAS DFE1 FELF ;
  86. LISCF = LISCEF ;
  87. *
  88. DFE1 = 50. ;
  89. FELI = FELF ;
  90. FELF = 2.75E3 ;
  91. *
  92. LISCEF = 'PROG' ( FELI + DFE1 ) PAS DFE1 FELF ;
  93. LISCF = LISCF ET LISCEF ;
  94. *
  95. * CHARGEMENT
  96. *
  97. FP = 'FORC' 'FY' 1. P1 ;
  98. EVOLU1 = 'EVOL' 'MANU' 'PAS' LISCF 'FORC' LISCF ;
  99. CHATOT = 'CHAR' MECA FP EVOLU1 ;
  100. *
  101. CARTOT = MATLIGN ET CARLIGN ;
  102. RIG1 = 'RIGI' modlign CARTOT ;
  103. RIGTOT = RIG1 ET rigblo1 ;
  104. *
  105. ZTAB = 'TABLE' ;
  106. ZTAB.blocages_mecaniques = rigblo1 ;
  107. ZTAB.caracteristiques = cartot;
  108. ztab.modele = modlign ;
  109. ZTAB.chargement = CHATOT ;
  110. ZTAB.temps_calcules = LISCF ;
  111. ZTAB.'CONVERGENCE_FORCEE' = FAUX ;
  112. ZTAB.'DELTAITER' = 1000 ;
  113. PASAPAS ZTAB ;
  114. *
  115. * POSTRAITEMENTS
  116. *
  117. NELCOU = 'NBEL' COUD1 ;
  118. EPA1 = ECOU1 ;
  119. REX1 = REXT1 ;
  120. RAC1 = RCO1 ;
  121. RIN1 = REX1 - EPA1 ;
  122. XIN1 = PI * ( ( REX1 ** 4. ) - ( RIN1 ** 4. ) ) * 0.25 ;
  123. XJN1 = 2. * XIN1 ;
  124. SEC1 = PI * ( ( REX1 ** 2. ) - ( RIN1 ** 2. ) ) ;
  125. *
  126. XLAM1 = EPA1 * RAC1 * 4. / ( ( REX1 + RIN1 ) ** 2. ) ;
  127. PI4 = PI / 4. ;
  128. R33 = ( 3. ** 0.5 ) ;
  129. RMOY1 = REX1 - ( EPA1 * 0.5 ) ;
  130. R89 = 8. / 9. ;
  131. PUI23 = -2. / 3. ;
  132. GAM1 = R89 * ( XLAM1 ** PUI23 ) ;
  133. KK1 = 1.65 / XLAM1 ;
  134. ALF2 = PI4 * GAM1 ;
  135. FACT1 = ( RMOY1 * ALF2 ) / KK1 ;
  136. LCOUD = RAC1 * 90. ;
  137. LELEM = LCOUD / NELCOU ;
  138. *
  139. depl1 = ztab . deplacements ;
  140. effo1 = ztab . contraintes ;
  141. *********coupl1 = ztab . deformations_inelastiques ;
  142. *
  143. ntpr1 = ( 'DIME' depl1 ) - 1 ;
  144. *
  145. * verif dernier pas elastique
  146. *
  147. ipas = 3 ;
  148. *
  149. * moments
  150. *
  151. fimp1 = extr liscf ( ipas + 1 ) ;
  152. mtheo1 = fimp1 * lh1 ;
  153. effor1 = redu ( effo1 . ipas ) coud1 ;
  154. mcalc1 = extr effor1 'MOMY' 1 nelcou 2 ;
  155. *
  156. * deformations
  157. *
  158. sigi1 = 'EXTR' sstbg 'ORDO' ;
  159. epsi1 = 'EXTR' sstbg 'ABSC' ;
  160. deptot = depl1 . ipas ;
  161. coutr1 = redu ( 'EPSI' modlign deptot cartot ) coud1 ;
  162. cycal1 = extr coutr1 'CY' 1 nelcou 2 ;
  163. defyt1 = ( alf2 * rmoy1 * cycal1 ) / kk1 ;
  164. vmitr1 = redu ( 'VMIS' modlign ( effo1 . ipas ) cartot ) coud1 ;
  165. vmise1 = extr vmitr1 'SCAL' 1 nelcou 2 ;
  166. epsit1 = ipol vmise1 sigi1 epsi1 ;
  167. mess ' ' ;
  168. mess 'moment theorique = ' mtheo1 'moment calcule = ' mcalc1 ;
  169. res1 = abs ( ( mtheo1 - mcalc1 ) / mcalc1 ) ;
  170. slim1 = ( extr sigi1 2 ) * 1.e-6 ;
  171. vmim1 = vmise1 * 1.e-6 ;
  172. mess 'Limite elastique = ' slim1 ' Vmis Max = ' vmim1 ;
  173. mess 'deftot theorique = ' epsit1 'deftot calcule = ' defyt1 ;
  174. res2 = abs ( ( epsit1 - defyt1 ) / epsit1 ) ;
  175. progelas = 'PROG' res1 res2 ;
  176. maxerr1 = 'MAXI' progelas ;
  177. 'SI' (maxerr1 <eg 1.e-4 ) ;
  178. erre 0 ;
  179. 'SINON' ;
  180. erre 5 ;
  181. 'FINSI' ;
  182. *
  183. * verif dernier pas plastique
  184. *
  185. ipas = ntpr1 ;
  186. *
  187. * moments
  188. *
  189. fimp1 = extr liscf ( dime liscf ) ;
  190. mtheo1 = fimp1 * lh1 ;
  191. effor1 = redu ( effo1 . ipas ) coud1 ;
  192. mcalc1 = extr effor1 'MOMY' 1 nelcou 2 ;
  193. *
  194. * deformations
  195. *
  196. deptot = depl1 . ipas ;
  197. coutr1 = redu ( 'EPSI' modlign deptot cartot ) coud1 ;
  198. cycal1 = extr coutr1 'CY' 1 nelcou 2 ;
  199. defyt1 = ( alf2 * rmoy1 * cycal1 ) / kk1 ;
  200. vmitr1 = redu ( 'VMIS' modlign ( effo1 . ipas ) cartot ) coud1 ;
  201. vmise1 = extr vmitr1 'SCAL' 1 nelcou 2 ;
  202. epsit1 = ipol vmise1 sigi1 epsi1 ;
  203. mess ' ' ;
  204. mess 'moment theorique = ' mtheo1 'moment calcule = ' mcalc1 ;
  205. res1 = abs ( ( mtheo1 - mcalc1 ) / mcalc1 ) ;
  206. slim1 = ( extr sigi1 2 ) * 1.e-6 ;
  207. vmim1 = vmise1 * 1.e-6 ;
  208. mess 'Limite elastique = ' slim1 ' Vmis Max = ' vmim1 ;
  209. mess 'deftot theorique = ' epsit1 'deftot calcule = ' defyt1 ;
  210. res2 = abs ( ( epsit1 - defyt1 ) / epsit1 ) ;
  211. progplas = 'PROG' res1 res2 ;
  212. maxerr2 = 'MAXI' progplas ;
  213. 'SI' (maxerr1 <eg 1.e-4 ) ;
  214. erre 0 ;
  215. 'SINON' ;
  216. erre 5 ;
  217. 'FINSI' ;
  218. fin ;
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  

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