Datasets examples

* fichier :  ex_proper.dgibi
GRAPH='N';                                                                      
SAUT PAGE ;    
*-----------------------------------------------------
* TEST  PLUS, MOINS, DEDU, TOUR
* OPERANDES MCHAML, MMODEL, CHPO, MAILLAGE (RIGIDITE)
*-----------------------------------------------------
*
opti dime 2 elem qua4 ;
*
p0 = 1. 0. ; e_x = p0 ; e_y = 0. 1. ; p_ori = 0. 0. ;
p1 = p0 plus e_x ;
s1 = (p0 d 1 p1) trans 1 e_y ;
lc1 = c (inve (s1 cote 4)) p0 p_ori 4 ;
cont = (d 6 p_ori p0)  et lc1 ;
su1 = surf cont plan ;
st1 = s1 et su1 ;
chp1 = manu chpo 1 s1  fy -1. ;
* tests PLUS, MOINS, DEDU, TOUR
m1 = mode st1 mecanique elastique ;
che1 = manu chml m1 sxx 1. ;
che3 = manu chml m1 smyy 2. smxy 1. ;
che5 = manu chml m1 smxx 5.66 eqeq 5. ;
m3 = mode (s1 plus e_x) mecanique elastique ;
m2 che2  = m3 che1 plus (2.*e_y) ;
m4 s3 = m1 s1 moins (2. 5.) ;
*list che2 ;
s2 = s1 plus (1. -1.) ;
chp10 = s2 moins s1 ;
*
chp21 = chp1 dedu trans s1 s2 ;
*
p3 = tourne p1 90. p_ori ;
m5 s4 = m1 s1 tourne 400. p0 ;
chp22 = chp1 dedu rota 400. p0 s1 s4 ;
*
chp5 = manu chpo 1 s1 scal 10. ;
chp6 = tour chp5 30. p1 ;
chp2 = tour chp1 180. p_ori ;
*
che4  = tour che3 90. p_ori ;
che6 = tour che5 45. p_ori ;
*
opti dime 3 elem cub8 ;
e_z = 0. 0. 1. ;
p5 = p1 tour 50. p_ori (p_ori plus e_z) ;
v1 = s1 volu 2 trans e_z ;
vt1 = st1 volu trans 1 e_z ;
vt2 = st1 volu 4 rota 45. p_ori (p_ori plus e_y) ;
mv1 = mode vt1 mecanique elastique ;
mv4 v3 = mv1 v1 moins (2. 5. 1.) ;
chp11 = v3 moins v1 ;
mv5 v4 = mv1 v1 tour 120. p0 (p0 plus e_z)  ;
vt6 = vt1 tour 120. p0 (p0 plus e_z)  ;
su7 = su1 dedu tran vt1 vt6 ;
su5 = su1 dedu 120. p0 (p0 plus e_z)  rota vt1 vt6 ;
vt5 = extr mv5 mail ;
su8 = su1 dedu 120. p0 (p0 plus e_z)  rota vt1 vt5  ;
*
mv6 che7 = mv1 che5 tour 50. p1 (p1 plus e_z) ;
*
t1 = extr che7 smyy 2 3 4 ;
tref = 3.3214 ;
resi = abs((t1 - tref)/tref) ;
si (resi