creep_rspec.h

Finite element program for mechanical problem. It can solve various problem in solid problem

#ifndef RSPEC_H
#define RSPEC_H

#include "genfile.h"
#include "alias.h"
#include "iotools.h"

/**
   class contains 
*/

class rspecmat
{
 public:
  rspecmat (void);
  ~rspecmat (void);
  void read (XFILE *in);
  void print (FILE *out);
  void give_ages (double &tb_age_dt,double &tb_age,double &tl_age,double &th_age,double &dt,double &maxtime,long &napptime,long ipp);
  long give_nceqother (long ipp);
  long give_nret_time (void);

  double give_L(double q3,double q,double tau);
  double give_J_E_mu(vector &e_mu,double t0, double tl, double t, long ipp,long ido);
  double give_q1();
  double give_C_const();
  double give_q4();
  double give_nonlin_func();

  double give_inv_v(double time_mid);

  void give_rettimes (vector &rettimes,long n_ret_times,long ipp);

  void store_emu(vector e_mu,long n_ret_times,long ipp,long ido);
  void give_emu(vector &e_mu,long n_ret_times,long ipp,long ido);

  void give_hidden_strains(matrix &gamma_mu,long ipp,long ido);
  void store_hidden_strains(matrix gamma_mu,long ipp,long ido);

  void give_stresses(vector &sigma,long ipp,long ido);
  void store_stresses(vector sigma,long ipp,long ido);
 
  void give_dstresses(vector &dsigma,long ipp,long ido);
  void store_dstresses(vector dsigma,long ipp,long ido);
  void store_dstresses_add(vector dsigma,long ipp,long ido);
 
  void give_strains(vector &eps,long ipp,long ido);
  void store_strains(vector eps,long ipp,long ido);
 
  void give_dstrains(vector &deps,long ipp,long ido);
  void store_dstrains(vector deps,long ipp,long ido);
  
  void give_creepdstrains(vector &deps_cr,long ipp,long ido);
  void store_creepdstrains(vector deps_cr,long ipp,long ido);
  
  void give_irrdstrains(vector &deps_sh,long ipp,long ido);
  void store_irrdstrains(vector deps_sh,long ipp,long ido);
  
  void give_stressirrdstrains(vector &deps_ss,long ipp,long ido);
  void store_stressirrdstrains(vector deps_ss,long ipp,long ido);

  double give_shrinkage(long ipp,long ido);
  void store_shrinkage(double eps_sh,long ipp,long ido);

  void store_hum(long ipp,long ido);
  void store_temp(long ipp,long ido);

  double give_tb_time(long ipp);
  double give_th_time(long ipp);
  long give_napproxtime(long ipp);
  void initvalues (long ipp, long ido);

  void nlstresses (long ipp,long im,long ido);
  void give_deps_free (vector &deps_sh, double t0, double t_dt, double t, double dt, long ipp,long im,long ido);
  void give_deps_stressinduced (vector &deps_ss, double t0, double t_dt, double t, vector sigma, long ipp,long im,long ido);
  void storeirrstrains (vector deps,long ipp, long ido);
  void giveirrstrains (long ipp, long ido, vector &epscr);

  double get_othervalue(long compother,long ipp);
  void print_othervalue_name(FILE *out,long compother);
  double give_actual_ym (long ipp,long im,long ido);
  double give_actual_ft (long ipp,long im,long ido);
  double give_actual_fc (long ipp,long im,long ido) ;

  // = 0 constant h, 
  long type_h;
  // = 0 constant temperature, 
  long type_temp;

  // = 1 measured Young's modulus E_28 (28 day) psi 
  long type_e;

  // environmental humidity and temperature (if they are not changing)
  double hum_env,temp_env;

  //according to Bazant's notation:
  // t  = t  ... age of concrete in days
  // tl = t' ... age at loading [days]
  // t0 = t0 ... age when drying begins [days] (only t0 <= tl is considered)

  //flags: drying shrinkage, drying shrinkage, thermal shrinkage
  int flag_drshr,flag_shr,flag_temp;

  //parameters
  double n_param,m_param,q1,q2,q3,q4,q5;
  double C_const,min_ret_time;//minimum retardation time used in calculation [days], truncated values will be hidden in C_{const}, should be about 1.0e-16 for exact q1 determination, otherwise can be like 1.e-3

  //  coefficient of thermal dilatancy
  double alpha;
  //  e28 is 28 day Young's modulus
  double e28;
  //  fc' is 28 day average cylinder strength fc' [psi] 1000psi=6.895 MPa(f.e.6.454=44.5MPa)     = 6381.0
  double fc;
  //  ft' is 28 day average cylinder strength ft' [psi] 1000psi=6.895 MPa(f.e.6.454=44.5MPa)     = 267.36
  double ft;
  // ft_ratio = actual tensile strength/actual young modulus
  double ft_ratio;
  //  w/c is water-cement ratio of the mix by weight                                             = 0.43
  double wc;
  //  s/c is sand-cement ratio of the mix by weight                                              = 3.4
  double sc;
  //  g/c is gravel-cement ratio of the mix by weight g/c=a/c-s/c                                = 1.98
  double gc;
  //  cs  cement content in m3  .. kg/m3
  double cs;
  //  coefficient of shape of structure
  double a1;
  //  coefficient for curing
  double a2;
  //  k_s shape factor slab=1.0, cylinder=1.15, square prism.=1.25, sphere=1.3, cube=1.55 
  double ks;
  //  k_d effective cross section thickness D=2*vs_s in inches (inch = 25.4mm)
  double kd;

  //time when structure is finished (concrete casting)
  double tb_time;
  // time when temperature and humidity start to change
  double th_time;
  //age of concrete when drying began [days]
  double t_0;
  //number of times of approximation
  long napproxtime;
  //number of retardation times
  long nRetTime;
  // type of reading of retardation times 0 = computing; 1 = reading
  long type_rt;

  //  asymptotic Young's modulus
  double e0;

  //actual time
  double actualtime;
  //time step
  double dt;

  //  stress/strain state
  strastrestate ss;

  double *retTime;
  double timeMax;

  //  total autogenous shrinkage for concrete
  double eps_ainf;
};

#endif