Method commitSensitivity

SRC/material/nD/BeamFiberMaterial2d.cpp:357–427

}

int
BeamFiberMaterial2d::commitSensitivity(const Vector &depsdh, int gradIndex,
				       int numGrads)
{
  static Vector dstraindh(6);

  const Matrix &threeDtangent = theMaterial->getTangent();

  static Matrix dd22(4,4);
  dd22(0,0) = threeDtangent(1,1);
  dd22(1,0) = threeDtangent(2,1);
  dd22(2,0) = threeDtangent(4,1);
  dd22(3,0) = threeDtangent(5,1);
  
  dd22(0,1) = threeDtangent(1,2);
  dd22(1,1) = threeDtangent(2,2);
  dd22(2,1) = threeDtangent(4,2);
  dd22(3,1) = threeDtangent(5,2);
  
  dd22(0,2) = threeDtangent(1,4);
  dd22(1,2) = threeDtangent(2,4);
  dd22(2,2) = threeDtangent(4,4);
  dd22(3,2) = threeDtangent(5,4);
  
  dd22(0,3) = threeDtangent(1,5);
  dd22(1,3) = threeDtangent(2,5);
  dd22(2,3) = threeDtangent(4,5);
  dd22(3,3) = threeDtangent(5,5);

  static Matrix dd21(4,2);
  dd21(0,0) = threeDtangent(1,0);
  dd21(1,0) = threeDtangent(2,0);
  dd21(2,0) = threeDtangent(4,0);
  dd21(3,0) = threeDtangent(5,0);
  
  dd21(0,1) = threeDtangent(1,3);
  dd21(1,1) = threeDtangent(2,3);
  dd21(2,1) = threeDtangent(4,3);
  dd21(3,1) = threeDtangent(5,3);
  
  static Vector sigma2(4);
  sigma2.addMatrixVector(0.0, dd21, depsdh, -1.0);

  const Vector &threeDstress = theMaterial->getStressSensitivity(gradIndex, true);
  //opserr << threeDstress;
  sigma2(0) -= threeDstress(1);
  sigma2(1) -= threeDstress(2);
  sigma2(2) -= threeDstress(4);
  sigma2(3) -= threeDstress(5);

  //const Vector &threeDstress2 = theMaterial->getStressSensitivity(gradIndex, false);
  //opserr << threeDstress2;
  //sigma2(0) += threeDstress2(1);
  //sigma2(1) += threeDstress2(2);
  //sigma2(2) += threeDstress2(4);
  //sigma2(3) += threeDstress2(5);