* fichier : joi1_coulomb.dgibi
'OPTI' 'ECHO' 1 ;
'OPTI' 'ELEM' 'SEG2' 'DIME' 3 'MODE' 'TRID' ;
'OPTI' 'EPSILON' 'LINEAIRE' ;
GRAPH = FAUX ;
P1 = 0 0 0 ;
P2 = 0 0 0 ;
joint = P1 P2 'DROIT' 1 ;
************************************************************************
* Début: MODELE *
************************************************************************
mo = 'MODL' joint 'MECANIQUE' 'ELASTIQUE' 'ORTHOTROPE'
'PLASTIQUE' 'COULOMB' 'JOI1' ;
************************************************************************
* Fin: MODELE *
************************************************************************
************************************************************************
* Début: MATERIAU *
************************************************************************
kn = 100. ;
ks1 = 50. ;
ks2 = ks1 ;
qn = 40. ;
qs1 = 20. ;
qs2 = 30. ;
*
ang = 45. ;
cosa = 'COS' ang ;
sina = 'SIN' ang ;
dir1 = cosa sina 0 ;
dir2 = (-1*sina) cosa 0 ;
*
phi = 45. ;
ma = 'MATE' mo 'DIRECTION' dir1 dir2
'KN ' kn 'KS1 ' ks1 'KS2 ' ks2
'QN ' qn 'QS1 ' qs1 'QS2 ' qs2
'TYPE' 1. 'FRIC' phi
'EF ' kn 'ECN ' 1E6 'COHE' 0. ;
************************************************************************
* Fin: MATERIAU *
************************************************************************
************************************************************************
* Début: CL & CHARGEMENT *
************************************************************************
cl1 = 'BLOQ' 'DEPL' 'DIRECTION' dir1 P1 ;
cl2 = 'BLOQ' 'DEPL' 'DIRECTION' dir2 P1 ;
cl3 = 'BLOQ' 'UZ' 'RX' 'RY' 'RZ' joint ;
cl4 = 'BLOQ' 'DEPL' 'DIRECTION' dir1 P2 ;
cl5 = 'BLOQ' 'DEPL' 'DIRECTION' dir2 P2 ;
cl = cl1 'ET' cl2 'ET' cl3 'ET' cl4 'ET' cl5 ;
vdpi = 0.1 ;
dpi1 = 'DEPI' cl4 vdpi ;
dpi2 = 'DEPI' cl5 vdpi ;
lt = 'PROG' 0 1 2 3 4 5 6 ;
lc1 = 'PROG' 0 -1 -1 -1 -0.5 -0.01 1 ;
ev1 = 'EVOL' 'MANU' 't' lt 'coef' lc1 ;
lc2 = 'PROG' 0 0 1 3 3 3 3 ;
ev2 = 'EVOL' 'MANU' 't' lt 'coef' lc2 ;
cha1 = 'CHAR' 'DIMP' dpi1 ev1 ;
cha2 = 'CHAR' 'DIMP' dpi2 ev2 ;
cha = cha1 'ET' cha2 ;
************************************************************************
* Fin: CL & CHARGEMENT *
************************************************************************
lis_t = 'PROG' 1 2 2.5 3 4 5 5.0098 6 ;
tab = 'TABLE' ;
tab.'MODELE' = mo ;
tab.'CARACTERISTIQUES' = ma ;
tab.'BLOCAGES_MECANIQUES' = cl ;
tab.'CHARGEMENT' = cha ;
tab.'TEMPS_CALCULES' = lis_t ;
tab.'TEMPS_SAUVES' = lis_t ;
tab.'FTOL' = 1E-3 ;
PASAPAS tab ;
time=tab.'TEMPS';
ntime = 'DIME' time ;
n=ntime - 1 ;
progt = 'PROG' ;
prN = 'PROG' ;
prFY = 'PROG' ;
prepse = 'PROG' ;
pretat = 'PROG' ;
prepou = 'PROG' ;
prlam1 = 'PROG' ;
prpla1 = 'PROG' ;
prpla3 = 'PROG' ;
************************************************************************
* Début: CONTROLE *
************************************************************************
REPETER contro ntime ;
i = &contro - 1 ;
t = time.i ;
progt = progt 'ET' ('PROG' t) ;
*-
sig = tab.'CONTRAINTES'.i ;
prN = prN 'ET' ('PROG' ('EXTR' sig 'EFFX' 1 1 1)) ;
prFY = prFY 'ET' ('PROG' ('EXTR' sig 'EFFY' 1 1 1)) ;
*-
vi = tab.'VARIABLES_INTERNES'.i ;
prepse = prepse 'ET' ('PROG' ('EXTR' vi 'EPSE' 1 1 1)) ;
pretat = pretat 'ET' ('PROG' ('EXTR' vi 'STAT' 1 1 1)) ;
prepou = prepou 'ET' ('PROG' ('EXTR' vi 'EPOU' 1 1 1)) ;
prlam1 = prlam1 'ET' ('PROG' ('EXTR' vi 'LAM1' 1 1 1)) ;
prpla1 = prpla1 'ET' ('PROG' ('EXTR' vi 'PLA1' 1 1 1)) ;
prpla3 = prpla3 'ET' ('PROG' ('EXTR' vi 'PLA3' 1 1 1)) ;
FIN contro ;
************************************************************************
* Début: CONTROLE *
************************************************************************
evN = evol manu 't' progt 'N' prN ;
evFY = evol manu 't' progt 'FY' prFY ;
evepse = evol manu 't' progt 'EPSE' prepse ;
evetat = evol manu 't' progt 'STAT' pretat ;
evepou = evol manu 't' progt 'EPOU' prepou ;
evpla1 = evol manu 't' progt 'PLA1' prpla1 ;
evpla3 = evol manu 't' progt 'PLA3' prpla3 ;
evlam1 = evol manu 't' progt 'LAM1' prlam1 ;
mu = ('SIN' phi) / ('COS' phi) ;
*-
* 0
© Cast3M 2003 - All rights reserved.
Disclaimer