Numérotation des lignes :

1. * fichier : fefp_powcap_bcn.dgibi
2. ************************************************************************
3. ************************************************************************
4. * COMPACTION OF A FLANGED COMPONENT
5. *
6. *------------- OPCIONES GENERALES --------------------------------
7. *
8. OPTION DIME 2 TRAC x ELEM qua8 MODE axis;
9. *
10. *------------- CREACION DE LA GEOMETRIA --------------------------
11. *
12. r1 = 6.3D0; r2 = r1 + 4.6D0; r3 = r2 + 15.6D0;
13. h1 = 0.D0; h2 = h1 + 13.7D0; h3 = h2 + 11.7D0;
14. p1 = r1 h1;p2 = r2 h1;p3 = r2 h2;
15. p4 = r3 h2;p5 = r3 h3;p6 = r2 h3;
16. p7 = r1 h3;p8 = r1 h2;paux = r3 h3;
17. den1 = 5; den2 = 10; den3 = 10; den4 = 8;
18. l12 = d den1 p1 p2;l23 = d den2 p2 p3;
19. l38 = d den1 p3 p8;l81 = d den2 p8 p1;
20. l34 = d den3 p3 p4;l45 = d den4 p4 p5;
21. l56 = d den3 p5 p6;l63 = d den4 p6 p3;
22. l83 = d den1 p8 p3;l36 = d den4 p3 p6;
23. l67 = d den1 p6 p7;l78 = d den4 p7 p8;
24. malla1 = daller l12 l23 l38 l81;
25. malla2 = daller l34 l45 l56 l63;
26. malla3 = daller l83 l36 l67 l78;
27. malla = malla1 et malla2 et malla3 ;
28. elim malla 1.D-4;
29. llbase = l12 et l34; llsupe = l56 et l67;
30. lllate = l23 et l45 et l78 et l81;
31. *
32. *------------- MODELO -------------------------------------------
33. *
34. E_ELAS = 2000.D0;
35. NUHOR = 0.37D0;
36. ROHOR = 2.5D3;
37. SIGINI = 90.D0;
38. EXPN1 = 1.D0;
39. EXPN2 = 2.7D0;
40. ETAINI = 0.489D0;
41. COHES = 15.D0;
42. ANGFR = 0.52359878D0 * 2.;
43. ANGFR0 = 0.52359878D0 * 1.5;
44. NNC = 2.D0;
45. mod1 = 'MODELI' malla mecanique elastique plastique
46. powdercap_fefp Q8RI;
47. mat1 = mater mod1 youn E_ELAS nu NUHOR rho ROHOR
48. sigy SIGINI nnn1 EXPN1 nnn2 EXPN2
49. eta0 ETAINI cohe COHES
50. phi ANGFR phi0 ANGFR0 nnnc NNC ;
51. *
52. *-------------- CONDICIONES DE CONTORNO -------------------------
53. *
54. rigcont = (bloq uz llbase) et (bloq ur lllate);
55. *
56. *-------------- MOVIMIENTO IMPUESTO -----------------------------
57. *
58. rigsupe = bloq uz llsupe;
59. valor = -6.06D0;
60. movyy = depi rigsupe valor;
61. evol1 = evol manu t (prog 0. 1.) level (prog 0. 1.);
62. evfut = char dimp evol1 movyy;
63. *
64. *--------------- RESOLUCION -------------------------------------
65. *
66. maxiter = 50;
67. t0 = table;
68. t0.MODELE = mod1;
69. t0.CARACTERISTIQUES = mat1;
70. t0.BLOCAGES_MECANIQUES = rigcont et rigsupe;
71. t0.CHARGEMENT = evfut;
72. t0.ACCELERATION = maxiter;
73. t0.MAXITERATION = maxiter;
74. t0.PRECISION = 1.E-10;
75. t0.CONVERGENCE_FORCEE = faux;
76. t0.K_TANGENT = vrai;
77. t0.NITER_KTANGENT = 1;
78. t0.PRECISINTER = 1.E-9;
79. t0.NITERINTER_MAX = 100;
80. t0.FEFP_FORMULATION = vrai;
81. t0.UPDATE_LAGRANGIAN = vrai;
82. t0.TEMPS_CALCULES = prog 0. pas 0.0001 0.001
83. pas 0.0005 0.01
84. pas 0.005 0.1
85. pas 0.010 0.2;
86. t0.TEMPS_SAUVES = t0.TEMPS_CALCULES;
87. t0.'PROCESSEURS'='MONO_PROCESSEUR';
88. pasapas t0;
89. *
90. *--------------- POSTPROCESO ------------------------------------
91. *
92. imax = (dime t0.TEMPS_CALCULES) - 1;
93. fuer = prog 0.;
94. i = 1;
95. geo0 = form;
96. repeter blocdefi imax;
97. aux = reac rigsupe t0.deplacements.i;
98. aux = resu aux;
99. pbas = extr aux MAIL ;
100. pbas = POIN 1 pbas ;
101. aux = extr aux fz pbas;
102. fuer = fuer et (prog aux);
103. i=i+1;
104. fin blocdefi;
105. des = valor * t0.TEMPS_CALCULES;
106. law = evol manu t des f(t) fuer;
107. err = (aux + 2.1063258E+04)/2.1063258E+04 ;
108. list err;
109. SI ((ABS err) < 0.5e-5) ;
110. ERRE 0 ;
111. SINON ;
112. ERRE 5 ;
113. FINSI ;
114. FIN;
115.
116.
117.
118.
119.
120.
121.
122.
123.
124.
125.
126.
127.

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