Difference between revisions of "Defining Cracks"

Revision as of 17:34, 28 September 2013

This section explains how to add explicit cracks to an MPM simulation.

Introduction

The extension of MPM to model explicit cracks is called CRAMP for CRAcks in the Material Point Method. The commands is this section are used to define the explicit cracks (and there can be any number of them). In addition, whenever a simulation has cracks you can use various crack settings commands to customize the features of the crack physics and set various material fracture properties.

The development of CRAMP is described first in a paper by Nairn (2003)^[1] (and its development was done using NairnMPM). Some other papers discuss calculation of J integral and stress intensity factor,^[2], propose energy balance propagation,^[3] handle cracks in 3D^[4] (although not currently available in NairnMPM), use cracks to model imperfect interfaces,^[5], and use traction laws with cracks.^[6] Some applications of using cracks in MPM include wood fracture,^[7] use of cracks to model glue bods in oriented strand board,^[8] simulation of dynamic fracture,^[9] and simulation of crack growth with fiber bridging.^[10]

Crack Definition Commmands

In CRAMP,^[1] a crack is defined as series of massless particles connected by crack segments. The first particle is the "start" tip and the last particle is the "end" tip. When CRAMP is active, nodes around the crack will be divided into two velocity fields to separately model motion of particle "above" the crack and "below" the crack. You define a crack by providing the points along the crack - (x₁,y₁) to (x_n,y_n). In addition, CRAMP fully accounts for crack contact, can model cracks with frictional contact, can use cracks to model imperfect interfaces, and can insert traction laws to model cohesive zones. When defining cracks, some additional are used to set crack contact mechanics, set crack tip materials, and insert traction laws.

In scripted files, a single crack is started with a NewCrack command, which is followed by one or more of the other commands to complete the crack definition:

NewCrack
GrowCrack
GrowCrackLine
GrowCrackArc
Friction
CrackInterface
CrackThickness

In XML files, each crack is defined in a CrackList element:

<CrackList friction='0.1' Dn='200' Dnc='-1' Dnt='200' Dt='5'>
  <pt units='mm' x='50.5' y='0' tip='1'/>
  <pt x='53' y='0'/>
     .
     .
     .
  <pt x='100.5' y='0'/>
  <Line xmin="102" ymin="0" xmax="130" ymax="10" resolution="5"
           start_tip="1" end_tip="1"/>
  <Circle xmin="102" ymin="0" xmax="130" ymax="10" resolution="10"
           start_tip="1" end_tip="1" start_angle="0" end_angle="90"/>
  <Thickness>1.0</Thickness>
</CrackList>

The commands and their functions are:

NewCrack command - this command starts a new crack (in XML files, the <CrackList> element starts the crack.
GrowCrack, GrowCrackLine, and GrowCrackArc commands - this commands extend the crack (the corresponding XML commands are <pt>, <Line>, and <Circle>).
Friction and CrackInterface commands - these commands set the crack contact properties.
CrackThickness command - this commands sets the crack thickness (the corresponding XML command is <Thickness>).

Starting the Crack

In scripted files, crack is started with the NewCrack command:

NewCrack (x1),(y1),<(mat)>,<(frict)>,<(tract)>

where

(x1),(y1) - defines the start of the crack with crack tip (the coordinates are (R,Z) if axisymmetric).
(mat) - this optional parameter can be one of three things:
1. A Material ID - It can define the material at the start of the crack. Crack propagation is only possible when one or both of its crack tips have a material ID. If the material ID is omitted, the MPM calculations will correctly model that crack, but no propagation or fracture parameter calculations will occur.
2. "exterior" - this text indicates a crack tip at the edge of the body.
3. "fixed" - this text indicates a fixed crack. A fixed crack will not translate with the object and will not track surfaces. It can calculate J Integral but not Stress Intensity factor. It is useful when using cracks to create surfaces in an object. The surface can change by propagating the crack at constant crack speed initiated by critical J or by initiation time (see steady state propagation).
(frict) - this optional parameter can specify a custom contact law for this crack. This option lets each crack have a different contact law such as different coefficients of friction. See the Friction Command for options in setting the contact law. If you want to set a contact law on a non-fixed crack, but do not want to define a tip material in parameter #3, you can enter parameter #3 as free, and then enter optional #4.
(frict) and (traction) - alternatively, the optional #4 and #5 parameters can assign a traction law to the initial crack tip. To use this option, set #4 to traction and #5 to a material ID for a traction law material. The traction law will applied only to the first crack particle. To continue the traction law along the crack, use the traction law options in subsequent GrowCrack, GrowCrackLine and GrowCrackArc commands.

Growing the Crack

Crack Contact Mechanics

Crack Thickness

Interacting Cracks

References

↑ ^1.0 ^1.1 J. A. Nairn, "Material Point Method Calculations with Explicit Cracks," Computer Modeling in Engineering & Sciences, 4, 649-664 (2003). (See PDF)
↑ Y. Guo and J. A. Nairn, "Calculation of J-Integral and Stress Intensity Factors using the Material Point Method," Computer Modeling in Engineering & Sciences, 6, 295-308 (2004). (See PDF)
↑ J. A. Nairn, "Simulation of Crack Growth in Ductile Materials,"; Engr. Fract. Mech., 72, 961-979 (2005). (See PDF)
↑ Y. Guo and J. A. Nairn, "Three-Dimensional Dynamic Fracture Analysis Using the Material Point Method," Computer Modeling in Eng. & Sci., 16, 141-156 (2006). (See PDF)
↑ J. A. Nairn, "Numerical Implementation of Imperfect Interfaces, Computational Materials Science, 40, 525-536 (2007). (See PDF)
↑ J. A. Nairn, "Analytical and Numerical Modeling of R Curves for Cracks with Bridging Zones," Int. J. Fracture, 155, 167-181 (2009). (See PDF)
↑ J. A. Nairn, "Material Point Method Simulations of Transverse Fracture in Wood with Realistic Morphologies," Holzforschung, 61, 375-381 (2007). (See PDF)
↑ J. A. Nairn and E. Le, "Numerical Modeling and Experiments on the Role of Strand-to-Strand Interface Quality on the Properties of Oriented Strand Board," Proc of 9th Int. Conf. on Wood Adhesives, Lake Tahoe, Neveda, USA, Sept. 28-30, 2009. (See PDF)
↑ S. G. Bardenhagen, J.A. Nairn, and H. Lu, "Simulation of dynamic fracture with the Material Point Method using a mixed J-integral and cohesive law approach," Int. J. Fracture, 170, 49-66 (2011).
↑ N. Matsumoto and J.A. Nairn, "Fracture Toughness of Wood and Wood Composites During Crack Propagation," Wood and Fiber Science, 44, 121-133 (2012). (See PDF)

[CRAMP-1] 1.0 ^1.1 J. A. Nairn, "Material Point Method Calculations with Explicit Cracks," Computer Modeling in Engineering & Sciences, 4, 649-664 (2003). (See PDF)

[GuoJ-2] Y. Guo and J. A. Nairn, "Calculation of J-Integral and Stress Intensity Factors using the Material Point Method," Computer Modeling in Engineering & Sciences, 6, 295-308 (2004). (See PDF)

[EB-3] J. A. Nairn, "Simulation of Crack Growth in Ductile Materials,"; Engr. Fract. Mech., 72, 961-979 (2005). (See PDF)

[Guo3D-4] Y. Guo and J. A. Nairn, "Three-Dimensional Dynamic Fracture Analysis Using the Material Point Method," Computer Modeling in Eng. & Sci., 16, 141-156 (2006). (See PDF)

[IIC-5] J. A. Nairn, "Numerical Implementation of Imperfect Interfaces, Computational Materials Science, 40, 525-536 (2007). (See PDF)

[RCurve-6] J. A. Nairn, "Analytical and Numerical Modeling of R Curves for Cracks with Bridging Zones," Int. J. Fracture, 155, 167-181 (2009). (See PDF)

[wood-7] J. A. Nairn, "Material Point Method Simulations of Transverse Fracture in Wood with Realistic Morphologies," Holzforschung, 61, 375-381 (2007). (See PDF)

[OSB-8] J. A. Nairn and E. Le, "Numerical Modeling and Experiments on the Role of Strand-to-Strand Interface Quality on the Properties of Oriented Strand Board," Proc of 9th Int. Conf. on Wood Adhesives, Lake Tahoe, Neveda, USA, Sept. 28-30, 2009. (See PDF)

[BardF-9] S. G. Bardenhagen, J.A. Nairn, and H. Lu, "Simulation of dynamic fracture with the Material Point Method using a mixed J-integral and cohesive law approach," Int. J. Fracture, 170, 49-66 (2011).

[MDF-10] N. Matsumoto and J.A. Nairn, "Fracture Toughness of Wood and Wood Composites During Crack Propagation," Wood and Fiber Science, 44, 121-133 (2012). (See PDF)

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

@@ Line 59: / Line 59: @@
 <li><tt>(mat)</tt> - this optional parameter can be one of three things:
 <ol>
-<li>It can be used to define a material ID for the material at the start of the crack. Crack propagation is only possible when one or both of its crack tips have a material ID. If the material ID is omitted, the MPM calculations will correctly model that crack, but no propagation or fracture parameter calculations will occur.
+<li>A Material ID - It can define the material at the start of the crack. Crack propagation is only possible when one or both of its crack tips have a material ID. If the material ID is omitted, the MPM calculations will correctly model that crack, but no propagation or fracture parameter calculations will occur.
-<li>It can be "exterior" to indicate a [[#Exterior Crack Tip|crack to the edge of the body]].
+<li><tt>"exterior"</tt> - this text indicates a [[#Exterior Crack Tip|crack tip at the edge of the body]].
-<li>It can be "fixed" to indicate a fixed crack. A fixed crack will not translate with the object and will not track surfaces. It can calculate J Integral but not Stress Intensity factor. It is useful when using cracks to create surfaces in an object. The surface can change by propagating the crack at constant crack speed initiated by critical J or by initiation time (see steady state propagation).
+<li><tt>"fixed"</tt> - this text indicates a fixed crack. A fixed crack will not translate with the object and will not track surfaces. It can calculate J Integral but not Stress Intensity factor. It is useful when using cracks to create surfaces in an object. The surface can change by propagating the crack at constant crack speed initiated by critical J or by initiation time (see steady state propagation).
 </ol>
 <li><tt>(frict)</tt> - this optional parameter can specify a custom contact law for this crack. This option lets each crack have a different contact law such as different coefficients of friction. See the Friction Command for options in setting the contact law. If you want to set a contact law on a non-fixed crack, but do not want to define a tip material in parameter #3, you can enter parameter #3 as free, and then enter optional #4.