Download full text - ELSA - Europa

S3=P33 PLUS VEC1;
*
N1=2;
N2=3;
N3=4;
*
P11_P12=P11 D N2 P12;
P11_P12P=P11_P12 PLUS VEC1;
P13_P14=P13 D N2 P14;
P13_P14P=P13_P14 PLUS VEC1;
*
DOM11=P11_P12 REGL N3 P11_P12P;
DOM12=P13_P14 REGL N3 P13_P14P;
*
P21_P22=P21 D N2 P22;
P21_P22P=P21_P22 PLUS VEC1;
P23_P24=P23 D N1 P24;
P23_P24P=P23_P24 PLUS VEC1;
*
DOM21=P21_P22 REGL N3 P21_P22P;
DOM22=P23_P24 REGL N3 P23_P24P;
*
P31_P32=P31 D N2 P32;
P31_P32P=P31_P32 PLUS VEC1;
P33_P34=P33 D N1 P34;
P33_P34P=P33_P34 PLUS VEC1;
*
DOM31=P31_P32 REGL N3 P31_P32P;
DOM32=P33_P34 REGL N3 P33_P34P;
*
TRAC (DOM11 ET DOM12 ET DOM21 ET DOM22 ET DOM31 ET DOM32);
*
ELIM 0.001 DOM11 (P1 ET R1);
ELIM 0.001 DOM12 (R1 ET S1);
ELIM 0.001 DOM21 (P2 ET R2);
ELIM 0.001 DOM22 (R2 ET S2);
ELIM 0.001 DOM31 (P3 ET R3);
ELIM 0.001 DOM32 (R3 ET S3);
*
TOUT=DOM11 ET DOM12 ET DOM21 ET DOM22 ET DOM31 ET DOM32;
*
OPTI SAUV FORMAT 'doma01.msh';
SAUV FORMAT TOUT;
*
*=======================================================================
* Reduction Craig-Bampton sur le domaine 2 du cas 3
*
NMOD1=10;
*
MOD1=MODELE DOM32 MECANIQUE ELASTIQUE ISOTROPE;
MAT1=MATERIAU MOD1 YOUN 210.E9 NU 0.3 RHO 7800.;
*
CL1=BLOQ DEPL P33_P34;
*
MAS1=MASSE MOD1 MAT1;
RIG1=(RIGIDITE MOD1 MAT1) ET CL1;
*
* Calcul des modes statiques
*
TDEPI1=TABLE;
TLREE1=TABLE;
NBN1=NBNO(P33_P34);
NSTA1=0;
*
REPETER BCL1 (NBN1);
SI (&BCL1 EGA 1);
LREEL1=PROG 1 * 1. (NBN1-1) * 0.;
LREEL2=PROG NBN1 * 0.;
SINON;
SI (&BCL1 EGA NBN1);
LREEL1=PROG (NBN1-1) * 0. 1 * 1.;
LREEL2=PROG NBN1 * 0.;
SINON;
LREEL1=PROG (&BCL1-1) * 0. 1. (NBN1-&BCL1) * 0.;
LREEL2=PROG NBN1 * 0.;
FINSI;
FINSI;
4