CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > ANSYS

crack growth modelling implementing Cohesive Zone Model (CZM)

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By imnull

Reply
 
LinkBack Thread Tools Display Modes
Old   May 6, 2011, 07:24
Question crack growth modelling implementing Cohesive Zone Model (CZM)
  #1
New Member
 
Denis
Join Date: May 2011
Posts: 17
Rep Power: 6
_Denis_ is on a distinguished road
Hello, everyone!
Could you, please, help me with the following issue?
Iím trying to simulate a crack growth using cohesive zone model (CZM), andthe results I get do not agree with my expectations, based on a common senseand simple estimations.



In short, the problem is: the cracking starts only after reaching the normal stress that is million times greater than the maximum normal contact stress parameter.



As I donít yet feel myself on a safe ground with APDL, I make use of AnsysMechanical, implementing command objects where needed.
At the moment the model is looking like a two beams (meshed with 20-nodebricks) bonded together with a surface-to-surface contact, beingseparated from each other.
In the Contacts branch I create a bonded Solid to Solid contact, and attacha command object to it. Here is the text from the command object:
ESEL,S,REAL,,4
EDELE,all
ESEL,ALL
tb,czm,3,,,cbde
tbdata,1,1e4,5e-8,,,1e-8

About the first 3 lines: if I get it right, because of the symmetry of mymodel Mechanical automatically creates two contact sets, with real constant ID3 and 4, respectively. Thus, I believe, it makes sense to delete one of thesesets.
The elastic properties of the material in use are: linear anisotropicelasticity with 11, 22, 33 components of the elasticity matrix around 1e11 Pa.All faces of the beams are perpendicular to the relevant coordinate axis, theseparation traction is applied along one of coordinate axis, normal to thecontact plane (to be precise, I impose on two faces of the different beams notthe force, but ramped displacement).
You can see that the maximum normal contact stress is set to be 1e4 Pa,which is not much, considering the elastic properties of the material. Iestimate that if I set this value to 1e10 Pa, the model should crack afterdeformation of less than 15 %. On the practice, if this value is greater than1e4, the body deforms until the deformation causes instabilities withoutcracking, while the normal stress result branch clearly shows that the actualstress at the interface is much greater than the current normal contact stressparameter value. How could these be? Where is my mistake?
Some additional information about my model:
The pinball radius is set to be program controlled, large deflection = 0n,Augmented Lagrangian method , Physics type: structural, Analysis type:transient.
Thank you very much for your attention!
_Denis_ is offline   Reply With Quote

Old   April 20, 2012, 15:51
Default
  #2
Member
 
Dmitry
Join Date: Nov 2011
Posts: 45
Rep Power: 5
imnull is on a distinguished road
Max_Stress=50 !MPa
Max_Udl=0.05 !/100%
t=0.1 !mm

Gap=Max_Udl*t

tb,czm,cid,,,cbdd
tbdata,1,Max_Stress,Gap,Max_Stress,Gap,0.001,1

Stiffness=-Max_Stress/(0.10*Gap)

rmodif,cid,3,Stiffness
rmodif,cid,12,Stiffness




PLUS


The pinball radius should be set less than your mesh size!!
and dot forget PURE PENALTY formulation in Contact Region
_Denis_ likes this.
imnull is offline   Reply With Quote

Old   April 20, 2012, 16:04
Default
  #3
New Member
 
Denis
Join Date: May 2011
Posts: 17
Rep Power: 6
_Denis_ is on a distinguished road
Thank you very much for answering.
The only problem is that my post is 1 year old, and I've alread got over that problem.
_Denis_ is offline   Reply With Quote

Old   April 20, 2012, 18:25
Default
  #4
Member
 
Dmitry
Join Date: Nov 2011
Posts: 45
Rep Power: 5
imnull is on a distinguished road
Quote:
Originally Posted by _Denis_ View Post
Thank you very much for answering.
The only problem is that my post is 1 year old, and I've alread got over that problem.
have u done a model validation with real experiment?
imnull is offline   Reply With Quote

Old   April 22, 2012, 03:24
Default
  #5
New Member
 
Denis
Join Date: May 2011
Posts: 17
Rep Power: 6
_Denis_ is on a distinguished road
I haven't conducted any experiments myself, but I've compared the speed of crack propagation resulting from my simulation with known experimental results.
_Denis_ is offline   Reply With Quote

Old   July 23, 2013, 10:41
Default cohesive mateial
  #6
New Member
 
ali
Join Date: Jul 2013
Posts: 1
Rep Power: 0
alizadeh is on a distinguished road
hi dear Denis
i m trying to model delamination via contact element and cohesive material which involve defining cohesive material between two planes meshed by shell99
would you let me know how can i define cohesive material between two planes

all the best
alizadeh
alizadeh is offline   Reply With Quote

Old   April 14, 2015, 08:41
Default
  #7
New Member
 
NithyaRamesh
Join Date: Apr 2015
Posts: 2
Rep Power: 0
NithyaRamesh is on a distinguished road
Hi
i have also one doubt regarding cohesive zone modeling. That, is it possible to give 1000mm shear separation for 100mm Double cantilever beam?
NithyaRamesh is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Implementing a new LES Model in OpenFoam fs82 OpenFOAM 6 October 13, 2009 09:58
Implementing new turbulence model sven OpenFOAM 13 October 12, 2009 13:18
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Implementing and running a new turbulence model svens OpenFOAM 3 August 21, 2009 02:59


All times are GMT -4. The time now is 04:45.