* * OPTI ECHO 1; * ******************* cette ligne fait 72 caractères ********************* * Essai de contact en 2D Axis * 2 cylindres en appui * Déplacement imposé (vers le bas) du cylindre supérieur * Essai de calcul direct avec RESO * * OPTI DIME 2 ELEM TRI3 MODE AXIS; * * maillage * b0 = 0. 0.; b1 = 5. 0.; b2 = 5. 20.; b3 = 0. 20.; * lb1 = D 10 b0 b1; lb2 = D 10 b1 b2; lb3 = D 10 b2 b3; lb4 = D 10 b3 b0; lingot = (lb1 ET lb2 ET lb3 ET lb4) SURF; lingot = lingot COUL BLEU ; * h0 = 0. 20.; h1 = 5. 20.; h2 = 5. 25.; h3 = 0. 25.; * lh1 = D 10 h0 h1; lh2 = D 5 h1 h2; lh3 = D 5 h2 h3; lh4 = D 5 h3 h0; poinceau = (lh1 ET lh2 ET lh3 ET lh4) SURF; poinceau = poinceau COUL ROUGE ; * TRAC (poinceau ET lingot); * ******************* cette ligne fait 72 caractères ********************* * calcul mo1 = MODE poinceau MECANIQUE ELASTIQUE ; mat1 = MATE mo1 YOUN 5.e5 NU 0.3 RHO 15.; * mo2 = MODE lingot MECANIQUE ELASTIQUE ; mat2 = MATE mo2 YOUN 4.e5 NU 0.3 RHO 1.; * cont1 = IMPO MAIL (INVE lb3) lh1 COUL VERT; tout = lingot ET poinceau ET cont1; * *TRAC (cont1 ET poinceau ET lingot); * * frottements mf1 = IMPF (cont1); * mof1 = MODE mf1 CONTACT FROTTANT COULOMB; * matf1 = MATE mof1 'MU' 0.6 'COHE' 0. 'ADHE' 0.; * motot1 = mo1 ET mo2 ET mof1; mattot1 = mat1 ET mat2 ; * *Conditions aux limites : cl1 = BLOQ lb1 UZ; cl2 = BLOQ lb1 UR; cl3 = BLOQ lb4 UR; cl4 = BLOQ lh4 UR; clpress = BLOQ lh3 UZ; * * *Chargement méca, déplacement imposé chapress = DEPI clpress (-1.) ; * *rig1 = RIGI mo1 mat1; *rig2 = RIGI mo2 mat2; *rigct = RIGI mof1 matf1; *dep1 = RESO (rig1 ET rig2 ET rigct ET clpress) chapress ; * rigtot = RIGI ( mo1 ET mo2 ET mof1) (mat1 ET mat2 ET matf1) ; dep1 = RESO (rigtot ET clpress) chapress ; * def0 = DEFORME ( tout) dep1 0. BLEU; def1 = DEFORME (tout) dep1 5. ROUGE; def11 = DEFORME (tout) dep1 2.; TRAC (def11 ) ; OPTIO DONN 5; FIN;