Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 725795, 11 pages

http://dx.doi.org/10.1155/2014/725795

Research Article

Dynamic Crack Analysis in Isotropic/Orthotropic Media via Extended Isogeometric Analysis

Tiantang Yu, Yongling Lai, and Shuohui Yin

Department of Engineering Mechanics, Hohai University, 1 Xikang Road, Nanjing 210098, China

Received 8 February 2014; Revised 13 March 2014; Accepted 13 March 2014; Published 8 April 2014

Academic Editor: Hung Nguyen-Xuan

Copyright © 2014 Tiantang Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The extended isogeometric analysis (X-IGA) is the combination of the extended finite element method (X-FEM) and the isogeometric analysis (IGA), so the X-IGA possesses the advantages of both methods. In this paper, the X-IGA is extended to investigate the dynamic stress intensity factors of cracked isotropic/orthotropic media under impact loading. For this purpose, a corresponding dynamic X-IGA model is developed, the Newmark time integration scheme is used to achieve a dynamic response, and the dynamic stress intensity factors are evaluated through the contour interaction integral technique. Numerical simulations show that the X-IGA results agree with other available reference solutions, and accurate results can be obtained by using the X-IGA with a relatively coarse mesh.

1. Introduction

Among existing numerical methods, the extended finite element method (X-FEM), which was developed in 1999 by Belytschko and his coworkers [1, 2], is the most effective method for solving discontinuity problems. In the X-FEM, the geometry of the discontinuity is independent of the computational mesh; therefore, the computational mesh does not need to be updated in the simulation of the propagation of discontinuities. In the past decades, the extensive literature on improving and applying the original X-FEM in modeling discontinuities has been published [3–15]. The X-FEM aims to enrich standard finite element approximation by using discontinuous basis functions in the framework of partition of unity.

The standard finite element approximation is element-based polynomial approximation, and it yields discretization errors in complex geometry. In order to overcome the drawback, the isogeometric analysis (IGA) was proposed by Hughes et al. [16] in 2005. The principle of the IGA is that nonuniform rational B-splines (NURBS) basis functions are employed as shape functions for geometric description and field approximation. The IGA has some unique advantages [16], so it has been developed and applied in many fields including structural mechanics [17–24], solid mechanics [25–27], fluid mechanics [28], and contact mechanics [29].

Recently, the IGA has been improved with enrichment functions to solve linear elastic fracture mechanics problems [30–32] and curved interface problems [33]. The extended isogeometric analysis (X-IGA) was proposed in [31] and contains the inherent advantages of both IGA and X-FEM. Currently, the X-IGA is further developed to analyze cracked orthotropic media [34]. In this study, the X-IGA is extended to investigate the dynamic fracture behavior of stationary cracks in isotropic/orthotropic media. A corresponding dynamic X-IGA model is developed, and the Newmark time integration scheme is used to achieve a dynamic response. The dynamic stress intensity factors (DSIFs) are evaluated by using the contour interaction integral technique.

The paper is organized as follows. Section 2 briefly reviews the fundamentals of NURBS-based IGA. The X-IGA for dynamic cracked isotropic/orthotropic media is described in Section 3. The DSIFs are derived by using the contour interaction integral technique in Section 4. Section 5 presents the numerical results obtained via X-IGA and compares the results with other solutions. Finally, conclusions and prospects are drawn in Section 6.

2. NURBS-Based Isogeometric Analysis

In the CAD, a two-dimensional NURBS surface can be constructed as [35]

where is the NURBS basis function; and are, respectively, the orders of basis functions in the direction and the direction; represents the coordinate of control point; and and are the numbers of basis functions in the direction and the direction, respectively.

In NURBS-based isogeometric analysis [16], the displacement field is approximated similarly to the isoparametric finite element method as follows:

where is the number of the control points; is the parametric coordinate; denotes the NURBS basis function at control point ; and is the displacements of control point .

3. X-IGA for Dynamic Cracked Media

3.1. Displacement Approximation

Similar to the X-FEM, the X-IGA aims to enrich the standard IGA approximation by using additional functions on the basis of the partition of unity to model discontinuities. For crack problems, the enriched displacement approximation can be expressed as

where , , and are the shape functions of standard IGA (NURBS basis functions); , , and are, respectively, the displacement and enrichment variable vectors at control point; , , and are, respectively, the set of all control points in the computational domain, the set of control points enriched with a modified Heaviside step function , and the set of control points enriched with the crack-tip branch enrichment functions ; ; and the basis function support of the control point in is completely split by the crack, whereas that of the control point in is partly split by the crack.

The modified Heaviside step function is given by

The crack-tip branch enrichment functions for orthotropic materials are defined as [36]

where and are the crack-tip local polar coordinates, and

where and are crack-tip material parameters and    are the roots of

where is the compliance coefficients.

Equation (5) cannot be directly used for isotropic materials because of the presence of 0/0 in the equation. The crack-tip branch enrichment functions for isotropic materials are defined as [1]

3.2. Discrete Equilibrium Equations

Considering a body with an initial traction-free crack in the state of dynamic equilibrium, the weak form of the momentum equation is

where is the displacement vector; and are the stress and strain tensors, respectively; and are the body force and external traction vectors, respectively; and is the mass density.

The discretized form of (9) using the X-IGA approximation (3) can be written as

where is the vector of the unknown variable at control points and , , and are, respectively, the global mass matrix, stiffness matrix, and external force vector at control points.

The element contribution to is

where

and the element contribution to is

where

with

Additionally, the element contribution to is

with

where and are the external traction and body force vectors, respectively.

3.3. Numerical Integration Scheme

The Gauss quadrature scheme is employed in the X-IGA. To obtain an accurate integration for crack tip elements and elements cut by crack, the triangular subdomain technique is used in the same way as that of the X-FEM [1]. For additional details, refer to [31].

3.4. Time Integration Scheme

The Newmark method is adopted for the time integration of (10). At time step , the discrete simultaneous equations are described as

where is the time step and the unconditionally stable conditions are and . In this study, and are chosen.

The implementation procedure of the Newmark time integration scheme applying into the X-IGA is outlined as follows.

(1)The initial velocities and initial displacements are set to be zero; then the initial accelerations are obtained with (10).(2)Compute the constants , , , , , , , and .(3)Compute .(4)Loop over the time steps as follows:

(a)compute ;(b)compute the displacements by solving ;(c)solve for the value , and the velocities are calculated by using .

(5)Repeat the loop for the next time step until the maximum time step is reached.

4. DSIF Calculations

The DSIFs are evaluated by using the domain form of the contour interaction integral technique. The following states are considered: state 1 (, , ), which corresponds to the actual state; state 2 (, , ), which is an auxiliary state that will be selected as the asymptotic field for model I or II. The interaction integral may be written as [37]

where is the interaction strain energy and is a weighing function. The second term of the right-hand side of the equation is the contribution of the inertia forces.

For isotropic materials, the relation between the interaction integral and the stress intensity factors is expressed as [37]

where is for plane stress and is for plane strain.

The relation between the interaction integral and the stress intensity factors for orthotropic materials is expressed as [36]

with

Choosing state 2 as mode I or mode II leads to mode-I or mode-II DSIF in terms of the interaction integral.

For isotropic materials,

For orthotropic materials,

5. Numerical Simulations

In this section, several examples of stationary cracks in isotropic/orthotropic media with available reference solutions are investigated to assess the accuracy of the proposed approach. The first two examples demonstrate the efficiency of the X-IGA for isotropic material problems, and the next two illustrate the application of the X-IGA to orthotropic material problems, whereas the last example demonstrates the capability and versatility of the X-IGA in modeling the complicated geometries. In all examples, the Heaviside step loading is considered and degree 3 X-IGA is adopted; that is, . Meshes are generated with linear parameterization method. In the case of investigating the effects of different meshes on the results, different uniform meshes are adopted; in other cases, the fine meshes are used around the crack, while the coarse meshes are used in other domains; thus excellent accuracy can be achieved at a low cost.

5.1. Semi-Infinite Crack in an Infinite Isotropic Plate

For the first example, we consider an infinite isotropic plate with a semi-infinite crack subjected to a tensile stress wave perpendicular to the crack face, as shown in Figure 1. The material properties of the media are as follows: Young’s modulus 210 GPa, Poisson’s ratio 0.3, and mass density 8000 kg/m³. The tensile stress 500 MPa. The theoretical solution of mode-I DSIF for the stationary crack is expressed as [38]

where   (: the dilatational wave speed).

Figure 1: Geometry and loading of a semi-infinite crack.

5.1.1. Convergence Study of the DSIF versus Meshes

In this study, five meshes with 27 × 9, 47 × 17, 67 × 25, 87 × 33, and 107 × 41 elements are considered. Figure 2 illustrates the typical regular mesh of 27 × 9 elements.

Figure 2: Control points and elements for degree 3 X-IGA with 27 × 9 elements (the blue cross represents the control points enriched by Heaviside function; the blue square represents the control points enriched by crack tip enrichment functions; the red line represents the crack).

Figure 3 presents the normalized mode-I DSIFs as a function of time for the five considered meshes. In the mesh refinement process, the mode-I DSIFs obtained by the X-IGA converge well to the theoretical solution. The maximum error increases beside for the coarse mesh and decreases significantly for the fine mesh. Similar observation was noted in the simulation of the singular edge-based smoothed finite element method (sES-FEM) [37]. The same simulations were conducted by other researches by using different methods. However, the same conclusion is recorded; that is, the error is lower in later stages than in earlier stages. The relative errors during are presented in Figure 4, in which the error is unchanged as the mesh is refined to a certain extent. In this study, the errors are the same for meshes with 67 × 25, 87 × 33, and 107 × 41 elements.

Figure 3: Convergence of the normalized mode-I DSIFs versus mesh densities.

Figure 4: Convergence in percentage errors for the normalized mode-I DSIFs versus mesh densities during .

5.1.2. Comparison with X-FEM Results

The X-IGA results using 67 × 25 mesh elements are compared with three available X-FEM solutions, including the X-FEM using mass lumping (120 × 60 quadrilateral elements) [39], the mesh-free enriched X-FEM with crack-tip enrichments (78 × 39 quadrilateral elements) [38], and X-FEM with a new enrichment function (120 × 59 quadrilateral elements) [40]. Figure 5 shows the comparative study of the normalized mode-I DSIFs. The results obtained with X-FEM using mass lumping [39] are unstable, and the accuracy of this method is the worst recorded among the four approaches. The results obtained with the X-IGA, the X-FEM using an enrichment function, and the mesh-free enriched X-FEM with crack-tip enrichments are very accurate during ; however, the results are relatively poor during . The X-IGA results are more accurate than other numerical results during .

Figure 5: Comparison of the normalized mode-I DSIFs among different methods.

5.2. Arbitrarily Oriented Central Crack in a Rectangular Isotropic Plate

We consider a rectangular isotropic plate with an arbitrarily oriented central crack subjected on the top and bottom of a uniform impact loading, as shown in Figure 6(a). The length of the crack is 4.8 mm. The material properties of the plate are as follows: Young’s modulus 200 GPa, Poisson’s ratio 0.3, and mass density 5000 kg/m³. Three crack inclination angles are investigated, and the total time of the simulation is 20 μs. In the X-IGA, the computational mesh is independent of the crack; therefore, the computation meshes for the three crack inclination angles are the same, as shown in Figure 6(b).

(a)

(b)

Figure 6: A rectangular plate with an arbitrarily oriented crack. (a) Geometry and loading, (b) computational meshes.

Figures 7 and 8, respectively, present the normalized mode-I and mode-II DSIFs for the three considered crack angles computed by the X-IGA, sES-FEM [37], and X-FEM. Figure 6(b) is also the mesh for the X-FEM, whereas the mesh of 60 × 120 triangular elements is used in the analysis of the sES-FEM.

Figure 7: Comparison of the normalized mode-I DSIFs obtained by the X-FEM, the sES-FEM, and the X-IGA for three crack inclination angles.

Figure 8: Comparison of the normalized mode-II DSIFs obtained by the X-FEM, the sES-FEM, and the X-IGA for three crack inclination angles.

A good agreement can be observed for the three considered crack inclination angles among the three numerical approaches.

5.3. Edge Crack in a Rectangular Orthotropic Plate

As a third example, a rectangular orthotropic plate with a 12 mm edge horizontal crack is considered (Figure 9). The plane stress condition is considered, and the material properties of the plate are as follows: 118.3 GPa, 54.8 GPa, 8.79 GPa, 0.083, and 1900 kg/m³. The time step is selected, where is the wave velocity along the material-axis.

Figure 9: Geometry and loads of an edge crack in a rectangular orthotropic plate.

Figure 10 shows the mesh with 75 × 57 elements used in the simulation. The normalized mode-I DSIFs achieved with the X-IGA (Figure 11) were compared with the reference boundary element method (BEM) (21 elements for the external boundary and eight elements for the crack) [41], the X-FEM (78 × 60 quadrilateral elements) [42], and the conventional finite element method (FEM) results (via ANSYS) [41]. Good agreement is observed among the methods. The X-IGA results more closely match the FEM and the BEM solution compared with the X-FEM results. In [42], the DSIFs are evaluated by using the domain separation integral method, which could be the reason that the X-FEM results are different.

Figure 10: The X-IGA discretization for an edge crack in a rectangular orthotropic plate.

Figure 11: Comparison of the normalized mode-I DSIFs among different methods.

To analyze the effects of different meshes on DSIFs, five different meshes with 27 × 19, 47 × 39, 75 × 57, 97 × 81, and 121 × 101 elements are adopted to simulate the problem. The results are shown in Figure 12. Accurate results can be obtained by using the X-IGA with a relatively coarse mesh.

Figure 12: Comparison of the normalized mode-I DSIFs for different meshes.

5.4. Central Crack with Different Orientations of the Axes of Orthotropy

A rectangular orthotropic plate with a single central horizontal crack was subjected to a Heaviside step tensile distributed load, as shown in Figure 13(a). The length of the crack is 4.8 mm and 20 mm. The material properties of the plate are the same as those in the example in Section 5.3. The plane stress state is assumed, and the time step is selected.

(a)

(b)

Figure 13: A single central crack in a rectangular orthotropic plate. (a) Geometry and loadings, (b) computational meshes.

Two different inclination angles between the material-axis and the crack-face are considered. Figures 14 and 15 show the comparison of the numerical results of BEM (24 elements for the external boundary and 10 for the crack) [41] and X-FEM (50 × 100 quadrilateral elements) [42]. General trends are in good agreement across all methods; the X-IGA results more closely match the BEM solution compared with the X-FEM; the reason for this phenomenon is the same as the above example.

Figure 14: Normalized mode-I DSIFs for a single central crack in a rectangular plate.

Figure 15: Normalized mode-II DSIFs for a single central crack in a rectangular plate.

5.5. Edge Crack in an Annular Isotropic Plate

To illustrate the capability and versatility of the X-IGA in modeling the complicated geometries, the last example deals with an annular isotropic plate under a uniform pressure. The geometry and loads of the plate are depicted in Figure 16. The length of the crack is 2 m; the material properties of the plate are as follows: Young’s modulus 210 GPa, Poisson’s ratio 0.3, and mass density 8000 kg/m³. The total time of the simulation is 0.004s. The plane strain state is assumed, and the time step 8 × 10⁻⁵ s is selected.

Figure 16: Geometry and loads of an edge crack in an annular isotropic plate.

The meshes are depicted in Figure 17 for 285 control points and 221 elements. To perform a comparison, this example is solved by the X-FEM, and the X-FEM meshes are shown in Figure 18 for 1292 nodes and 1221 quadrilateral elements.

Figure 17: Meshes and control points. (The cross represents the control points enriched by Heaviside function; the square represents the control points enriched by crack tip enrichment functions.)

Figure 18: The X-FEM meshes.

The DSIFs are computed with the X-FEM and the X-IGA, and the results are shown in Figure 19. It is found that the results from the X-IGA with a relatively coarse mesh are in good agreement with those from the X-FEM with a very fine mesh. The amplitude of the normalized mode-I DSIF is negative in some time ranges, which may be due to the contact between two crack-faces which is not taken into consideration in the present study.

Figure 19: Comparison of the normalized mode-I DSIFs between the X-FEM and the X-IGA.

6. Conclusions and Prospects

The NURBS-based extended isogeometric analysis (X-IGA) was extended to investigate the dynamic fracture behavior of stationary cracks in isotropic/orthotropic media. A corresponding dynamic X-IGA model was developed, the Newmark time integration scheme was used to achieve the dynamic response, and the DSIFs were evaluated with the contour interaction integral technique. Numerical results indicate that accurate results can be obtained by using the X-IGA with a relatively coarse mesh, and the X-IGA is suitable to model the cracked complicated geometries.

The NURBS-based isogeometric analysis has several drawbacks including the handling trimmed geometries and the local refinement. The isogeometric analysis based on T-splines can effectively overcome these drawbacks [43]. A key advantage of the extended isogeometric analysis over isogeometric analysis is the ability to model crack propagation. Therefore, the application of extended isogeometric analysis based on T-splines to dynamic crack propagation analysis in isotropic/orthotropic media is an area of considerable promise. On the other hand, the smoothed finite element method (SFEM) developed by combining the smoothing technique with the finite element method is an efficient and accurate numerical simulation tool for the dynamic fracture problems [44]; therefore combining the smoothing technique with the X-IGA for fracture analysis is another area of considerable promise.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (51179063). The financial support is gratefully acknowledged.

References

T. Belytschko and T. Black, “Elastic crack growth in finite elements with minimal remeshing,” International Journal for Numerical Methods in Engineering, vol. 45, no. 5, pp. 601–620, 1999.

N. Moës, J. Dolbow, and T. Belytschko, “A finite element method for crack growth without remeshing,” International Journal for Numerical Methods in Engineering, vol. 46, no. 1, pp. 131–150, 1999.

X. Y. Liu, Q. Z. Xiao, and B. L. Karihaloo, “XFEM for direct evaluation of mixed mode SIFs in homogeneous and bi-materials,” International Journal for Numerical Methods in Engineering, vol. 59, no. 8, pp. 1103–1118, 2004.

T. T. Yu and P. Liu, “Improved implementation of the extended finite element method for stress analysis around cracks,” Archives of Civil and Mechanical Engineering, vol. 11, no. 3, pp. 787–805, 2011.

J. V. Cox, “An extended finite element method with analytical enrichment for cohesive crack modeling,” International Journal for Numerical Methods in Engineering, vol. 78, no. 1, pp. 48–83, 2009.

A. R. Khoei and K. Karimi, “An enriched-FEM model for simulation of localization phenomenon in Cosserat continuum theory,” Computational Materials Science, vol. 44, no. 2, pp. 733–749, 2008.

R. Duddu, S. Bordas, D. Chopp, and B. A. Moran, “A combined extended finite element and level set method for biofilm growth,” International Journal for Numerical Methods in Engineering, vol. 74, no. 5, pp. 848–870, 2008.

J. Chessa, P. Smolinski, and T. Belytschko, “The extended finite element method (XFEM) for solidification problems,” International Journal for Numerical Methods in Engineering, vol. 53, no. 8, pp. 1959–1977, 2002.

S. Groß and A. Reusken, “An extended pressure finite element space for two-phase incompressible flows with surface tension,” Journal of Computational Physics, vol. 224, no. 1, pp. 40–58, 2007.

P. Liu, T. T. Yu, T. Q. Bui, and C. Zhang, “Transient dynamic crack analysis in non-homogeneous functionally graded piezoelectric materials by the X-FEM,” Computational Materials Science, vol. 69, pp. 542–558, 2013.

T. T. Yu and L. Y. Shi, “Determination of sharp V-notch stress intensity factors using the extended finite element method,” The Journal of Strain Analysis for Engineering Design, vol. 47, no. 2, pp. 95–103, 2012.

T. T. Yu, “The extended finite element method (XFEM) for discontinuous rock masses,” Engineering Computations, vol. 28, no. 3, pp. 340–369, 2011.

L. Chen, T. Rabczuk, G. R. Liu, S. Bordas, K. Y. Zeng, and P. Kerfriden, “Extended finite element method with edge-based strain smoothing (ESm-XFEM) for linear elastic crack growth,” Computer Methods in Applied Mechanics and Engineering, vol. 209–212, pp. 250–265, 2012.

N. Sukumar, D. L. Chopp, N. Moës, and T. Belytschko, “Modeling holes and inclusions by level sets in the extended finite-element method,” Computer Methods in Applied Mechanics and Engineering, vol. 190, no. 46-47, pp. 6183–6200, 2001.

J. Réthoré, A. Gravouil, and A. Combescure, “An energy-conserving scheme for dynamic crack growth using the extended finite element method,” International Journal for Numerical Methods in Engineering, vol. 63, no. 5, pp. 631–659, 2005.

T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, “Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement,” Computer Methods in Applied Mechanics and Engineering, vol. 194, no. 39–41, pp. 4135–4195, 2005.

R. Bouclier, T. Elguedj, and A. Combescure, “Locking free isogeometric formulations of curved thick beams,” Computer Methods in Applied Mechanics and Engineering, vol. 245-246, pp. 144–162, 2012.

D. J. Benson, Y. Bazilevs, M. C. Hsu, and T. J. R. Hughes, “Isogeometric shell analysis: the Reissner-Mindlin shell,” Computer Methods in Applied Mechanics and Engineering, vol. 199, no. 5–8, pp. 276–289, 2010.

C. H. Thai, H. Nguyen-Xuan, N. Nguyen-Thanh, T.-H. Le, T. Nguyen-Thoi, and T. Rabczuk, “Static, free vibration, and buckling analysis of laminated composite Reissner–Mindlin plates using NURBS-based isogeometric approach,” International Journal for Numerical Methods in Engineering, vol. 91, no. 6, pp. 571–603, 2012.

H. Nguyen-Xuan, C. H. Thai, and T. Nguyen-Thoi, “Isogeometric finite element analysis of composite sandwich plates using a higher order shear deformation theory,” Composites B: Engineering, vol. 55, pp. 558–574, 2013.

C. H. Thai, A. J. Ferreira, E. Carrera, and H. Nguyen-Xuan, “Isogeometric analysis of laminated composite and sandwich plates using a layerwise deformation theory,” Composite Structures, vol. 104, pp. 196–214, 2013.

C. H. Thai, A. J. M. Ferreira, S. P. A. Bordas, T. Rabczuk, and H. Nguyen-Xuan, “Isogeometric analysis of laminated composite and sandwich plates using a new inverse trigonometric shear deformation theory,” European Journal of Mechanics: A/Solids, vol. 43, pp. 89–108, 2014.

L. V. Tran, C. H. Thai, and H. Nguyen-Xuan, “An isogeometric finite element formulation for thermal buckling analysis of functionally graded plates,” Finite Elements in Analysis and Design, vol. 73, pp. 65–76, 2013.

C. H. Thai, T. Rabczuk, and H. Nguyen-Xuan, “A rotation-free isogeometric analysis for composite sandwich thin plates,” International Journal of Composite Materials, vol. 3, no. 6A, pp. 10–18, 2013.

C. V. Verhoosel, M. A. Scott, T. J. R. Hughes, and R. de Borst, “An isogeometric analysis approach to gradient damage models,” International Journal for Numerical Methods in Engineering, vol. 86, no. 1, pp. 115–134, 2011.

P. Fischer, M. Klassen, J. Mergheim, P. Steinmann, and R. Müller, “Isogeometric analysis of 2D gradient elasticity,” Computational Mechanics, vol. 47, no. 3, pp. 325–334, 2011.

W. A. Wall, M. A. Frenzel, and C. Cyron, “Isogeometric structural shape optimization,” Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 33–40, pp. 2976–2988, 2008.

P. N. Nielsen, A. R. Gersborg, J. Gravesen, and N. L. Pedersen, “Discretizations in isogeometric analysis of Navier–Stokes flow,” Computer Methods in Applied Mechanics and Engineering, vol. 200, no. 45-46, pp. 3242–3253, 2011.

T. Temizer, P. Wriggers, and T. J. R. Hughes, “Contact treatment in isogeometric analysis with NURBS,” Computer Methods in Applied Mechanics and Engineering, vol. 200, no. 9–12, pp. 1100–1112, 2011.

E. D. Luycker, D. J. Benson, T. Belytschko, Y. Bazilevs, and M. C. Hsu, “X-FEM in isogeometric analysis for linear fracture mechanics,” International Journal for Numerical Methods in Engineering, vol. 87, no. 6, pp. 541–565, 2011.

S. S. Ghorashi, N. Valizadeh, and S. Mohammadi, “Extended isogeometric analysis for simulation of stationary and propagating cracks,” International Journal for Numerical Methods in Engineering, vol. 89, no. 9, pp. 1069–1101, 2012.

D. J. Benson, Y. Bazilevs, E. de Luycker et al., “A generalized finite element formulation for arbitrary basis functions: from isogeometric analysis to XFEM,” International Journal for Numerical Methods in Engineering, vol. 83, no. 6, pp. 765–785, 2010.

G. Haasemann, M. Kästner, S. Prüger, and V. Ulbricht, “Development of a quadratic finite element formulation based on the XFEM and NURBS,” International Journal for Numerical Methods in Engineering, vol. 86, no. 4-5, pp. 598–617, 2011.

S. S. Ghorashi, N. Valizadeh, and S. Mohammadi, “Analysis of cracked orthotropic media using the extended isogeometric analysis (XIGA),” in Proceedings of the 2nd European Conference on Extended Finite Element (XFEM '11), Cardiff University, Cardiff, United Kingdom, June 2011.

S.-M. Hu, Y.-F. Li, T. Ju, and X. Zhu, “Modifying the shape of NURBS surfaces with geometric constraints,” Computer-Aided Design, vol. 33, no. 12, pp. 903–912, 2001.

A. Asadpoure, S. Mohammadi, and A. Vafai, “Crack analysis in orthotropic media using the extended finite element method,” Thin-Walled Structures, vol. 44, no. 9, pp. 1031–1038, 2006.

P. Liu, T. Q. Bui, Ch. Zhang, T. T. Yu, G. R. Liu, and M. V. Golub, “The singular edge-based smoothed finite element method for stationary dynamic crack problems in 2D elastic solids,” Computer Methods in Applied Mechanics and Engineering, vol. 233–236, pp. 68–80, 2012.

T. Menouillard and T. Belytschko, “Dynamic fracture with meshfree enriched XFEM,” Acta Mechanica, vol. 213, no. 1-2, pp. 53–69, 2010.

T. Elguedj, A. Gravouil, and H. Maigre, “An explicit dynamics extended finite element method. Part 1: mass lumping for arbitrary enrichment functions,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 30–32, pp. 2297–2317, 2009.

T. Menouillard, S.-H. Song, Q. Duan, and T. Belytschko, “Time dependent crack tip enrichment for dynamic crack propagation,” International Journal of Fracture, vol. 162, no. 1-2, pp. 33–49, 2010.

F. García-Sánchez, C. Zhang, and A. Sáez, “A two-dimensional time-domain boundary element method for dynamic crack problems in anisotropic solids,” Engineering Fracture Mechanics, vol. 75, no. 6, pp. 1412–1430, 2008.

D. Motamedi and S. Mohammadi, “Dynamic analysis of fixed cracks in composites by the extended finite element method,” Engineering Fracture Mechanics, vol. 77, no. 17, pp. 3373–3393, 2010.

N. Nguyen-Thanh, H. Nguyen-Xuan, S. Bordas, and T. Rabczuk, “Isogeometric analysis using polynomial splines over hierarchical T-meshes for two-dimensional elastic solids,” Computer Methods in Applied Mechanics and Engineering, vol. 200, no. 21-22, pp. 1892–1908, 2011.

L. Chen, G. R. Liu, Y. Jiang, K. Y. Zeng, and J. Zhang, “A singular edge-based smoothed finite element method (ES-FEM) for crack analyses in anisotropic media,” Engineering Fracture Mechanics, vol. 78, no. 1, pp. 85–109, 2011.