00001 #ifndef CREEPB_H 00002 #define CREEPB_H 00003 00004 #include "genfile.h" 00005 #include "alias.h" 00006 00007 /** 00008 class contains B3 model of creep 00009 00010 from t0 to t 00011 K_s shape factor slab=1.0, cylinder=1.15, sguare prism.=1.25, sphere=1.3, cube=1.55 00012 tb from concrete starts 00013 t_w age when drying begins 00014 (fc') is 28 day average cilinder strenght fc' [ksi] ksi=1000psi=6.895 MPa(f.e.6.454=44.5MPa)***6.381 00015 (w/c) is water-cement ratio of the mix by weight ***0.43 00016 (s/c) is send-cement ratio of the mix by weight ***3.4 00017 (g/c) is gravel-cement ratio of the mix by weight g/c=a/c-s/c ***1.98 00018 (a/c) is aggregate-cement ratio of the mix by weight a/c=g/c+s/c 00019 (a1) is coef. for cements of type I,II a1=1.00, III a1=0.93, IV a1=1.05 ***1.05 00020 (ro) is mass of concrete in [lb/ft3] =16.03 kg/m3 ***156 00021 (k_d) effective cross section thickness D=2*vs_s 00022 cs cement content in m3 .. kg/m3 00023 E0=(0.09+1/(1.7*(0.5*ro*ro*fc*1e-4)*(0.5*ro*ro*fc*1e-4))) 00024 Et=E0*sqrt(t/(4+0.85*t)) podle ACI Commite 209/II 00025 00026 components of other array: 00027 previous total strains (nc components) 00028 internal variables describing history (7 x nc components for this model) 00029 shrinkage and thermal strain (1 component) 00030 previous moisture (1 component) 00031 previous temperature (1 component) 00032 00033 */ 00034 class creepb 00035 { 00036 public: 00037 creepb (void); 00038 ~creepb (void); 00039 void creepinit (long ipp,double val,nonmechquant nmq); 00040 void read (FILE *in); 00041 double approx (vector &areacoord,vector &nodval); 00042 void inv_sym (matrix &a); 00043 void updateval(); 00044 void nlstresses (long ipp); 00045 void phase1 (long ipp); 00046 void phase2 (long ipp); 00047 void get_h (long ipp); 00048 void get_temp (long ipp); 00049 void matstiff (matrix &d, long ipp); 00050 void seps_time (matrix &screep,vector &sig); 00051 void get_desht (double &des_hn, double t0, double t); 00052 void b3_law (double &jt, double t0, double t); 00053 00054 00055 // =1 constant h, 00056 long type_h; 00057 // =1 constant temperature, 00058 long type_temp; 00059 // array containing numbers of components of stress and strain tensors 00060 long *cncomp; 00061 // total number of components of stress and strain tensors 00062 long tncomp; 00063 // number of approximated functions on the element 00064 long napfun; 00065 // stress/strain state 00066 strastrestate ss; 00067 // number of stresses on the element 00068 long nc; 00069 00070 vector retTime; 00071 vector ert; 00072 double desht; 00073 double ccTime; 00074 double ddTime; 00075 double t0; 00076 double timemat; 00077 double timeMax; 00078 long nRetTime; 00079 long imat; 00080 // Young's modulus of elasticity 00081 double e0; 00082 // Poisson's ratio 00083 double mi; 00084 // coefficient of thermal dilatancy 00085 double alfa; 00086 // time of casting 00087 double tb; 00088 00089 double t_w; 00090 // compression strength in MPa 00091 double fc; 00092 // water-cement ratio 00093 double wc; 00094 // sound/cement ratio 00095 double sc; 00096 // gravel-cement ratio 00097 double gc; 00098 double c_s; 00099 // coefficient of shape of structure 00100 double a1; 00101 // relative moisture 00102 double h_s; 00103 double h_slast; 00104 // temperature 00105 double temp_s; 00106 double temp_slast; 00107 // 00108 double k_s; 00109 // 00110 double r_s; 00111 // 00112 double ts; 00113 // 00114 double k_d; 00115 00116 matrix apom; 00117 }; 00118 00119 #endif