Datasets examples

* fichier :  tufi_relax.dgibi
OPTION ECHO 0;
opti dime 3 mode TRID elem seg2 ;
rap=10.;T1=10.;RMOY=T1 * rap;REXT1=RMOY + (T1 / 2.);
P1=0. 0. 0.;P2=0. 0. 0.;                                                        
TUFSS=MANU SEG2 P1 P2;                                                          
OBJTUFI=MODL TUFSS  MECANIQUE ELASTIQUE FLUAGE NORTON TUFI;
CARTOT=mate OBJTUFI YOUN 2.E5 NU 0.3
AF1 3.e-15 AF2 5. AF3 1. SMAX 200. RAYO REXT1 EPAI T1
VX 1. VY 0. VZ 0. VXF 0. VYF 0. VZF 1. ANGL 120;
CDL=BLOQ depla rota p1;
RIG1= RIGI OBJTUFI CARTOT;RIG = RIG1 ET CDL;
DD = RESO RIG (MOME 100.E5 'MY' P2);
cdli=bloq 'RY' P2;CDL=CDL et cdli;
RIG = RIG1 ET CDL;
vdepli=extr DD 'RY' P2;
depli=depi cdli vdepli;
DD = RESO RIG depli;
M0=react DD cdli;vM0=extr M0 'MY' P2;
CMM=vdepli / vM0;
SS = SIGMA CARTOT OBJTUFI DD;
EVT = EVOL MANU 'T' (PROG 0. 1.E10) 'F(T)' (PROG 1. 1.);
dt=10.;LIS = PROG 0. pas dt 110. ;
FDT=CHAR 'DIMP' depli EVT;                                                                
TENTR=TABLE;
TENTR.CARACTERISTIQUES=CARTOT;
TENTR.MODELE=OBJTUFI;
TENTR.CHARGEMENT=FDT;
TENTR.'CONTRAINTES'=tabl;
TENTR.'DEPLACEMENTS'=tabl;
TENTR.'CONTRAINTES' . 0 =SS;
TENTR.'DEPLACEMENTS' . 0 = DD;
TENTR.BLOCAGES_MECANIQUES=CDL;
TENTR.TEMPS_CALCULES= LIS;
pasapas TENTR;
t10=TENTR.TEMPS . 10;
U10=TENTR.DEPLACEMENTS . 10;
VAR10=TENTR.VARIABLES_INTERNES . 10;
M10=react U10 cdli;vM10=extr M10 'MY' P2;phi_el=CMM * vM10;
phi_f=vdepli - phi_el;
vC10=extr VAR10 'JPOI' 1 1 1;
phi_fb=extr VAR10 'EPSE' 1 1 1;
U9=TENTR.DEPLACEMENTS . 9;
Mi=react U9 cdli;vMi=extr Mi 'MY' P2;phi_el=CMM * vMi;
phi_f9=vdepli - phi_el;
U10=TENTR.DEPLACEMENTS . 10;
Mi=react U10 cdli;vMi=extr Mi 'MY' P2;phi_el=CMM * vMi;
phi_f10=vdepli - phi_el;
U11=TENTR.DEPLACEMENTS . 11;
Mi=react U11 cdli;vMi=extr Mi 'MY' P2;phi_el=CMM * vMi;
phi_f11=vdepli - phi_el;
phi_fpt1=(phi_f10 - phi_f9) / dt;
phi_fpt2=(phi_f11 - phi_f10) / dt;
phi_fpt=(phi_fpt1 + phi_fpt2) / 2.;

* COMPARAISON A LA SOLUTION ANALYTIQUE
Mth10=7.747E6;phi_f10=1.850E-4;vCth10=4.576E-3;
ph_fpt10=1.020E-6;

errM=abs ((vM10 - Mth10) / Mth10);
errphif=abs ((phi_f - phi_f10) / phi_f10);
errphifb=abs ((phi_fb - phi_f10) / phi_f10);
errC=abs ((vC10 - vCth10) / vCth10);
errphfp=abs ((phi_fpt - ph_fpt10) / ph_fpt10);                    
err=errM * errphif * errphifb * errC * errphfp;
SI (ERR < (1.E-3 ** 5))ERRE 0;SINO;ERRE 5;FINSI;
fin;