00001 /* 00002 00003 Copyright (C) University of Oxford, 2005-2011 00004 00005 University of Oxford means the Chancellor, Masters and Scholars of the 00006 University of Oxford, having an administrative office at Wellington 00007 Square, Oxford OX1 2JD, UK. 00008 00009 This file is part of Chaste. 00010 00011 Chaste is free software: you can redistribute it and/or modify it 00012 under the terms of the GNU Lesser General Public License as published 00013 by the Free Software Foundation, either version 2.1 of the License, or 00014 (at your option) any later version. 00015 00016 Chaste is distributed in the hope that it will be useful, but WITHOUT 00017 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00018 FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 00019 License for more details. The offer of Chaste under the terms of the 00020 License is subject to the License being interpreted in accordance with 00021 English Law and subject to any action against the University of Oxford 00022 being under the jurisdiction of the English Courts. 00023 00024 You should have received a copy of the GNU Lesser General Public License 00025 along with Chaste. If not, see <http://www.gnu.org/licenses/>. 00026 00027 */ 00028 00029 #include "MooneyRivlinMaterialLaw.hpp" 00030 00031 template<unsigned DIM> 00032 double MooneyRivlinMaterialLaw<DIM>::Get_dW_dI1(double I1, double I2) 00033 { 00034 return mC1; 00035 } 00036 00037 template<unsigned DIM> 00038 double MooneyRivlinMaterialLaw<DIM>::Get_dW_dI2(double I1, double I2) 00039 { 00040 /* 00041 * This is covered, but gcov doesn't see this as being covered 00042 * for some reason, maybe because of optimisations. 00043 */ 00044 #define COVERAGE_IGNORE 00045 assert(DIM == 3); 00046 #undef COVERAGE_IGNORE 00047 return mC2; 00048 } 00049 00050 template<unsigned DIM> 00051 double MooneyRivlinMaterialLaw<DIM>::Get_d2W_dI1(double I1, double I2) 00052 { 00053 return 0.0; 00054 } 00055 00056 template<unsigned DIM> 00057 double MooneyRivlinMaterialLaw<DIM>::Get_d2W_dI2(double I1, double I2) 00058 { 00059 /* 00060 * This is covered, but gcov doesn't see this as being covered 00061 * for some reason, maybe because of optimisations. 00062 */ 00063 #define COVERAGE_IGNORE 00064 assert(DIM == 3); 00065 #undef COVERAGE_IGNORE 00066 return 0.0; 00067 } 00068 00069 template<unsigned DIM> 00070 double MooneyRivlinMaterialLaw<DIM>::Get_d2W_dI1I2(double I1, double I2) 00071 { 00072 /* 00073 * This is covered, but gcov doesn't see this as being covered 00074 * for some reason, maybe because of optimisations. 00075 */ 00076 #define COVERAGE_IGNORE 00077 assert(DIM == 3); 00078 #undef COVERAGE_IGNORE 00079 return 0.0; 00080 } 00081 00082 template<unsigned DIM> 00083 double MooneyRivlinMaterialLaw<DIM>::GetC1() 00084 { 00085 return mC1; 00086 } 00087 00088 template<unsigned DIM> 00089 double MooneyRivlinMaterialLaw<DIM>::GetC2() 00090 { 00091 assert(DIM==3); 00092 return mC2; 00093 } 00094 00095 template<unsigned DIM> 00096 MooneyRivlinMaterialLaw<DIM>::MooneyRivlinMaterialLaw(double c1, double c2) 00097 : mC1(c1), 00098 mC2(c2) 00099 { 00100 assert(DIM==2 || DIM ==3); 00101 00102 // If dim==3, check that c2 was passed in, ie c2 isn't the default value 00103 if ((DIM==3) && (c2<MINUS_LARGE+1)) 00104 { 00105 EXCEPTION("Two parameters needed for 3d Mooney-Rivlin"); 00106 } 00107 00108 if (c1 < 0.0) 00109 { 00110 EXCEPTION("c1 must be positive in mooney-rivlin"); // is this correct? 00111 } 00112 } 00113 00114 template<unsigned DIM> 00115 void MooneyRivlinMaterialLaw<DIM>::ScaleMaterialParameters(double scaleFactor) 00116 { 00117 assert(scaleFactor > 0.0); 00118 mC1 /= scaleFactor; 00119 mC2 /= scaleFactor; 00120 } 00121 00123 // Explicit instantiation 00125 00126 template class MooneyRivlinMaterialLaw<2>; 00127 template class MooneyRivlinMaterialLaw<3>;