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