PARADISCYL:Cylinder-Benchmark2

Benchmarks to check traction boundary condition on the surface

ILL RYU

May 1, 2010

From these benchmarks, we like to confirm if the traction boundary condition holds on the surface, which means that the traction force should be zero on the cylinder boundary. To do that, two models are prepared , varying the geometry of dislocation loop, as following figures.

Model geometries

• 

  [flat loop]

• 

  [Inclined loop ]

To satisfy traction zero boundary, we need to check if $\sigma _{rr}\,\!$ ,$\sigma _{r\theta }\,\!$, and $\sigma _{rz}\,\!$ can be zero on the surface.

---

Total stress on the cylinder surface

• 

  [flat loop]

• 

  [Inclined loop ]

Conclusion

Based on the plot, we can see that $\sigma _{rr}\,\!$ ,$\sigma _{r\theta }\,\!$, and $\sigma _{rz}\,\!$ can be zero on the surface.Therefore, the traction boundary condition is satisfied well.

Input & plot files

ParaDiS input
Matlab file for plot