full constructor
| 140 | |
| 141 | // full constructor |
| 142 | SSPquadUP::SSPquadUP(int tag, int Nd1, int Nd2, int Nd3, int Nd4, NDMaterial &theMat, |
| 143 | double thick, double Kf, double Rf, double k1, double k2, |
| 144 | double eVoid, double alpha, double b1, double b2, |
| 145 | double Pup, double Plow, double Pleft, double Pright) |
| 146 | :Element(tag,ELE_TAG_SSPquadUP), |
| 147 | theMaterial(0), |
| 148 | mExternalNodes(SQUP_NUM_NODE), |
| 149 | mTangentStiffness(SQUP_NUM_DOF,SQUP_NUM_DOF), |
| 150 | mInternalForces(SQUP_NUM_DOF), |
| 151 | Q(SQUP_NUM_DOF), |
| 152 | mMass(SQUP_NUM_DOF,SQUP_NUM_DOF), |
| 153 | mDamp(SQUP_NUM_DOF,SQUP_NUM_DOF), |
| 154 | mNodeCrd(2,4), |
| 155 | dN(4,2), |
| 156 | Mmem(3,8), |
| 157 | Kstab(8,8), |
| 158 | mSolidK(8,8), |
| 159 | mSolidM(8,8), |
| 160 | mPerm(4,4), |
| 161 | mThickness(thick), |
| 162 | fBulk(Kf), |
| 163 | fDens(Rf), |
| 164 | mAlpha(alpha), |
| 165 | mPorosity(0), |
| 166 | applyLoad(0), |
| 167 | pressureLoad(12), |
| 168 | pressureUpperSide(Pup), |
| 169 | pressureLowerSide(Plow), |
| 170 | pressureLeftSide(Pleft), |
| 171 | pressureRightSide(Pright) |
| 172 | { |
| 173 | mExternalNodes(0) = Nd1; |
| 174 | mExternalNodes(1) = Nd2; |
| 175 | mExternalNodes(2) = Nd3; |
| 176 | mExternalNodes(3) = Nd4; |
| 177 | |
| 178 | mThickness = thick; |
| 179 | fBulk = Kf; |
| 180 | fDens = Rf; |
| 181 | mAlpha = alpha; |
| 182 | |
| 183 | b[0] = b1; |
| 184 | b[1] = b2; |
| 185 | |
| 186 | appliedB[0] = 0.0; |
| 187 | appliedB[1] = 0.0; |
| 188 | |
| 189 | perm[0] = k1; |
| 190 | perm[1] = k2; |
| 191 | |
| 192 | // LM change |
| 193 | P[0] = Pup; |
| 194 | P[1] = Plow; |
| 195 | P[2] = Pleft; |
| 196 | P[3] = Pright; |
| 197 | // LM change |
| 198 | |
| 199 | mPorosity = eVoid/(1.0 + eVoid); |