set transformation matrix from global to local system
| 1252 | |
| 1253 | // set transformation matrix from global to local system |
| 1254 | void TwoNodeLinkSection::setTranGlobalLocal() |
| 1255 | { |
| 1256 | // resize transformation matrix and zero it |
| 1257 | Tgl.resize(numDOF,numDOF); |
| 1258 | Tgl.Zero(); |
| 1259 | |
| 1260 | // switch on dimensionality of element |
| 1261 | switch (elemType) { |
| 1262 | case D1N2: |
| 1263 | Tgl(0,0) = Tgl(1,1) = trans(0,0); |
| 1264 | break; |
| 1265 | case D2N4: |
| 1266 | Tgl(0,0) = Tgl(2,2) = trans(0,0); |
| 1267 | Tgl(0,1) = Tgl(2,3) = trans(0,1); |
| 1268 | Tgl(1,0) = Tgl(3,2) = trans(1,0); |
| 1269 | Tgl(1,1) = Tgl(3,3) = trans(1,1); |
| 1270 | break; |
| 1271 | case D2N6: |
| 1272 | Tgl(0,0) = Tgl(3,3) = trans(0,0); |
| 1273 | Tgl(0,1) = Tgl(3,4) = trans(0,1); |
| 1274 | Tgl(1,0) = Tgl(4,3) = trans(1,0); |
| 1275 | Tgl(1,1) = Tgl(4,4) = trans(1,1); |
| 1276 | Tgl(2,2) = Tgl(5,5) = trans(2,2); |
| 1277 | break; |
| 1278 | case D3N6: |
| 1279 | Tgl(0,0) = Tgl(3,3) = trans(0,0); |
| 1280 | Tgl(0,1) = Tgl(3,4) = trans(0,1); |
| 1281 | Tgl(0,2) = Tgl(3,5) = trans(0,2); |
| 1282 | Tgl(1,0) = Tgl(4,3) = trans(1,0); |
| 1283 | Tgl(1,1) = Tgl(4,4) = trans(1,1); |
| 1284 | Tgl(1,2) = Tgl(4,5) = trans(1,2); |
| 1285 | Tgl(2,0) = Tgl(5,3) = trans(2,0); |
| 1286 | Tgl(2,1) = Tgl(5,4) = trans(2,1); |
| 1287 | Tgl(2,2) = Tgl(5,5) = trans(2,2); |
| 1288 | break; |
| 1289 | case D3N12: |
| 1290 | Tgl(0,0) = Tgl(3,3) = Tgl(6,6) = Tgl(9,9) = trans(0,0); |
| 1291 | Tgl(0,1) = Tgl(3,4) = Tgl(6,7) = Tgl(9,10) = trans(0,1); |
| 1292 | Tgl(0,2) = Tgl(3,5) = Tgl(6,8) = Tgl(9,11) = trans(0,2); |
| 1293 | Tgl(1,0) = Tgl(4,3) = Tgl(7,6) = Tgl(10,9) = trans(1,0); |
| 1294 | Tgl(1,1) = Tgl(4,4) = Tgl(7,7) = Tgl(10,10) = trans(1,1); |
| 1295 | Tgl(1,2) = Tgl(4,5) = Tgl(7,8) = Tgl(10,11) = trans(1,2); |
| 1296 | Tgl(2,0) = Tgl(5,3) = Tgl(8,6) = Tgl(11,9) = trans(2,0); |
| 1297 | Tgl(2,1) = Tgl(5,4) = Tgl(8,7) = Tgl(11,10) = trans(2,1); |
| 1298 | Tgl(2,2) = Tgl(5,5) = Tgl(8,8) = Tgl(11,11) = trans(2,2); |
| 1299 | break; |
| 1300 | } |
| 1301 | } |
| 1302 | |
| 1303 | |
| 1304 | // set transformation matrix from local to basic system |