07.12.2012 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!