CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions > OpenFOAM CC Toolkits for Fluid-Structure Interaction

Defining pre-existing Fracture as function of time and Position

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 13, 2022, 19:36
Default Defining pre-existing Fracture as function of time and Position
  #1
New Member
 
Axel Dorian
Join Date: Nov 2018
Posts: 12
Rep Power: 6
Axello is on a distinguished road
Using the basics model of Solid mechanics or Solids4Foam on FE4.0/4.1 for Solid Tutorials, on fractures, the pre-exisiting crack is predefined :

in the fvSolution file, we have :

crackLimitingBoxes
1(
(0 -1e-5 0)(0.06 1e-5 1) // -1e-5 ----- 1e5
);

and in the case file 0/U (D) for displacement, the crack is defined as :

crack
{
type solidCohesive;
relaxationFactor 0.01;
contact no;
penaltyScale 1;
frictionCoeff 0.1;
explicitSeparationDistance no;
value uniform ( 0 0 0 );
}


This assumes the crack height and width are fixed.

1. Please, how can I do to set the crack height/width as a function of time and position?

We have that G = - dΩ/ds , where :
Ω is the strain energy only affected by the change in fracture surface area: dΩ = (dΩ/ds)dΩ
G is the surface Energy release
s : crack surface area

in the file .....\src\solidModels\arbitraryCrack\solidCohesive \solidCohesiveFvPatchVectorField.C, the formula :

if ( ((currentGI_[i]/GIc[i]) + (currentGII_[i]/GIIc[i])) >= 1 )
{
//Pout << "GIc[i] is " << GIc[i] << ", curG is " << currentG << endl;
if (!cracked_[i])
{
Pout << "Face " << i << " is fully cracked" << endl;
}

cracked_[i] = true;

gives the condition criteria for propagation.

In the syntax below
if (dict.found("oldGI"))
{
oldGI_ = scalarField("oldGI", dict, p.size());
}
a declaration and definition of the variable oldGI ?

Or

In FE4.1/tut/solidMechanics/elasticAcpSolidFoam/crackingBiMatDcbDugdale/0 is it proper to define oldGI as follow:

crack
{
type solidCohesive;
relaxationFactor 0.01;
contact no;
penaltyScale 1;
frictionCoeff 0.1;
explicitSeparationDistance no;
oldGI 220; for instance
value uniform ( 0 0 0 );
}


2. Regarding the Surface Energy release, how can one record or extract the data file of Energy release variation w.r.t. crack surface area, to time, and Temperature as well in thermal stress is included?

Thank you.
Axello is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Defining pre-existing Fracture as function of time and Position Axello OpenFOAM CC Toolkits for Fluid-Structure Interaction 6 October 18, 2022 08:39
Defining pre-existing Fracture as function of time and Position Axello OpenFOAM Running, Solving & CFD 0 August 13, 2022 19:35
[General] Extracting ParaView Data into Python Arrays Jeffzda ParaView 19 June 11, 2020 16:54
calculation diverge after continue to run zhajingjing OpenFOAM 0 April 28, 2010 05:35
DPM UDF particle position using the macro P_POS(p)[i] dm2747 FLUENT 0 April 17, 2009 02:29


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