tscam.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 "util/array.h"
#include "util/exception.h"
#include "models/subcam.h"

using std::cout;
using std::endl;

int main(int argc, char **argv) try
{
    REAL OCR    = 1.0;
    REAL DSigX  = 10;
    REAL DSigY  = 10;
    REAL DSigZ  = 10;
    int  FEndiv = 1;
    if (argc>=2) OCR    = atof(argv[1]);
    if (argc>=3) DSigX  = atof(argv[2]);
    if (argc>=4) DSigY  = atof(argv[3]);
    if (argc>=5) DSigZ  = atof(argv[4]);
    if (argc>=6) FEndiv = atoi(argv[5]);

    cout << "OCR    = " << OCR    << endl;
    cout << "DSigX  = " << DSigX  << endl;
    cout << "DSigY  = " << DSigY  << endl;
    cout << "DSigZ  = " << DSigZ  << endl;
    cout << "FEndiv = " << FEndiv << endl;
    cout << endl << endl;

    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;
    ini_data.push_back(2);     // Sx
    ini_data.push_back(2);     // Sy
    ini_data.push_back(2);     // Sz
    ini_data.push_back(1.691); // v_ini (specific volume)
    ini_data.push_back(1);     // OCR

    SubCam SC(sub_cam_prms, ini_data);

    // Set Integration Schemes Constants Structure
    IntegSchemesCtes isc;
    isc.Type     (IntegSchemesCtes::FE);
    isc.FE_ndiv  (FEndiv);
    isc.ME_maxSS (20);
    isc.ME_STOL  (1e-1);
    isc.ME_dTini (0.1);
    isc.EdMax    (20);
    EquilibModel::SetIntegSchemesCtes(isc);

    // Actualize
    Tensor2 DSig; DSig=DSigX,DSigY,DSigZ, 0,0,0;
    Tensor2 DEps;
    cout << ".... Actualize ...\n";
    SC.BackupState();
    SC.Actualize(DSig, DEps);
    SC.RestoreState();

    // StressUpdate
    Tensor2 dsig;
    cout << ".... StressUpdate ...\n";
    SC.StressUpdate(DEps, dsig);

    // Output
    cout << "DSig     = " << DSig*1     << endl;
    cout << "DEps (%) = " << DEps*100.0 << endl;
    cout << "dsig     = " << dsig*1     << endl;

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

---