Crack PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 23

Crack Growth Analysis Using BEASY

S.M. Niku, R.A. Adey Computational Mechanics BEASY, Ashurst Lodge, Ashurst, Southampton, SO4 2AA

ABSTRACT
Catastrophic failure of engineering structures is caused by cracks that extend beyond a safe size. Cracks are present to some extent in all structures and can lead to failure or decrease in structural strength if they grow. This paper described an automatic procedure for predicting the growth of cracks in two dimensional structures and the key data for fracture mechanics based design. The stress intensity factors. The procedure is implemented in a general purpose engineering analysis system BEASY. The system uses a dual boundary element technique to simplify the modeling of the cracks and to provide highly accurate stress intensity factors. Applications are presented for single and multiple cracks.

1.0 INTRODUCTION
The cost to industry of fracture was recently estimated in a report of the US Department of Commerce entitled "The economic effect of fracture in the United States". In this report it estimated that the cost of fracture was $119 billion dollars per year (4% of gross national product). It further estimated that approximately one third could be saved through the use of current fracture control technology and a further 25% could be saved through fracture related research. Therefore the annual cost of fracture could be halved by the application of better design tools based on fracture mechanics technology. The combination of the easy to use and apply boundary element approach in BEASY and the new fracture growth tools incorporated in the system are aimed to satisfy this goal. Catastrophic fracture failure of engineering structures is caused by cracks that extend beyond a

Crack Growth Analysis Using BEASY

Page 1

safe size. Cracks, present to some extent in all structures, either as a result of manufacturing fabrication defects or localised damage in service, may grow by mechanisms such as fatigue, stresscorrosion or creep. The crack growth leads to a decrease in the structural strength. Thus, when the service loading cannot be sustained by the current residual strength, fracture occurs leading to the failure of the structure. Fracture, the final catastrophic event that takes place very rapidly, is preceded by crack growth which develops slowly during normal service conditions, mainly by fatigue due to cyclic loading. Damage tolerance assessment is a procedure that defines whether a crack can be sustained safely during the projected service life of the structure. Damage tolerance assessment is, therefore, required as a basis for any fracture control plan, generating the following information, upon which fracture control decisions can be made: The effect of cracks on the structural residual strength, leading to the evaluation of their maximum permissible size. The cracks growth as a function of time, leading to the evaluation of the life of the cracks to reach their maximum permissible size, from which the safe operational life of the structure is defined. Linear elastic fracture mechanics can be used in damage tolerance analyses to describe the behaviour of cracks. The fundamental assumption of linear elastic fracture mechanics is that the crack behaviour is determined solely by the values of the stress intensity factors which are a function of the applied load and the geometry of the cracked structure. The stress intensity factors, thus play a fundamental role in linear elastic fracture mechanics applications. Crack-growth processes are simulated with an incremental crack-extension analysis. For each increment of the crack extension, a stress analysis is carried out and the stress intensity factors are evaluated. The crack path, predicted on an incremental basis, is computed by a criterion defined in terms of the stress intensity factors. The boundary element method is well established as a powerful solution tool for fracture mechanics. See Aliabadi & Rooke [2]. The reason for its success is the boundary only representation, the high accuracy and the methods ability to represent the high stress fields near the crack front. In this paper the implementation within the BEASY boundary element system is described and application presented showing the wide range of applications possible.

2.0 THEORETICAL FOUNDATION


The theoretical foundation of the dual boundary elements and crack growth algorithm follows closely the work of Portela and Aliabadi [1]. The authors will not repeat the theoretical basis but would refer the reader to reference [1] for the detailed basis of the crack growth and dual boundary elements and to Brebbia and Dominguez [3] for a description of the boundary element

Crack Growth Analysis Using BEASY

Page 2

method.

Zone I

Zone II

a) Dual Boundary Standard Boundar

b)

FIGURE 1. Boundary element models of edge crack. a)Shows the crack modeled using dual boundary elements. b) Shows the crack modeled using two zones.

2.1 CRACK MODELING


A crack can be represented in a boundary element model using two main approaches. The traditional approach requires the user to define a zone boundary along the crack surface and continue this through the body of the components being studied. This can be clearly seen in Fig. 1b where the problem is split into two zones and the edge crack is extended by a zone interface (shown as a dotted line) across to another external boundary. This technique while not requiring any special theoretical development places on extra burden on the user. The second approach recently developed is to use dual boundary elements to represent the crack. Fig. 1a shows the model of the edge crack using this approach. In this case the modeling is extremely simple and economical. The crack is represented by two elements occupying the same physical location, each element representing a face of the crack.

Crack

FIGURE 2. Dual boundary element representation of the crack.

Crack Growth Analysis Using BEASY

Page 3

2.2 STRESS INTENSITY FACTORS


The system automatically computes the stress intensity factor using both J integral and displacement formula. This can be used for computing both Mode I and Mode II intensity factors for edge cracks and embedded cracks. (Note: The system currently uses displacement and stress formulas for computing stress intensity factors for three dimensional problems). The J integral path is automatically computed by the software and follows the crack tip as it grows into the material.

Crack

Integral path

FIGURE 3. Contour for J integral.

2.3 INCREMENTAL CRACK-EXTENSION ANALYSIS


The incremental crack-extension analysis assumes a piece-wise linear discretisation of the unknown crack path. For each increment of the crack extension, the dual boundary element method is applied to carry out a stress analysis of the cracked structure and the J-integral is the technique used for the evaluation of the stress intensity factors. The steps of this basic computational cycle, repeatedly executed for any number of crack-extension increments, are summarised as follows: Carry out a stress analysis of the structure Compute the stress intensity factors Compute the direction of the crack-extension increment Extend the crack one increment along the direction computed in the previous step Repeat all the above steps sequentially until a specified number of crack-extension increments is reached For the sake of simplicity, the increment of crack extension is discretised with a fixed number of new boundary elements. In order to avoid numerical problems, concerned with the relative size of neighbouring elements, the crack increment length is kept between convenient limiting bounds, defined in terms of the size of the crack-tip element. Apart from this constraint, the length of the crack extension increment may be defined as the result of a compromise between accuracy and

Crack Growth Analysis Using BEASY

Page 4

computational cost; the smaller the crack increment the more accurate and expensive is the analysis. The results obtained from an incremental crack-extension analysis are a crack path diagram and diagrams of the stress intensity factor variation along the crack path.

2.4 CRACK-EXTENSION CRITERION


Several criteria have been proposed to describe the mixed-mode crack growth. Among them, one of the most commonly used is based on the maximum principal stress at the crack tip, Erdogan & Sih. The maximum principal stress criterion postulates that the growth of the crack will occur in a direction perpendicular to the maximum principal stress. Thus, the local crack-growth direction is determined by the condition that the local shear stress is zero. In practice this requirement gives a unique direction irrespective of the length of the crack extension increment. Therefore the procedure adopted in the system is to use a predictor corrector technique to ensure the crack path is unique and independent of the crack extension increment used.

2.5 SOLUTION EFFICIENCY


The crack growth simulation process requires the prediction of a new solution to the whole problem each time the model geometry is changed. Therefore during the incremental growth of the boundary elements representing the crack hundreds of new solutions are required. While this may not produce a high computational demand for small problems, it is essential for larger models that an efficient solution strategy is adopted. The tasks for the analysis can be summarised as follows: Perform stress analysis of model. Compute stress intensity factors and crack growth direction. Modify the model geometry to extend the crack. Repeat starting at step 1. The bulk of the computation is involved in step 1 as the boundary element equation system has to be formed and solved. By careful choice of the equation solution technique, the stages can be modified as follows: Perform stress analysis of model. Compute stress intensity factors and crack growth direction. Modify the model geometry to extend the crack. Add to the system of equation the coefficient corresponding to the elements added to the crack increment. Solve the modified system of equation and compute stresses and displacements. Repeat from step 2.

Crack Growth Analysis Using BEASY

Page 5

Using this procedure the computation effort is dramatically reduced as only a small part of the system of equation has to be created at each step. Fig. 4. Also only the additional equations have to be factual at each step.

Initial Geometry Ist Crack 2nd Crack 3rd Crack

FIGURE 4. Schematic view of the system of equation for difference crack lengths.

3.0 BEASY
The new crack growth software forms part of a family of application under the BEASY name. The system is integrated with major CAD systems and pre and post processors.
Mechanical Design

Fracture Mechanics

Stress Analysis Product Simulation Thermal Analysis Corrosion Simulation Based Design Analysis Acoustic Analysis

Contact Analysis

Electrostatic Analysis

Field Analysis

Crack Growth Analysis Using BEASY

Page 6

4.0 NUMERICAL VALIDATION


The accurate prediction of the crack growth is highly dependent upon the accuracy of the stress intensity factors. Any minor innaccuracy in these values is accumulated during crack growth process which can lead to a completely unrealistic prediction. Therefore considerable care has been taken in selecting and implementing the techniques used to compute the SIF values. A large number of tests have been carried out to verify the accuracy and reliability of the software for simple modes and mixed mode cracks. Comparisons of the stress intensity factors have been made with a number of published test solutions and the accuracy found to be very high. For single mode and mixed mode cracks the solution has been found to be almost exact.

4.1 Edge crack under Pure Mode I Loading

h a h

FIGURE 5. Edge crack under Pure Mode I loading.

In this problem Fig. 5 the growth of straight edge crack is investigated. Under the pure mode I load it is expected that the crack will grow along the direction of the initial crack. The problem

Crack Growth Analysis Using BEASY

Page 7

was analysed using the BEASY fracture software and the results are shown in Fig. 6 and Fig. 7.

FIGURE 6. Predicted crack propagation under Pure Mode I

FIGURE 7. Predicted deflected shape after cracks propagation. Note effectiveness of the dual BEM is representing the crack faces.

Crack Growth Analysis Using BEASY

Page 8

It can be clearly seen that the prediction follows the direction expected.

4.2 Edge crack under pure mode II loading.


In this case the model was the same as that used in the Mode I case but the model was subject to Pure Mode II. Fig. 8 shows the predicted crack paths and clearly shows that the crack path follows the line expected. The deformed shape is shown in Fig 9

FIGURE 8. Predicted crack path due to Pure Mode II loading on plate with edge crack

FIGURE 9. Deformed shape of plate following crack growth under Pure Mode II loading.

Crack Growth Analysis Using BEASY

Page 9

The variation of Mode I & Mode II SIF with respect to crack increments is shown in Figure 10. As expected KI increases as crack grows whileKII is reduceds to zero.

3 2.5 2 1.5 1 0.5 0 -0.5

Stress Intensity Factors

KI

KII

No of crack extensions
FIGURE 10. Stress Intensity factors along the crack.

4.3 Gear Tooth


The growth of a crack from the root of a gear is investigated in this case.

Crack Growth Analysis Using BEASY

Page 10

The problem was originally studied by Flasker and Pehan [5] and Fig. 11 shows the mesh used

FIGURE 11. Model of gear tooth showing initial crack.

and the loading. The crack growth was predicted using the software and the crack path is shown in Fig. 12. This compares well with experimental results and manual finite element calculations reported by [5].

FIGURE 12. Predicted crack path.

Crack Growth Analysis Using BEASY

Page 11

The predicted deflected shape is shown in Fig. 12.

FIGURE 13. Predicted deflected shape.

5.0 Further Applications


This section contains a series of applications solved using BEASY Crack Growth.

5.1 Cross Plate


In this application a crack is initiated from the corner of the plate and the crack path predicted.

Crack

FIGURE 14. Diagram of cross plate showing loading

Under the applied loading the crack is expected to grow towards the corner vertically above the

Crack Growth Analysis Using BEASY

Page 12

crack.

FIGURE 15. Model of cross plate showing initial crack

FIGURE 16. Cross plate showing the crack path

The crack path has been automatically computed and as shown in figure 16. It grows in the direction expected.

Crack Growth Analysis Using BEASY

Page 13

5.2 Multiple cracks


BEASY can be used to solve problems with multiple cracks and the problem can be split into multiple zones representing different parts or materials. In this example a series of cracks initiating from circular holes in a plate are analysed. The stress intensity factors are automatically computed for each crack. A further study to compute the crack propagation can also be made.

FIGURE 17. BEASY model showing model geometry and mesh

Crack Growth Analysis Using BEASY

Page 14

FIGURE 18. Deformed shape of multiple cracks

5.3 Embedded crack propagating at both ends


This example demonstrates the propagation of multiple cracks. The plate is subjected to a vertical load. applied at both ends of the plate.

FIGURE 19. Initial model of embedded crack

Crack Growth Analysis Using BEASY

Page 15

The crack can be seen more clearly in Figure 20. Note the small number of elements required to solve the problem.

FIGURE 20. Close up view of the crack

FIGURE 21. Predicted paths of the two cracks

Crack Growth Analysis Using BEASY

Page 16

The complete stress distribution and the deformed shape is also predicted by BEASY.

FIGURE 22. Deformed shape of the model

5.4 Multi hole plate

FIGURE 23. BEASY model showing the initial crack position

Crack Growth Analysis Using BEASY

Page 17

In this application the path of a crack propagating between a number of holes in a plate is studied. The complete history of the stress intensity factors is computed along the path as well as the crack tip.

FIGURE 24. Predicted crack path

FIGURE 25. Exaggerated deformed shape display

Crack Growth Analysis Using BEASY

Page 18

A very useful feature of BEASY is the accuracy report. This describes the convergence of the solution and the equilibrium. The accuracy reported describes the overall accuracy of the stress solution as a percentage. This should not be confused with the accuracy of the stress intensity factors which are computed with very high accuracy.

FIGURE 26. Display of the accuracy report

5.5 Cracks growing towards each other


A very interesting example is the case of two cracks propagating towards each other. This can be solved very easily using the system as can be seen in the Figure 26 where the mesh is shown.

Crack Growth Analysis Using BEASY

Page 19

FIGURE 27. Initial boundary element model

The initial deformed shape of the plate before the crack grows can be predicted and is shown in Figure 28.

FIGURE 28. Initial deformed shape

Crack Growth Analysis Using BEASY

Page 20

The predicted crack path is computed.

FIGURE 29. Predicted crack path

The crack path becomes very complicated as the cracks become very long as can be seen in the deformed shape display.

FIGURE 30. Deformed shape

Crack Growth Analysis Using BEASY

Page 21

The deformed shape shows a complex stress distribution as the cracks have approached each other and then passed. The cracks have not maintained complete symmetry as they approach the edge of the plate. A clearer view is shown by zooming into the crack area.

FIGURE 31. Enlarged view of the cracks deformed shape

6.0 CONCLUSION
The paper has described a powerful system for the prediction of fracture data and crack growth in two dimensional structures. The techniques have shown high accuracy and the ability to correctly predict the growth of cracks. The system is particularly powerful because of the very simple modeling and complete lack of meshing problems associated with FEM solution of these problems. The reported work is part of an ongoing project which is currently investigating three dimensional crack growth.

ACKNOWLEDGEMENT
This work was partly funded under Project MASS of the European Community. Partners in the project include Computational Mechanics BEASY, Microparts, Technical University of Berlin and Wessex Institute of Technology. Additional support was provided by the European Space Agency. The authors acknowledge the support of Dr Aliabadi during this work.

REFERENCES
1. Portela A., Aliabadi M.H. & Rooke D.P. Dual Boundary Element incremental analysis of

Crack Growth Analysis Using BEASY

Page 22

crack propagation. Int Journal Computers and Structures Vol. 46 1993. pp 237-247. 2. Aliabadi M.H., Rooke D.R., Numerical Fracture Mechanics, Computational Mechanics Publications, SOUTHAMPTON, U.K. 1991. 3. Brebbia C.A., Dominguez J., Boundary Elements - An Introductory Course. Computational Mechanics Publications, SOUTHAMPTON, U.K. 1989. 4. Erdogan, F and Sih, G.C. On the crack extension in plates under plane loading and transverse shear. J. Basic Eng., 1963, 85, 519-527. 5. Flasker J., Pehan S., "Crack propagation in tooth root with variable loading". Comm in Num Meth in Eng. Vol. 9. 103-110 (1993).

Crack Growth Analysis Using BEASY

Page 23

You might also like