Prbléme sur castem

dfhbbbb Messages postés 1 Date d'inscription samedi 23 mars 2024 Statut Membre Dernière intervention 23 mars 2024 - Modifié le 23 mars 2024 à 20:02
 dfhbbbb - 25 mars 2024 à 13:46

Bonjour,

Je travaille sur un projet sur castem , le but est de calculer les modes propre d'un portique en forme de triangle rectangle. Mon probléme est du a une erreur " INDETERMINATION DETECTEE AU NOEUD     24 INCONNUE RY
 INDETERMINATION LEVEE PAR LA MISE A ZERO DE LA SUSDITE dans mondes
 5 corrections de residu dans la resolution. Precision obtenue: 1.49335E-07 Modes d'ensemble: 19"

OPTI DIME 3 ELEM SEG2 ;

OEIL = 0. 0. 100. ; 

*DONNEES

LONG1 = 10. ;

NELEM=20;

RHOP = 7800. ;

NUP = 0.3 ;

YOP = 2.1E11 ;

*GEOMETRIE

PA1 = 0. 0. 0. ;

PC1 = LONG1 0. 0. ;

PB1= 0. LONG1. 0.;

L1 = PA1 DROI NELEM  PC1;

L2 = PB1 DROI NELEM  PC1;

L3 = PA1 DROI NELEM  PB1 ; 

ST1= L1 ET L2 ET L3;

TRAC OEIL QUAL ST1;

*CALCUL DES DONNEES

SEC1 = 1.982E-4; 

IY = 2.973E-10 ;

IZ =2.673E-10  ;

TOR = 1.E10 ;

*DEFINITION DU MODELE ET DU MATERIAU

MOD1 = MODEL ST1 MECANIQUE ELASTIQUE ISOTROPE POUT ;

MAT1 = MATE MOD1 YOUN YOP NU NUP RHO RHOP ;

CAR1 = CARAC MOD1 SECT SEC1 INRY IY INRZ IZ TORS TOR ;

MATTOT = MAT1 ET CAR1 ; 

*CONDITION LIMITE

AMPFY = 500. ; 

F1= FORCE ( AMPFY. 0. 0.) PC1 ; 

CL1 = BLOQ DEPLA ROTA PA1;

CL2 = BLOQ DEPLA ROTA PB1;

CLT= CL1 ET CL2;

*MATRICE DE RIGIDITE

RIG1 = RIGI MOD1 MATTOT ;

RIGCL1 = RIG1 ET CLT ;

MAS1 = MASS MOD1 MATTOT ;

FREQINF=0;

FREQSUP=500;



*CALCUL DES MODES

TAB1 =  VIBR INTER 0 500 RIGCL1 MAS1 TBAS  IMPR;

TAB2 = TAB1.modes;

*RESULTATS

F1= (TAB2.1).frequence ;

F2= (TAB2.2).frequence ;

F3= (TAB2.3).frequence ;

F4= (TAB2.4).frequence ;

F5= (TAB2.5).frequence ;

F6= (TAB2.6).frequence ;

MESS 'FREQUENCE 1 (HZ):' F1;

MESS 'FREQUENCE 2 (HZ):' F2;

MESS 'FREQUENCE 3 (HZ):' F3;

MESS 'FREQUENCE 4 (HZ):' F4;

MESS 'FREQUENCE 5 (HZ):' F5;

MESS 'FREQUENCE 6 (HZ):' F6;

*TRACE DES DEFORMEES MODALES

DEP1=(TAB2.1).deformee_modale ;

DEF1= DEFO DEP1 ST1 1. ROUGE;

DEP2=(TAB2.2).deformee_modale ;

DEF2= DEFO DEP2 ST1 1. VERT ;

DEP3=(TAB2.3).deformee_modale ;

DEF3= DEFO DEP3 ST1 1. BLEU;

DEP4=(TAB2.4).deformee_modale ;

DEF4= DEFO DEP4 ST1 1. JAUNE;

DEP5=(TAB2.5).deformee_modale ;

DEF5= DEFO DEP5 ST1 1. ROSE;

DEP6=(TAB2.6).deformee_modale ;

DEF6= DEFO DEP6 ST1 1. ORANGE;

DEF0=DEFO DEP1 ST1 0.;



TITRE 'Deformees modales (MODE 1)' ;

TRAC (DEF0 ET DEF1);

TITRE 'Deformees modales (MODE 2)' ;

TRAC (DEF0 ET DEF2);

TITRE 'Deformees modales (MODE 3)' ;

TRAC (DEF0 ET DEF3);

TITRE 'Deformees modales (MODE 4)' ;

TRAC (DEF0 ET DEF4);

TITRE 'Deformees modales (MODE 5)' ;

TRAC (DEF0 ET DEF5);

TITRE 'Deformees modales (MODE 6)' ;

TRAC (DEF0 ET DEF6);

FIN;

1 réponse

Je dirais qu'il y a un problème dans vos conditions aux limites.

0

bonjour 

Merci pour votre réponse , Sachant que je veux encastrer la ligne 3, j'ai donc encastré les point A et B, donc il ne devrait plus avoir de problème de rotation ?   

0
Rejoignez-nous