int
FiberSection3dThermal::commitSensitivity(const Vector& defSens, int gradIndex, int numGrads)
{
// here add SHVs to store the strain sensitivity.
if (SHVs == 0) {
SHVs = new Matrix(3,numGrads);
}
(*SHVs)(0,gradIndex) = defSens(0);
(*SHVs)(1,gradIndex) = defSens(1);
(*SHVs)(2,gradIndex) = defSens(2);
int loc = 0;
double d0 = defSens(0);
double d1 = defSens(1);
double d2 = defSens(2);
for (int i = 0; i < numFibers; i++) {
UniaxialMaterial *theMat = theMaterials[i];
double y = matData[loc++] - yBar;
double z = matData[loc++] - zBar;
loc++; // skip A data.
double strainSens = d0 + y*d1 + z*d2;
theMat->commitSensitivity(strainSens,gradIndex,numGrads);
}
return 0;
}
int
FiberSection3d::commitSensitivity(const Vector& defSens, int gradIndex, int numGrads)
{
double d0 = defSens(0);
double d1 = defSens(1);
double d2 = defSens(2);
//dedh = defSens;
static double yLocs[10000];
static double zLocs[10000];
if (sectionIntegr != 0)
sectionIntegr->getFiberLocations(numFibers, yLocs, zLocs);
else {
for (int i = 0; i < numFibers; i++) {
yLocs[i] = matData[3*i];
zLocs[i] = matData[3*i+1];
}
}
static double dydh[10000];
static double dzdh[10000];
if (sectionIntegr != 0)
sectionIntegr->getLocationsDeriv(numFibers, dydh, dzdh);
else {
for (int i = 0; i < numFibers; i++) {
dydh[i] = 0.0;
dzdh[i] = 0.0;
}
}
double y, z;
double depsdh = 0;
for (int i = 0; i < numFibers; i++) {
UniaxialMaterial *theMat = theMaterials[i];
y = yLocs[i] - yBar;
z = zLocs[i] - zBar;
// determine material strain and set it
depsdh = d0 - y*d1 + z*d2 - dydh[i]*e(1) + dzdh[i]*e(2);
theMat->commitSensitivity(depsdh,gradIndex,numGrads);
}
return 0;
}
int
FiberSection2d::commitSensitivity(const Vector& defSens,
int gradIndex, int numGrads)
{
double d0 = defSens(0);
double d1 = defSens(1);
dedh = defSens;
static double fiberLocs[10000];
if (sectionIntegr != 0)
sectionIntegr->getFiberLocations(numFibers, fiberLocs);
else {
for (int i = 0; i < numFibers; i++)
fiberLocs[i] = matData[2*i];
}
static double locsDeriv[10000];
static double areaDeriv[10000];
if (sectionIntegr != 0) {
sectionIntegr->getLocationsDeriv(numFibers, locsDeriv);
sectionIntegr->getWeightsDeriv(numFibers, areaDeriv);
}
else {
for (int i = 0; i < numFibers; i++) {
locsDeriv[i] = 0.0;
areaDeriv[i] = 0.0;
}
}
double y;
double kappa = e(1);
for (int i = 0; i < numFibers; i++) {
UniaxialMaterial *theMat = theMaterials[i];
y = fiberLocs[i] - yBar;
// determine material strain and set it
double strainSens = d0 - y*d1 - locsDeriv[i]*kappa;
theMat->commitSensitivity(strainSens,gradIndex,numGrads);
}
return 0;
}
