Extractors

"Extractors" have been added to the model. They allow the user to read particular results such as the equivalent Von-Mises stress at given positions in the mesh after the FE analysis. The values are displayed in the console at the end of the calculation (if only 1 node is concerned). They are also saved to disc

• in a .tsv (tab-separated values) file when fossils is used,
• or in a .ascii file when Metafor is used.

The input syntax is documented [here in the wiki].(https://gitlab.uliege.be/rboman/fossils/-/wikis/input-file-for-boneload-simulations#extractors).

Some extractors have been added to models/dolicorhynchops/dolicorhynchops_10k.py:

    p['extractors'] = [
        {
            'name': 'tooth1',
            'nodes': [-10.20362, -17.46838, -229.9061],
            'variables': [ 'sigvm' ]
        },
        {
            'name': 'tooth2',
            'nodes': [-11.92466, 26.3042, -229.5354],
            'variables': [ 'sigxx', 'sigyy', 'sigzz', 'sigvm' ]
        },
        {
            'name': 'Lmuscle',
            'nodes': f'{path}/Lmuscle.stl',
            'variables': [ 'x0', 'y0', 'z0', 'sigvm' ]
        }
    ]

The values are calculated at the classical extractor list:

extracting results...
        internal energy = 95.219 N.mm
        contact_pts_Fx = 17.915 N
        [...]
        tooth1_sigvm = 4.507 MPa
        tooth2_sigxx = -1.000 MPa
        tooth2_sigyy = -1.025 MPa
        tooth2_sigzz = -0.687 MPa
        tooth2_sigvm = 3.803 MPa
        Lmuscle_x0 = (95 values in "Lmuscle_x0.tsv")
        Lmuscle_y0 = (95 values in "Lmuscle_y0.tsv")
        Lmuscle_z0 = (95 values in "Lmuscle_z0.tsv")
        Lmuscle_sigvm = (95 values in "Lmuscle_sigvm.tsv")
Total Execution time: 2.2s

cxxfem/src/femPost.cpp

@@ -24,11 +24,45 @@ Post::~Post()
         std::cout << "~Post()\n";
 }
 
-/// rebuild solution vectors for gmsh and add views
+/// rebuild solution vectors for gmsh and add views:
+///
+///  X
+///  Y
+///  Z
+///  displacement_vector
+///  force_vector
+///  stress_tensor
+///      stress_xx
+///      stress_yy
+///      stress_zz
+///      stress_xy
+///      stress_xz
+///      stress_yz
+///  strain_tensor
+///      strain_xx
+///      strain_yy
+///      strain_zz
+///      strain_xy
+///      strain_xz
+///      strain_yz
+///  smooth_stress_tensor
+///      smooth_stress_xx
+///      smooth_stress_yy
+///      smooth_stress_zz
+///      smooth_stress_xy
+///      smooth_stress_xz
+///      smooth_stress_yz
+///  smooth_strain_tensor
+///      smooth_strain_xx
+///      smooth_strain_yy
+///      smooth_strain_zz
+///      smooth_strain_xy
+///      smooth_strain_xz
+///      smooth_strain_yz
 
 void Post::build_views()
 {
-    std::cout << "creating views...\n";
+    std::cout << "creating views..." << std::endl;
     // nodal dofs
     add_dofs_view(solver.pbl, solver.dofs, solver.X, solver.dofPart, solver.dofIdx);