Télécharger la plaquette Cast3M


fr en

Télécharger ex_proper.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ex_proper.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. GRAPH='N';                                                                      
  5. SAUT PAGE ;    
  6. *-----------------------------------------------------
  7. * TEST  PLUS, MOINS, DEDU, TOUR
  8. * OPERANDES MCHAML, MMODEL, CHPO, MAILLAGE (RIGIDITE)
  9. *-----------------------------------------------------
  10. *
  11. opti dime 2 elem qua4 ;
  12. *
  13. p0 = 1. 0. ; e_x = p0 ; e_y = 0. 1. ; p_ori = 0. 0. ;
  14. p1 = p0 plus e_x ;
  15. s1 = (p0 d 1 p1) trans 1 e_y ;
  16. lc1 = c (inve (s1 cote 4)) p0 p_ori 4 ;
  17. cont = (d 6 p_ori p0)  et lc1 ;
  18. su1 = surf cont plan ;
  19. st1 = s1 et su1 ;
  20. chp1 = manu chpo 1 s1  fy -1. ;
  21. * tests PLUS, MOINS, DEDU, TOUR
  22. m1 = mode st1 mecanique elastique ;
  23. che1 = manu chml m1 sxx 1. ;
  24. che3 = manu chml m1 smyy 2. smxy 1. ;
  25. che5 = manu chml m1 smxx 5.66 eqeq 5. ;
  26. m3 = mode (s1 plus e_x) mecanique elastique ;
  27. m2 che2  = m3 che1 plus (2.*e_y) ;
  28. m4 s3 = m1 s1 moins (2. 5.) ;
  29. *list che2 ;
  30. s2 = s1 plus (1. -1.) ;
  31. chp10 = s2 moins s1 ;
  32. *
  33. chp21 = chp1 dedu trans s1 s2 ;
  34. *
  35. p3 = tourne p1 90. p_ori ;
  36. m5 s4 = m1 s1 tourne 400. p0 ;
  37. chp22 = chp1 dedu rota 400. p0 s1 s4 ;
  38. *
  39. chp5 = manu chpo 1 s1 scal 10. ;
  40. chp6 = tour chp5 30. p1 ;
  41. chp2 = tour chp1 180. p_ori ;
  42. *
  43. che4  = tour che3 90. p_ori ;
  44. che6 = tour che5 45. p_ori ;
  45. *
  46. opti dime 3 elem cub8 ;
  47. e_z = 0. 0. 1. ;
  48. p5 = p1 tour 50. p_ori (p_ori plus e_z) ;
  49. v1 = s1 volu 2 trans e_z ;
  50. vt1 = st1 volu trans 1 e_z ;
  51. vt2 = st1 volu 4 rota 45. p_ori (p_ori plus e_y) ;
  52. mv1 = mode vt1 mecanique elastique ;
  53. mv4 v3 = mv1 v1 moins (2. 5. 1.) ;
  54. chp11 = v3 moins v1 ;
  55. mv5 v4 = mv1 v1 tour 120. p0 (p0 plus e_z)  ;
  56. vt6 = vt1 tour 120. p0 (p0 plus e_z)  ;
  57. su7 = su1 dedu tran vt1 vt6 ;
  58. su5 = su1 dedu 120. p0 (p0 plus e_z)  rota vt1 vt6 ;
  59. vt5 = extr mv5 mail ;
  60. su8 = su1 dedu 120. p0 (p0 plus e_z)  rota vt1 vt5  ;
  61. *
  62. mv6 che7 = mv1 che5 tour 50. p1 (p1 plus e_z) ;
  63. *
  64. t1 = extr che7 smyy 2 3 4 ;
  65. tref = 3.3214 ;
  66. resi = abs((t1 - tref)/tref) ;
  67. si (resi <eg 5.e-2) ;
  68.  erre 0 ;
  69. sinon ;
  70.  erre 5 ;
  71. finsi ;
  72.  
  73. fin ;
  74.  
  75.  
  76.  
  77.  
  78.  

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