ABAQUS user subroutines for the simulation of viscoplastic - loicz
ABAQUS user subroutines for the simulation of viscoplastic - loicz
ABAQUS user subroutines for the simulation of viscoplastic - loicz
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
5.2 Plate containing a hole — BODNER-PARTOM model coupled with damage<br />
The material parameters <strong>of</strong> <strong>the</strong> BODNER-PARTOM model <strong>of</strong> IN 738 LC at 850 °C have also been<br />
determined by OLSCHEWSKI et al. [1990], and <strong>the</strong>ir values at 850 °C are shown in Table 4. The<br />
material parameters <strong>of</strong> <strong>the</strong> damage model are <strong>the</strong> same as listed in Table 3.<br />
E 149650 MPa ν 0.33 K0 4.18 10 5 MPa<br />
D0<br />
8.82 10 9 h -1<br />
n 0.289 K1 3.76 10 5 MPa<br />
A1=A2 1.65 10 -7 MPa/h m1 0.581 K2 3.07 10 5 MPa<br />
r1=r2 5.4 m2 0.344 K3 1.54 10 5 MPa<br />
Table 4: Material parameters <strong>of</strong> <strong>the</strong> BODNER-PARTOM model <strong>for</strong> IN 738 LC at 850 °C.<br />
In gas turbine blades with cooling channels, stress concentration occurs due at <strong>the</strong>se channels. A<br />
square plate with a central circular hole is <strong>the</strong>re<strong>for</strong>e chosen as a model representation <strong>of</strong> <strong>the</strong> area <strong>of</strong><br />
blades where <strong>the</strong> air cooling channels are located. The FE model used <strong>for</strong> <strong>the</strong> calculation is shown in<br />
Fig. 6. First, <strong>the</strong> plate is subjected to a creep load <strong>of</strong> σ3 = 180 MPa. After 40000 hours <strong>the</strong> maximum<br />
damage reaches a value <strong>of</strong> about 0.1. A second load <strong>of</strong> σ2 = 180 MPa is <strong>the</strong>n applied. There is only<br />
one element in <strong>the</strong> thickness direction so that any gradient over <strong>the</strong> thickness can not be captured. The<br />
three-dimensional 8-node linear brick continuum element with reduced integration, C3D8R, from <strong>the</strong><br />
element library <strong>of</strong> <strong>ABAQUS</strong> is used, and geometric non-linearity has been considered. Figs. 7 and 8<br />
show <strong>the</strong> contour plot <strong>of</strong> <strong>the</strong> maximum principal damage after 40000 h and 98000 h, respectively. β is<br />
assumed to be 0.5. Distribution <strong>of</strong> <strong>the</strong> maximum principal value <strong>of</strong> <strong>the</strong> strain and stress, after 40000 h<br />
and 98000 h, are shown in <strong>the</strong> Figs. 9-12, respectively. The <strong>ABAQUS</strong> input-file used <strong>for</strong> <strong>the</strong><br />
computation is given in <strong>the</strong> Appendix 2.<br />
25