Accueil | Cast3M

* fichier :  dyna_nl3.dgibi
************************************************************************
************************************************************************
 
 
*         Test Dyna_nl3.dgibi: Jeux de données        *
*         ------------------------------------        *
*                                                     *
complet = vrai;
*
*************************************************
*  pour calcul complet mettre complet à : vrai;
***************************************************
*======================================================
*       TEST DES LA PRESSION SUIVEUSE EN DYNAMIQUE
*
*       FLOTTEMENT D'UNE POUTRE ENCASTREE-LIBRE
*
* On etudie la reponse de la poutre à un lacher quand 
* elle est soumise à une force de compression axiale
* legerement superieure à la valeur theorique de la 
* force critique de flottement (Fcr = 20.05*E*I/l**2)
* On doit obtenir une reponse d'amplitude croissante.
*
*******************************************************
*                                                    ==   
*    /|aa___________________bb                       ==
*    /|______________________|<----- F               ==
*    /|                        cc                    ==
*           EQUIVALENT A     /|                      ==  
*                        dd / |                      ==  
*    /|aa_________________bb|<--                     ==
*    /|_____________________|<--  PRESSION           ==
*    /|                     |<--                     ==
*                        ee |/ ff                    ==
*                                                    ==
*                                                    ==
*====================================================== 
graph = faux;
option dime 3 elem seg2 ;
*tri3 ;
*======================================================
* definition de la geometrie
* la force suiveuse est modelisee par une pression équi
* valente appliquée sur une petite plaque fixée à l'ex
* trémité de la poutre
*======================================================
aa = 0 0 0 ;
bb = 12 0 0;
cc = 12 0.01 0.01 ;
dd = 12 0.01 -0.01 ;
ee = 12 -0.01 -0.01 ;
ff = 12 -0.01 0.01 ;
dens1 = 10 ;
lig1 = aa d dens1 bb ;
lig2 = cc d 2 dd d 2 ee d 2 ff d 2 cc ;
surf1 = surf plan lig2 ;
geo = lig1 et surf1 ;
elim 1.e-4 geo ;
*=====================================================
*definition des conditions aux limites
* encastrement à la base de la poutre et mouvement 
* dans le plan xOy
*=====================================================
cl1 = bloq rota depla aa ;
cl2 = bloq uz rx ry surf1 ;
cl3 = bloq bb uy ;
*======================================================
*      definition du modele, du materiau et 
*               des caracteristiques
*======================================================
mo1 = mode lig1 mecanique pout ;
mo2 = mode surf1 mecanique coq3 ;
ma1 = mate mo1 youn 2.1e7 nu 0.3 rho 7.8e-3 ;
ma2 = mate mo2 youn 2.1e10 nu 0.3 rho 7.8e-10 ;
ca1 = cara mo1 sect 0.15 inry 2.8125e-4 inrz 2.8125e-4
                                            tors 1.e-5;
ca2 = cara mo2 epai 0.01 ;
mat1 = ma1 et ca1 ;
mat2 = ma2 et ca2 ;
*
* Chargement de pression
  MOPRE = 'MODE' SURF1 'CHARGEMENT' 'PRESSION' 'COQ3' ;
  CHMZ = pres MOPRE pres (2050781.25*1.05) ;
*
  TTPS = 'TABL' ;
  TTPS. 0 = 0. ;
  TTPS. 1 = 0.00001 ;
  TTPS. 2 = 100. ;
  TPRE = 'TABL' ;
  TPRE. 0 = 0.*CHMZ ;
  TPRE. 1 = 1.*CHMZ ;
  TPRE. 2 = TPRE. 1 ;
  CHAA = 'CHAR' 'PRES' TTPS TPRE ;
*
mattotal = mat1 et mat2 ;
mototal = mo1 et mo2 et mopre ;
*======================================================
* calcul des matrices de rigidite et de masse
*======================================================
ri1 = rigi mototal mattotal ;
ritotal = ri1 et cl1 et cl2 ;
masst = masse mototal mattotal ;
*======================================================
* chargement par pression equivalente a la force 
* critique de flottement pondérée par un coefficient 
* de 1.05 pour se placer legerement au dessus de 
* l'instabilité.
*
* on impose un petit deplacement transverse pour 
* perturber le systeme.
*======================================================
fo2 = depi cl3 .01 ;
fo0 = manu chpo surf1 1 p (2050781.25*1.05) 
      nature discret ;
*======================================================
*               calcul non lineaire
* calcul statique sous la force de compression et la
* force transversale de perturbation
*======================================================
dt = 0.00001 ;
li1 = prog  0. dt ;
li2 = prog  0. 1 ;
ev = evol manu t (li1 ) f(t) (li2 )  ;
fco = char dimp  fo2 ev ;
lis1a = prog 0. pas (dt/5) dt ;
 
 
tab3 = table ;
tab3.grands_deplacements = vrai ;
tab3.chargement = fco et chaa ;
tab3.blocages_mecaniques = cl1  et cl2 et cl3 ;
tab3.modele =  mototal ;
tab3.caracteristiques = mattotal ;
tab3.temps_calcules =  lis1a ;
tab3.precision      =  1e-6  ;
tab3.deltaiter = 20 ;
tab3.maxdefor=1e-3;
tab3.k_sigma=faux;
*  raideur negative ==>
*tab3.'STABILITE' = faux;
*tab3.hypothese_deformations='LINEAIRE';
 
*optio veri 1;
pasapas tab3 ;
 
oeil = 0 0 10000 ;
dim1 = dime (tab3.deplacements) ;
lisx = prog  ; lisy = prog ;
 repeter bou1 dim1 ;
   n = &bou1 - 1 ;
   uy1 = extr tab3.deplacements.n 'UY' bb ;
   lisy = lisy et (prog uy1) ;
   lisx = lisx et (prog tab3.temps. n);
   def1 = defo tab3.deplacements.n  geo 100 ;
     si (&bou1 ega 1) ; deft = def1 ; sinon ;
     deft = deft et def1 ; finsi ;
 fin bou1 ;
ev1 = evol manu 'temps' lisx 'depl' lisy ;
si graph ; trac deft oeil cach anime ; finsi ;
 
 
*======================================================
*           calcul non lineaire dynamique
*======================================================
* definition de l'évolution du chargement 
* durée de la phase de chargement statique :dt
 
ev2 = evol manu (prog dt 100) (prog 1 1 ) ;
*
*fncons = char psui fo0 ev2 ;
 
si complet;
lis1c = prog dt pas .0003 .045 ;
sinon;
lis1c = prog dt  pas .001 0.005;
finsi;
ffbi= force (0 0 0) dd;
fbidon = 'CHAR' 'MECA' ffbi ev2;
tab3.'CHARGEMENT'= fbidon et chaa ;
vv = manu chpo ( surf1 et lig1 ) 1 'UX' 0.;
tab3.'VITESSES' = table;
nso= dime tab3.'TEMPS' - 1; 
tab3.'VITESSES'. nso = vv;
tab3.'ACCELERATIONS' = table;
tab3.'ACCELERATIONS'. nso = vv;
tab3.blocages_mecaniques = cl1  et cl2 ;
tab3.dynamique = vrai ;
tab3.temps_calcules = lis1c ;
tab3.hypothese_deformations='LINEAIRE';
 
pasapas tab3 ;
 
*====================================================
* exploitation des resultats : deplacement de l'
* extrémité de la poutre (point bb)
*====================================================
 
dim1 = dime (tab3.deplacements) ;
lisx = prog  ; lisy = prog ;
 repeter bou1 dim1 ;
   n = &bou1 - 1 ;
   uy1 = extr tab3.deplacements.n 'UY' bb ;
   lisy = lisy et (prog uy1) ;
   lisx = lisx et (prog tab3.temps. n);
 fin bou1 ;
ev1 = evol manu 'temps' lisx 'depl' lisy ; 
list lisy;
 si graph ; dess ev1  ; finsi ;
 
* la reference pour le test est le resultat d'une 
*             certaine version de C2000
dmax = maxi abs lisy ; list dmax ;
si complet ;
valtes = 5.8565;
sinon;
valtes =  .082;
finsi;
err1 = abs ((dmax - valtes)/valtes);
list err1 ;
 si (err1 > .05) ; erreur 5 ;
 sinon ; erre 0 ; finsi ;  
fin ;