void
FourNodeQuadWithSensitivity::setPressureLoadAtNodes(void)
{
pressureLoad.Zero();
if (pressure == 0.0)
return;
const Vector &node1 = theNodes[0]->getCrds();
const Vector &node2 = theNodes[1]->getCrds();
const Vector &node3 = theNodes[2]->getCrds();
const Vector &node4 = theNodes[3]->getCrds();
double x1 = node1(0);
double y1 = node1(1);
double x2 = node2(0);
double y2 = node2(1);
double x3 = node3(0);
double y3 = node3(1);
double x4 = node4(0);
double y4 = node4(1);
double dx12 = x2-x1;
double dy12 = y2-y1;
double dx23 = x3-x2;
double dy23 = y3-y2;
double dx34 = x4-x3;
double dy34 = y4-y3;
double dx41 = x1-x4;
double dy41 = y1-y4;
double pressureOver2 = pressure/2.0;
// Contribution from side 12
pressureLoad(0) += pressureOver2*dy12;
pressureLoad(2) += pressureOver2*dy12;
pressureLoad(1) += pressureOver2*-dx12;
pressureLoad(3) += pressureOver2*-dx12;
// Contribution from side 23
pressureLoad(2) += pressureOver2*dy23;
pressureLoad(4) += pressureOver2*dy23;
pressureLoad(3) += pressureOver2*-dx23;
pressureLoad(5) += pressureOver2*-dx23;
// Contribution from side 34
pressureLoad(4) += pressureOver2*dy34;
pressureLoad(6) += pressureOver2*dy34;
pressureLoad(5) += pressureOver2*-dx34;
pressureLoad(7) += pressureOver2*-dx34;
// Contribution from side 41
pressureLoad(6) += pressureOver2*dy41;
pressureLoad(0) += pressureOver2*dy41;
pressureLoad(7) += pressureOver2*-dx41;
pressureLoad(1) += pressureOver2*-dx41;
//pressureLoad = pressureLoad*thickness;
}
void
Tri31::setPressureLoadAtNodes(void)
{
pressureLoad.Zero();
if (pressure == 0.0) return;
const Vector &node1 = theNodes[0]->getCrds();
const Vector &node2 = theNodes[1]->getCrds();
const Vector &node3 = theNodes[2]->getCrds();
double x1 = node1(0);
double y1 = node1(1);
double x2 = node2(0);
double y2 = node2(1);
double x3 = node3(0);
double y3 = node3(1);
double dx12 = x2-x1;
double dy12 = y2-y1;
double dx23 = x3-x2;
double dy23 = y3-y2;
double dx31 = x1-x3;
double dy31 = y1-y3;
double pressureOver2 = pressure/2.0;
// Contribution from side 12
pressureLoad(0) += pressureOver2*dy12;
pressureLoad(2) += pressureOver2*dy12;
pressureLoad(1) += pressureOver2*-dx12;
pressureLoad(3) += pressureOver2*-dx12;
// Contribution from side 23
pressureLoad(2) += pressureOver2*dy23;
pressureLoad(4) += pressureOver2*dy23;
pressureLoad(3) += pressureOver2*-dx23;
pressureLoad(5) += pressureOver2*-dx23;
// Contribution from side 31
pressureLoad(4) += pressureOver2*dy31;
pressureLoad(0) += pressureOver2*dy31;
pressureLoad(5) += pressureOver2*-dx31;
pressureLoad(1) += pressureOver2*-dx31;
//pressureLoad = pressureLoad*thickness;
}
