Télécharger q8ri_bcn.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : q8ri_bcn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 1 ;
  5. ******************************************************************
  6. * TEST Q8RI BCN
  7. *
  8. * Verification of element Q8RI (QUA8 with 2x2 Gauss points)
  9. *
  10. * Elastic analysis of a square subjected to biaxial extension
  11. *
  12. ******************************************************************
  13. TITRE 'VERIFICATION OF ELEMENT Q8RI' ;
  14. OPTI DIME 2 ELEM QUA8 MODE PLAN DEFO ;
  15.  
  16. GRAPH='N';
  17. *
  18. *------------ Geometry definition --------------------------------
  19. *
  20. pa = 0. 0. ; pb = 10. 0. ; pc = 10. 10. ; pd = 0. 10 ;
  21. lab = pa DROIT 10 pb ;lbc = pb DROIT 10 pc ;
  22. lcd = pc DROIT 10 pd ;lda = pd DROIT 10 pa ;
  23. mesh = DALLER lab lbc lcd lda ;
  24. *
  25. SI (NEG GRAPH 'N');
  26. TRACE mesh ;
  27. FINSI ;
  28. *
  29. *------------ Model and material ---------------------------------
  30. *
  31. mod1 = MODE mesh MECANIQUE ELASTIQUE Q8RI ;
  32. mat1 = MATE mod1 YOUN 1.E7 NU 0 ;
  33. *
  34. *------------ Stiffness and boundary conditions ------------------
  35. *
  36. rig1 = RIGI mod1 mat1 ;
  37. bc_lab = BLOQ UY lab ; bc_lbc = BLOQ UX lbc ;
  38. bc_lcd = BLOQ UY lcd ; bc_lda = BLOQ UX lda ;
  39. bc_tot = bc_lab ET bc_lbc ET bc_lcd ET bc_lda ;
  40. rigtot = rig1 ET bc_tot ;
  41. *
  42. *------------ Prescribed displacements: dilatation ---------------
  43. *
  44. di_lbc = DEPI bc_lbc 0.01 ;
  45. di_lcd = DEPI bc_lcd 0.01 ;
  46. ftot = di_lbc ET di_lcd ;
  47. *
  48. *------------ Linear system --------------------------------------
  49. *
  50. desp = RESO rigtot ftot ;
  51. *=================================================================
  52. * Postprocess
  53. *=================================================================
  54. *
  55. * Deformed mesh
  56. *
  57. SI (NEG GRAPH 'N') ;
  58. def0 = DEFO mesh desp 0. BLAN ;
  59. def1 = DEFO mesh desp 10. ROUG ;
  60. TRACE (def0 ET def1) ;
  61. FINSI ;
  62. *
  63. * Stresses computed with SIGM
  64. *
  65. stress1 = SIGM mod1 mat1 desp 'LINE' ;
  66. *
  67. * Stresses computed with EPSI and ELAS
  68. *
  69. strain = EPSI mod1 desp 'LINE' ;
  70. stress2 = ELAS mod1 strain mat1 ;
  71. *
  72. * Gradient of displacements computed with GRAD
  73. *
  74. grad1 = GRAD mod1 desp ;
  75. *
  76. * Verifications
  77. ERRE 0 ;
  78. * 1) stress1 = stress2
  79. *
  80. difstres = MAXI (stress1-stress2) ABS ;
  81. SI (difstres > 1.E-4) ; ERRE 5 ; FINSI ;
  82. *
  83. * 2) smxx = smyy = 1.e4 , smzz = smxy = 0
  84. *
  85. smxxmin = MINI (EXCO stress1 SMXX) ; dxxmin = ABS (smxxmin-1.E4) ;
  86. smxxmax = MAXI (EXCO stress1 SMXX) ; dxxmax = ABS (smxxmax-1.E4) ;
  87. smyymin = MINI (EXCO stress1 SMYY) ; dyymin = ABS (smyymin-1.E4) ;
  88. smyymax = MAXI (EXCO stress1 SMYY) ; dyymax = ABS (smyymax-1.E4) ;
  89. smzzmin = MINI (EXCO stress1 SMZZ) ; dzzmin = ABS smzzmin ;
  90. smzzmax = MAXI (EXCO stress1 SMZZ) ; dzzmax = ABS smzzmax ;
  91. smxymin = MINI (EXCO stress1 SMXY) ; dxymin = ABS smxymin ;
  92. smxymax = MAXI (EXCO stress1 SMXY) ; dxymax = ABS smxymax ;
  93. dif = PROG dxxmin dxxmax dyymin dyymax
  94. dzzmin dzzmax dxymin dxymax ;
  95. SI ( (MAXI dif ABS)> 1.E-4) ;
  96. ERRE 5 ;
  97. FINSI ;
  98. FIN ;
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  

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