Télécharger la plaquette Cast3M


fr en

Télécharger liai_ar1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  liai_ar1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *     Test liai_ar1.dgibi: Jeux de donnees            
  6. *     ------------------------------------
  7. **********************************************************         
  8. *                TEST DE VALIDATION                      *
  9. *                ------------------                      *
  10. *              MODELE PLASTIQUE BILIN_EFFX               *
  11. *       SIMULATION DE LA LIAISON CRAYON/GRILLE           *
  12. *                                                        *
  13. *  MAILLAGE:                                             *
  14. *    CRAYON                   = POUTRES ELASTIQUES       *
  15. *    1 CELLULE DE GRILLE      = PLAQUE EQUIVALENTE       *
  16. *                               ELASTIQUE EN COQ4        *
  17. *    1 LIAISON CRAYON/CELLULE = 1 POUTRE AVEC MODELE     *
  18. *                               PLASTIQUE BILIN_EFFX     *
  19. *                                                        *
  20. *  CHARGEMENT:                                           *
  21. *     FORCE DE TRACTION AXIALE SUR L'EXTREMITE           *
  22. *     SUPERIEURE DU CRAYON: CHARGE PUIS DECHARGE         *
  23. *                                                        *
  24. *  CONDITIONS AUX LIMITES:                               *
  25. *     CELLULE DE GRILLE ENCASTREE SUR SES 4 COTES        *
  26. *                                                        *
  27. *  COMPARAISON / SOLUTION ANALYTIQUE:                    *
  28. *     RAIDEUR AU CENTRE EN FLEXION DE LA                 *
  29. *     PLAQUE EQUIVALENTE PAR CALCUL ANNEXE               *
  30. ********************************************************** 
  31.  
  32. * ::::::::::::::::::::::::::::::::::::::::::::::::
  33. * DEFINITION DES PROCEDURES POUR POST TRAITEMENT :
  34. * ::::::::::::::::::::::::::::::::::::::::::::::::
  35. 'DEBP' FEPS   mo*'MMODEL'    mat*'MCHAML'    n*'ENTIER'
  36.              dep*'TABLE'     sig*'TABLE'
  37.             mot1*'MOT'      mot2*'MOT'       P*'POINT'                 ;
  38. ************************************************************************
  39. prx = prog                                                             ;
  40. pry = prog                                                             ;
  41.  
  42. repe stru n                                                            ;
  43. *-
  44.  i   = &stru-1                                                         ; 
  45.  eps = 'EPSI' mo mat dep.i 'LINE'                                      ;
  46.  eps = 'CHAN' 'CHPO' mo eps                                            ;
  47.  x   = 'EXTR' eps mot1 P                                               ;
  48.  prx = prx 'ET' (prog x)                                               ;
  49.  sg  = sig.i                                                           ;
  50.  sg  = 'CHAN' 'CHPO' mo sg                                             ;
  51.  y   = 'EXTR' sg mot2 P                                                ;
  52.  pry = pry 'ET' (prog y)                                               ;
  53.  ev  = 'EVOL' 'MANU' mot1 prx mot2 pry                                 ;
  54. *-
  55. fin stru                                                               ;
  56. ************************************************************************
  57. 'FINP' ev                                                              ;
  58. ************************************************************************
  59.  
  60. ************************************************************************
  61. 'DEBP' FTEMPS   mo*'MMODEL'    mat*'MCHAML'    n*'ENTIER'
  62.               tab1*'TABLE'    tab2*'TABLE'
  63.               mot1*'MOT'                       P*'POINT'               ;
  64. ************************************************************************
  65. prx = prog                                                             ;
  66. pry = prog                                                             ;
  67.  
  68. repe stru n                                                            ;
  69. *-
  70.  i   = &stru-1                                                         ;
  71.  x   = tab1.i                                                          ;
  72.  prx = prx 'ET' (prog x)                                               ;
  73.  y   = 'EXTR' tab2.i  mot1  P                                          ;
  74.  pry = pry 'ET' (prog y)                                               ;
  75.  ev  = 'EVOL' 'MANU' 'temps' prx mot1 pry                              ;
  76. *-
  77. fin stru                                                               ;
  78. ************************************************************************
  79. 'FINP' ev                                                              ;
  80. ************************************************************************
  81.  
  82. ************************************************************************
  83. 'DEBP' SOLEFFX
  84.          F*'FLOTTANT' Fseuil*'FLOTTANT'    k*'FLOTTANT' 
  85.         l1*'FLOTTANT'     A1*'FLOTTANT'   E1*'FLOTTANT'  E1T*'FLOTTANT'
  86.         l2*'FLOTTANT'     A2*'FLOTTANT'   E2*'FLOTTANT'   Fa/'FLOTTANT';
  87. ************************************************************************
  88. k1  = E1  * A1 ;
  89. k1T = E1T * A1 ;
  90. k2  = E2  * A2 ;
  91.  
  92. pr1 = prog 0   l1     ; pr1un = prog (dime pr1)*1 ;
  93. pr2 = prog l1 (l1+l2) ; pr2un = prog (dime pr2)*1 ;
  94.  
  95. 'SI' ('NON' ('EXIS' Fa))                                               ;
  96. *- Charge
  97.    'SI' (F '<EG' Fseuil)                                               ;
  98. *- Solution élastique pour la liaison
  99.    epsg = F/k1                                                         ;
  100.    u2b  = (F*l1)/k1                                                    ;
  101.    'SINON'                                                             ;
  102. *- Solution plastique pour la liaison
  103.    eps1 = Fseuil*(1-(E1T/E1))                                          ;
  104.    epsg = (F - eps1) / k1T                                             ;
  105.    u2b  = epsg*l1                                                      ;
  106.    'FINSI'                                                             ;
  107. 'SINON'                                                                ;
  108. *- Décharge
  109.    eps1 = (Fa-Fseuil)*(1-(E1/E1T))                                     ;
  110.    epsg = (F - eps1) / k1                                              ;
  111.    u2b  = epsg*l1                                                      ;
  112. 'FINSI'                                                                ;
  113.  
  114. u1  = epsg*pr1 + (F *(pr1un/k))                                        ;
  115. ev1 = 'EVOL' 'MANU' 'ABSC' pr1 'DEPL' u1                               ;
  116.  
  117. u2a = (F/k2) * (pr2 - (l1*pr2un))                                      ;
  118. u2  = u2a + ((u2b +(F/k))*pr2un)                                       ;
  119. ev2 = 'EVOL' 'MANU' 'ABSC' pr2 'DEPL' u2 ;
  120.  
  121. ev = 'CONCAT' ev1 ev2 ;
  122.  
  123. ************************************************************************
  124. 'FINP' ev                                                              ;
  125. ************************************************************************
  126. * :::::::::::::
  127. * DEBUT DU PB :
  128. * :::::::::::::
  129.  
  130. 'OPTI' 'DIME' 3     ;
  131.  
  132. *'OPTI' 'TRAC' 'PSC' ;
  133. 'TEMPS' 'ZERO'      ;
  134.  
  135. DESSI = faux ;
  136.  
  137. lcel = 12.595E-3 ;
  138. long = 0.15      ;
  139. l1   = long / 10 ;
  140.  
  141. ************
  142. *- CRAYON -*
  143. ************
  144. 'OPTI' 'ELEM' 'SEG2' ;
  145.  
  146. p1 = 0  0     0      ;
  147. p3 = 0  0    l1      ;
  148. p2 = 0  0  long      ;
  149. q3 = 0  0  (-1*l1)   ;
  150. q2 = 0  0  (-1*long) ;
  151.  
  152. nelt = 10      ;
  153. n1   =  1      ;
  154. n2   = nelt-n1 ;
  155.  
  156. geo1 = p1 p3  'DROI' n1 ;
  157. geo2 = p3 p2  'DROI' n2 ;
  158. geo3 = p1 q3  'DROI' n1 ;
  159. geo4 = q3 q2  'DROI' n2 ;
  160. geo  = geo1 'ET' geo2 'ET' geo3 'ET' geo4   ;
  161.  
  162. ************
  163. *- GRILLE -*
  164. ************
  165. 'OPTI' 'ELEM' 'QUA4'         ;
  166.  
  167. r1 = (-1*lcel)  (-1*lcel)  0 ;
  168. r2 = lcel       (-1*lcel)  0 ;
  169. r3 = lcel       lcel       0 ;
  170. r4 = (-1*lcel)  lcel       0 ;
  171.  
  172. n3 = 4 ;
  173.  
  174. lplac1 = r1 r2 'DROIT' n3 ;
  175. lplac2 = r2 r3 'DROIT' n3 ;
  176. lplac3 = r3 r4 'DROIT' n3 ;
  177. lplac4 = r4 r1 'DROIT' n3 ;
  178. plac   = lplac1 lplac2 lplac3 lplac4 'DALLER' 'PLAN' ;
  179. r0     = plac 'POIN' 'PROC' p1                       ;
  180.  
  181. *************
  182. *- LIAISON -*
  183. *************
  184. 'OPTI' 'ELEM' 'SEG2'   ;
  185.  
  186. lien = 'COUL' (r0 p3 'DROIT' 1) 'ROUG' ;
  187.  
  188. ****************
  189. *- ASSEMBLAGE -*
  190. ****************
  191. mailtot = geo 'ET' plac 'ET' lien ;
  192.  
  193. *************
  194. *- MODELES -*
  195. *************
  196. moc = 'MODE' geo  'MECANIQUE' 'ELASTIQUE'                      'POUTRE';
  197. mop = 'MODE' plac 'MECANIQUE' 'ELASTIQUE'                      'COQ4'  ;
  198. mol = 'MODE' lien 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE'
  199.                                     'BILIN_EFFX' 'POUTRE';
  200. mo  = moc 'ET' mop 'ET' mol                                            ;
  201.  
  202.  
  203. youn_M5 = 0.984E11        ;
  204. nu_M5   = 0.368           ;
  205. yp      = youn_M5 / 10000.;
  206. fseuilx = 6.              ;
  207.  
  208. mac = 'MATE' moc 'YOUN' youn_M5  'NU' nu_M5                           ;
  209. map = 'MATE' mop 'YOUN' youn_M5  'NU' nu_M5                           ;
  210. mal = 'MATE' mol 'YOUN' youn_M5  'NU' nu_M5  'ETAN' yp 'YIEF' fseuilx ;
  211. ma  = mac 'ET' map 'ET' mal                                           ;
  212.  
  213. *ep_plaq = 0.425E-3  ;
  214. ep_plaq = 2.879E-3  ;
  215. ep_cray = 0.57E-3   ;
  216. Dext    = 9.5E-3    ;
  217. Dint    = 8.36E-3   ;
  218.  
  219. air = (pi/4)  * ((Dext**2)-(Dint**2)) ;
  220. iy  = (pi/64) * ((Dext**4)-(Dint**4)) ;
  221. iz  = iy                              ;
  222. ix  = 2*iy                            ;
  223. sr1 = (Dint/Dext)/(1+((Dint/Dext)**2));
  224. sr2 = 6./(7.+(20.*(sr1**2)))          ;
  225. sry = sr2*air                         ;
  226. srz = sry                             ;
  227.  
  228. cac = 'CARA' moc 'TORS' ix 'INRY' iy 'INRZ' iz 'SECT' air 
  229.                                                'SECY' sry 'SECZ' srz ;
  230. cap = 'CARA' mop 'EPAI' ep_plaq                                      ;
  231. cal = 'CARA' mol 'TORS' ix 'INRY' iy 'INRZ' iz 'SECT' air 
  232.                                                'SECY' sry 'SECZ' srz ;
  233. ca  = cac 'ET' cap 'ET' cal                                          ;
  234.  
  235. mat = ma 'ET' ca                                                     ;
  236.  
  237.  
  238. cl = bloq 'DEPL' 'ROTA' (lplac1 et lplac2 et lplac3 et lplac4) ;
  239.  
  240. valfo = 1.                                             ;
  241. fo    = 'FORCE' 'FZ' valfo  p2                         ;
  242. evfo  = 'EVOL' 'MANU' ('PROG' 0 12 24) ('PROG' 0 12 0) ; 
  243. cha   = 'CHAR' fo evfo                                 ;
  244.  
  245. lisc = 'PROG' 1 pas 1 24;
  246. *lisc = 'PROG' 1 pas 1 6 pas 1 7 pas 1 12 pas 1 13 pas 1  24 ;
  247. *pas 0.1 13 pas 1 24 ; yp = youn_M5 / 100. ;
  248.  
  249. tab = 'TABLE' ;
  250. tab.'MODELE'              = mo                                        ;
  251. tab.'CARACTERISTIQUES'    = mat                                       ;
  252. tab.'CHARGEMENT'          = cha                                       ;
  253. tab.'BLOCAGES_MECANIQUES' = cl                                        ;
  254. *tab.deltaiter= 20;
  255. tab.maxdefor=5.e-5;
  256. tab.'TEMPS_CALCULES'      = lisc                                      ;
  257. tab.'TEMPS_SAUVES'        = lisc                                      ;
  258. tab.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  259.  
  260. PASAPAS tab ;
  261.  
  262. time  = tab.TEMPS   ;
  263. ntime = 'DIME' time ;
  264. n     = ntime-1     ;
  265.  
  266. dep = tab.'DEPLACEMENTS'       ;
  267. sig = tab.'CONTRAINTES'        ;
  268. vi  = tab.'VARIABLES_INTERNES' ;
  269.  
  270. titr 'P2'                                               ;
  271. FXEPS_P2 = FEPS   mo mat ntime dep  sig 'EPS' 'EFFX' P2 ;
  272. FXt_P2   = FTEMPS mo mat ntime time dep 'UZ'         P2 ;
  273. 'SI' DESSI; dess FXEPS_P2 ; dess FXt_P2 ; 'FINSI'       ;
  274.  
  275. titr 'P3'                   ;
  276. FXEPS_P3 = FEPS   mo mat ntime dep  sig 'EPS' 'EFFX' P3 ;
  277. FXt_P3   = FTEMPS mo mat ntime time dep 'UZ'         P3 ;
  278. 'SI' DESSI; dess FXEPS_P3 ; dess FXt_P3 ; 'FINSI'       ;
  279.  
  280. titr 'R0'                   ;
  281. FXEPS_R0 = FEPS   mo mat ntime dep  sig 'EPS' 'EFFX' R0 ;
  282. FXt_R0   = FTEMPS mo mat ntime time dep 'UZ'         R0 ;
  283. 'SI' DESSI; dess FXEPS_R0 ; dess FXt_R0 ; 'FINSI'       ;
  284.  
  285. titr 'Q2'                                               ;
  286. FXEPS_Q2 = FEPS   mo mat ntime dep  sig 'EPS' 'EFFX' Q2 ;
  287. FXt_Q2   = FTEMPS mo mat ntime time dep 'UZ'         Q2 ;
  288. 'SI' DESSI; dess FXEPS_Q2 ; dess FXt_Q2 ; 'FINSI'       ;
  289.  
  290. titr 'Q3'                   ;
  291. FXEPS_Q3 = FEPS   mo mat ntime dep  sig 'EPS' 'EFFX' Q3 ;
  292. FXt_Q3   = FTEMPS mo mat ntime time dep 'UZ'         Q3 ;
  293. 'SI' DESSI; dess FXEPS_Q3 ; dess FXt_Q3 ; 'FINSI'       ;
  294.  
  295. tdes         = 'TABLE'        ;
  296. tdes.'TITRE' = 'TABLE'        ;
  297. tdes . 1     = 'MARQ ETOI'    ;
  298. tdes . 2     = 'MARQ CARR'    ;
  299. tdes . 3     = 'MARQ TRIA'    ;
  300. tdes . 4     = 'MARQ TRIB'    ;
  301. tdes.'TITRE' . 1 = 'MOT' 'P2' ;
  302. tdes.'TITRE' . 2 = 'MOT' 'P3' ;
  303. tdes.'TITRE' . 3 = 'MOT' 'Q2' ;
  304. tdes.'TITRE' . 4 = 'MOT' 'Q3' ;
  305.  
  306. titr 'UX(t)'                      ;
  307. 'SI' DESSI ; dess (FXt_P2 et FXt_P3 et FXt_Q2 et FXt_Q3) LEGE tdes ;
  308. 'FINSI'    ;
  309.  
  310. ni = 6                                         ;
  311. titr 'UX le long de la poutre t=' time.ni      ;
  312. evUXgeo = 'EVOL' 'CHPO' dep.ni 'UZ' (inve geo) ;
  313. 'SI' DESSI; dess evUXgeo ; 'FINSI'             ;
  314.  
  315. geolien = lien et geo2                                        ;
  316. titr 'UX le long de la liaison t=' time.ni                    ;
  317. evUXlien = 'EVOL' 'CHPO' dep.ni 'UZ' geolien                  ;
  318. evUXliey = extr evUXlien 'ORDO'                               ;
  319. evUXliex = (extr evUXlien 'ABSC')+(prog (dime evUXliey)*long) ;
  320. evUXlie = 'EVOL' 'MANU' evUXliex evUXliey                     ;
  321. 'SI' DESSI; dess (evUXlie et evUXgeo) ; 'FINSI'               ;
  322.  
  323. def  = 'DEFO' dep.ni mailtot 100. 'ROUG' ;
  324. def0 = 'DEFO' dep.ni mailtot   0.        ;
  325. 'SI' DESSI; trac (def0 et def) ; 'FINSI' ;
  326.  
  327. l2 = long-l1              ;
  328. k1 = (youn_M5 * air) / l1 ;
  329. k2 = (youn_M5 * air) / l2 ;
  330. k  = k1 / 2.              ;
  331.  
  332. *- pour ep_plaq = 0.425E-3
  333. *k  = 1./(4.26698E-6) ;
  334.  
  335. PREXABSO = 1E-6 ; 
  336.  
  337. pinstant = 'PROG' 6. 7. 12. 12.5 18. 24.         ;
  338. instants = 'LECT'                                ;
  339. 'REPE' boc00 ('DIME' pinstant)                   ;
  340.  tinstant = 'EXTR' pinstant &boc00               ; 
  341.  'REPE' boc0 ntime                               ;
  342.     nt = &boc0 - 1                               ;
  343.     t  = time.nt                                 ;
  344.     'SI' (t 'NEG' 0.)                            ;
  345.        test = (time.nt - tinstant) / tinstant    ;
  346.        test = 'ABS' test                         ;
  347.      'SI' (test '<EG' 1E-4)                      ;
  348.         instants = instants 'ET' ('LECT' nt)     ;       
  349.      'FINSI'                                     ;
  350.     'FINSI'                                      ;
  351.  'FIN' boc0                         ;
  352. 'FIN' boc00                         ;
  353.  
  354. 'REPE' boc1 ('DIME' instants)                   ;
  355.  
  356.    ti = 'EXTR' instants &boc1                   ;
  357.    t  = time.ti                                 ;
  358.    ft = 'MAXI' ('TIRE' cha t)                   ;
  359.    'SI' (t '<EG' 12.)                           ;
  360.    evUXsol = SOLEFFX     ft    fseuilx    k  
  361.                   l1     air   youn_M5   yp
  362.                   l2     air   youn_M5          ;
  363.    fa=ft                                        ;
  364.    'SINON'                                      ;
  365.    evUXsol = SOLEFFX     ft    fseuilx    k  
  366.                   l1     air   youn_M5   yp
  367.                   l2     air   youn_M5   fa     ;
  368.    'FINSI'                                      ;
  369.    evUXlien = 'EVOL' 'CHPO' dep.ti 'UZ' geolien ;
  370.  
  371.    evUXsolx = 'EXTR' evUXsol  'ABSC'            ;
  372.    evUXsoly = 'EXTR' evUXsol  'ORDO'            ;
  373.    evUXlenx = 'EXTR' evUXlien 'ABSC'            ;
  374.    evUXleny = 'EXTR' evUXlien 'ORDO'            ;
  375.    evUXleny = 'IPOL' evUXsolx evUXlenx evUXleny ;
  376.    test     = 'ABS' (evUXsoly - evUXleny)       ;
  377.    test     = 'MAXI' test                       ;
  378.    'SI' ('<' test PREXABSO)                     ;
  379.       'ERRE' 0                                  ;
  380.    'SINON'                                      ;
  381.    'MESS' 'ERREUR DEPLACEMENT AXIAL a t= ' t 
  382.         test   ' > ' PREXabso                   ;
  383. mess ' evuxsoly';       list evuxsoly;
  384. mess ' evuxleny';       list evuxleny;
  385. mess ' evuxsolx';       list evuxsolx;
  386. *      'ERRE' 5                                  ;
  387.    'FINSI'                                      ;
  388.  
  389.  
  390.    tdes         = 'TABLE'                       ;
  391.    tdes.'TITRE' = 'TABLE'                       ;
  392.    tdes . 1     = 'MARQ ETOI'                   ;
  393.    tdes . 2     = 'MARQ CARR'                   ;
  394.    tdes.'TITRE' . 1 = 'MOT' 'Solution'          ;
  395.    tdes.'TITRE' . 2 = 'MOT' 'Calcul'            ;        
  396.    titr 'UX le long de la poutre t=' t          ;
  397.    'SI' DESSI;
  398.       dess (evUXsol et evUXlien) tdes LEGE      ;
  399.    'FINSI'                                      ;
  400.  
  401. 'FIN'  boc1                                     ;
  402.  
  403. 'FIN' ;
  404.  
  405.  
  406.  
  407.  
  408.  
  409.  
  410.  
  411.  
  412.  
  413.  
  414.  
  415.  
  416.  
  417.  
  418.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales