Télécharger la plaquette Cast3M


fr en

Télécharger Cast_test_RelaCoq.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : cast_test_relacoq.dgibi
  2. *
  3. **********************************************************************
  4. ************************************************************************
  5. ************************************************************************
  6. *
  7. * Tranfert of load for shell elements - load known for a solid model
  8. *
  9. *
  10. * D. COMBESCURE
  11. *
  12. * F4E - 6th January 2014
  13. ********************************************************************
  14. OPTI DIME 3 ELEM CUB8;
  15. *
  16. * Modele 1: Forces calculated with another software/here created manually
  17. *
  18. FLAGVISU = FAUX ;
  19. *
  20. L1 = 0.1;
  21. L2 = 0.25;
  22. L3 = 3.;
  23. *
  24. fac1 = 1;
  25. n1 = 2*fac1;
  26. n2 = 3*fac1;
  27. n3 = 10*fac1;
  28. *
  29. p0 = (-0.5*L1) 0. 0.;
  30. p1 = (0.5*L1) 0. 0.;
  31. vy = 0. L2 0.;
  32. vz = 0. 0. L3;
  33. p0c = 0. 0. 0.;
  34. p1c = 0. L2 0.;
  35. p0top = p0 plus (vz);
  36. p1top = p1 plus (vz);
  37. *
  38. Lig0 = D n1 p0 p1;
  39. Sur0 = Lig0 tran n2 vy;
  40. Vol0 = Sur0 volu n3 tran vz;
  41. elim 0.001 (vol0 et p0top et p1top);
  42. *
  43. Mod0 = mode vol0 mecanique elastique CUB8;
  44. Emat0 = 200.D9;
  45. Nmat0 = 0.3;
  46. romat0 = 7800.;
  47. mat0 = MATE mod0 youn Emat0 NU Nmat0 RHO romat0;
  48. *
  49. rig0 = rigi Mod0 mat0;
  50. mas0 = mass Mod0 Mat0;
  51. *
  52. * Uniform acceleration
  53. *for0 = mas0*accuni;
  54. *
  55. * Triangular acceleration
  56. *for0 = mas0 * acc0;
  57. *
  58. * Lumped force
  59. for0 = FORC  FX 1. p0top;
  60. *
  61. SI FLAGVISU;
  62.   TRAC VOL0 (VECTEUR for0 FX FY FZ) cach;
  63. FINSI;
  64. *
  65. *
  66. * Modele 1 - Mechanical model - Shell model
  67. *
  68. n1b = 1*fac1;
  69. n2b = 3*fac1;
  70. n3b = 8*fac1;
  71. ep1 = 0.00000000001;
  72. ep0 = (-1.)*ep1;
  73. ep2 = ep1;
  74. ep3 = ep1;
  75. *
  76. p0c = 0. 0. 0.;
  77. p1c = 0. L2 0.;
  78. vz = 0. 0. L3;
  79. *
  80. opti elem tri3;
  81. Lig1c = D n2b p0c p1c;
  82. vol1 = Lig1c tran n3b vz;
  83. *
  84. *
  85. Mod1 = mode vol1 mecanique elastique DKT;
  86. Emat1 = 200.D9;
  87. Nmat1 = 0.3;
  88. romat1 = 7800.;
  89. mat1 = MATE mod1 youn Emat1 NU Nmat1 RHO romat1 EPAI L1;
  90. *
  91. rig1 = rigi Mod1 mat1;
  92. mas1 = mass Mod1 Mat1;
  93. *
  94. *
  95. bl0 = bloq depl vol1;
  96. bl00 = bloq depl rota lig1C;
  97. *
  98. rigmod = rig1 et bl00;
  99. *
  100. tab11 = table;
  101. tab11.'RELATION_3D' = VRAI;
  102. vol1_3D = coq2mas mod1 mat1 tab11;
  103. blcoq = tab11.'RELATION_3D';
  104. *
  105. *
  106. * Transfert through the Lagrangian multiplier
  107. *
  108. *
  109. frdep0 = for0;
  110. *
  111. * Relation maitre-esclave
  112. *
  113. * Vol0 suit Vol1_3D qui suit Vol1
  114. *
  115. rel2 = rela vol0 'ACCRO' vol1_3D;
  116. rigtot = rig1 et bl00 et rel2 et blcoq;
  117. depj1 = reso rigtot frdep0;
  118. *
  119. * Calcul des forces nodales a partir des deplacements sur le modele 1
  120. frdep1 = (rig1 et bl00) * depj1;
  121. *
  122. SI FLAGVISU;
  123.  oeil = 2. 2. 4.;
  124.  TITRE 'Transfer through Lagrangian muliplier';
  125.  trac oeil (defo vol1_3D depj1) CACH;
  126.  trac oeil (defo vol1 depj1) CACH;
  127.  trac oeil vol1 (vecteur frdep1 FX FY FZ) cach;
  128. FINSI;
  129. *
  130. *
  131. res0 = maxi (resu frdep0);
  132. res1 = maxi (resu frdep1);
  133. *
  134. ERR = (res0 - res1)/res0;
  135. *                                                                      *
  136. 'SI' (ERR '>' 1.E-6) ; 'ERRE' 5 ; 'FINS' ;
  137. *
  138. 'FIN' ;
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  

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