crsecplstr.cpp

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

#include "crsecplstr.h"
#include "global.h"
#include "stochdriver.h"

crsecplstr::crsecplstr (void)
{
  t=0.0;  rho=0.0;  m=0.0;
}
crsecplstr::~crsecplstr (void)
{
  
}

void crsecplstr::read (XFILE *in)
{
  switch (Mp->tprob){
  case linear_statics:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case eigen_dynamics:{
    xfscanf (in,"%lf %lf %lf",&t,&rho,&m);
    break;
  }
  case forced_dynamics:{
    xfscanf (in,"%lf %lf %lf",&t,&rho,&m);
    break;
  }
  case mat_nonlinear_statics:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case geom_nonlinear_statics:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case earth_pressure:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case mech_timedependent_prob:
  case growing_mech_structure:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case layered_linear_statics:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case lin_floating_subdomain:{
    xfscanf (in,"%lf",&t);
    break;
  }
  case nonlin_floating_subdomain:{
    xfscanf (in,"%lf",&t);
    break;
  }
  default:{
    fprintf (stderr,"\n\n unknown problem type is required in function read (file %s, line %d)",__FILE__,__LINE__);
  }
  }
}

void crsecplstr::changeparam (atsel &atcs,vector &val)
{
  long i;
  
  for (i=0;i<atcs.num;i++){
    switch (atcs.atrib[i]){
    case 0:{
      t=val[i];
      break;
    }
    case 1:{
      rho=val[i];
      break;
    }
    default:{
      fprintf (stderr,"\n\n wrong number of atribute in function changeparam (file %s, line %d).\n",__FILE__,__LINE__);
    }
    }
  }
}