Accueil | Cast3M

* fichier :  testICQ4.dgibi
************************************************************************
************************************************************************
*           Test TestICQ4.dgibi: Jeux de données      *
*           ------------------------------------      *
*                                                     *
*
* FICHIER GIBIANE POUR TESTER LES ELEMENTS
* INCOMPRESSIBLES ICQ4
*
* Maillage d'une plaque de forme carrée avec les 
* éléments ICQ4.
* 
* Cas déformation plane.
*
* On soumet la plaque à un champ de contrainte 
* hydrostatique (contribution: alpha=1) et à un champ 
* de contrainte déviatorique (contribution: beta=0.001)
* 
* Les résultats théoriques sont déterminés dans le cas 
* incompressible (coefficient de Poisson proche de 0.5)
*
* cf rapport DMT sur ces elements
*
*
*-----MAILLAGE-----------
*
opti dime 2 mode plan defo elem qua4 ;
*
p1 = 0. 0. ;
p2 = 100. 0. ;
p3 = 100. 100. ;
p4 = 0. 100. ;
*
li1 = d 10 p1 p2 ;
li2 = d 10 p2 p3 ;
li3 = d 10 p3 p4 ;
li4 = d 10 p4 p1 ;
*
sur1 = dall li1 li2 li3 li4 'PLAN' ;
*
*
*-----MODELE MECANIQUE------------
*
* Valeur du coefficient de Poisson
*
nu1 = 0.499 ;
*
mod1 = modele sur1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 
       icq4 ;
mat1 = mate mod1 'YOUNG' 2.e11 'NU' nu1 ;
rigi1 = rigi mod1 mat1 ;
*
*
*-----CONDITIONS AUX LIMITES------------
*
bl1 = bloq li1 'UY' ;
bl2 = bloq li4 'UY' ;
bl3 = bloq p1 'UX' ;
*
*
*----------------CHARGEMENT-----------------
*
* Coefficients de contribution de chargement
* (alpha-->hydrostatique)
* (beta-->déviatorique)
*
alpha = 1. ;
beta = 1.e-3 ;
*
*-----Chargement déviatorique------------
*
*-----sur li1-----
*
chpx1 = coor 1 li1 ;
chpresy1 = pres 'MASS' mod1 (chpx1 * (-2.)*1.e5) ;
chpres1 = beta * chpresy1 ;
vecfor1 = vect chpres1 1.E-7 'FX' 'FY' rouge ;
*
*-----sur li2-----
*
chpy2 = coor 2 li2 ;
chpresx2 = pres 'MASS' mod1 (200.*1.e5) li2 ;
chpresy2 = pres 'MASS' mod1 (chpy2 * (-2.)*1.e5) ;
chforx2 = chpresx2 ;
chfory2 = exco 'FX' chpresy2 'FY';
chpres2 = beta * ( chforx2 et chfory2 ) ;
vecfor2 = vect chpres2 'FX' 'FY' 1.E-7 rouge ; 
*
*-----sur li3-----
*
chpx3 = coor 1 li3 ;
chpresx3 = pres 'MASS' mod1 (-200.*1.e5) li3 ;
chpresy3 = pres 'MASS' mod1 (chpx3 * (-2.)*1.e5) ;
chforx3 = exco 'FY' chpresx3 'FX';
chfory3 = chpresy3 ;
chpres3 = beta * ( chforx3 et chfory3 ) ;
vecfor3 = vect chpres3 'FX' 'FY' 1.E-7 rouge ; 
*
*-----sur li4-------
*
chpy4 = coor 2 li4 ;
chpresy4 = pres 'MASS' mod1 (chpy4 * (-2.)*1.e5) ;
chfory4 = exco 'FX' chpresy4 'FY';
chpres4 =  beta * chfory4 ;
vecfor4 = vect chpres4 1.E-7 'FX' 'FY' rouge ;
*
*
*-----Chargement hydrostatique------------
*
*-----sur li1-----
*
presy1 = pres 'MASS' mod1 (1.e5) li1 ;
pres1 = alpha * presy1 ;
vfor1 = vect pres1 1.E-5 'FX' 'FY' vert ;
*
*-----sur li2-----
*
presx2 = pres 'MASS' mod1 (1.e5) li2 ;
pres2 = alpha * presx2 ;
vfor2 = vect pres2 'FX' 'FY' 1.E-5 vert ; 
*
*-----sur li3-----
*
presy3 = pres 'MASS' mod1 (1.e5) li3 ;
pres3 = alpha * presy3 ;
vfor3 = vect pres3 'FX' 'FY' 1.E-5 vert ; 
*
*-----sur li4-------
*
presx4 = pres 'MASS' mod1 (1.e5) li4 ;
pres4 = alpha * presx4 ;
vfor4 = vect pres4 1.E-5 'FX' 'FY' vert ;
*
*
*---------CALCUL DES DEPLACEMENTS --------------------
*
*
* Solution analytique au point P3: uyan.
*
x3 = coor 1 p3 ;
y3 = coor 2 p3 ;
mu = 1.e11 / ( 1 + nu1 ) ;
uyan = ( 2 * x3 * y3 ) / ( 2 * mu ) * 1.e5 * beta ;
*
* Solution Castem au point P3: uycas.
*
deptot = reso (rigi1 et bl1 et bl2 et bl3)
(chpres2 et chpres3 et chpres4 et chpres1 
et pres2 et pres3 et pres4 et pres1) ;
uycas = extr deptot 'UY' p3 ;
*
* Calcul de l'erreur sur le déplacement.
*
erry = abs (( uyan - uycas ) / ( uyan )) * 100 ;
*
* Test de précision sur le calcul du déplacement uy 
* (<1%).
*
si ( erry >eg 1. ) ;
     erreur 5 ;
finsi ;
*
*
*-----------CALCUL CONTRAINTES-------------------------
*
*
* Solution analytique: pression = alpha*1.E5.
*
* Solution Castem : pression maximale = maxpres1.
*
sig1 = sigma mod1 deptot mat1 ;
sxx = exco 'SMXX' sig1 'P' ;
syy = exco 'SMYY' sig1 'P' ;
szz = exco 'SMZZ' sig1 'P' ;
sp = (sxx + syy + szz) / 3. ;
maxpres1 = maxi sp ;
*
* Calcul de l'erreur sur la pression.
*
errpress = abs ((maxpres1 + (alpha * 1.E5)) / 
           (alpha * 1.E5)) * 100. ;
*
* Test de précision sur le calcul de la pression (<1%).
*
si ( errpress >eg 1. ) ;
     erreur 5 ;
finsi ;
*
*
fin;