tcscam.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/coupled/csubcam.h"

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

int main(int argc, char **argv) try
{
    Array<REAL> prms;
    prms.push_back(0.0891); // lam
    prms.push_back(0.0196); // kap
    prms.push_back(0.20);   // nu
    prms.push_back(3.2);    // Rcs
    prms.push_back(5000);   // c
    prms.push_back(1.0e-4); // k

    Array<REAL> idat;
    idat.push_back(2);     // Sx
    idat.push_back(2);     // Sy
    idat.push_back(2);     // Sz
    idat.push_back(1.691); // v_ini (specific volume)
    idat.push_back(1);     // OCR
    idat.push_back(0);     // Pp (pore-pressure)

    CSubCam csc(prms, idat);
    cout << "Sig_ini  = " << csc.Sig() << endl;

    
    // Actualize
    Tensor2 DSig; DSig=10,10,5, 0,0,0;
    REAL    DPp;  DPp=-2;
    Tensor2 DEps;
    REAL    DnSr;
    csc.Actualize(DSig,DPp, DEps,DnSr);

    // Output
    cout << "Sig_fin  = " << csc.Sig()  << endl;
    cout << "DEps (%) = " << DEps*100.0 << endl;
    cout << "DnSr     = " << DnSr       << 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);
} // }}}

// vim:fdm=marker

---