00001 /************************************************************************************* 00002 * MechSys - A C++ library to simulate (Continuum) Mechanical Systems * 00003 * Copyright (C) 2005 Dorival de Moraes Pedroso <dorival.pedroso at gmail.com> * 00004 * Copyright (C) 2005 Raul Dario Durand Farfan <raul.durand at gmail.com> * 00005 * * 00006 * This file is part of MechSys. * 00007 * * 00008 * MechSys is free software; you can redistribute it and/or modify it under the * 00009 * terms of the GNU General Public License as published by the Free Software * 00010 * Foundation; either version 2 of the License, or (at your option) any later * 00011 * version. * 00012 * * 00013 * MechSys is distributed in the hope that it will be useful, but WITHOUT ANY * 00014 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A * 00015 * PARTICULAR PURPOSE. See the GNU General Public License for more details. * 00016 * * 00017 * You should have received a copy of the GNU General Public License along with * 00018 * MechSys; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, * 00019 * Fifth Floor, Boston, MA 02110-1301, USA * 00020 *************************************************************************************/ 00021 00022 #include <iostream> 00023 00024 #include "util/array.h" 00025 #include "util/exception.h" 00026 #include "models/subtij.h" 00027 00028 using std::cout; 00029 using std::endl; 00030 00031 int main(int argc, char **argv) try 00032 { 00034 Array<REAL> sub_tij_prms; 00035 sub_tij_prms.push_back(0.0891); // lam 00036 sub_tij_prms.push_back(0.0196); // kap 00037 sub_tij_prms.push_back(0.20); // nu 00038 sub_tij_prms.push_back(3.5); // Rcs 00039 sub_tij_prms.push_back(500); // c 00040 sub_tij_prms.push_back(1.5); // bet 00041 00043 Array<REAL> ini_data; 00044 ini_data.push_back(2); // Sx 00045 ini_data.push_back(2); // Sy 00046 ini_data.push_back(2); // Sz 00047 ini_data.push_back(1.691); // v_ini (specific volume) 00048 ini_data.push_back(1); // OCR 00049 00051 SubTij ST(sub_tij_prms, ini_data); 00052 00053 return 0; 00054 } 00055 catch (Exception * e) 00056 { 00057 e->Cout(); 00058 if (e->IsFatal()) exit(1); 00059 delete e; 00060 } 00061 catch (char const * s) 00062 { 00063 cout << "[1;31m Fatal:" << s << "[0m\n"; 00064 exit(1); 00065 }