Nana Jefferson
Inscrit le: 04 Juil 2016 |
Messages: 4 |
|
|
|
Posté le: 30 Aoû , 2016 14:38 |
|
|
|
|
|
Bonjour, j'ai commencé à utiliser castem il y a quelques semaines et je n'arrive pas à appliquer une pression dans un reservoir cônique qui décroit avec la hauteur d'eau dans ce dernier. voici ce que j'ai essayé de faire, mais je compte utilisé pas à pas et je ne sais pas si cela conviendra.
OPTI DIME 2 MODE AXIS ELEM QUA4;
*** CREATION D_UN TABLEAU CONTENANT TOUTE LES VALEURS DE LA HAUTEUR D_EAU;
HAUT = TABL;
HAUT.1 = 0.23;
HAUT.2 = 0.194;
HAUT.3 = 0.1663;
HAUT.4 = 0.147;
HAUT.5 = 0.1315;
HAUT.6 = 0.1178;
HAUT.7 = 0.104;
HAUT.8 = 0.09324;
HAUT.9 = 0.0828;
HAUT.10 = 0.0769;
HAUT.11 = 0.068;
HAUT.12 = 0.0642;
HAUT.13 = 0.0577;
HAUT.14 = 0.0472;
HAUT.15 = 0.0431;
HAUT.16 = 0.0401;
HAUT.17 = 0.0359;
HAUT.18 = 0.0323;
HAUT.19 = 0.029;
HAUT.20 = 0.0255;
HAUT.21 = 0.0228;
HAUT.22 = 0.0198;
HAUT.23 = 0.0172;
HAUT.24 = 0.0166;
HAUT.25 = 0.0148;
HAUT.26 = 0.0138;
HAUT.27 = 0.0096;
HAUT.28 = 0.0077;
HAUT.29 = 0.0053;
HAUT.30 = 0.;
*** DEFINITION DU PAS D'INCREMENTATION
I = 0;
REPETER BOUCLE 29;
*** DEFINITION DES LIGNES
I = I+1;
P0 = POIN 0. 0.;
P1 = POIN 0. 0.01;
P2 = POIN 0.10388 0.;
P3 = POIN 0.09500 0.01;
P4 = POIN 0.13140 0.23;
P5 = POIN 0.12250 0.24;
P6 = POIN 0.15250 0.23;
P7 = POIN 0.15250 0.24;
P5P = POIN (0.095 +(0.1196*HAUT. I)) HAUT. I;
*** DEFINITION DES LIGNES
L1 = D 20 P1 P3;
L2 = D 20 P3 P5P;
L2P = D 20 P5P P5;
L3 = D 20 P5 P7;
L4 = D 20 P7 P6;
L5 = D 20 P6 P4;
L6 = D 20 P4 P2;
L7 = D 20 P2 P0;
L8 = D 20 P0 P1;
LT = L1 ET L2 ET L2P ET L3 ET L4 ET L5 ET L6 ET L7 ET L8;
*** DEFINITION DES SURFACES
S1 = LT SURF PLAN;
MOD1 = MODE S1 MECANIQUE ELASTIQUE ISOTROPE;
MAT1 = MATE MOD1 YOUN 1.5E7 NU 0.2 RHO 2000.;
*** CONDITIONS AUX LIMITES
CL1 = BLOQ UZ P4;
*** CHARGEMENT
* POIDS PROPRE
NOMKG = MASS MOD1 MAT1;
NOMG = MANU CHPO S1 1 UZ -9.81;
PP = NOMKG*NOMG;
VECP = VECT PP 1.E-1 FR FZ BLEU;
* POIDS DE L'EAU AGISSANT SUR LA PAROIS INTERNE
CHCO = COOR 2 L2;
H = MANU CHPO L2 1 'SCAL' HAUT.I;
CHPR1 = (H- CHCO);
CHEAU = MANU CHPO L2 1 'SCAL' 1.E4;
CHPR2 = CHEAU*CHPR1;
CHPRES1 = PRES MASS MOD1 CHPR2;
VEC1 = VECT CHPRES1 1.E-2 FR FZ ROUG;
* POIDS DE L'EAU AGISSANT SUR LA BASE
CHPRES2 = PRES MASS MOD1 L1 ((HAUT.I)*1.E4);
VEC2 = VECT CHPRES2 1.E-2 FR FZ BLEU;
* CHARGEMENT TOTAL : POIDS PROPRE + POIDS DE L'EAU
CHAG1 = PP ET CHPRES1 ET CHPRES2;
FIN BOUCLE;
|
_________________ @wanaje
|