Télécharger pilotage_indirect_1.dgibi

Retour à la liste

Numérotation des lignes :

  1. *+++++++++++++++++++++++++++++++++++++++++
  2. *
  3. * Exemple: perforated plate under traction
  4. * Several elements with MAZARS model in the middle
  5. * Testing new version of PASAPAS/PILOINDI
  6. * -> Discplacement control at the right edge <-
  7. * -> Path equation: control maximum strain increment <-
  8. *
  9. *+++++++++++++++++++++++++++++++++++++++++
  10. 'TITR' 'Displacement applied at the right edge';
  11. 'OPTI' 'DIME' 2 'ELEM' 'TRI3' 'MODE' 'PLAN CONTR';
  12. 'OPTI' 'EPSI' 'LINEAIRE';
  13. 'OPTI' 'TRAC' 'PSC';
  14. *+++++++++++++++++++++++++++++++++++++++++
  15. *
  16. * USER DEFINED PROCEDURES
  17. *
  18. *+++++++++++++++++++++++++++++++++++++++++
  19. 'DEBP' PILOINDI PRECED*'TABLE' U*'CHPOINT' DU*'CHPOINT' DUI*'CHPOINT' DUII*'CHPOINT' DTAU*'FLOTTANT';
  20.  
  21. WTAB=PRECED.'WTABLE';
  22. MODTOT = WTAB.'MO_TOT';
  23.  
  24. * --- equation de pilotage indirecte --------
  25. CHATAU = 'MANU' 'CHML' MODTOT 'EPXX' DTAU 'TYPE' 'DEFORMATIONS' 'STRESSES';
  26.  
  27. DEF1 = 'EPSI' MODTOT duI ;
  28. DEF2 = 'EPSI' MODTOT duII ;
  29. DEFT = 'EPSI' MODTOT dU;
  30.  
  31. EPS_TOT = 'EXCO' 'EPXX' DEFT;
  32. def_I = 'EXCO' 'EPXX' DEF1;
  33. CH_A_1 = 'EXCO' 'EPXX' DEF2;
  34.  
  35. CH_A_0 = EPS_TOT + def_I;
  36. CH_AUX2 = CH_A_1**(-1.0);
  37. CH_ETA = (CHATAU - CH_A_0)*CH_AUX2;
  38.  
  39. MA1 = 'MASQ' CH_A_1 'EGSUPE' 0. ;
  40. CH_LAMB_MAX = (MA1*CH_ETA);
  41. MA2 = 'MASQ' CH_A_1 'INFERIEUR' 0. ;
  42. CH_LAMB_MIN = (MA2*CH_ETA);
  43.  
  44. big_value = 1.0E10;
  45. CH_LAMB_MAX = CH_LAMB_MAX + (big_value*MA2);
  46. CH_LAMB_MIN = CH_LAMB_MIN - (big_value*MA1);
  47.  
  48. LIM_INF = 'MAXI' CH_LAMB_MIN 'AVEC' (MOTS 'EPXX');
  49. LIM_SUP = 'MINI' CH_LAMB_MAX 'AVEC' (MOTS 'EPXX');
  50.  
  51. D_eta= LIM_SUP;
  52.  
  53. 'FINPROC' D_ETA ;
  54.  
  55. 'DEBPROC' RESFORCE CH1*CHPOINT GEO*MAILLAGE COMP1*MOT ;
  56. CH2 = REDU CH1 GEO ;
  57. CHTOT = RESU CH2 ;
  58. P1 = (EXTR CHTOT MAIL) POIN 1 ;
  59. TOT1 = EXTR CHTOT COMP1 P1 ;
  60. 'FINPROC' TOT1 ;
  61.  
  62. *+++++++++++++++++++++++++++++++++++++++++
  63. *
  64. * END OF USER DEFINED PROCEDURES
  65. *
  66. *+++++++++++++++++++++++++++++++++++++++++
  67.  
  68. *+++++++++++++++++++++++++++++++++++++++++
  69. *
  70. * PARAMETRIC GEOMETRY
  71. *
  72. *+++++++++++++++++++++++++++++++++++++++++
  73. pig = 3.141592654;
  74. R = 0.5;
  75. e = (pig/80.)*R;
  76. e1 = (pig/80.)*R;
  77. CC = 1. 0.85;
  78.  
  79. DTETA = ((e/R)/pig)*180.;
  80. SINDT = SIN(DTETA);
  81. COSDT = COS(DTETA);
  82.  
  83. DTETA1 = ((e1/R)/pig)*180.;
  84. SINDT1 = SIN(DTETA1);
  85. COSDT1 = COS(DTETA1);
  86.  
  87. dens1 = 1./5.;
  88. dens2 = 1./8.;
  89. dens3 = 2.*R*(SINDT);
  90. dens4 = 2.*R*(SINDT1);
  91.  
  92. PA = 0. 0.;
  93. PB = 2. 0.;
  94. PC = 2. 2.;
  95. PD = 0. 2.;
  96.  
  97. PE = (1.-(R*SINDT)) 0.;
  98. PF = (1.+(R*SINDT)) 0.;
  99. PG = (1.+(R*SINDT1)) 2.;
  100. PH = (1.-(R*SINDT1)) 2.;
  101.  
  102. PI = 2. 0.85;
  103. PL = 0. 0.85;
  104.  
  105. PC1 = (1.+(R*SINDT)) (0.35+(R*(1-COSDT)));
  106. PC2 = 1.5 0.85;
  107. PC3 = (1.+(R*SINDT1)) (1.35-(R*(1-COSDT1)));
  108. PC4 = (1.-(R*SINDT1)) (1.35-(R*(1-COSDT1)));
  109. PC5 = 0.5 0.85;
  110. PC6 = (1.-(R*SINDT)) (0.35+(R*(1-COSDT)));
  111.  
  112. LAE = DROI PA PE 'DINI' dens1 'DFIN' dens3;
  113. LC6E = DROI PE PC6 'DINI' dens3 'DFIN' dens3;
  114. LC5C6 = CERC PC6 CC PC5 'DINI' dens3 'DFIN' dens2;
  115. LC5L = DROI PC5 PL 'DINI' dens2 'DFIN' dens1;
  116. LLA = DROI PL PA 'DINI' dens1 'DFIN' dens1;
  117. LSW = LAE ET LC6E ET LC5C6 ET LC5L ET LLA;
  118. SSW = SURF LSW;
  119.  
  120. LEF = DROI 1 PE PF;
  121. LC1F = DROI PF PC1 'DINI' dens3 'DFIN' dens3;
  122. LC6C1 = CERC PC1 CC PC6 'DINI' dens3 'DFIN' dens3;
  123. LS = LEF ET LC1F ET LC6C1 ET (INVE LC6E);
  124. SS = SURF LS;
  125.  
  126. LFB = DROI PF PB 'DINI' dens3 'DFIN' dens1;
  127. LBI = DROI PB PI 'DINI' dens1 'DFIN' dens1;
  128. LC2I = DROI PI PC2 'DINI' dens1 'DFIN' dens2;
  129. LC2C1 = CERC PC2 CC PC1 'DINI' dens2 'DFIN' dens3;
  130. LSE = LFB ET LBI ET LC2I ET LC2C1 ET (INVE LC1F);
  131. SSE = SURF LSE;
  132.  
  133. LIC = DROI PI PC 'DINI' dens1 'DFIN' dens1;
  134. LCG = DROI PC PG 'DINI' dens1 'DFIN' dens4;
  135. LGC3 = DROI PG PC3 'DINI' dens4 'DFIN' dens4;
  136. LC3C2 = CERC PC3 CC PC2 'DINI' dens4 'DFIN' dens2;
  137. LNE = (INVE LC2I) ET LIC ET LCG ET LGC3 ET LC3C2;
  138. SNE = SURF LNE;
  139.  
  140. LC3C4 = CERC PC3 CC PC4 'DINI' dens4 'DFIN' dens4;
  141. LGH = DROI 1 PG PH;
  142. LHC4 = DROI PH PC4 'DINI' dens4 'DFIN' dens4;
  143. LN = LC3C4 ET (INVE LGC3) ET LGH ET LHC4;
  144. SN = SURF LN;
  145.  
  146. LC4C5 = CERC PC5 CC PC4 'DINI' dens2 'DFIN' dens4;
  147. LHD = DROI PH PD 'DINI' dens4 'DFIN' dens1;
  148. LDL = DROI PD PL 'DINI' dens1 'DFIN' dens1;
  149. LNW = (INVE LC5L) ET LC4C5 ET (INVE LHC4) ET LHD ET LDL;
  150. SNW = SURF LNW;
  151.  
  152. DAMMPART = SS ET SN;
  153. ELASPART = SSW ET SSE ET SNE ET SNW;
  154.  
  155. LAT_E = LDL ET LLA;
  156. LAT_D = LIC ET LBI;
  157.  
  158. ESTRUT = (DAMMPART 'ET' ELASPART);
  159. *TRAC ESTRUT;
  160. *+++++++++++++++++++++++++++++++++++++++++
  161. *
  162. * END OF PARAMETRIC GEOMETRY
  163. *
  164. *+++++++++++++++++++++++++++++++++++++++++
  165.  
  166. *+++++++++++++++++++++++++++++++++++++++++
  167. *
  168. * MODEL AND MATERIAL
  169. *
  170. *+++++++++++++++++++++++++++++++++++++++++
  171. AT=1.0e0;
  172. AC=1.2e0;
  173. BC=1.5e+03;
  174. BT=1.0E4;
  175. epsd0=1.0E-4;
  176. YG=1.0E+9;
  177. BBT=1.06e0;
  178. NNU=0.0;
  179. RO =2.3E+03;
  180.  
  181. MOD1 = 'MODE' ELASPART 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';
  182. MAT1 = 'MATE' MOD1 'YOUN' YG 'NU' NNU 'RHO' RO;
  183.  
  184. MOD2 = 'MODE' DAMMPART 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'ENDOMMAGEMENT' 'MAZARS';
  185. MAT2 = 'MATE' MOD2 'YOUN' YG 'NU' NNU 'RHO' RO 'BTRA' BT 'KTR0' epsd0 'ATRA' AT 'ACOM' AC 'BCOM' BC 'BETA' BBT;
  186.  
  187. MODTOT = MOD1 'ET' MOD2;
  188. MATTOT = MAT1 'ET' MAT2;
  189.  
  190. *+++++++++++++++++++++++++++++++++++++++++
  191. *
  192. * END OF MODEL AND MATERIAL
  193. *
  194. *+++++++++++++++++++++++++++++++++++++++++
  195.  
  196. *+++++++++++++++++++++++++++++++++++++++++
  197. *
  198. * BOUNDARY CONDITIONS
  199. *
  200. *+++++++++++++++++++++++++++++++++++++++++
  201. CL1 = 'BLOQ' 'UX' LAT_E;
  202. CL2 = 'BLOQ' 'UY' PA;
  203. CLT = CL1 'ET' CL2;
  204.  
  205. LA = ESTRUT 'POIN' 'DROIT' PA PB 0.00001;
  206. CL4 = 'BLOQ' 'UY' LA;
  207. CLT = CLT ET CL4;
  208.  
  209. *+++++++++++++++++++++++++++++++++++++++++
  210. *
  211. * END OF BOUNDARY CONDITIONS
  212. *
  213. *+++++++++++++++++++++++++++++++++++++++++
  214.  
  215. *+++++++++++++++++++++++++++++++++++++++++
  216. *
  217. * SOLUTION
  218. *
  219. *+++++++++++++++++++++++++++++++++++++++++
  220. * -> some initializations
  221. *number of time increments
  222. Nmaxincre=132;
  223. * control parameter
  224. TAU = 7.0E-5;
  225. *u_hat (direction of the applied displacement)
  226. CL3 = 'BLOQ' 'UX' LAT_D;
  227. CLT = CLT 'ET' CL3;
  228. uch = 'DEPI' CL3 1.0 ;
  229.  
  230. lis1 = PROG 0. pas TAU (Nmaxincre * TAU) ;
  231. EV1 = evol manu 't' lis1 'tau' lis1 ;
  232.  
  233. **** preparing PASAPAS input table
  234. tab1=table;
  235. tab1.mova='D ';
  236. tab1.blocages_mecaniques=CLT ;
  237. tab1.caracteristiques=mattot ;
  238. tab1.modele=modtot;
  239. tab1.temps_calcules= lis1 ;
  240.  
  241. **** new lines to be included for ativating path-control *******
  242.  
  243. tab1 . 'PILOTAGE_INDIRECT' = VRAI ;
  244. tab1 . 'DEPLACEMENTS_PILOTES' = uch ;
  245. tab1 . 'PARAMETRE_DE_PILOTAGE' = EV1;
  246.  
  247. pasapas tab1;
  248. *+++++++++++++++++++++++++++++++++++++++++
  249. *
  250. * END OF SOLUTION
  251. *
  252. *+++++++++++++++++++++++++++++++++++++++++
  253.  
  254. *+++++++++++++++++++++++++++++++++++++++++
  255. *
  256. * POST-PROCESSING
  257. *
  258. *+++++++++++++++++++++++++++++++++++++++++
  259. nn = DIME tab1 . temps;
  260.  
  261. PROF = PROG 0.;
  262. PROT = PROG 0.;
  263. PROD = PROG 0.;
  264.  
  265. REPE LAB (nn - 1 ) ;
  266. ii = &LAB ;
  267. TT = tab1 . TEMPS . ii ;
  268. deplac = tab1 . 'DEPLACEMENTS' . ii ;
  269. force_indirect = 'REAC' CL1 deplac ;
  270. FF = RESFORCE force_indirect LAT_E 'FX' ;
  271. DD = 'EXTR' deplac 'UX' PB ;
  272.  
  273. PROT = PROT et ( PROG TT) ;
  274. PROF = PROF et ( PROG (-1.0*FF)) ;
  275. PROD = PROD et ( PROG DD) ;
  276.  
  277. Fin LAB ;
  278.  
  279. PROETA = TAB1 . 'COEFFICIENT_DE_PILOTAGE' ;
  280. evo11 = evol turq manu 'd' PROD 'f' PROF ; dess evo11 ;
  281. evo12 = evol jaun manu 't' PROT 'eta' PROETA ; dess evo12 ;
  282.  
  283. *+++++++++++++++++++++++++++++++++++++++++
  284. *
  285. * END OF POST-PROCESSING
  286. *
  287. *+++++++++++++++++++++++++++++++++++++++++
  288. fin;
  289.  
  290.  

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