regular tables removed, loading scenarios adapted

7b5457c3 · Rostislav Chudoba · 71b67507 · 7b5457c3 · 7b5457c3 · 7b5457c3
Commit 7b5457c3 authored 2 years ago by Rostislav Chudoba
--- a/environment.yml
+++ b/environment.yml
@@ -10,7 +10,7 @@ dependencies:
 - k3d # ==2.12.0
 - ipywidgets # ==7.6.5
 - ipympl # ==0.8.8
- ipyregulartable # ==0.2.0
+# - ipyregulartable # ==0.2.0
 - ipytree # ==0.2.1
 #- pyface==7.3.0
 #- traitsui==7.1.0

--- a/tour4_plastic_bond/4_1_PO_multilinear_unloading.ipynb
+++ b/tour4_plastic_bond/4_1_PO_multilinear_unloading.ipynb
--- a/tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb
+++ b/tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb
--- a/tour4_plastic_bond/4_3_PO_trc_cfrp_cyclic.ipynb
+++ b/tour4_plastic_bond/4_3_PO_trc_cfrp_cyclic.ipynb
@@ -134,7 +134,7 @@
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "b685da8c49be4f8fbc15b9d36f9edabb",
+       "model_id": "0d185daa73fa4e6f8c227f7406fddfda",
       "version_major": 2,
       "version_minor": 0
      },
@@ -166,13 +166,20 @@
  },
  {
   "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 4,
   "metadata": {},
   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "xxxxx\n"
+     ]
+    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "c0500e76cc7748b1b1953833102698e9",
+       "model_id": "34857243edc54fe680b6fc142ea3effc",
       "version_major": 2,
       "version_minor": 0
      },
@@ -185,13 +192,12 @@
    }
   ],
   "source": [
-    "po_trc.loading_scenario = 'cyclic'\n",
+    "po_trc.loading_scenario.profile = 'cyclic-nonsym-const'\n",
-    "po_trc.loading_scenario_.trait_set(number_of_cycles=2,\n",
+    "po_trc.loading_scenario.profile_.trait_set(number_of_cycles=2,\n",
    "                              unloading_ratio=0.8,\n",
    "                            number_of_increments=200,\n",
-    "                              amplitude_type='constant',\n",
+    "                              );\n",
-    "                              loading_range='non-symmetric');\n",
+    "po_trc.loading_scenario.interact()"
-    "po_trc.loading_scenario_.interact()"
   ]
  },
  {
@@ -205,13 +211,20 @@
  },
  {
   "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 5,
   "metadata": {},
   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Warning: Update progress bar from model - not implemented yet\n"
+     ]
+    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "10fe2e84cd274cc5baa7b3d69b757625",
+       "model_id": "c5a1d36bcd8c45c78566e2c4e0c5abda",
       "version_major": 2,
       "version_minor": 0
      },
@@ -255,7 +268,7 @@
  },
  {
   "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
@@ -295,13 +308,13 @@
  },
  {
   "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "c99b02c3148a4abb94cc07a68078ad18",
+       "model_id": "bac2295e25b244d4adf0ac01cea34083",
       "version_major": 2,
       "version_minor": 0
      },
@@ -361,13 +374,20 @@
  },
  {
   "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 16,
   "metadata": {},
   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Warning: Update progress bar from model - not implemented yet\n"
+     ]
+    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "e4fbee9827424bf7bbe9dc2d6da695a6",
+       "model_id": "717e68a2e6774d2f96acdec9dbc5d997",
       "version_major": 2,
       "version_minor": 0
      },
@@ -380,12 +400,19 @@
    }
   ],
   "source": [
-    "pm.loading_scenario = 'cyclic'\n",
+    "pm.loading_scenario.profile = 'cyclic-sym-incr'\n",
-    "pm.loading_scenario_.trait_set(number_of_cycles=4,\n",
+    "pm.loading_scenario.profile_.trait_set(number_of_cycles=2,\n",
    "                              unloading_ratio=0.5,\n",
-    "                               number_of_increments=200,\n",
+    "                            number_of_increments=200,\n",
-    "                              amplitude_type='increasing',\n",
+    "                              );\n",
-    "                              loading_range='symmetric')\n",
+    "\n",
+    "# pm.loading_scenario = 'cyclic'\n",
+    "# pm.loading_scenario_.trait_set(number_of_cycles=4,\n",
+    "#                               unloading_ratio=0.5,\n",
+    "#                                number_of_increments=200,\n",
+    "#                               amplitude_type='increasing',\n",
+    "#                               loading_range='symmetric')\n",
+    "pm.reset()\n",
    "pm.run()\n",
    "pm.interact()"
   ]
@@ -396,7 +423,7 @@
   "source": [
    "## Discussion of the study\n",
    "\n",
-    " - By checking the `history` showing the field variables along the bond zone, we can observe the the alternating tension and compression profiles in the reinforcement and in the matrix.\n",
+    " - By checking the `history` showing the field variables along the bond zone, we can observe the alternating tension and compression profiles in the reinforcement and in the matrix.\n",
    " - The shear transfer is again taking place only on a short length due to the softening nature of the bond behavior.\n",
    " - The yielding of the bond within the process zone occurs both under pullout and under push-in load as apparent from the existence of the hysteretic loop.\n",
    " - In contrast to the [study of TRC test](#trc_results), the unloading stiffness does not correspond to the initial stiffness. This question is discussed in a more detail below.  "
@@ -461,7 +488,7 @@
 ],
 "metadata": {
  "kernelspec": {
-   "display_name": "Python 3",
+   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
@@ -475,7 +502,7 @@
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
-   "version": "3.8.10"
+   "version": "3.9.0"
  },
  "toc": {
   "base_numbering": 1,

 %% Cell type:markdown id: tags:
 <a id="top"></a>
 %% Cell type:markdown id: tags:
 # **4.3: Cyclic pullout of textile fabrics and CFRP sheets**
 %% Cell type:markdown id: tags:
 [![title](../fig/bmcs_video.png)](https://moodle.rwth-aachen.de/mod/page/view.php?id=619022)
 %% Cell type:markdown id: tags:
 <div style="background-color:lightgray;text-align:left"> <img src="../icons/start_flag.png" alt="Previous trip" width="40" height="40">
    &nbsp; &nbsp; <b>Starting point</b> </div>
 %% Cell type:markdown id: tags:
 Once we have seen that there is a relation between the shape of the constitutive law and the  stress redistribution process within a bond zone which has an immediate consequence on the shape of the pullout curve, we extend our horizon to the case of **non-monotonic loading**.
 %% Cell type:markdown id: tags:
 <div style="background-color:lightgray;text-align:left"> <img src="../icons/destination.png" alt="Previous trip" width="40" height="40">
    &nbsp; &nbsp; <b>Where are we heading</b> </div>
 %% Cell type:markdown id: tags:
 To motivate a more general development of the model for the bond-slip behavior let us consider the case of non-monotonic loading. What happens in the material structure of the bond if the load is reduced and than it grows again?
 %% Cell type:markdown id: tags:
 # **TRC pullout test modeled using plastic material model**
 %% Cell type:markdown id: tags:
 ## Test results
 %% Cell type:markdown id: tags:
 <a id="trc_test_results"></a>
 Consider again the test double-sided pullout test that was used to motivate the development of plastic material model in [notebook 4.1](4_1_PO_multilinear_unloading.ipynb#trc_study_monotonic) which demonstrated that no unloading branch could be reproduced. We reprint the test results for convenience here with the goal to examine the ability of the plastic model to reproduce the unloading stiffness of the specimen with the total length of 300 mm, i.e. $L_\mathrm{b} = 150$ mm.
 ![image](../fig/test_unloading.png)
 %% Cell type:markdown id: tags:
 ## Pullout model setup
 <a id="trc_test"></a>
 Let us first construct the model with the same geometry as specified in [notebook 3.2](../tour3_nonlinear_bond/3_2_anchorage_length.ipynb#trc_parameters). Note again that the maximum control displacement $w_\mathrm{max}$ is
 the half of the crack opening displacement measured in the test. Therefore `w_max` has been set to 3 mm.
 %% Cell type:code id: tags:
 ``` python
 %matplotlib widget
 import matplotlib.pyplot as plt
 import numpy as np
 from bmcs_cross_section.pullout import PullOutModel1D
 po_trc = PullOutModel1D(n_e_x=100, w_max=3) # mm
 po_trc.geometry.L_x=150 # [mm]
 po_trc.time_line.step = 0.02
 po_trc.cross_section.trait_set(A_m=1543, A_f=16.7, P_b=10);
 ```
 %% Cell type:markdown id: tags:
 ## Material model with kinematic hardening
 To apply the elastic-plastic model in this pullout model, set the the `elasto-plasticity` option for the `material_model` attribute. Then, the plastic model with the parameters $\tau_Y = \bar{\tau}, K, \gamma$ introduced in the [notebook 3.2](../tour3_nonlinear_bond/3_2_anchorage_length.ipynb#trc_pullout_study) can be set to the attribute `material_model_` using the `trait_set` method as follows
 | parameter | name | value | unit |
 | - | - | -: | - |
 | bond stiffness | `E_b` | 6.4 | MPa |
 | elastic limit stress | `tau_bar` | 4 | MPa |
 | isotropic hardening modulus | `K` | 0 | MPa |
 | kinematic hardening modulus | `gamma` | 1.0 | MPa |
 %% Cell type:markdown id: tags:
 Using these parameters, we aim to reproduce the same shape of the [bond slip law that was used previously in notebook 3.2](../tour3_nonlinear_bond/3_2_anchorage_length.ipynb#trc_parameters), represented as a tri-linear function for the simulation of the monotonically increasing load.
 %% Cell type:code id: tags:
 ``` python
 po_trc.material_model='elasto-plasticity'
 po_trc.material_model_.trait_set(
    E_m=28000, E_f=170000, E_b=6.4, K = 0, gamma=1.0, tau_bar=4
 );
 po_trc.material_model_.interact()
 ```
 %% Output
 %% Cell type:markdown id: tags:
 ## Load scenario
 Finally, following the same mimic, we change the loading scenario from `monotonic` to `cyclic` and can set the parameters of the cyclic loading through the attribute with the trailing underscore `loading_scenario_`.
 To verify the loading scenario we directly render the model component using the `interact` method.
 %% Cell type:code id: tags:
 ``` python
-po_trc.loading_scenario = 'cyclic'
+po_trc.loading_scenario.profile = 'cyclic-nonsym-const'
-po_trc.loading_scenario_.trait_set(number_of_cycles=2,
+po_trc.loading_scenario.profile_.trait_set(number_of_cycles=2,
                              unloading_ratio=0.8,
                            number_of_increments=200,
-                              amplitude_type='constant',
+                              );
-                              loading_range='non-symmetric');
+po_trc.loading_scenario.interact()
-po_trc.loading_scenario_.interact()
 ```
 %% Output
+    xxxxx
 %% Cell type:markdown id: tags:
 ## Run the simulation
 Apparently, we apply directly the full control displacement, then reduce it to 80% and reload again to the maximum. Our intention is not to simulate the whole history displayed in the [test results](trc_test_results) but to compare the unloading stiffness reproduced by the model.
 The model is now ready for execution so that let us run it and render the user interface to inspect the results.
 %% Cell type:code id: tags:
 ``` python
 po_trc.reset()
 po_trc.run()
 po_trc.interact()
 ```
 %% Output
+    Warning: Update progress bar from model - not implemented yet
 %% Cell type:markdown id: tags:
 ## Discussion of the study
 - The predicted response reproduces the experimental response reasonably well. For the **short bond length** used in the test, the slip and shear profiles are nearly uniform along the bond zone.  - **Task:** Change the bond length to 200 mm and compare the response with the pull-out curve depicted in the [test results](#trc_test_results)
 - The main focus of the study at hand was the question, if the **unloading stiffness** corresponds to the initial stiffness, similarly to the test response. Since this is the case, the conclusion can be made, that the bond behavior in the tested TRC specimens is primarily governed by plasticity.
 - Since the model uses **linear hardening**, the pull-out response has a linear inelastic branch up to the point of unloading. The non-linear shape of the experimentally obtained curve cannot be reproduced in a better way. **Non-linear hardening** function instead of a single hardening modulus $\gamma$ can be theoretically included into the framework of plasticity. However, such enhancement requires a non-trivial systematic calibration procedure which would be able to uniquely associate the inelastic mechanisms governing the response to the hardening process, which is a non-trivial task. Such effort would only be justified for a material combination with large amount of practical applications.
 %% Cell type:markdown id: tags:
 <a id="cfrp_test"></a>
 # **CFRP sheet test with softening plastic material**
 It is interesting to examine the ability of the plastic model with negative hardening, i.e. softening, to simulate the cyclic response of the CFRP sheets.
 ## Pullout model setup
 Let us construct the model with these parameters and check the shape of the bond-slip relation
 %% Cell type:code id: tags:
 ``` python
 %matplotlib widget
 import matplotlib.pyplot as plt
 import numpy as np
 from bmcs_cross_section.pullout import PullOutModel1D
 A_f = 16.67 # [mm^2]
 A_m = 1540.0 # [mm^2]
 p_b = 100.0 #
 E_f = 170000 # [MPa]
 E_m = 28000 # [MPa]
 pm = PullOutModel1D()
 pm.sim.tline.step = 0.01 # 100 time increments
 pm.cross_section.trait_set(A_f=A_f, P_b=p_b, A_m=A_m)
 pm.geometry.L_x = 300 # length of the specimen [mm]
 pm.w_max = 2.8 # maximum control displacement [mm]
 pm.n_e_x = 100 # number of finite elements
 ```
 %% Cell type:markdown id: tags:
 ## Material model with isotropic softening
 The bond-slip law identified during the [**Tour 3:** (CFRP sheet test)](../tour3_nonlinear_bond/3_1_nonlinear_bond.ipynb#cfrp_bond_slip) can be reproduced by setting the parameters of
 the plastic material model with isotropic hardening to:
 | parameter | name | value | unit |
 | - | - | -: | - |
 | bond stiffness | `E_b` | 80 | MPa |
 | elastic limit stress | `tau_bar` | 8 | MPa |
 | isotropic hardening modulus | `K` | -20 | MPa |
 | kinematic hardening modulus | `gamma` | 0 | MPa |
 %% Cell type:code id: tags:
 ``` python
 pm.material_model='elasto-plasticity'
 pm.material_model_.trait_set(
    E_m=28000, E_f=170000, E_b=80, K = -20, gamma=0, tau_bar=8
 )
 pm.material_model_.interact()
 ```
 %% Output
 %% Cell type:markdown id: tags:
 <div style="background-color:lightgray;text-align:left"> <img src="../icons/remark.png" alt="Previous trip" width="40" height="40">
    &nbsp; &nbsp; <b>What is the difference compared to a trilinear bond-slip law used before?</b> </div>
 %% Cell type:markdown id: tags:
 Apparently, the shape is identical to the [trilinear bond slip law](../tour3_nonlinear_bond/3_1_nonlinear_bond.ipynb#cfrp_trilinear_bond) law applied in the first analysis of this test setup. However, now unloading in inelastic regime would be governed by initial stiffness.
 It is important to remark that the model has been augmented with an additional feature that prevents softening into negative stresses. To explain this minor but essential extension recal the shape of the yield condition for plastic model with isotropic hardening
 \begin{align}
 f = | \tau | - (\tau_Y - Z)
 \end{align}
 where the hardening stress is given as $Z = K z$. With the negative value of hardening modulus $K$ and implicitly positive value of the hardening state variable $z$ we can notice that the size of the elastic domain $(\tau_Y - Z)$ shrinks with increasing $z$. Without any treatment, the elastic domain could reach a negative size, which is nonphysical. Thus, the criterion of the positive size of the elastic domain is used to define a third branch in the plastic model.
 **In fact,** with this necessary extension, we introduce **damage into a plastic model**, because we explicitly set the material stiffness at a predefined value of deformation, namely at the value of slip inducing the loss of elastic domain. As we will lear in Tour 5, relating the loss of stiffness to the measure of strain is the general concept behind damage models. :
 In this context, where we introduced an argument of physical admissibility, it is worth reminding that inelastic behavior implicitly induces **dissipation of energy**. A sound formulation of material models is only possible in the framework of thermodynamic laws that can guarantee an admissible behavior for any loading history. The evaluation of energy dissipation and its importance for objectivity and robustness of material models will be illuminated in **Tour 6**.
 With the material model at hand, let us explore the response of the CFRP sheets in a cyclic loading scenario with symmetric, increasing amplitudes
 %% Cell type:markdown id: tags:
 ## Load scenario
 To see the effect of propagating softening plasticity upon unloading, let us apply four loading cycles oscillating applying the pull-out and push-in load in each cycle. Note that with this,
 a compressive stress is transferred by the FRP sheet upon push-in loading, which is unrealistic.
 Such response might be expected from a test on CFRP bars. Still the purpose of the study is
 to show the consequences of the applied assumptions on a broad scale of loading configuration and check the plausibility of the model response for a wider range of loading scenarios.
 %% Cell type:code id: tags:
 ``` python
-pm.loading_scenario = 'cyclic'
+pm.loading_scenario.profile = 'cyclic-sym-incr'
-pm.loading_scenario_.trait_set(number_of_cycles=4,
+pm.loading_scenario.profile_.trait_set(number_of_cycles=2,
                              unloading_ratio=0.5,
-                               number_of_increments=200,
+                            number_of_increments=200,
-                              amplitude_type='increasing',
+                              );
-                              loading_range='symmetric')
+# pm.loading_scenario = 'cyclic'
+# pm.loading_scenario_.trait_set(number_of_cycles=4,
+#                               unloading_ratio=0.5,
+#                                number_of_increments=200,
+#                               amplitude_type='increasing',
+#                               loading_range='symmetric')
+pm.reset()
 pm.run()
 pm.interact()
 ```
 %% Output
+    Warning: Update progress bar from model - not implemented yet
 %% Cell type:markdown id: tags:
 ## Discussion of the study
- - By checking the `history` showing the field variables along the bond zone, we can observe the the alternating tension and compression profiles in the reinforcement and in the matrix.
+ - By checking the `history` showing the field variables along the bond zone, we can observe the alternating tension and compression profiles in the reinforcement and in the matrix.
 - The shear transfer is again taking place only on a short length due to the softening nature of the bond behavior.
 - The yielding of the bond within the process zone occurs both under pullout and under push-in load as apparent from the existence of the hysteretic loop.
 - In contrast to the [study of TRC test](#trc_results), the unloading stiffness does not correspond to the initial stiffness. This question is discussed in a more detail below.
 %% Cell type:markdown id: tags:
 <div style="background-color:lightgray;text-align:left"> <img src="../icons/question_puzzle.png" alt="Question" width="50" height="50">
    &nbsp; &nbsp; <b>Stimulating question</b> </div>
 %% Cell type:markdown id: tags:
 # **Qualitative comparison of the two studies**
 %% Cell type:markdown id: tags:
 ## **How to interpret unloading stiffness at structural level?**
 In the classification of plasticity and damage models explained in  [4.1](4_1_PO_multilinear_unloading.ipynb#damage_and_plasticity) we have learned that in a single material point, plastic behavior unloads with initial stiffness and induces a persistent deformation $s_\mathrm{pl}$ at zero level of load. On the other hand, unloading of material point governed by damage behavior would recover all deformation and reach zero deformation at zero stress. **Can we assume such behavior of both class of models at the structural level?**
 Comparing the slopes of the first unloading branch obtained for the [TRC test](trc_test) and of the [CFRP test](cfrp_test), we can notice that in the former case, the unloading stiffness is indeed equal to the initial stiffness. In this case, governed by the hardening bond-slip law, the shear stress increases and decreases upon loading and unloading in all material points simultaneously. On the other hand, in the case of [CFRP test](cfrp_test) governed by softening, the debonded zone reaches zero stress level and does not contribute to stress transfer any more. As a result, the effective length of the bond zone becomes shorter, so that the structural stiffness upon unloading must be smaller than the initial stiffness.
 **The conclusion** is, that an application of a plastic model with softening,
 is only valid up to the point at which the elastic domain vanished. From that point on,
 a full damage must be assumed. Therefore, the kind of model used in the simulation of the CFRP test can be regarded as a first version of a combined damage-plasticity model.
 %% Cell type:markdown id: tags:
 ## **When to apply isotropic and when kinematic hardening?**
 In the two examples above, both kinematic and isotropic hardening has been used. In the first yielding step and upon the first unloading, there is no difference between the two parameters. They will only show a different behavior upon cyclic loading scenario. A correct distinction is only possible based on a test delivering several unloading and reloading cycles. Correct identification of the parameters is a demanding task which calls for a systematic and automated support in terms of calibration algorithms and validation procedures.
 %% Cell type:markdown id: tags:
 <div style="background-color:lightgray;text-align:left;width:45%;display:inline-table;"> <img src="../icons/previous.png" alt="Previous trip" width="50" height="50">
    &nbsp; <a href="../tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb#top">4.2 Basic concept of plasticity</a>
 </div><div style="background-color:lightgray;text-align:center;width:10%;display:inline-table;"> <a href="#top"><img src="../icons/compass.png" alt="Compass" width="50" height="50"></a></div><div style="background-color:lightgray;text-align:right;width:45%;display:inline-table;">
    <a href="../tour5_damage_bond/5_1_Introspect_Damage_Evolution_Damage_initiation.ipynb#top">5.1 Pullout simulation using damage model</a>&nbsp; <img src="../icons/next.png" alt="Previous trip" width="50" height="50"> </div>
 %% Cell type:code id: tags:
 ``` python
 ```

--- a/tour5_damage_bond/5_1_Introspect_Damage_Evolution_Damage_initiation.ipynb
+++ b/tour5_damage_bond/5_1_Introspect_Damage_Evolution_Damage_initiation.ipynb