Télécharger Cast_test_RelaCoq.dgibi

Retour à la liste

Numérotation des lignes :

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

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