Issue 37

I.Llavori et alii, Frattura ed Integrità Strutturale, 37 (2016) 87-93; DOI: 10.3221/IGF-ESIS.37.12

MPC

P

δ app

U x

= 0

F

= 0

U y

Figure 2 : Boundary condition of the numerical model.

Start analysis: Read data and generate mesh

FEM simulation

1.- Compute SWT and Miner damage 2.- Compute wear and set new geometry 3.- Interp. of Miner’s damage to new node location

No

ω =1?

Yes

Generate crack segment and LSM analysis

X-FEM simulation

1.- Compute Stress Intensity Factor 2.- Compute Crack propagation rate 3.- Compute wear and set new geometry

No

Failure?

Yes

End Analysis

Figure 3: Flow chart showing computational sequence for numerical analysis.

R ESULTS

ig. 7 illustrates the evolution of the Von Mises (VM) stress field at different stages of the δ = 36 µm amplitude displacement simulation (Test 9A), when the maximum axial load is applied. The top left image corresponds to the unworn specimen at the start of the simulation. At this stage, the maximum VM peak stress (1348 MPa) lies at the right end of the contact, that is, on the side where the axial load is applied. The top right figure corresponds to the stress state prior to crack initiation, at 10,000 cycles. It is observed that the contact area is slightly bigger than the initial one, while the reduction of the maximum VM value is around 25% due to the contact stress redistribution originated by wear. F

90