ttij.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 <cmath>

#include "tensors/tensors.h"
#include "tensors/tijtensor.h"

using Tensors::Tensor2;
using Tensors::TijTensor;
using namespace std;

int main(int argc, char **argv) try
{
    //Tensor2 Sig; Sig = 2.01, 2.0, 2.0, 0.0, 0.0, 0.0;
    Tensor2 Sig; Sig = 10, 20, 30, sqrt(2.0), 2.0*sqrt(2.0), 3.0*sqrt(2);
    //Tensor2   Sig; Sig = 2.5,1.0,1.0, 0.0,0.0,-1.5*sqrt(2.0);
    REAL    tn;
    REAL    ts;
    REAL    Is[3];
    REAL    sqrt_rz;
    Tensor2 Tau;
    Tensor2 inv_Tau;
    Tensor2 a;
    Tensor2 t;
    Tensor2 dtn_ds;
    Tensor2 dts_ds;
    Tensor2 dts_dt;

    TijTensor Tij(0.00001);

    Tij.Derivs(Sig,
               tn,
               ts,
               Is,
               sqrt_rz,
               Tau,
               inv_Tau,
               a,
               t,
               dtn_ds,
               dts_ds,
               dts_dt);

    cout << "Sig     = " << Sig     << endl;
    cout << "tn      = " << tn      << endl;
    cout << "ts      = " << ts      << endl;
    cout << "Is[0]   = " << Is[0]   << endl;
    cout << "Is[1]   = " << Is[1]   << endl;
    cout << "Is[2]   = " << Is[2]   << endl;
    cout << "sqrt_rz = " << sqrt_rz << endl;
    cout << "Tau     = " << Tau     << endl;
    cout << "inv_Tau = " << inv_Tau << endl;
    cout << "a       = " << a       << endl;
    cout << "t       = " << t       << endl;
    cout << "dtn_ds  = " << dtn_ds  << endl;
    cout << "dts_ds  = " << dts_ds  << endl;
    cout << "dts_dt  = " << dts_dt  << endl;
    
    return 0;
}
catch (Exception * e)
{
    e->Cout();
    if (e->IsFatal()) exit(1);
    delete e;
}

---