| 1246 | } |
| 1247 | |
| 1248 | void |
| 1249 | SSPquadUP::GetPermeabilityMatrix(void) |
| 1250 | // this function computes the permeability matrix for the element |
| 1251 | { |
| 1252 | mPerm.Zero(); |
| 1253 | Matrix k(2,2); |
| 1254 | Matrix dNp(2,4); |
| 1255 | |
| 1256 | // permeability tensor |
| 1257 | k(0,0) = perm[0]; |
| 1258 | k(1,1) = perm[1]; |
| 1259 | |
| 1260 | // transpose the shape function derivative array |
| 1261 | dNp(0,0) = dN(0,0); dNp(0,1) = dN(1,0); dNp(0,2) = dN(2,0); dNp(0,3) = dN(3,0); |
| 1262 | dNp(1,0) = dN(0,1); dNp(1,1) = dN(1,1); dNp(1,2) = dN(2,1); dNp(1,3) = dN(3,1); |
| 1263 | |
| 1264 | // compute permeability matrix |
| 1265 | mPerm.addMatrixTripleProduct(1.0, dNp, k, 4.0*J0*mThickness); |
| 1266 | |
| 1267 | return; |
| 1268 | } |
| 1269 | |
| 1270 | // LM change |
| 1271 | // nodes numbering |
nothing calls this directly
no test coverage detected