tbar.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/barcelona.h"

using Tensors::Tensor2;
using std::cout;
using std::endl;

int main(int argc, char **argv) try
{
    // Parameters
    Array<REAL> prms; prms.resize(13);
    prms[0]  = 0.25  ; // lam0
    prms[1]  = 0.05  ; // kap
    prms[2]  = 10000 ; // G (kPa)
    prms[3]  = 20.0  ; // phics (deg)
    prms[4]  = 0.9   ; // r
    prms[5]  = 0.1   ; // bet
    prms[6]  = 0.5   ; // k
    prms[7]  = 0.002 ; // lams
    prms[8]  = 0.001 ; // kaps
    prms[9]  = 101.3 ; // patm (kPa)
    prms[10] = 1     ; // pref
    prms[11] = 10    ; // A
    prms[12] = 10    ; // B

    // Initial Data
    Array<REAL> inidata; inidata.resize(7); //{ SigX SigY SigZ vini OCR z1 Pp }
    inidata[0] = 24;  // =z0 (kPa)
    inidata[1] = 24;  // =z0 (kPa)
    inidata[2] = 24;  // =z0 (kPa)
    inidata[3] = 1.8; // v_ini
    inidata[4] = 1;   // OCR
    inidata[5] = 20;  // z1 (kPa)
    inidata[6] = -10; // Pp (kPa);  ==> Suc = 0-(-10) = +10 (kPa)

    // Barcelona Model
    Barcelona bar(prms, inidata);

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

    // Output
    cout << "DEps (%) = " << DEps*100.0 << 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

---