Method commitSensitivity

SRC/material/nD/BeamFiberMaterial2dPS.cpp:277–311

}

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

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

  static Matrix dd22(1,1);
  dd22(0,0) = PStangent(1,1);

  static Matrix dd21(1,2);
  dd21(0,0) = PStangent(1,0);
  dd21(0,1) = PStangent(1,2);
  
  static Vector sigma2(1);
  sigma2.addMatrixVector(0.0, dd21, depsdh, -1.0);

  const Vector &PSstress = theMaterial->getStressSensitivity(gradIndex, true);
  //opserr << PSstress;
  sigma2(0) -= PSstress(1);

  //const Vector &PSstress2 = theMaterial->getStressSensitivity(gradIndex, false);
  //opserr << PSstress2;
  //sigma2(0) += PSstress2(1);

  static Vector strain2(1);
  dd22.Solve(sigma2,strain2);

  dstraindh(0) = depsdh(0);
  dstraindh(1) = strain2(0);
  dstraindh(2) = depsdh(1);

  return theMaterial->commitSensitivity(dstraindh, gradIndex, numGrads);
}

const Matrix&  
BeamFiberMaterial2dPS::getTangent()