tcamclay.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/camclay.h"
#include "tensors/tensors.h"
#include "tensors/functions.h"
#include "fem/solver/intsolverdata.h"

using Tensors::Tensor2;
using Tensors::Tensor4;
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;

    // Parameters
    Array<REAL> prms; prms.resize(4);
    prms[0] = 0.25  ; // lam
    prms[1] = 0.05  ; // kap
    prms[2] = 20.0  ; // phics (deg)
    prms[3] = 10000 ; // G (kPa)

    // Initial Data
    Array<REAL> inidata; inidata.resize(5); //{ SigX SigY SigZ vini OCR }
    inidata[0] = 24;  // =z0 (kPa)
    inidata[1] = 24;  // =z0 (kPa)
    inidata[2] = 24;  // =z0 (kPa)
    inidata[3] = 1.8; // v_ini
    inidata[4] = OCR; // OCR

    // Cam-clay Model
    CamClay cam(prms, inidata);

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

    // TgStiffness
    Tensor4 Dep;
    cam.TgStiffness(Dep);
    cout << "Dep =\n" << Dep << endl;
    cout << endl << endl;

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

    // StressUpdate
    Tensor2 dsig;
    cout << ".... StressUpdate ...\n";
    cam.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);
} // }}}

// vim:fdm=marker

---