... | @@ -7,7 +7,7 @@ damageinfo.dat is read by DamageInfoParser.h |
... | @@ -7,7 +7,7 @@ damageinfo.dat is read by DamageInfoParser.h |
|
## initial control:
|
|
## initial control:
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| initial_defect_width = 0.2 | here, the initial defect will be **0.2** times the characteristic length `lc` of the material. |
|
|
| initial_defect_width = 0.2 | here, the initial defect will be **0.2** times the characteristic length `lc` of the material. |
|
|
| restart_file = ls_control_timestep_0730.bin | (no warranty on restart procedure) |
|
|
| restart_file = ls_control_timestep_0730.bin | (no warranty on restart procedure) |
|
|
|
|
|
... | @@ -19,7 +19,7 @@ damageinfo.dat is read by DamageInfoParser.h |
... | @@ -19,7 +19,7 @@ damageinfo.dat is read by DamageInfoParser.h |
|
## job control :
|
|
## job control :
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| main_debug = 0 | **1** for verbose prints, else **0** |
|
|
| main_debug = 0 | **1** for verbose prints, else **0** |
|
|
| max_number_iterations = 200 | maximum number of iteration for crack propagation.|
|
|
| max_number_iterations = 200 | maximum number of iteration for crack propagation.|
|
|
| maximum_elastic_nl_iteration = 20 | maximum number of iteration for Newton-Raphson solver. |
|
|
| maximum_elastic_nl_iteration = 20 | maximum number of iteration for Newton-Raphson solver. |
|
... | @@ -43,7 +43,7 @@ damageinfo.dat is read by DamageInfoParser.h |
... | @@ -43,7 +43,7 @@ damageinfo.dat is read by DamageInfoParser.h |
|
## physical caracteristique:
|
|
## physical caracteristique:
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| lc = 0.01 | Characteristic length `lc` of the material |
|
|
| lc = 0.01 | Characteristic length `lc` of the material |
|
|
| Yc = 100 | Characteristic Yield `Yc` of the material |
|
|
| Yc = 100 | Characteristic Yield `Yc` of the material |
|
|
| YPCoef = 0. | Generally **0**, except for Cohesive Zone Model equivalence where it is ???|
|
|
| YPCoef = 0. | Generally **0**, except for Cohesive Zone Model equivalence where it is ???|
|
... | @@ -54,10 +54,10 @@ See article **Fig. 3** in **Comparison between thick level set (TLS) and cohesiv |
... | @@ -54,10 +54,10 @@ See article **Fig. 3** in **Comparison between thick level set (TLS) and cohesiv |
|
Nicolas Moës and Claude Stolz, **Advanced Modeling and Simulation in Engineering Sciences**, 2015, pp 2:18,
|
|
Nicolas Moës and Claude Stolz, **Advanced Modeling and Simulation in Engineering Sciences**, 2015, pp 2:18,
|
|
DOI: 10.1186/s40323-015-0041-9
|
|
DOI: 10.1186/s40323-015-0041-9
|
|
|
|
|
|
[[https://https://static-content.springer.com/image/art%3A10.1186%2Fs40323-015-0041-9/MediaObjects/40323_2015_41_Fig3_HTML.gif]]
|
|
[[https://static-content.springer.com/image/art%3A10.1186%2Fs40323-015-0041-9/MediaObjects/40323_2015_41_Fig3_HTML.gif]]
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| hardeningCZM = 1 | **1**: use Cohesive Zone Model equivalence |
|
|
| hardeningCZM = 1 | **1**: use Cohesive Zone Model equivalence |
|
|
| CZM_wf = 7.80e-5 | *wf* maximun opening value (stress is 0). |
|
|
| CZM_wf = 7.80e-5 | *wf* maximun opening value (stress is 0). |
|
|
| CZM_w1 = 3.00e-5 | *w1* opening value. |
|
|
| CZM_w1 = 3.00e-5 | *w1* opening value. |
|
... | @@ -70,7 +70,7 @@ DOI: 10.1186/s40323-015-0041-9 |
... | @@ -70,7 +70,7 @@ DOI: 10.1186/s40323-015-0041-9 |
|
## numerical methode control:
|
|
## numerical methode control:
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| polynomial_order = 2 | **2** (but, level set order is 1) |
|
|
| polynomial_order = 2 | **2** (but, level set order is 1) |
|
|
| space_type = 1 | ... |
|
|
| space_type = 1 | ... |
|
|
| use_enrichment = 5 | **5** to use Heavyside function enrichment of the crack, else **3** |
|
|
| use_enrichment = 5 | **5** to use Heavyside function enrichment of the crack, else **3** |
|
... | @@ -83,7 +83,7 @@ DOI: 10.1186/s40323-015-0041-9 |
... | @@ -83,7 +83,7 @@ DOI: 10.1186/s40323-015-0041-9 |
|
## geometric representation of the tls
|
|
## geometric representation of the tls
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| tls_geom_interface = 2 | **2** for Double Cut to use with `use_enrichment` equal to **5** or **0** , **0** for Simple Cut to use with `use_enrichment` equal to **3** or **0** (**1** for experimental double cut implemented by Nicolas Chevaugeon for 2D problems). |
|
|
| tls_geom_interface = 2 | **2** for Double Cut to use with `use_enrichment` equal to **5** or **0** , **0** for Simple Cut to use with `use_enrichment` equal to **3** or **0** (**1** for experimental double cut implemented by Nicolas Chevaugeon for 2D problems). |
|
|
|
|
|
|
|
|
|
... | @@ -96,7 +96,7 @@ DOI: 10.1186/s40323-015-0041-9 |
... | @@ -96,7 +96,7 @@ DOI: 10.1186/s40323-015-0041-9 |
|
damage front (level set) displacement control
|
|
damage front (level set) displacement control
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| front_advance_algo = 3 | generally **3**, except for Cohesive Zone Model equivalence where it is **4**. |
|
|
| front_advance_algo = 3 | generally **3**, except for Cohesive Zone Model equivalence where it is **4**. |
|
|
| k_max = 100. | maximum value of `kmax`. |
|
|
| k_max = 100. | maximum value of `kmax`. |
|
|
| power_law_coef = 100. | initial value of `kmax`, which may evolve during the TLS process but will be limited by `k_max`(above). Small value induices a global expansion of the damage zone (level set displacment), big value induices propagation arround th process zone (crack tip). Must be big compared to `lc`. See P-E **A level set based model for damage growth: The thick level set approach**, N. Moës, C. Stolz, P.-E. Bernard, N. Chevaugeon, **International Journal For Numerical Methods in Engineering**, 86,pp 358–380, 2011 |
|
|
| power_law_coef = 100. | initial value of `kmax`, which may evolve during the TLS process but will be limited by `k_max`(above). Small value induices a global expansion of the damage zone (level set displacment), big value induices propagation arround th process zone (crack tip). Must be big compared to `lc`. See P-E **A level set based model for damage growth: The thick level set approach**, N. Moës, C. Stolz, P.-E. Bernard, N. Chevaugeon, **International Journal For Numerical Methods in Engineering**, 86,pp 358–380, 2011 |
|
... | @@ -109,7 +109,7 @@ damage front (level set) displacement control |
... | @@ -109,7 +109,7 @@ damage front (level set) displacement control |
|
## damage law control:
|
|
## damage law control:
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| damage_shape = Poly3 | the shape of the function Phi(d). Generally `Poly3 ` except for Cohesive Zone Model equivalence where it is `Poly2revert`.|
|
|
| damage_shape = Poly3 | the shape of the function Phi(d). Generally `Poly3 ` except for Cohesive Zone Model equivalence where it is `Poly2revert`.|
|
|
| hcoef = 0. | coefficient before `alpha_i` variables of the asymmetric traction/compression behavior of the material : **0** means that asymmetric is taken into account, **1** means that the material is symmetric (standard elasticity). Generally **0** or **1** but can be fixed numerically in [0,1]. |
|
|
| hcoef = 0. | coefficient before `alpha_i` variables of the asymmetric traction/compression behavior of the material : **0** means that asymmetric is taken into account, **1** means that the material is symmetric (standard elasticity). Generally **0** or **1** but can be fixed numerically in [0,1]. |
|
|
| penalty_coef = 1. | involved in the computation of the field averaging (see Eq.9.a in the paper of A. Salzman 2016). |
|
|
| penalty_coef = 1. | involved in the computation of the field averaging (see Eq.9.a in the paper of A. Salzman 2016). |
|
... | @@ -124,7 +124,7 @@ damage front (level set) displacement control |
... | @@ -124,7 +124,7 @@ damage front (level set) displacement control |
|
Only used withDamageTestContext <DamageTestOctree> template parameter.
|
|
Only used withDamageTestContext <DamageTestOctree> template parameter.
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| band_width_octree = 1.5 | ratio of `lc`. Distance arround the damage level set where the mesh is fine. |
|
|
| band_width_octree = 1.5 | ratio of `lc`. Distance arround the damage level set where the mesh is fine. |
|
|
| band_width_geo = 1. | default value is **1** |
|
|
| band_width_geo = 1. | default value is **1** |
|
|
| band_width_fr = 0.25 | ratio of `lc`. Define mesh refinement where boundary conditions are applied. |
|
|
| band_width_fr = 0.25 | ratio of `lc`. Define mesh refinement where boundary conditions are applied. |
|
... | @@ -136,7 +136,7 @@ Only used withDamageTestContext <DamageTestOctree> template parameter. |
... | @@ -136,7 +136,7 @@ Only used withDamageTestContext <DamageTestOctree> template parameter. |
|
## search for initiation control:
|
|
## search for initiation control:
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| search_for_initiations = -10000 | when smaller than **0**, initation zone is automatically computed at the beginning of the job. When greater than **0** the initiation zone is evaluated every `search_for_initiations` iterations. |
|
|
| search_for_initiations = -10000 | when smaller than **0**, initation zone is automatically computed at the beginning of the job. When greater than **0** the initiation zone is evaluated every `search_for_initiations` iterations. |
|
|
| start_with_search_for_initiations = 1 | **1**: start the job with a virgin part and `search_for_initiations` must be less than **0**. **0**: the initiation zone is defined by the level set given in the main file. |
|
|
| start_with_search_for_initiations = 1 | **1**: start the job with a virgin part and `search_for_initiations` must be less than **0**. **0**: the initiation zone is defined by the level set given in the main file. |
|
|
| sphere_radius_search_for_initiations = 0.03 | ratio of `lc` for level set initiation. |
|
|
| sphere_radius_search_for_initiations = 0.03 | ratio of `lc` for level set initiation. |
|
... | @@ -148,7 +148,7 @@ Only used withDamageTestContext <DamageTestOctree> template parameter. |
... | @@ -148,7 +148,7 @@ Only used withDamageTestContext <DamageTestOctree> template parameter. |
|
## export control:
|
|
## export control:
|
|
|
|
|
|
| **command**| **description**|
|
|
| **command**| **description**|
|
|
|:---------------------------------|:-------------------------------|
|
|
|:------------------------------------------------|:-------------------------------|
|
|
| full_elastic_out = 0 | **0**: export mesh without the full-damaged zone (**d=1**). **1**: the zone where **d=1** is meshed for output. |
|
|
| full_elastic_out = 0 | **0**: export mesh without the full-damaged zone (**d=1**). **1**: the zone where **d=1** is meshed for output. |
|
|
| load_factor_out = 1 | **1**: fields are ouput by taken into account the load factor. **0** fields are output for the load given in the **main.dat** file |
|
|
| load_factor_out = 1 | **1**: fields are ouput by taken into account the load factor. **0** fields are output for the load given in the **main.dat** file |
|
|
| output_name = output | ... |
|
|
| output_name = output | ... |
|
... | | ... | |