| 1572 | } |
| 1573 | |
| 1574 | Matrix MixedBeamColumn2d::getMd(int sec, Vector dShapeFcn, Vector dFibers, double L) { |
| 1575 | double xi[MAX_NUM_SECTIONS]; |
| 1576 | beamIntegr->getSectionLocations(numSections, L, xi); |
| 1577 | |
| 1578 | double x, A, B; |
| 1579 | x = L*xi[sec]; |
| 1580 | A = ( x/L - 2*pow(x/L,2) + pow(x/L,3) )*L; |
| 1581 | B = (-pow(x/L,2) + pow(x/L,3) )*L; |
| 1582 | |
| 1583 | theNaturalMatrix.Zero(); |
| 1584 | theNaturalMatrix(0,1) = A * ( dShapeFcn(1) - dFibers(1) ); |
| 1585 | theNaturalMatrix(0,2) = B * ( dShapeFcn(1) - dFibers(1) ); |
| 1586 | |
| 1587 | return theNaturalMatrix; |
| 1588 | } |
| 1589 | |
| 1590 | Matrix MixedBeamColumn2d::getNld_hat(int sec, const Vector &v, double L, bool geomLinear) { |
| 1591 | double xi[MAX_NUM_SECTIONS]; |
no test coverage detected