Method commitSensitivity

SRC/element/XMUelements/AC3D8HexWithSensitivity.cpp:1762–1838

int
AC3D8HexWithSensitivity::commitSensitivity(int gradNumber, int numGrads)
{
  double r  = 0.0;
  double rw = 0.0;
  double s  = 0.0;
  double sw = 0.0;
  double t  = 0.0;
  double tw = 0.0;
  double weight = 0.0;
  double det_of_Jacobian = 0.0;
  static const int nstress = 3 ;
  int success;

  static Vector stress (nstress);
  short where = 0;
  
  Matrix sigma(1,3);
  Matrix NF(1,nodes_in_elem);
  
  this->computeDiff();
  NF.Zero();

  static Matrix ul(1,3);

  for(short GP_c_r = 1; GP_c_r <= r_integration_order; GP_c_r++) {
    r = get_Gauss_p_c(r_integration_order, GP_c_r);
    rw = get_Gauss_p_w(r_integration_order, GP_c_r);
    for( short GP_c_s = 1; GP_c_s <= s_integration_order; GP_c_s++ ) {
      s = get_Gauss_p_c(s_integration_order, GP_c_s);
      sw = get_Gauss_p_w(s_integration_order, GP_c_s);
      for(short GP_c_t = 1; GP_c_t <= t_integration_order; GP_c_t++) {
        t = get_Gauss_p_c(t_integration_order, GP_c_t);
        tw = get_Gauss_p_w(t_integration_order, GP_c_t);
        
        det_of_Jacobian = detJ[where];
        Matrix &dhGlobal = *L[where];
        
        weight = sw * rw * tw * det_of_Jacobian;

        where++;

      }
    }
  }
  
 
  int i;

  Vector epsilon(3); 
  Matrix sstrain(3,1);
  
  ul(0,0)= theNodes[0]->getDispSensitivity(1,gradNumber);
  ul(0,1)= theNodes[1]->getDispSensitivity(2,gradNumber);
  ul(0,2)= theNodes[2]->getDispSensitivity(3,gradNumber);
 // opserr<<"ul"<<ul<<endln;

  for(i = 0; i < numGP; i++) {