Method calculateAll

SRC/element/shell/ASDShellT3.cpp:1572–1848

}

int ASDShellT3::calculateAll(Matrix& LHS, Vector& RHS, int options)
{
    // Check options
    if (!m_transformation->isLinear()) {
        // corotational calculation of the tangent LHS requires the RHS!
        if (options & OPT_LHS)
            options |= OPT_RHS;
    }

    // Zero output
    int result = 0;
    if (options & OPT_RHS)
        RHS.Zero();
    if (options & OPT_LHS)
        LHS.Zero();

    // Global displacements
    auto& UG = ASDShellT3Globals::instance().UG;
    m_transformation->computeGlobalDisplacements(UG);

    // update transformation
    if (options & OPT_UPDATE)
        m_transformation->update(UG);

    // Compute the reference coordinate system
    ASDShellT3LocalCoordinateSystem reference_cs =
        m_transformation->createReferenceCoordinateSystem();

    // Compute the local coordinate system.
    ASDShellT3LocalCoordinateSystem local_cs =
        m_transformation->createLocalCoordinateSystem(UG);

    // Some matrices/vectors
    auto& N = ASDShellT3Globals::instance().N;
    auto& dN = ASDShellT3Globals::instance().dN;
    auto& dNdX = ASDShellT3Globals::instance().dNdX;
    auto& jac = ASDShellT3Globals::instance().jac;
    auto& B = ASDShellT3Globals::instance().B;
    auto& Bd = ASDShellT3Globals::instance().Bd;
    auto& Bhx = ASDShellT3Globals::instance().Bhx;
    auto& Bhy = ASDShellT3Globals::instance().Bhy;
    auto& B1 = ASDShellT3Globals::instance().B1;
    auto& B1TD = ASDShellT3Globals::instance().B1TD;
    auto& E = ASDShellT3Globals::instance().E;
    auto& S = ASDShellT3Globals::instance().S;
    auto& D = ASDShellT3Globals::instance().D;
    auto& Dsection = ASDShellT3Globals::instance().Dsection;

    // matrices for orienting strains in section coordinate system
    auto& Re = ASDShellT3Globals::instance().Re;
    auto& Rs = ASDShellT3Globals::instance().Rs;
    if (m_angle != 0.0) {
        if (options & OPT_UPDATE)
            getRotationMatrixForGeneralizedStrains(-m_angle, Re);
        if ((options & OPT_RHS) || (options & OPT_LHS))
            getRotationMatrixForGeneralizedStresses(m_angle, Rs);
    }