Fracture in Abaqus
Submitted by Michael Borowicz on Tue, 2007-03-27 09:57.
Good morning
I'm working in Abaqus on modelling crack propagation. Regards this i have two questions.
1)I would like to ask if it is possible for Abaqus to create automatic propagation in 2D. What should I do to receive the effect. I've defined the crack and in the history Output i've defined the crack initiation criterion. There are no warnings, no errors, and in visualization module the effects are miserable.
2)When I'm trying to plot Crack prop. direction I'm receiving only points in vertical for time=1s, and I would like to receive crack propagation for different state of times, different moments.
P.s. Sorry for my english, I hope everything is clear enough to understand this post. I hope someone can help me.
Regards
Michael Borowicz


prediction of crack propagation using ABAQUS
Hi Michael,
I am not sure whether I understand your question. ABAQUS surely has capability of predicting crack propagation. Particularly in ABAQUS/Explicit, element deletion technique is provided such that the element,where the failure criterion is locally reached, will be removed from the calculation. Thereby dynamic element deletion can be visualized as a crack progration.
It actually provides a variety of failure models for different materials and different problems. But for doing it correctly, you must do some benchmark tests to make sure the failure criterion that you used is reasonable. Except that, the numercial results very depend on mesh design.
I suggest you to check ABAQUS manual (ver6.5) section 11.2.8 and chapter 11.6 for details. Also in the ABAQUS/Example problems manual, problems 2.13 and 2.14 are both good examples that you should perform for being familiar with the method.
Hope this helpful for you.
Zhenyu
another way of handling crack in Abaqus
In addition to Zhenyu's reply, I want to add that Abaqus can handle crack using cohesive zone model. In old versions, one has to write an UEL to handle cohesive elements. In the new version (6.5), it allows the user to define the traction-separation behavior in the contact elements, so you can regard this as cohesive elements, to my understanding.
For the ad hoc problem, if you want to implement your own fracture criterion, the cohesive element can be easily implement using a template that Prof. Bower and I wrote. Please see paper #7 in http://web.utk.edu/~ygao7/publication.htm
two cracks reaching their strength peak
Dear Yanfei,
I am reading your paper
Y.F. Gao, A.F. Bower, "A simple technique for avoiding convergence problems in finite element simulations of crack nucleation and growth on cohesive interfaces." Modelling and Simulation in Materials Science and Engineering 12, 453-463, 2004
with great interest and want to raise some questions for discussion.
In your paper, you applied your technique (by introducing a small viscosity in the interface cohesive law) to a system having only one crack. Can this technique be applied to a system with many cracks, or say just two cracks, in the case when all those cracks are reaching their strength peak?
about cohesive zone model
Yanfei,
In my knowledge, cohesive zone model is very powerful for the problem where the crack path is known, such as the delamination of a bilayer (thinfilm on substrate ect.) where the crack is always propagating along the interface between two materials. While for the problem that the crack paths are too complicate to be predicted such as petalling process of a plate subjected to blast loads, the cohesive zone model cannot handle it in an efficient way.
Please correct me if I misunderstood something.
Zhenyu
advantages and disadvantages of coheisve zone model
Zhenyu,
You are right about the advantages and disadvantages of cohesive zone model. To my observation, here are two major advantages:
(1) crack path is predefined
(2) the system will be artificially more compliant with the addition of cohesive elements
The first one can be avoided by introducing cohesive elements everywhere, but for highly symmetric problem, the path will be mesh dependent. With cohesive elements everywhere, the second problem is more severe. So we need to make sure the initial response is almost rigid. But doing so will introduce snap-open instability. That's why Prof. Bower and I modified the traditional cohesive elements along this line.
thanks
yanfei
compliance of the system
Dear Yanfei,
You modified the traditional cohesive elements by introducing a small viscosity in the interface cohesive law.
Will this affect the compliance of the system?
cohesive zone formulations avoiding mesh bias
Hiya,
Nice and informative discussion.
Cohesive zone modelling always have some disadvantages as discussed by Yanfei and others.
One way of removing the mesh bias is to use cohesive elements everywhere but it is a tedious job and many problems regarding stability and convergence.:)
Simulation of cup–cone fracture using the cohesive model Engineering Fracture Mechanics, Volume 70, Issue 14, September 2003, Pages 1943-1961
I. Scheider and W. Brocks
To remove mesh bias and predict a crack path with cohesive zone model, following reference may be of interest to people using cohesive models and stuck due to mesh bias.
Mesh-independent discrete numerical representations of cohesive-zone models , Engineering Fracture Mechanics, Volume 73, Issue 2, January 2006, Pages 160-177
René de Borst, Joris J.C. Remmers and Alan Needleman
I still remember those days when i was using cohesive model in simulating interface fracture. Breath taking and sometimes frustrating, i would say.:)
Best Regards,
Amir
*******************************************************
Nature does not care whether you can describe a phenomenon by an equation. It just manifests them.
*******************************************************
XFEM for interface debonding
In the eXtended Finite Element Method (XFEM), the interface discontinuity can be well solved.
X-FEM can be embedded in ABAQUS using UEL
Dear Henry,
Maybe the method you metioned above is a good choice for problems about moving discontinuities. I am using X-FEM to simulate fracture process in ABAQUS. The user subroutine of UEL is used to construct special element with enriched DOFs. The X-FEM can make use of the powerful ability of solving nonlinear equations of ABAQUS or other FE commercial software.
Any open-source X-FEM subroutine for ABAQUS
Hi Xiujun,
Some of us have realized how powerful ABAQUS can can be, e.g. http://imechanica.org/node/470. It's such a good idea to take advantage of ABAQUS to perform X-FEM calculation via UEL. Is there any open-source X-FEM subroutine for ABAQUS? I noticed Prof. Y. Gao has posted some UEL subroutines on line (refer to publication 7). Is there more?
-XH
Joint-work is needed in developing program
Dear XH,
I have ever tried to search the open-source X-FEM subroutines for ABAQUS but can't find any more than you mentioned so far( Eugenio's work is expected). According to my experiences of writing UEL I find it is difficult to develop a robust program individually. A lot of work which have little relation to X-FEM itself takes much time. If there are a few ones to do the similar thing, maybe joint-work is a choice to avoid the repeated labor. I will be very glad to share my program with others.
Can XFEM-LSM be embedded in ABAQUS?
Dear Xiujun,
You propose a very appealing approach to simulate crack propagation by embedding XFEM into ABAQUS commercial code. As far as I know, XFEM can also be implemented by a combination of XFEM and Level-Set method, eg. the work by Ji and Dolbow (IJNME, 2002). This new version of XFEM-LSM can be used to simulate a variety of time-dependent moving boundary problems. The key of the strategy is to solve a time-dependent mechanical equation as well as a level-set update equation. My question concerns whether it is possible to implement solution of these time-dependent equations in ABAQUS. If that is a case, then simulation of dynamic crack propagation is possible! Can I have your comments?
Dear Jinxiong The question
Dear Jinxiong
The question you mentioned is interesting and forsighting. I am not sure I understand you correctly .The dynamic crack propagation problem may be understood in two ways: 1) considering the velocity of crack propagation itself;2) only taking account of the inertial force without considering the v of crack propagation. The latter is relatively simple and there is no foundamental differences but the dynamic loadings between the static problem and dynamic problem. XFEM-LSM embedded into ABAQUS can be used to simulate the problem in this case. I am trying to simulate the dynamic response of gravity dam with crack propagation under earthquake. As for the former I am not sure. But I think it is hopeful to realize it.
Crack propagation with high speed
Dear Xiujun,
My question is actually the first case you mentioned, i.e., unsteady crack motion with fast speed. In this case, crack velocity is an important factor and may be comparable with shear wave speed. One of Karma's phase field method simulation results show that the crack velocity,v, can reach 0.4c, where c is the shear wave speed. The conventional steady-state crack propagation criterion will be inaccurate in this case and crack branching and instability may observed. This problem is really challenging!
Jinxiong
Interface modelling to predict delamination pattern.
Hi Amir,
I am
doing the similar project like what you did. My project is to predict
the delamination pattern in the interface of a bi-layer system
(substrate and overlay) using ABAQUS. Substrate is made of normal
concrete and overlay is made of High-Performance concrete. Delamination
is induced by maintaining high temperature (say 95 F) in the enclosed
room due to which different materials used will absorb different
temperatures and because of this difference of temperature in two
layers, shear stresses are developed at the interface due to which the
delamination occurs.
My questions are: What are the properties required to input in ABAQUS?
How to model the interface?
How to input the temperature loads?
Any suggestions will be appreciated. Thank you very much Sir.
Sathish konduru.
West Virginia university
Re: Interface modelling in ABAQUS
Hi Sathish,
You might know already about them but as starting point, basic ingredients of the cracking behaviour would be:
- construction of an optimized finite element model with proper assumptions of loading, boundary, symmetry, etc conditions and size, type of elements,...
- selection of approperiate material model for the material, and identify or obtain the corresponding material parameters (most likely from experiments or literature search).
- There are different material models already available in ABAQUS based on concrete cracking, damage plasticity, etc. Please, have a look in the documentation and see if you find anyone of them useful. Although, I am not sure if they are based on thermal loading or not. Otherwise there is always a possibility to write a user material subroutine.
- In most of the cases you must know the crack path apriori to perform crack propagation analyses in ABAQUS. Do you already know?
- Another tedious way would be to find the potential direction of crack propagation and separate nodes in that direction but it is a tedious job to do. You might need many steps/increments etc.
- There are many criterion present, based on your explaination a criterion based on shear stresses might serve the purpose but not sure until further details of the experimental results are known.
I hope this helps. Please let me know if you require further clarifications.
All the best. ENJOY THE BUMPY RIDE.:) and never lose hope.
Amir
Crack propagation
Thank You very much for all replies.
I just wasn't sure whether Abaqus can handle with remeshing while predicting the crack propagation, but now I will try to look at the problem using "Dynamic element dylatation" suggested in first somment. I will study the documentation advised by Zhenyu Xue and I will try to get more familiar with the problem.
Thank You one more time
Best regards
Michael Borowicz
How I can get XFEM-LSM?
Dear All:
I am very interested to hear about the possibility of embedding X-FEM in ABAQUS. I agreed with Jinxiong who has mentioned about this possibility. It is a very appealing approach to simulate crack propagation. How to do that? Any ideas? Or maybe someone already did.
Right now, I am using ABAQUS to get the stress intensity factors (KI, KII, and KIII). To get crack growth, I am using a simplified numerical procedure starting with the Paris's law. I would like to use X-FEM to continue the work from ABAQUS.
I was wondering how X-FEM can be used from here to simulate crack propagation. Also, could you please tell me how I can get a version of XFEM-LSM?
Any advice will be greatly appreciated. Thanks.
Regards,
David.
Open XFEM++
Hi David,
Your research is very interesting. The key idea of XFEM is an enrichment of additional DOFs which may include some information of problems, e.g., the singularity of crack tip, the topology description of the level set function. Enrichment of XFEM for static problems is possible and actually some researchers have succeeded. My previous comment was wondering the difficulties involved in dynamic problems. Maybe this also has been realized by someone in ABAQUS. Can anyone tell me some related works?
There is an open XFEM++ website by Prof. Stephane Bordas. Hopefully it can help you!
Jinxiong
X-FEM implementation in ABAQUS
Dear all,
I'm also of the opinion that ABAQUS capabilities are very powerful and are worth using. We're modelling 2D crack propagation in fretting fatigue problems (i.e. under contact effects such as those due to square ended or cylindrical indenters). We use X-FEM embedded in ABAQUS via UEL, enabling the introduction of extra dof (this must be something very similar to what Xiujun Fang mentions above). We pre-process and post-process the information using Matlab as an interface. In this way, we benefit from the contact procedures and non-linear solver capabilities of ABAQUS in combination with the great advantages of X-FEM. We presented this work at a computational mechanics conference in Lisbon last year. In a similar way, we have also used the more general approach of partition of unity to enrich the singularity that arises at the end of a squared ended contact zone under a sliding condition (also combining X-FEM+ABAQUS). This was presented at WCCM in LA last summer. Some journal papers are now under review. This is joint-work with N. Sukumar at UC Davis.
Since we have realized that the ABAQUS implementation of X-FEM may be of interest to many researchers, we are preparing a paper focusing especially on the implementation aspects and LEFM applications. Once we "clean and tidy up" the code, we'll be glad to make it available as an open-source code.
Regards,
Eugenio.
Some papers about VCCT and DCZM using ABAUQS/UEL
You might be interested and thanks.
1. Xie D and Biggers, Jr. SB, Progressive crack growth analysis using interface element based on the virtual crack closure technique, Finite Elements in Analysis and Design, 42(2006): 977-984. (Science Direct Top 25 Hottest Articles, #4, Apr-Jun, 2006, #11, Jul-Sep, 2006)
2. Xie D and Waas AM, Discrete cohesive zone model for mixed-mode fracture using finite element analysis, Engineering Fracture Mechanics, 73(2006): 1783-1796. (Science Direct Top 25 Hottest Articles, #2, Apr-Jun, 2006, #8, Jul-Sep, 2006)
3. Xie D, Waas AM, Shahwan KW, Schroeder JA and Boeman RG, Computation of strain energy release rate for kinking cracks based on virtual crack closure technique, CMES: Computer Modeling in Engineering & Sciences, 6(2004): 515-524.
Some other papers with application of VCCT and DCZM are
1. Xie D and Biggers, Jr. SB, Delamination growth and residual strength of compressively loaded sandwich panels with stiffness tailored face sheets, Journal of Sandwich Structures and Materials, in press
2. Qian Q and Xie D, Analysis of mixed-mode dynamic crack propagation by interface element based on virtual crack closure technique, Engineering Fracture Mechanics, 74(2007): 807-814.
3. Xie D and Biggers, Jr. SB, Calculation of transient strain energy release rate under impact loading based on virtual crack closure technique, International Journal of Impact Engineering, 34(2007): 1047-1060.
4. Xie D, Salvi AG, Sun C, Waas AM, Caliskan A, Discrete cohesive zone model to study static fracture in a plate of textile braided carbon fiber composites, Journal of Composite Materials, 40(2006): 2025-2046.
5. Xie D and Biggers, Jr. SB, Strain energy release rate calculation for a moving delamination front of arbitrary shape based on virtual crack closure technique, Part I: Formulation and Validation, Engineering Fracture Mechanics, 73(2006): 771-785. (Science Direct Top 25 Hottest Articles, #6, Jan-Mar, 2006).
6. Xie D and Biggers, Jr. SB, Strain energy release rate calculation for a moving delamination front of arbitrary shape based on virtual crack closure technique, Part II: Sensitivity Study on Modeling Details, Engineering Fracture Mechanics, 73(2006): 786-801. (Science Direct Top 25 Hottest Articles, #20, Jan-Mar, 2006)
7. Xie D, Chung J, Waas AM, Shahwan KW, Schroeder JA, Boeman RG, Kunc V, and Klett LB, Failure analysis of adhesively bonded structures: from coupon level data to structure level predictions and verification, International Journal of Fracture, 134(2005): 231-250.
8. Xie D, Waas AM, Shahwan KW Schroeder JA, and Boeman RG, Fracture criteria for kinking cracks in triple material bonded joints, Engineering Fracture Mechanics, 72(2005): 2487-2504. (Science Direct Top 25 Hottest Articles, #3, Jul-Sep, 2005, #13, Oct-Dec, 2005)
Element Splitting
Hello Dear all.
I am working on my project , encountering with one problem about how i can use element spilliting methode with my desired criteria in abaqus whether using abaqus toolbars or with introducing subroutine
is any one can help me?
thanck you in advance
what kind of problem you
what kind of problem you are going to solve? Do you mean element split or node split?
information regarding fracture mechaincs
i am particularly new to fracture mechanics. and new in solving ansys and abaqus software does any body have the idea how the software like ansys and abaqus calculates stress intensity factors and j-integral
Regards
Zoeb Lakdawala
Calculation of SIFs with ANSYS
Hi Zoeb,
I have some experience of calculation of stress intensity factors by using ANSYS software. Can you describe what type of problem do you want to solve? And also if you go to verification problems in ANSYS help, in verification problem 143, there is an example input file for you to calculate stress intensity factors in a plate having a crack at the center of the plate.
Regards,
Erkan.
description of the problem
a shell having an internal pressure of 10.54 bar having a crack of elliptical shape.
Simulation Residual stress With ABAQUS
Iwant to simulate residual stress in inner portion of a pipe in ABAQUS Please HELP ME!!!
Open source code for Embedded Strong Discontinuity (EFEM)
Hi all,
I am now learning strain localization in solids and geometerial by using Embedded Strong Discontinuity theory (EFEM or elemental enrichment). I am looking for an open source code to learn how to understand and how to reproduce some examples by using EFEM. Did anyone know that?
Any advice will be appreciated.
Cheers.
crack propagation using Zencrack with abaqus
hi
any one have idea about how to use zencrack for cracking ,, how i get this software ( free download link) , please let me know
i have to simulate the crack progation in annular disk for my project work
tnks
Zencrack for evaluation
Hello Prashant:
For details, you can use this link:
http://www.zentech.co.uk/zencrack_download_evaluation.htm
Cheers.
hi David
thanks for the link , but i was wondering, that we can use only up 40 element or some thing , with this limitation i am not able complete my project , as my model will be having some 1 Lakhs element.
so could u tell me how can i can solve my problem.
thanks & regards
Metal Foam Failure critrion
Dear All
I'm working on studing the use of metal foam in sandwich structre containing Al or steel plate to overcome an impact shells or plates at differnet impact velocity but the prblem that I face how to define the failure critrion of Meat foam except difine the crushable foam properties, and I've read a paper on using cohesive element but i don'nt know how to use it.
Hopefully some one can help me solving this prblem, by the way I'm using ABAQUS v6.6
Regards
plastic zones arround the crack tip
Hi everyone,
I have a lot of doubts:
1) i wonder if it is posible to find the plastic zones trough abaqus arround the crack tipi.
2) What do i need exactly to define a crack in a viscoelastic step?
3) When i ask the hystory output for the J-integral how can i interpretate the graphics and how can i know if its working because in the model, in visualization nothing is happening is almost the same figure that in the asembly.
Iknow is alot but is really important, so thank you in advance for helping me.
Violeta
let's see the plastic zone arround the tip
Hi,
To your question 1):
You can define the field output PEEQ (Equivalent plastic strian at integration points) , then you can check the plastic zone afterwards.
Xiaobo
element and node splitting
hi guys,
i am going to simulate material failure by using Johnson-cook pasticity and damage model, but i need to split my elements when the material fail, so i have some questions:
1- how can i split element ( or nodes) in my model to show the effect of material separation
2- is it going to work with ALE formulation implemented in Abaqus?
i did some simulations but when it starts for plasticity calculation it exit with ERROR CODE 5 or 114? do you know whats that?
Regards
Hassan
Sheffield
crack propagation
hello
I read the above intersting discussion. I have a doubt, whether I can give my own fracture criterion in ABAQUS. I am trying to simulate crack propagation using ANSYS right now I am using Max normal stress criterion now I want change fracture criteron. I want to include fracture toughness along with max normal stress can I able to implement in ANSYS or is there any provision in ABAQUA
Regards
Gouse
Elliptical Surface crack Modeling
Dear All,
can anybody help me to Model an Elliptical surface crack.Please give me an instruction manual kind of thing . Otherwise mention where to start whether in ABAQUS/CAE or in ABAQUS/Explicit.
Thanks in advance. please respond.........
Regards,
MVN
Modeling of Crack in thru-hole of a plate/Lamb wave propagation
Hi All,
I'm presently working on modelling of guided ultrasonic wave scattering at defect using Abaqus/CAE. Can anyone help me with some tutorial materials on modelling of crack or notch in a through-hole of a plate as well as loading the model with Lamb Wave?
Warm Rgrds.
Jo Ukpai
JIU
Modelling Elliptical Surface Crack(Quarter Plate model)
Hai Zhenyu Xue,
Hope you are having rich experience with ABAQUS. Can you please help me in " Modelling a Surface Elliptical Crack ( Quarter Plate)". If you have any kind of notes on the same please forward it. Otherwise atleast tell some source from where I can get the help in this regard.
Thanking you
MVN