00001 #include "relaxeuroc.h"
00002 #include "stochdriver.h"
00003 #include "global.h"
00004
00005
00006
00007
00008
00009
00010
00011 relaxeuroc::relaxeuroc (void)
00012 {
00013
00014 siginit=0.0;
00015
00016 fpk=0.0;
00017
00018 a=0.0;
00019
00020 b=0.0;
00021
00022 e=0.0;
00023 }
00024
00025
00026
00027
00028
00029
00030
00031
00032 relaxeuroc::~relaxeuroc (void)
00033 {
00034
00035 }
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048 void relaxeuroc::read (XFILE *in)
00049 {
00050
00051
00052
00053
00054
00055 xfscanf (in,"%lf %lf %lf %lf %lf",&siginit,&fpk,&a,&b,&e);
00056 }
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067 void relaxeuroc::print (FILE *out)
00068 {
00069
00070
00071
00072
00073
00074 fprintf (out,"%le %le %le %le %le",siginit,fpk,a,b,e);
00075 }
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085 double relaxeuroc::stress_decrement (void)
00086 {
00087 double c,d,f,time;
00088
00089
00090 time = Mp->time/3600;
00091
00092
00093 c=b*siginit/fpk;
00094 d=0.75*(1.0-siginit/fpk);
00095
00096 f = siginit*a*exp(c)*pow(time/1000.0,d)/100000.0;
00097
00098 return f;
00099 }
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109 void relaxeuroc::stress (vector &sig,vector &eps,strastrestate ssst)
00110 {
00111 double c,d,f,time;
00112
00113
00114 time = Mp->time/3600.0;
00115
00116
00117 c=b*siginit/fpk;
00118 d=0.75*(1.0-siginit/fpk);
00119
00120 f = -(siginit - siginit*a*exp(c)*pow(time/1000.0,d)/100000.0);
00121
00122 switch (ssst){
00123 case bar:{
00124 sig[0]=f;
00125 eps[0]=f/e;
00126 break;
00127 }
00128 case spacestress:{
00129 sig[0]=f;
00130 sig[1]=f;
00131 sig[2]=f;
00132 sig[3]=0.0;
00133 sig[4]=0.0;
00134 sig[5]=0.0;
00135
00136
00137 eps[0]=f/e;
00138 eps[1]=f/e;
00139 eps[2]=f/e;
00140 eps[3]=0.0;
00141 eps[4]=0.0;
00142 eps[5]=0.0;
00143 break;
00144 }
00145 default:{
00146 print_err("unknown number of components of stress tensor is required", __FILE__, __LINE__, __func__);
00147 }
00148 }
00149
00150 }