testtet10.cpp

/*************************************************************************************
 * MechSys - A C++ library to simulate (Continuum) Mechanical Systems                *
 * Copyright (C) 2005 Dorival de Moraes Pedroso <dorival.pedroso at gmail.com>       *
 * Copyright (C) 2005 Raul Dario Durand Farfan  <raul.durand at gmail.com>           *
 *                                                                                   *
 * This file is part of MechSys.                                                     *
 *                                                                                   *
 * MechSys is free software; you can redistribute it and/or modify it under the      *
 * terms of the GNU General Public License as published by the Free Software         *
 * Foundation; either version 2 of the License, or (at your option) any later        *
 * version.                                                                          *
 *                                                                                   *
 * MechSys is distributed in the hope that it will be useful, but WITHOUT ANY        *
 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A   *
 * PARTICULAR PURPOSE. See the GNU General Public License for more details.          *
 *                                                                                   *
 * You should have received a copy of the GNU General Public License along with      *
 * MechSys; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, *
 * Fifth Floor, Boston, MA 02110-1301, USA                                           *
 *************************************************************************************/

#include <iostream>

#include "fem/fem.h"
#include "fem/ele/tet10equilib.h"
#include "linalg/vector.h"

using std::cout;
using std::endl;
using LinAlg::Vector;

int main(int argc, char **argv) try
{
    /* {{{                        
     *                         t
     *                         |
     *                         |                                           
     *                         | 3                                         
     *                         @,                                          
     *                        /|`                                          
     *                        ||  `,                                       
     *                       / |    ',                                     
     *                       | |      \                                    
     *                      /  |       `.                                  
     *                      |  |         `,  9                             
     *                     /   @ 7         `@                              
     *                     |   |             \                             
     *                    /    |              `.                           
     *                    |    |                ',                         
     *                 8 @     |                  \                        
     *                   |     @.,,_       6       `.                      
     *                  |     / 0   ``'-.,,@_        `.                    
     *                  |    /              ``''-.,,_  ', 2                
     *                 |    /                        ``'@.,,,              
     *                 |   '                       ,.-``     ``''- s       
     *                |  ,@ 4                 _,-'`                        
     *                ' /                 ,.'`                             
     *               | /             _.@``                                 
     *               '/          ,-'`   5                                  
     *              |/      ,.-``                                          
     *              /  _,-``                                               
     *            .@ '`                                                    
     *           / 1                                                       
     *          /                                                          
     *         /                                                           
     *        r                                                            
     }}} */
    Array<REAL> sub_cam_prms;
    sub_cam_prms.push_back(0.0891); // lam
    sub_cam_prms.push_back(0.0196); // kap
    sub_cam_prms.push_back(0.20);   // nu
    sub_cam_prms.push_back(3.2);    // Rcs
    sub_cam_prms.push_back(5000);   // c

    Array<REAL>         ini_data;
    Array<Array<REAL> > ini_vals;
    ini_data.push_back(2);     // Svert
    ini_data.push_back(2);     // Shorz
    ini_data.push_back(2);     // Shorz
    ini_data.push_back(1.691); // v_ini (specific volume)
    ini_data.push_back(1);     // OCR
    ini_vals.push_back(ini_data);

    Array<FEM::Node> N;
    N.resize(10);
    N[0].Initialize(0,  0, 0, 0);
    N[1].Initialize(1, 10, 0, 0);
    N[2].Initialize(2,  0,10, 0);
    N[3].Initialize(3,  0, 0,10);
    N[4].Initialize(4,  5, 0, 0);
    N[5].Initialize(5,  5, 5, 0);
    N[6].Initialize(6,  0, 5, 0);
    N[7].Initialize(7,  0, 0, 5);
    N[8].Initialize(8,  5, 0, 5);
    N[9].Initialize(9,  0, 5, 5);

    FEM::Tet10Equilib E;

    // Setup element connectivity and nodal vars
    for (int i=0; i<10; ++i)
        E.SetNode(i, &N[i]);

    // Setup element model
    E.ReAllocateModel("SubCam", sub_cam_prms, ini_vals);

    
    // Temporary array of pointers to the nodes of the face {1,2,3, 5,9,8}
    Array<FEM::Node*> a_ptr_face_nodes;
    a_ptr_face_nodes.push_back( &N[1] );
    a_ptr_face_nodes.push_back( &N[2] );
    a_ptr_face_nodes.push_back( &N[3] );
    a_ptr_face_nodes.push_back( &N[5] );
    a_ptr_face_nodes.push_back( &N[9] );
    a_ptr_face_nodes.push_back( &N[8] );

    // Calculate the nodal values for face {1,2,3,5,9,8} according to the type of the correspondent element
    String       str_dof_name;
    Vector<REAL> a_nodal_values;
    E.CalcFaceNodalValues("tz"   ,   // In:  Face DOF Name, ex.: "tx, ty, tz"
                          -10    ,   // In:  Face DOF Value, ex.: 10 kPa, 20 kPa
                          a_ptr_face_nodes,   // In:  Array of ptrs to face nodes
                          str_dof_name    ,   // Out: Nodal DOF Name, ex.: "fx, fy, fz"
                          a_nodal_values  );  // Out: Vector of nodal values

    // Loop along the nodes of face {1,2,3,5,9,8}
    for (int m=0; m<6; ++m)
    {
        // Get a reference to a DOFVarsStruct
        FEM::Node::DOFVarsStruct & dof_var = a_ptr_face_nodes[m]->DOFVar(str_dof_name);

        // Setup Natural nodal variable
          dof_var.NaturalBry += a_nodal_values(m); // Upgrade (add) existent values
        dof_var.IsEssenPresc  = false;
    }
    
    
    cout << "str_dof_name = " << str_dof_name << endl;
    cout << "Vector of nodal values=\n" << a_nodal_values << endl;

    cout << "Nodes=\n";
    for (int i=0; i<10; ++i) cout << N[i];

    return 0;
}
catch (Exception * e)
{
    e->Cout();
    if (e->IsFatal()) exit(1);
    delete e;
}
catch (char const * s)
{
    cout << " Fatal:" << s << "\n";
    exit(1);
}

// vim:fdm=marker

---