Method setDomain

SRC/element/UWelements/SSPquadUP.cpp:281–332

}

void
SSPquadUP::setDomain(Domain *theDomain)
{
    theNodes[0] = theDomain->getNode(mExternalNodes(0));
    theNodes[1] = theDomain->getNode(mExternalNodes(1));
    theNodes[2] = theDomain->getNode(mExternalNodes(2));
    theNodes[3] = theDomain->getNode(mExternalNodes(3));

    for (int i = 0; i < 4; i++) {
        if (theNodes[i] == 0) {
            return;  // don't go any further - otherwise segmentation fault
        }
    }

    // initialize coordinate vectors
    const Vector &mIcrd_1 = theNodes[0]->getCrds();
    const Vector &mIcrd_2 = theNodes[1]->getCrds();
    const Vector &mIcrd_3 = theNodes[2]->getCrds();
    const Vector &mIcrd_4 = theNodes[3]->getCrds();

    // coordinate matrix
    mNodeCrd(0,0) = mIcrd_1(0);
    mNodeCrd(1,0) = mIcrd_1(1);
    mNodeCrd(0,1) = mIcrd_2(0);
    mNodeCrd(1,1) = mIcrd_2(1);
    mNodeCrd(0,2) = mIcrd_3(0);
    mNodeCrd(1,2) = mIcrd_3(1);
    mNodeCrd(0,3) = mIcrd_4(0);
    mNodeCrd(1,3) = mIcrd_4(1);

    // establish jacobian terms 
    J0 = ((mNodeCrd(0,1)-mNodeCrd(0,3))*(mNodeCrd(1,2)-mNodeCrd(1,0))+(mNodeCrd(0,2)-mNodeCrd(0,0))*(mNodeCrd(1,3)-mNodeCrd(1,1)))/8;
    J1 = ((mNodeCrd(0,1)-mNodeCrd(0,0))*(mNodeCrd(1,2)-mNodeCrd(1,3))+(mNodeCrd(0,2)-mNodeCrd(0,3))*(mNodeCrd(1,0)-mNodeCrd(1,1)))/24;
    J2 = ((mNodeCrd(0,0)-mNodeCrd(0,3))*(mNodeCrd(1,2)-mNodeCrd(1,1))+(mNodeCrd(0,2)-mNodeCrd(0,1))*(mNodeCrd(1,3)-mNodeCrd(1,0)))/24;

    // establish stabilization terms (based on initial material tangent, only need to compute once)
    GetStab();

    // establish mass matrix for solid phase (constant, only need to compute once)
    GetSolidMass();

    // establish permeability matrix (constant, only need to compute once)
    GetPermeabilityMatrix();

    //LM change
	// Compute consistent nodal loads due to surface pressure (at any side)
    this->setPressureLoadAtNodes();
	//LM change
		
    // call the base-class method
    this->DomainComponent::setDomain(theDomain);
}

int
SSPquadUP::commitState(void)