computes the B matrix for drilling DOF according to Huges and Brezzi used only in case of reduced integration
| 286 | // computes the B matrix for drilling DOF according to Huges and Brezzi |
| 287 | // used only in case of reduced integration |
| 288 | void computeBdrilling( |
| 289 | const Matrix& dNdX, const Vector& N, |
| 290 | Vector& Bd, Vector& Bhx, Vector& Bhy) |
| 291 | { |
| 292 | Bd.Zero(); |
| 293 | Bhx.Zero(); |
| 294 | Bhy.Zero(); |
| 295 | for (int i = 0; i < 3; ++i) { |
| 296 | int ii = i * 6; |
| 297 | Bd(ii) = -0.5 * dNdX(i, 1); |
| 298 | Bd(ii + 1) = 0.5 * dNdX(i, 0); |
| 299 | Bd(ii + 5) = -N(i); |
| 300 | Bhx(ii + 5) = dNdX(i, 0); |
| 301 | Bhy(ii + 5) = dNdX(i, 1); |
| 302 | } |
| 303 | } |
| 304 | |
| 305 | // computes the complete B matrix (membrane, bending and shear) |
| 306 | void computeBMatrix( |