# Blog

With the ESIS blog we want to promote and intensify relevant scientific discussions on recent publications in Engineering Fracture Mechanics. The blog is hosted by IMechanica on:

http://imechanica.org/node/9794

The editors, Professors Karl-Heinz Schwalbe and Tony Ingraffea, support this initiative. ESIS hopes that this blog will achieve the following objectives:

• To initiate scientific discussions on relevant topics by highlighting, leaving comments, suggestions, questions etc. related to recent publications;

• To suggest re-reading, re-examination, comparison of results from the past, that may be overlooked by the authors or have fallen into general oblivion;

• To promote and give reference to groups with similar or related scientific goals and to promote collaboration, as well as bridging gaps between different disciplines;

• To focus attention on new ideas that may face a risk of drowning in the noise of today's extensive scientific production;

Per Ståhle

## Latest posts from ESIS Blog on IMechanica

#### Discussion of fracture paper #19 - Fracture mechanical properties of graphene

Extreme thermal and electrical conductivity, blocks out almost all gases, stiff as diamond and stronger than anything else. The list of extreme properties seems never ending. The paper Growth speed of single edge pre-crack in graphene sheet under tension, Jun Hua et al., Engineering Fracture Mechanics 182 (2017) 337–355 , deals with the fracture mechanical properties of graphene. A sheet of armchair graphene can be stretched up to 15 percent which is much for a crystalline material but not so much when compared with many polymers. The ultimate load, on the other hand, becomes huge almost 100 GPa or more. Under the circumstances, it is problematic to say the least, that the fracture toughness is that of a ceramic, only a few MPam^(1/2). Obviously cracks must be avoided if the high ultimate strength should be useful. Already a few microns deep scratches will bring the strength down to a a few hundred MPa. The research group consisting of Jun Hua, Qinlong Liu, Yan Hou, Xiaxia Wu and Yuhui Zhang from the dept. of engineering mechanics, school of science, Xi’an University of Architecture and Technology, Xi’an, China, has studied fast crack growth in a single atomic layer graphene sheet with a pre-crack. They are able to use molecular dynamics simulations to study the kinetics of a quasi-static process. They pair the result with continuum mechanical relations to find crack growth rates. A result that provide confidence is that the fracture toughness obtained from molecular primitives agrees well with what is obtained at experiments. The highlighted results are that the crack growth rate increases with increasing loading rate and decreasing crack length. The tendencies are expected and should be obtained also by continuum mechanical simulations, however then not be first principle and requiring a fracture criterion. Another major loss would be the possibility to directly observe the details of the fracture process. According to the simulation results the crack runs nicely between two rows of atoms without branching or much disturbances of the ordered lattice. The fracture process itself would not be too exciting if it was not for some occasional minor disorder that is trapped along the crack surfaces. The event does not seem to occur periodically but around one of ten atoms suffers from what the authors call abnormal failure. Remaining at the crack surface are dislocated atoms with increased bond orders. All dislocated atoms are located at the crack surface. The distorted regions surrounding solitary dislocated carbon atoms are small. A motivated question would be if the dissipated energy is of the same order of magnitude as the energy required to break the bonds the connects the upper and lower half planes before fracture. Can this be made larger by forcing the crack to grow not along a symmetry plane as in the present study. Without knowing much about the technical possibilities I assume that if two graphene sheets connected to each other rotated so that the symmetry planes do not coincide, the crack would be forced to select a less comfortable path in at least one of the sheets. Everyone with comments or questions is cordially invited raise their voice. Per Ståhle

31.Dec.2017

#### Discussion of fracture paper #18 - A crack tip energy release rate caused by T-stress

A T-stress is generally not expected to contribute to the stress intensity factor because its contribution to the free energy is the same before and after crack growth. Nothing lost, nothing gained. Some time ago I came across a situation when a T-stress, violates this statement. The scene is the atomic level. As the crack is producing new crack surfaces the elastic stiffness in the few atomic layers closest to the crack plane are modified. This changes the elastic energy which could provide, contribute to or at least modify the energy release rate. If the energy is sufficient depends on the magnitude of the T-stress, the change of the elastic modulus and how many atomic layers that are involved. If I should make an estimate it would be that the energy release rate is the change of the T-stress times the fraction of change of the elastic modulus times the square root of the thickness on the affected layer. Assuming that the T-stress is a couple of GPa, the change of the fraction of change of the elastic modulus is 10% and the affected layer is around ten atomic layers one ends up with 100kPa m^(1/2). Fairly small and the stress and its change are taken at its upper limits but still it is there. The only crystalline material I could find is ice with a toughness of the same level. Other materials are affected but require some additional remote load. Interestingly enough I came across a paper describing a different mechanism leading to a T-stress contribution to the energy release rate. The paper is: Zi-Cheng Jiang, Guo-Jin Tang, Xian-Fang Li, Effect of initial T-stress on stress intensity factor for a crack in a thin pre-stressed layer, Engineering Fracture Mechanics, pp. 19-27. This is a really read worthy paper. The reasons for the coupling between the T-stress and the stress intensity factor is made clear by their analysis. The authors have an admirable taste for simple but accurate solutions. The paper describes a crack with a layer of residual stress, that gives a T-stress in the crack tip vicinity. As the crack advances increasing more material end up behind the crack tip rather than in front of it. The elastic energy density caused by the T-stress is larger in front of the tip than it is behind it. The energy released on the way and can only disappear at the singular crack tip, not anywhere else in the elastic material. The reason for the energy release is the assumed buckling in the direction perpendicular to the crack plane. An Euler-Bernoulli beam theory is used to calculate the contribution to the energy release rate. Having read the paper I realise that in a thin sheet buckling out of its own plane in the presence of a crack and a compressive T-stress there will be energy released that should contribute to crack growth. The buckling will give a more seriously distorted stress state around the crack tip, but never the less. In this case the buckling area would be proportional to the squared crack length in stead of crack length times the height of the layer as in the Jiang et al. paper. The consequence is that the contribution to the stress intensity factor should scale with the T-stress times square root of the crack length. Suddenly I feel that it would be very interesting to hear if anyone, maybe the authors themselves, know of other mechanisms that could lead to this kind of surprising addition to the energy release rate caused by T-stresses. It would be great if we could add more to the picture. Anyone with information is cordially invited to contribute. Per Ståhle

01.Nov.2017

#### Discussion of fracture paper #17 - What is the second most important quantity at fracture?

No doubt the energy release rate comes first. What comes next is proposed in a recently published study that describes a method based on a new constraint parameter Ap . The paper is: Fracture assessment based on unified constraint parameter for pressurized pipes with circumferential surface cracks, M.Y. Mu, G.Z. Wang, F.Z. Xuan, S.T. Tu, Engineering Fracture Mechanics 175 (2017), 201–218 The parameter Ap is compared with established parameters like T , Q etc. The application is to pipes with edge cracks. I would guess that it should also apply to other large structures with low crack tip constraint. As everyone knows, linear fracture mechanics works safely only at small scales of yielding. Despite this, the approach to predict fracture by studying the energy loss at crack growth, using the stress intensity factor K I and its critical limit, the fracture toughness, has been an engineering success story. K I captures the energy release rate at crack growth. This is a well-founded concept that works for technical applications that meet the necessary requirements. The problem is that many or possibly most technical applications hardly do that. The autonomy concept in combination with J -integral calculations, which gives a measure of the potential energy release rate of a stationary crack, widens the range of applications. However, it is an ironin that the J -integral predicts the initiation of crack growth which is an event that is very difficult to observe, while global instability, which is the major concern and surely easy to detect, lacks a basic single parameter theory. For a working concept, geometry and load case must be classified with a second parameter in addition to K I or J . The most important quantity is no doubt the energy release rate, but what is the second most important. Several successful parameters have been proposed. Most of them describe some type of crack tip constraint, such as the T -stress, Q , the stress triaxiality factor h , etc. A recent suggestion that, as it seems to me, have great potential is a measure of the volume exposed to high effective stress, Ap . It was earlier proposed by the present group GZ Wang and co-authors. Ap is defined as the relative size of the region in which the effective stress exceeds a certain level. As pointed out by the authors, defects in large engineering structures such as pressure pipes and vessels are often subjected to a significantly lower level of crack tip constraint than what is obtained in laboratory test specimens. The load and geometry belong to an autonomy class to speak the language of KB Broberg in his book "Fracture and Cracks". The lack of a suitable classifying parameter is covered by Ap . The supporting idea is that K I or J describe the same series of events that lead to fracture both in the lab and in the application if the situations meet the same class requirements, i.e. in this case have the same Ap . The geometry and external loads are of course not the same, while a simpler and usually smaller geometry is the very idea of the lab test. The study goes a step further and proposes a one-parameter criterion that combines the K I or J with Ap by correlation with data. The method is reinforced by several experiments that show that the method remains conservative, while still avoiding too conservative predictions. The latter of course makes it possible to avoid unnecessary disposal and replacement or repair of components. The authors' conclusions are based on experience of a particular type of application. I like the use of the parameter. I guess more needs to be done extensively map of the autonomy classes that is covered by the method. I am sure the story does not end here. A few questions could be sent along: Like "Is it possible to describe or give name to the second most important quantity after the energy release rate?" The paper mentions that statistical size effects and loss of constraint could affect Ap . Would it be possible to do experiments that separates the statistical effect from the loss of constraint? Is it required or even interesting? It would be interesting to hear from the authors or anyone else who would like to discuss or comment the paper, the proposed method, the parameter or anything related. Per Ståhle

18.Oct.2017

#### Discussion of fracture paper #16 - What is wrong with pure mode I and II? A lot it seems.

It is common practice when solving boundary value problems to split the solution into a symmetric and an antisymmetric part to temporarily reduce the number of variables and the mathematical administration. As soon as the symmetric problem is solved, the antisymmetric problem, or vice versa, is almost solving itself. Any problem can be split into a symmetric and an antisymmetric part which is a relief for anyone who analyses mixed cases. It gives a clearer view but it is an academic exercise while nature usually doesn't have any comprehension of symmetry and antisymmetry. Fracture is no exception. The fracture processes will be activated when sufficient conditions are fulfilled. Even the smallest deviation from the pure mode I or II caused by geometry or load will not affect the conditions at the crack tip in any decisive way. Everything is almost pure mode I or II and it may be convenient ignore the small deviation and still treat the problem as a pure case. This seems simple enough but the paper reviewed tells that it has been a tripwire for many. The selected paper is the recently published: "An improved definition for mode I and mode II crack problems" by M.R. Ayatollahi, M. Zakeri in Engineering Fracture Mechanics 175 (2017) 235–246. The authors examine a power series expansion for an Airy stress function about the crack tip. The series give stress as a sum of powers r-1/2, 1, r1/2, r, etc. of the distance to the crack tip. Each term has an known angular dependence. The application is to a plane crack with any in-plane load. The series starts with a square root singular term while it is assumed that the crack tip is sharp and the material is linear elastic. The assumption requires that the geometrical features of the crack tip and the nonlinear region is not visible from where the expansion with some accuracy describes the stress field. The problem that the authors emphasise is that the splitting in symmetric and antisymmetric modes that leads to two similar expansions of the radial power functions with symmetric and antisymmetric angular functions. The representations so far has been called pure if the solution is strictly symmetric or strictly antisymmetric, i.e. the notation has been pure mode I and pure mode II. The problem is that not only seldom, has a vanishing mode I stress intensity factor misled investigators to drop all symmetric terms of the series expansion. Also mode II has been unfairly treated in the same way. The most striking problem is of course when the constant stress acting along the crack plane, the T-stress, by mistake is neglected. The authors are doing a nice work sorting this out. They describe a range of cases where one stress intensity factor vanishes but for sure the crack tip stress state is neither strictly symmetric nor strictly antisymmetric. They also provide quite many examples to demonstrate the necessity to consider the T-stress even if the mode I singular stress term is absent. I commend the authors for doing a conscientious work. I f I should bring up something where different positions may be assumed it would be the selection of the series. The powers of r-1, r-3/2, r-2 etc are never mentioned and I agree that it is not always necessary. It should be commonly known that a sharp crack, a linear elastic material and traction free crack surfaces says it. There cannot be any stronger singularities than r-1/2. However, isn't one consequence that close enough to the crack tip any constant stress should be insignificant as compared to the singular stress terms. If so, it should not have any significant effect on the stresses closest to the crack tip and neither affect the fracture processes nor the selection of crack path. On the other hand, if the constant term has a real influence on the course of events, that would as far as I understand mean that the nonlinear region has to have a substantial extent so that its state is given by both singular terms and the T-stress. The contradiction is then that the stronger singular terms r-1, r-3/2, etc. cannot be neglected. These terms are there. Already the r-1 term seems obvious if the crack has grown because of the residual stress caused by plastic strain along the crack surface that in the wake region behind the crack tip. Also, the region of convergence, which is at most the length of the crack, is another pothole. Outside the convergence region a different series or an analytical continuation, may be used. For the series expansions the symmetric and antisymmetric solutions have to be treated as well, with the difference that there are constant stresses in both symmetric and antisymmetric modes that have to be included. It would be interesting to hear if there are any thoughts regarding this. Per Ståhle P.S. On the courtesy of Elsevier there is a 3 month promotional access to the latest article in the blog, meaning the articles are freely available to everyone. Now everyone who wishes to comment or discuss the paper here can do so. (Dr. Kumar, I hope you are reading this).

20.Jul.2017

#### Discussion of fracture paper #15 - Designing for crack arrest

Everyone loves an elegant engineering solution. It is particularly true when the alternatives are terrifying. In the paper: ”Brittle crack propagation/arrest behaviour in steel plate – Part I: Model formulation” by Kazuki Shibanuma, Fuminori Yanagimoto, Tetsuya Namegawa, Katsuyuki Suzuki, Shuji Aihara in Engineering Fracture Mechanics, 162 (2016) 324-340 . a team from University of Tokyo proposes a model for prediction of the arrest of propagating brittle cracks in steel plates. The approach, in spite of its simplicity, captures the physics of the fracture process. The model formulates the energy release rate in simple and comprehensible terms and gives accurate predictions. The theory is validated on several experiments described in a subsequent paper, a ”Part II: Experiments and model validation” also in Engineering Fracture Mechanics. The characteristics are those of a pilot study with the goal to provide a design tool for predicting crack arrest in steel plates. In the model, the energy to complete the fracture process is at most what is left of the released energy when the work of plastic deformation and the part of the kinetic energy that is reflected away from crack tip region have been covered. The energy dissipation at plastic deformation is reduced at increasing crack tip velocity while the opposite applies to the dissipated kinetic energy. The energy required for the fracture processes is supposed to be constant. If it at some velocity is more than required then the energy is in balance only at a single stable higher crack tip velocity. If crack growth is initiated then the crack accelerates until the energy balance is obtained. When the crack subsequently loose driving force or require additional work, caused, e.g., by elevated temperature which decreases the material viscosity or by whatever, the crack decelerates until zero velocity or until the minimum energy release rate is obtained and the crack arrest comes abruptly. Surface-ligaments are assumed to consume a serious part of the available energy. The slower the crack grows the wider these ligaments become which rapidly increases the plastic dissipation. Finally the energy balance and the stability of the crack tip velocity cannot be maintained and the crack will come to a stop. Considering that one has to keep trac k of the complicated sequence of processes that keep the crack growing, it seemed obvious to me that this would end up in a horrible and time consuming analysis. Then, to my surprise, the investigators present an ingenious solution that simplifies the analysis a lot. It is based on three assumptions: 1) that the crack front is assumed to be straight through the plate, 2) that the unbroken side-ligaments are regarded as integrated parts of the crack front, and 3) that the evaluation of the state of the crack front in done at the plate mid-plane. In the subsequent part II the functionality of the model is verified. The validation is performed on different grades of steel that are exposed to different load levels. The authors believe that this model can be used to establish a design strategy for steel plates. I too believe that, even if more possibly needs to be done to qualify the method as a design standard. I understand that the authors are familiar with the series of wide plate experiments on crack arrest in very large specimens (around 11x1x0.1 m3) reported by Naus et al., NUREG/CR-4930 ORNL-6388, Oakridge Laboratories, USA, 1987. In the aftermath of the experiments a variety of models where proposed. An interesting observation made by D. Alexander and I.B. Johansson at Oakridge Labs when they examined the crack surfaces was that remains of plastic deformation framed the cloven grains. The guess was that this was remaining parts of broken ligaments between the crack surfaces and that these ligaments were ripped apart during the fracture process. The area covered by these remains was clearly increasing with decreasing crack tip velocity. Just before crack arrest they could cover as much as 10 to 20 % of the ”brittle” part of the crack surface. I have a feeling that this may mean something. The plastic ligaments per se consume large amounts of energy and with increasing fractions they might influence the crack tip velocity at arrest. Only 10% may seem as small or even insignificant, but considering that the plastic ligaments that bridge a crack may consume many times more energy than the pure cleavage of the remaining 90%, even 10% must be important. It would be interesting to know if the authors observed any remains of plastic ligaments. If so, did the fraction of them change in any systematic way during crack growth? Any contribution to this blog is gratefully acknowledged. Per Ståhle

20.Nov.2016