creep_dpl.h

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

#ifndef DPL_H
#define DPL_H

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

/**
   class contains Double power law model of creep
*/

class dplmat
{
 public:
  dplmat (void);
  ~dplmat (void);
  void read (XFILE *in);
  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);
  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 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);
  double double_power_law (double tt, double t,long ipp);
  void give_deps_free (long ipp);
  void give_deps_stressinduced (long ipp);
  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);

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

  //according to Bazant's notation:
  // t  = t  ... age of concrete in days
  // tt = t' ... age at loading [days]

  //  density .. kg/m3
  double ro;
  //  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;
  //  w/c is water-cement ratio of the mix by weight                                             = 0.43
  double wc;
  //  s/c is send-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;
  //  coefficient of shape of structure
  double a1;

  //time when structure is finished (concrete casting)
  double tb_time;
  // time when temperature and humidity start to change
  double th_time;
  //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;
  //  number of stresses on the element
  long nc;

  double *retTime;
  double *emu;
  double timeMax;
};

#endif