Salome HOME
Update of CheckDone
[modules/smesh.git] / src / Tools / blocFissure / doc / principles.rst
1 .. _general_principles:\r
2 \r
3 ###################\r
4 General principles\r
5 ###################\r
6 \r
7 « Bloc Fissure » is based on GEOM module geometrical Boolean operations. The initial structure being a mesh, a conversion from mesh to geometry becomes necessary. This operation is called extraction and reconstruction because it is only applied to a small part of the mesh, which is around the crack. The extracted mesh is called the « Box » and only external faces of this mesh are kept and converted into several geometrical surfaces. This operation implies some limitations on the input mesh. When all the Boolean operations are done, the geometry that contains the crack is meshed again with a ruled mesh in the torus and a free mesh elsewhere.\r
8 \r
9 In order to illustrate « Bloc Fissure » principle, the simple case of a crack insertion in a parallelepipedic specimen is detailed step by step:\r
10 \r
11 1. The first step consists in loading the structure mesh (a) as well as the crack surface geometry (b).\r
12 \r
13 2. The crack is then meshed (c). A length criterion defines the size of the extracted « Box ». This length is called the length of influence. Every element of the structure mesh having a node within this zone is included in the « Box ». A second operation adds those elements in the Box in order to have continuous faces (d).\r
14 \r
15 3. A geometrical Box is reconstructed from the extracted Box mesh. The reconstruction is limited to faces which intersect the crack (e). A torus is created following the crack front (f).\r
16 \r
17 4. The geometrical Box is then cut by the torus and the crack (g). Several plans are created in order to partition the box and the torus into radiuses for the future mesh (h).\r
18 \r
19 5. The Box, the crack and the torus are meshed on their external surface (i) and then filled with volumetric elements (j). Crack nodes are doubled to « open » crack lips.\r
20 \r
21 6. Finally the cracked box mesh is reinserted in the initial mesh ensuring the connectivity (k).\r
22 \r
23 +---------------------------------------+---------------------------------------+---------------------------------------+\r
24 | .. image:: images/general_princ_a.png | .. image:: images/general_princ_b.png | .. image:: images/general_princ_c.png |\r
25 |   :height: 600                        |    :height: 184                       |   :height: 184                        |\r
26 |   :align: center                      |    :align: center                     |   :align: center                      |\r
27 |                                       |                                       |                                       |\r
28 | (a)                                   | (b)                                   | (c)                                   |\r
29 +                                       +---------------------------------------+---------------------------------------+\r
30 |                                       | .. image:: images/general_princ_d.png | .. image:: images/general_princ_e.png |\r
31 |                                       |    :height: 184                       |   :height: 184                        |\r
32 |                                       |    :align: center                     |   :align: center                      |\r
33 |                                       |                                       |                                       |\r
34 |                                       | (d)                                   | (e)                                   |\r
35 +                                       +---------------------------------------+---------------------------------------+\r
36 |                                       | .. image:: images/general_princ_f.png | .. image:: images/general_princ_g.png |\r
37 |                                       |    :height: 184                       |   :height: 184                        |\r
38 |                                       |    :align: center                     |   :align: center                      |\r
39 |                                       |                                       |                                       |\r
40 |                                       | (f)                                   | (g)                                   |\r
41 +---------------------------------------+---------------------------------------+---------------------------------------+\r
42 | .. image:: images/general_princ_h.png | .. image:: images/general_princ_k.png |                                       |\r
43 |   :height: 184                        |    :height: 600                       |                                       |\r
44 |   :align: center                      |    :align: center                     |                                       |\r
45 |                                       |                                       |                                       |\r
46 | (h)                                   | (k)                                   |                                       |\r
47 +---------------------------------------+                                       +                                       +\r
48 | .. image:: images/general_princ_i.png |                                       |                                       |\r
49 |   :height: 184                        |                                       |                                       |\r
50 |   :align: center                      |                                       |                                       |\r
51 |                                       |                                       |                                       |\r
52 | (i)                                   |                                       |                                       |\r
53 +---------------------------------------+                                       +                                       +\r
54 | .. image:: images/general_princ_j.png |                                       |                                       |\r
55 |   :height: 184                        |                                       |                                       |\r
56 |   :align: center                      |                                       |                                       |\r
57 |                                       |                                       |                                       |\r
58 | (j)                                   |                                       |                                       |\r
59 +---------------------------------------+---------------------------------------+---------------------------------------+\r
60 \r
61 The length of influence is important. It defines the size of the extracted Box. Here is a 2D example showing how are selected the elements of the Box :\r
62 \r
63 .. image:: images/longueur_influence.png\r
64    :width: 600\r
65    :align: center\r
66    :alt: The length of influence\r
67 \r
68 Every element having a node at a smaller distance to the crack than the length of influence is selected. Then a filling algorithm fulfill the Box with elements to get a Box. The Box is not limited to rectangular shapes. See the section on :ref:`test cases <test_cases>` to see examples.\r
69 \r
70 .. _recommendations:\r
71 \r
72 ######################################\r
73 Recommendations and limitations\r
74 ######################################\r
75 \r
76 1) **Linear input mesh:**\r
77 \r
78 The input structure mesh must have linear elements in the extracted Box. Elements elsewhere do not have this limitation. The user can easily convert the mesh from quadratic to linear in MESH module in SALOME. Note that the element curvature of quadratic element is lost in the conversion.\r
79 \r
80 2) **Hexahedral mesh:**\r
81 \r
82 The extracted Box must have ruled quadrangles on the external faces which are intersected by the crack. The algorithm for geometrical reconstruction of the box needs 4 nodes elements on the surface of the box mesh. It means that free tetrahedral elements meshes are incompatible with « Bloc Fissure».\r
83 \r
84 3) **Crack shall exceed from the structure:**\r
85 \r
86 Surface crack geometry shall exceed from the structure mesh. Boolean operation can be problematic if crack surface edges are strictly on the Box faces. On the other hand, surface crack shall not exceed too much. A good rule would be to ensure that the area of the crack in the structure is greater than the area of the crack outside of the structure. It also means that internal crack can't be handled.\r
87 \r
88 4) **Crack front edges must exceed from the structure:**\r
89 \r
90 For similar reasons, crack front edges must exceed from the structure mesh. The user shall be really careful when fusing crack front edges within the structure with edges outside of the structure because junction mustn’t be on the box external face. For example the following figure shows the bad and the good practice. In grey a 2D view of a structure to cut and in red the crack surface. Line 1 is the edge declared as the crack front. On the left case, Line 1 stops on the box boundary. Even if Line 1 is extended with Line 2 and 5, « Bloc Fissure» will fail. The good practice is to extend the Line 1 with the same shape. See how to extend the front edges in the :ref:`tutorials section <tutorials>`\r
91 \r
92 .. image:: images/schema_lignes1.png\r
93    :scale: 80\r
94    :alt: Crack front edges\r
95 \r
96 5) **Loss of element groups:**\r
97 \r
98 Crack insertion with « Bloc Fissure» doesn’t ensure elements groups conservation. All the groups that are partially or totally in the Box mesh are impacted. Groups are not deleted but newly created elements are not added to the initial groups. Users must be really careful and verify their groups after crack insertion.\r
99 \r
100 +-------------------------------------------+------------------------------------------+\r
101 |.. image:: images/recom_groupe_avant.png   | .. image:: images/recom_groupe_apres.png |\r
102 |   :height: 400                            |    :height: 400                          |\r
103 |   :align: center                          |    :align: center                        |\r
104 +-------------------------------------------+------------------------------------------+\r
105 \r
106 In this example, the initial group contains all the external faces of the mesh (left). After insertion, the group still exists but all elements of the Box are not in the group anymore (right).\r
107 \r
108 6) **Loss of interfaces:**\r
109 \r
110 Interfaces between elements are not kept in the Box zone during insertion. For example, in the case of a crack inserted at an interface between two materials, there is no way to rebuild exactly the original groups for affecting materials properties.\r
111 \r
112 +--------------------------------------------+---------------------------------------------+\r
113 |.. image:: images/recom_interface_avant.png | .. image:: images/recom_interface_apres.png |\r
114 |   :height: 300                             |    :height: 300                             |\r
115 |   :align: center                           |    :align: center                           |\r
116 +--------------------------------------------+---------------------------------------------+\r
117 \r
118 7) **Element normal:**\r
119 \r
120 Elements normal orientation doesn’t always follow the outgoing rule. User must verify faces orientation if fluxes boundary conditions are to be applied to the mesh.\r
121 \r
122 \r
123 \r
124 \r