Exemples de jeux de données Cast3M

* fichier :  grandksi.dgibi
************************************************************************
************************************************************************
complet = faux;
*
*************************************************
*  pour calcul complet mettre complet à : vrai;
***************************************************
opti echo 0;
********************************************************
*                                                      *
*   calcul fleche d'une plaque sous son propre poids   *
*                                                      *
* la plaque est horizontale : on joue sur inclinaison  *
*                                                      *
*            de G et de la Force de Traction           *
*                                                      *
*         TRAITE par METHODE des GRANDS DEPLACEMEMTS   *
*                                                      *
*                                                      *
*                                                      *
*                                                      * 
********************************************************
*
*
*                    geometrie
*
*
opti dime 3 elem qua4 ;
*
inclin = 0. ;
*
teta = 90. ;
*
temp = 30. ;
mess ' effort traction en Newton  0 si pas d''  efforts   sinon  15 ' ;
  eft = 15;
*obte eft ;
*
epu = 0. ;
ru=19520.;
epc=7E-04 ;
rc=1400.;
*
rcf=2.*epu*ru/(2*epu+epc)+rc;
dens 1 ;
vu2 = 1000. 0. 0. ;
vu3 = 0 100. 100. ;
*
long = 0.92 ;
larg = 0.04 ;
*
*     coordonnees des 6 sommets initiaux 
*
p1 =  0.  0.   0. ;
p2 =  0. .24   0. ;
p3 =  0. .68   0. ;
p4 =  0. .92   0. ;
*
p5 =  0.  0.   larg ;
p6 =  0. .24   larg ;
p7 =  0. .68   larg ;
p8 =  0. .92   larg ;
*
*
*     creation du point de traction du bec
*
 af =  p4 plus p8 ; af = af /2 ;
*
*
*     creation lignes s'appuyant sur les sommets
*
*
l1 = droi 2 p1 p5 ;
l2 = droi 2 p2 p6 ;
l3 = droi 2 p3 p7 ;
l4 = droi 2 p4 p8 ;
*
*     creation des milieux des lignes l2 et l3
*
m2 =  0. .24   ( larg/2 ) ;
m3 =  0. .35   ( larg/2 ) ;
m4 =  0. .46   ( larg/2 ) ;
m5 =  0. .57   ( larg/2 ) ;
m6 =  0. .68   ( larg/2 ) ;
*
*
*
si complet ;
na = 9;nb = 12; nc = 6;
sinon;
aa = 7; nb = 8; nc = 4;
finsi;
c1b = droi  9 p1 p2 ;
c2b = droi 12 p2 p3 ;
c3b = droi  6 p3 p4 ;
*
c1h = droi  9 p5 p6 ;
c2h = droi 12 p6 p7 ;
c3h = droi  6 p7 p8 ;
*
*      creation maillage
*
s1 = dall l1 c1b l2 c1h ;
s2 = dall l2 c2b l3 c2h ;
s3 = dall l3 c3b l4 c3h ;
*
sbec = s1 et s2 et s3 ;
*trac sbec;
sbec = sbec coul jaun ;
pf = sbec poin proc af ;
*
*********************************************************************
*                                                                   *
*                       impression des donnees                      *
*                                                                   *
*********************************************************************
*
mess ' angle inclinaison inclin  ' inclin 'degres' ;
mess ' angle rotation    teta    ' teta 'degres' ;
mess ' angle = 90  plaque a plat ' ;
mess ' temperature       temp    ' temp ;
mess ' force de traction en N eft   ' eft ;
mess ' epaisseur U       epu     ' epu ;
mess ' longueur plaque   long    ' long ;
mess ' largeur plaque    larg    ' larg ;
mess ' epaisseur bec     epc     ' epc ;
mess ' ro equivalent     rcf     ' rcf ;
*
*
*********************************************************************
*********************************************************************
*
**********   modele du materiau et  caracteristiques
*
*
mo = mode sbec MECANIQUE ELASTIQUE coq4 ;
matx = mate mo 'YOUN' 35E+09 'RHO' rcf 'ALPHA' 2.5E-06 'NU' 0.3 ;
ca = cara mo EPAI epc ;
mat1 = matx et ca ;
*
**********   matrice de masse ; acceleration pesanteur ; poids
*
m = mass mo mat1 ;
*
stet = sin teta  ; ctet = cos teta ;
cic = cos inclin ; sic = sin inclin ;
gx = +9.81 * stet * cic ;
gy = +9.81 * sic ;
gz = -9.81 * ctet * cic ;
accel = manu chpo sbec 6 ux gx uy gy uz gz mx 0. my 0. mz 0. ;
*
poids = m * accel  ;
*
**********   condition aux limites
*
*
*    encastrement = bloq rota et depl
*
*
*    appui = bloq depl ;  variante  = bloq  ux et-ou uy et-ou uz
*
**********   blocage rotation et deplacement l1
*
bl1rd = bloq rota depl l1  ;
*
**********   blocage deplacement [ ux uz ]  l4
*
bl4d = bloq l4 ux uz ;
*
blqtot = bl1rd et bl4d ;
*
**********   creation force de traction
*
fo1 = force (0. eft  0.) pf ;
*
**********    RESOLUTION PAR RESO
*
   sc = poids + fo1 ;
********************************    RESOLUTION PAR NONLIN
*
*
**********   creation chargement
*
scf = fo1 ; scp = poids ;
ev1 = evol manu xx1 ( prog 0. 1.1 ) ff1 ( prog 0. 1.1 ) ;
ch1 = char meca scf ev1 ;
ch2 = char meca scp ev1 ;
cht = ch1 et ch2 ;
sctot = scf + scp;
*
**********
*
tgd = table ;
tgd.grands_deplacements = vrai ;
tgd.k_sigma = vrai ;
tgd.'REAC_GRANDS'=100.;
*  liste des temps a calculer
si complet ;
tfin = 0.37;
sinon;
tfin = 0.07;
finsi;
 tcal = prog 0.  0.0001 0.0003 0.001 0.003 0.01  0.02         
 0.04   pas 0.03   tfin;                                                                 
 
tgd.blocages_mecaniques = blqtot; 
tgd.caracteristiques = mat1;
tgd.modele = mo;                             
tgd.temps_calcules = tcal;
tgd.chargement = cht;  
tgd.hypothese_deformations = 'LINEAIRE' ;
 
 
pasapas tgd ;
 
**********   resultats
*
dep1 = peche tgd deplacements  tfin;
aa = maxi ( exco dep1 UX);
si complet;
valcomp = 2.0938e-03;
sinon;
valcomp = 1.24824E-03;
finsi;
err = ( aa - valcomp ) / valcomp  abs;
err100 = err * 100;
mess '       ';
mess '       ';
mess '       ';
mess '       ';
mess ' la fleche maximum attendue est  ' valcomp;
mess '       ';
mess '       ';
mess ' valeur de la fleche  trouvée : ' aa ;
mess '       ';
mess '       ';
mess  ' erreur commise en % ' err100;
mess '       ';
mess '       ';
mess '       ';
mess '       ';
si ( err > 0.005) ;  erreur 5; finsi;
fin;