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

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