User:Cbonneville/sandbox
Computational Methods In Fracture Mechanics
[edit]Many methods exist for calculating G with finite elements. Although a direct calculation of the J-integral is possible (using the strains and stresses outputted by finite element analysis), approximate approaches for some types of crack growth exist and provide reasonable accuracy with straightforward calculations.
Nodal Release Method
[edit]![](http://upload.wikimedia.org/wikipedia/commons/thumb/e/ed/Nodal_Release_Debounded.png/540px-Nodal_Release_Debounded.png)
If the crack is growing straight, the energy release rate can be decomposed as a sum of 3 terms associated with the energy in each 3 modes. As a result, the Nodal Release method (NR) can be used to determine from FEA results. The energy release rate is calculated at the nodes of the finite element mesh for the crack at an initial length and extended by a small distance . First, we calculate the displacement variation at the node of interest (before and after the crack tip node is released). Secondly, we keep track of the nodal force outputted by FEA. Finally, we can find each components of using the following formulas:
Where is the width of the element bounding the crack tip. The accuracy of the method highly depends on the mesh refinement, both because the displacement and forces depend on it, and because . Note that the equations above are derived using the crack closure integral.
If the energy release rate exceed a critical value, the crack will grow. In this case, a new FEA simulation is perfomed (for the next time step) where the node at the crack tip is released. For a bounded substrate, we may simply stop enforcing fixed Dirichlet boundary conditions at the crack tip node of the previous time step (i.e. displacements are no longer restrained) . For a symmetric crack, we would need to update the geometry of the domain with a longer crack opening (and therefore generate a new mesh[1]).
MCCI Method
[edit]Terrence is doing this part
Domain Integral Approach for J
[edit]![](http://upload.wikimedia.org/wikipedia/commons/thumb/d/d2/Contour_Domain_Integral_J-Integral.png/309px-Contour_Domain_Integral_J-Integral.png)
The J-integral may be calculated directly using the finite element mesh and shape functions[2]. We consider a domain countour as shown in figure 2 and choose an arbitrary smooth function such that on and on .
For linear elastic cracks growing straight ahead, . The energy release rate can then be calculated over the area bounded by the contour using an updated formulation:
The formula above may be applied to any annular area surrounding the crack tip (in particular, a set of neighboring elements can be used). This method is very accurate, even with a coarse mesh arround the crack tip (one may choose an integration domain located far away, with stresses and displacement less sensitive to mesh refinment)
Derivation of the J-integral for domain integral approach The J-intregral may be expressed over the full contour as follow:
With . , on and the work and stresses cancel out on and , hence by application of the divergence theorem this leads to:
Finally, by noting that and using the equilibrium equation:
Crack Tip Singular Elements
[edit]Terrence is doing this part
![]() | This is a user sandbox of Cbonneville. You can use it for testing or practicing edits. This is not the sandbox where you should draft your assigned article for a dashboard.wikiedu.org course. To find the right sandbox for your assignment, visit your Dashboard course page and follow the Sandbox Draft link for your assigned article in the My Articles section. |
- ^ "ScienceDirect". www.sciencedirect.com. doi:10.1016/s0045-7825(97)00287-9. Retrieved 2019-05-10.
- ^ Zehnder, Alan T. (2012). Fracture Mechanics. Lecture Notes in Applied and Computational Mechanics. Vol. 62. Dordrecht: Springer Netherlands. doi:10.1007/978-94-007-2595-9. ISBN 9789400725942.