diff --git a/geometry/3d_beam_nx/beam_3d.jl b/geometry/3d_beam_nx/beam_3d.jl index 53b92b4..a6e8ec2 100644 --- a/geometry/3d_beam_nx/beam_3d.jl +++ b/geometry/3d_beam_nx/beam_3d.jl @@ -64,8 +64,24 @@ function create_interface(mesh, slave::Problem, master::Problem) return create_interface(mesh, slave_surface, master_surface) end +""" Convert Mesh object from quadratic to linear. """ +function to_linear!(mesh) + mapping = Dict(:Tet10 => :Tet4, :Tri6 => :Tri3) + nnodes = Dict(:Tet4 => 4, :Tri3 => 3) + for elid in keys(mesh.elements) + eltype = mesh.element_types[elid] + if haskey(mapping, eltype) + mesh.element_types[elid] = mapping[eltype] + nnodes_new = nnodes[mapping[eltype]] + mesh.elements[elid] = mesh.elements[elid][1:nnodes_new] + end + end +end + # start of simulation -mesh = abaqus_read_mesh("beam_3d_2nd_order_tetra_30mm.inp") +mesh = abaqus_read_mesh("beam_3d_1st_order_tetra_30mm.inp") +#mesh = abaqus_read_mesh("beam_3d_2nd_order_tetra_30mm.inp") +#to_linear!(mesh) info("element sets = ", collect(keys(mesh.element_sets))) info("surface sets = ", collect(keys(mesh.surface_sets))) @@ -95,7 +111,7 @@ update!(load, "surface pressure", 20.0) # solution isfile("results.h5") && rm("results.h5") isfile("results.xmf") && rm("results.xmf") -solver = Solver(Linear, beam, load) +solver = Solver(Linear, beam, bc1, load) solver.xdmf = Xdmf("results") solver() @@ -116,35 +132,40 @@ end function lsq_fit(elements, field) A = SparseMatrixCOO() b = SparseMatrixCOO() + volume = 0.0 for element in elements - gdofs = get_connectivity(element) - for ip in get_integration_points(element) - detJ = element(ip, time, Val{:detJ}) - w = ip.weight*detJ - N = element(ip, time) - f = field(element, ip, time) - add!(A, gdofs, gdofs, w*kron(N', N)) - for i=1:length(f) - add!(b, gdofs, w*f[i]*N, i) + gdofs = get_connectivity(element) + # increase integration order by 1 from default + for ip in get_integration_points(element, 1) + detJ = element(ip, time, Val{:detJ}) + w = ip.weight*detJ + N = element(ip, time) + f = field(element, ip, time) + add!(A, gdofs, gdofs, w*kron(N', N)) + for i=1:length(f) + add!(b, gdofs, w*f[i]*N, i) + end + volume += w end - end end + info("Mass matrix for least-squares fit is assembled. Total volume to fit: $volume") A = sparse(A) b = sparse(b) A = 1/2*(A + A') SparseArrays.droptol!(A, 1.0e-6) + SparseArrays.dropzeros!(A) nz = get_nonzero_rows(A) - dropzeros!(A) - F = ldltfact(A) + F = ldltfact(A[nz,nz]) x = zeros(size(b)...) x[nz, :] = F \ b[nz, :] - nodal_values = Dict(i => vec(x[i,:]) for i=1:size(x,1)) + nodal_values = Dict(i => vec(x[i,:]) for i in nz) return nodal_values end +xdmf = get(solver.xdmf) S = lsq_fit(beam.elements, get_stress) # store values to xml/h5 @@ -161,7 +182,6 @@ end using JuliaFEM: new_dataitem, new_child, set_attribute, add_child, save! -xdmf = get(solver.xdmf) dataitem = new_dataitem(xdmf, "/Results/Time $(solver.time)/Nodal Fields/Stress", stress) frame = read(xdmf, "/Domain/Grid/Grid") attribute = new_child(frame, "Attribute") @@ -169,7 +189,7 @@ set_attribute(attribute, "Name", "Stress") set_attribute(attribute, "Center", "Node") set_attribute(attribute, "AttributeType", "Tensor6") add_child(attribute, dataitem) - save!(xdmf) + close(xdmf.hdf) diff --git a/geometry/3d_beam_nx/beam_3d_1st_order_tetra_30mm.inp b/geometry/3d_beam_nx/beam_3d_1st_order_tetra_30mm.inp new file mode 100644 index 0000000..7f140a7 --- /dev/null +++ b/geometry/3d_beam_nx/beam_3d_1st_order_tetra_30mm.inp @@ -0,0 +1,1002 @@ +*NODE , Nset = nset_csys0 + 1, 2.649428 , -21.93735 , 217.2934 + 2, 27.54531 , 1.108443 , 228.8077 + 3, 22.96453 , -8.66196 , 268.0618 + 4, -21.40007 , 25.22895 , 273.0769 + 5, 12.89412 , 25.35329 , 270.296 + 6, -26.46135 , 9.751407 , 240.604 + 7, -23.60307 , -0.5336556 , 132.7367 + 8, 0.8954823 , -16.16293 , 164.7006 + 9, -24.94313 , -25.07641 , 145.2563 + 10, 0.4479259 , -16.8575 , 135.2276 + 11, 0.02428401 , -16.84218 , 106.2713 + 12, -24.98255 , -0.06760177 , 102.1092 + 13, 0.0443055 , 16.39546 , 127.902 + 14, -24.45267 , 24.77173 , 148.6451 + 15, -24.99471 , 24.96161 , 118.372 + 16, -0.8997279 , 18.81353 , 158.0369 + 17, 25.03252 , -0.07771355 , 102.112 + 18, 24.95259 , -25.0321 , 115.7944 + 19, 27.72841 , -0.293104 , 71.95084 + 20, 25.11274 , 24.96942 , 88.37872 + 21, 25.00076 , -25.0199 , 86.33241 + 22, -24.97938 , -25.02628 , 115.7949 + 23, -25.72485 , 0.06487084 , 72.6856 + 24, -24.92338 , -25.0053 , 86.32683 + 25, 24.48854 , -0.5154036 , 132.3097 + 26, 25.93127 , 21.28833 , 150.2057 + 27, 27.64374 , -19.95454 , 170.0083 + 28, 23.88761 , -25.35627 , 145.2014 + 29, 25.11609 , 24.9708 , 58.08149 + 30, -0.2447172 , 14.5861 , 65.53416 + 31, 0.07591873 , 16.49115 , 97.88775 + 32, 22.75688 , -25.2836 , 57.18617 + 33, 27.63651 , -12.01023 , 32.6123 + 34, -13.12758 , -22.63027 , 28.33912 + 35, 14.68236 , 23.23247 , 23.33398 + 36, -19.17568 , 5.446418 , 27.5244 + 37, -8.809165 , -0.2430688 , 270.0499 + 38, 0.2661249 , -17.80944 , 76.81774 + 39, -24.49397 , -24.98426 , 56.02147 + 40, -24.98573 , 24.98295 , 88.37307 + 41, -20.69963 , 26.99919 , 54.47268 + 42, -24.72781 , 25.86675 , 178.8007 + 43, -24.81836 , 25.39212 , 209.6436 + 44, 6.460481 , 27.61566 , 215.3756 + 45, -16.78045 , 0.3982796 , 186.3691 + 46, 25.09822 , 24.95513 , 118.3715 + 47, 25.0783 , 23.82875 , 177.486 + 48, -25.0212 , -24.87815 , 174.2111 + 49, 29.75592 , 6.467265 , 198.0093 + 651, -50.0 , 25.22785 , 193.2267 + 652, -50.0 , 0.7430228 , 206.6289 + 653, -50.0 , -23.99774 , 220.0183 + 654, -50.0 , -22.07244 , 253.9352 + 655, -50.0 , 2.436702 , 237.9985 + 656, -50.0 , 25.99597 , 223.6004 + 657, -50.0 , 27.49457 , 253.4913 + 658, -50.0 , 29.9268 , 278.449 + 659, -50.0 , 6.497106 , 270.4035 + 660, -50.0 , -24.94024 , 100.9857 + 661, -50.0 , 0.06480544 , 87.29038 + 662, -50.0 , 24.81696 , 73.43247 + 663, -50.0 , 0.06294202 , 57.36694 + 664, -50.0 , -23.72317 , 41.31783 + 665, -50.0 , -24.7577 , 71.664 + 666, -50.0 , 23.77166 , 42.37573 + 667, -50.0 , 0.04352627 , 28.56641 + 668, -50.0 , 24.9761 , 103.515 + 669, -50.0 , 25.00222 , 133.3688 + 670, -50.0 , 25.06035 , 163.2248 + 671, -50.0 , -24.99111 , 130.1574 + 672, -50.0 , 0.01580934 , 116.9823 + 673, -50.0 , -24.94941 , 159.341 + 674, -50.0 , 0.03374581 , 146.5468 + 675, -50.0 , 0.2419424 , 176.3142 + 676, -50.0 , -24.65243 , 188.9462 + 778, -24.92087 , -50.0 , 100.0034 + 779, -0.02707883 , -50.0 , 114.6107 + 780, 24.7223 , -50.0 , 129.145 + 781, -0.07242509 , -50.0 , 143.9754 + 782, -23.78324 , -50.0 , 159.1944 + 783, -24.7903 , -50.0 , 129.1356 + 784, 23.70813 , -50.0 , 159.1897 + 785, -0.04064292 , -50.0 , 172.1079 + 786, 0.07782174 , -50.0 , 85.36618 + 787, 24.83621 , -50.0 , 70.84835 + 788, 24.93575 , -50.0 , 99.99269 + 789, -23.6763 , -50.0 , 40.79128 + 790, 0.08447849 , -50.0 , 56.03104 + 791, -24.73715 , -50.0 , 70.86278 + 792, 23.7887 , -50.0 , 40.80619 + 793, 0.05147156 , -50.0 , 27.8916 + 859, -2.844261 , -50.0 , 274.4003 + 860, 3.067045 , -50.0 , 225.8414 + 861, 17.79876 , -50.0 , 254.1803 + 862, 26.13106 , -50.0 , 278.7815 + 863, -17.83471 , -50.0 , 245.9473 + 864, -26.11139 , -50.0 , 221.3009 + 892, 24.4003 , -2.844261 , 0.0 + 893, -24.15862 , 3.067045 , 0.0 + 894, 4.180344 , 17.79876 , 0.0 + 895, 28.78155 , 26.13106 , 0.0 + 896, -4.052692 , -17.83471 , 0.0 + 897, -28.69907 , -26.11139 , 0.0 + 925, -24.15862 , 3.067045 , 300.0 + 926, 24.4003 , -2.844261 , 300.0 + 927, 4.180344 , 17.79876 , 300.0 + 928, 28.78155 , 26.13106 , 300.0 + 929, -4.052692 , -17.83471 , 300.0 + 930, -28.69907 , -26.11139 , 300.0 + 958, 25.03691 , 50.0 , 134.9316 + 959, 0.06483927 , 50.0 , 149.9933 + 960, -24.94041 , 50.0 , 165.0632 + 961, 0.1047811 , 50.0 , 180.2299 + 962, 25.04272 , 50.0 , 195.2832 + 963, 25.04349 , 50.0 , 165.0721 + 964, -24.65016 , 50.0 , 225.6887 + 965, 0.126205 , 50.0 , 210.5384 + 966, -24.8927 , 50.0 , 195.2433 + 967, 0.1264049 , 50.0 , 241.2084 + 968, 23.81343 , 50.0 , 257.0779 + 969, 24.83341 , 50.0 , 225.6572 + 970, -23.68546 , 50.0 , 257.0829 + 971, 0.06618169 , 50.0 , 270.7359 + 972, 0.103386 , 50.0 , 119.7668 + 973, -24.87592 , 50.0 , 104.7186 + 974, -24.94078 , 50.0 , 134.9249 + 975, 24.88604 , 50.0 , 74.31196 + 976, 0.1278253 , 50.0 , 89.46007 + 977, 25.02198 , 50.0 , 104.7554 + 978, 0.1271092 , 50.0 , 58.79149 + 979, -23.68474 , 50.0 , 42.91735 + 980, -24.70657 , 50.0 , 74.34253 + 981, 23.81315 , 50.0 , 42.9215 + 982, 0.06645916 , 50.0 , 29.26398 + 1081, 50.0 , -24.04267 , 219.9838 + 1082, 50.0 , 0.7125605 , 206.5824 + 1083, 50.0 , 25.24923 , 193.2424 + 1084, 50.0 , 25.9456 , 223.5992 + 1085, 50.0 , 2.424309 , 237.9667 + 1086, 50.0 , -22.01633 , 253.9517 + 1087, 50.0 , 27.48433 , 253.4867 + 1088, 50.0 , 29.9289 , 278.4534 + 1089, 50.0 , 6.500372 , 270.3959 + 1090, 50.0 , 24.81494 , 73.40613 + 1091, 50.0 , 0.03667517 , 87.30159 + 1092, 50.0 , -24.95728 , 101.01 + 1093, 50.0 , -23.71272 , 41.31017 + 1094, 50.0 , 0.04684983 , 57.37447 + 1095, 50.0 , -24.76602 , 71.67199 + 1096, 50.0 , 23.76693 , 42.36515 + 1097, 50.0 , 0.04084962 , 28.56608 + 1098, 50.0 , 25.05472 , 163.2151 + 1099, 50.0 , 25.00519 , 133.3589 + 1100, 50.0 , 24.96884 , 103.5075 + 1101, 50.0 , 0.003502569 , 116.9813 + 1102, 50.0 , -24.99949 , 130.165 + 1103, 50.0 , 0.01926796 , 146.5365 + 1104, 50.0 , -24.95747 , 159.3337 + 1105, 50.0 , 0.2317527 , 176.3 + 1106, 50.0 , -24.66696 , 188.9299 + 1208, -50.0 , -50.0 , 171.4286 + 1209, -50.0 , -50.0 , 142.8571 + 1210, -50.0 , -50.0 , 114.2857 + 1211, -50.0 , -50.0 , 85.71429 + 1212, -50.0 , -50.0 , 57.14286 + 1213, -50.0 , -50.0 , 28.57143 + 1221, -50.0 , -50.0 , 266.6667 + 1222, -50.0 , -50.0 , 233.3333 + 1226, -50.0 , 16.66667 , 0.0 + 1227, -50.0 , -16.66667 , 0.0 + 1231, -50.0 , 16.66667 , 300.0 + 1232, -50.0 , -16.66667 , 300.0 + 1236, -50.0 , 50.0 , 270.0 + 1237, -50.0 , 50.0 , 240.0 + 1238, -50.0 , 50.0 , 210.0 + 1239, -50.0 , 50.0 , 180.0 + 1240, -50.0 , 50.0 , 150.0 + 1241, -50.0 , 50.0 , 120.0 + 1242, -50.0 , 50.0 , 90.0 + 1243, -50.0 , 50.0 , 60.0 + 1244, -50.0 , 50.0 , 30.0 + 1255, -16.66667 , -50.0 , 0.0 + 1256, 16.66667 , -50.0 , 0.0 + 1260, -16.66667 , -50.0 , 200.0 + 1261, 16.66667 , -50.0 , 200.0 + 1265, -16.66667 , -50.0 , 300.0 + 1266, 16.66667 , -50.0 , 300.0 + 1270, 16.66667 , 50.0 , 0.0 + 1271, -16.66667 , 50.0 , 0.0 + 1275, 16.66667 , 50.0 , 300.0 + 1276, -16.66667 , 50.0 , 300.0 + 1280, 50.0 , -50.0 , 171.4286 + 1281, 50.0 , -50.0 , 142.8571 + 1282, 50.0 , -50.0 , 114.2857 + 1283, 50.0 , -50.0 , 85.71429 + 1284, 50.0 , -50.0 , 57.14286 + 1285, 50.0 , -50.0 , 28.57143 + 1293, 50.0 , -50.0 , 266.6667 + 1294, 50.0 , -50.0 , 233.3333 + 1298, 50.0 , -16.66667 , 0.0 + 1299, 50.0 , 16.66667 , 0.0 + 1303, 50.0 , -16.66667 , 300.0 + 1304, 50.0 , 16.66667 , 300.0 + 1308, 50.0 , 50.0 , 270.0 + 1309, 50.0 , 50.0 , 240.0 + 1310, 50.0 , 50.0 , 210.0 + 1311, 50.0 , 50.0 , 180.0 + 1312, 50.0 , 50.0 , 150.0 + 1313, 50.0 , 50.0 , 120.0 + 1314, 50.0 , 50.0 , 90.0 + 1315, 50.0 , 50.0 , 60.0 + 1316, 50.0 , 50.0 , 30.0 + 1327, -50.0 , -50.0 , 200.0 + 1328, -50.0 , -50.0 , 0.0 + 1329, -50.0 , -50.0 , 300.0 + 1330, -50.0 , 50.0 , 0.0 + 1331, -50.0 , 50.0 , 300.0 + 1332, 50.0 , -50.0 , 0.0 + 1333, 50.0 , -50.0 , 200.0 + 1334, 50.0 , -50.0 , 300.0 + 1335, 50.0 , 50.0 , 0.0 + 1336, 50.0 , 50.0 , 300.0 +**HWCOLOR COMP 1 3 +*ELEMENT,TYPE=C3D4,ELSET=beam + 642, 654, 930, 659, 1232 + 665, 1221, 1232, 1329, 930 + 1, 1, 1086, 2, 3 + 2, 4, 5, 6, 967 + 3, 4, 5, 967, 971 + 4, 7, 8, 9, 10 + 5, 7, 11, 12, 13 + 6, 7, 669, 14, 15 + 7, 7, 14, 16, 13 + 8, 17, 18, 1092, 1101 + 9, 17, 1091, 19, 20 + 10, 17, 21, 1092, 18 + 11, 12, 22, 672, 660 + 12, 12, 661, 23, 24 + 13, 12, 22, 660, 24 + 14, 25, 1103, 26, 27 + 15, 25, 8, 10, 28 + 16, 25, 11, 18, 10 + 17, 25, 16, 26, 13 + 18, 19, 1094, 29, 1090 + 19, 19, 20, 30, 31 + 20, 19, 1095, 21, 32 + 21, 19, 1094, 33, 29 + 22, 19, 30, 33, 32 + 23, 34, 35, 30, 33 + 24, 35, 896, 34, 36 + 25, 35, 978, 29, 981 + 26, 35, 981, 1270, 982 + 27, 892, 35, 896, 33 + 28, 676, 653, 1, 652 + 29, 863, 3, 1, 37 + 30, 35, 896, 36, 894 + 31, 35, 894, 982, 1270 + 32, 12, 23, 38, 24 + 33, 12, 31, 11, 38 + 34, 23, 36, 30, 39 + 35, 23, 40, 41, 662 + 36, 23, 30, 31, 38 + 37, 4, 925, 37, 927 + 38, 37, 654, 930, 659 + 39, 37, 5, 3, 2 + 40, 859, 3, 863, 37 + 41, 654, 859, 1221, 863 + 42, 930, 859, 1221, 654 + 43, 30, 32, 38, 39 + 44, 37, 929, 926, 927 + 45, 36, 897, 893, 896 + 46, 968, 1308, 1087, 1309 + 47, 34, 896, 793, 1255 + 48, 39, 32, 38, 790 + 49, 3, 1303, 1089, 926 + 50, 34, 30, 36, 39 + 51, 34, 39, 789, 790 + 52, 5, 3, 1087, 1089 + 53, 5, 926, 37, 3 + 54, 926, 5, 1089, 3 + 55, 5, 1087, 3, 2 + 56, 5, 968, 2, 967 + 57, 42, 960, 966, 961 + 58, 5, 1275, 971, 927 + 59, 5, 971, 1275, 968 + 60, 43, 964, 6, 44 + 61, 43, 42, 44, 45 + 62, 43, 966, 1238, 964 + 63, 43, 966, 651, 1238 + 64, 6, 652, 656, 43 + 65, 45, 14, 674, 42 + 66, 42, 16, 14, 960 + 67, 42, 45, 651, 675 + 68, 42, 960, 670, 1239 + 69, 43, 42, 45, 651 + 70, 42, 960, 1239, 966 + 71, 17, 25, 18, 1101 + 72, 20, 31, 976, 30 + 73, 17, 1100, 46, 1101 + 74, 20, 46, 977, 31 + 75, 20, 31, 977, 976 + 76, 20, 977, 1314, 975 + 77, 20, 975, 1314, 1090 + 78, 26, 1105, 27, 1103 + 79, 26, 16, 959, 13 + 80, 26, 47, 1105, 1098 + 81, 26, 16, 963, 959 + 82, 26, 47, 963, 16 + 83, 26, 958, 963, 1312 + 84, 26, 958, 1312, 1099 + 85, 14, 959, 16, 13 + 86, 674, 42, 675, 45 + 87, 42, 14, 674, 670 + 88, 14, 960, 16, 959 + 89, 42, 45, 14, 16 + 90, 14, 974, 1240, 960 + 91, 14, 1240, 974, 669 + 92, 17, 11, 25, 13 + 93, 17, 11, 21, 18 + 94, 25, 26, 46, 13 + 95, 25, 27, 8, 28 + 96, 26, 46, 958, 1099 + 97, 46, 13, 958, 972 + 98, 46, 972, 31, 13 + 99, 46, 958, 1313, 977 + 100, 46, 977, 1313, 1100 + 101, 19, 32, 1094, 1095 + 102, 21, 786, 11, 38 + 103, 21, 787, 32, 1095 + 104, 21, 38, 787, 786 + 105, 21, 788, 787, 1283 + 106, 21, 788, 1283, 1092 + 107, 17, 31, 38, 11 + 108, 12, 15, 31, 40 + 109, 12, 15, 40, 668 + 110, 31, 15, 972, 973 + 111, 31, 15, 973, 40 + 112, 46, 977, 31, 972 + 113, 20, 30, 976, 975 + 114, 7, 672, 12, 22 + 115, 7, 9, 674, 671 + 116, 17, 38, 19, 21 + 117, 31, 976, 40, 973 + 118, 14, 15, 974, 13 + 119, 15, 13, 972, 974 + 120, 4, 1276, 1236, 658 + 121, 15, 1241, 974, 973 + 122, 15, 1241, 973, 668 + 123, 34, 1255, 897, 896 + 124, 1327, 676, 1260, 1208 + 125, 48, 45, 675, 676 + 126, 48, 1261, 1260, 785 + 127, 48, 782, 8, 785 + 128, 48, 9, 8, 782 + 129, 48, 673, 782, 1208 + 130, 27, 8, 1, 49 + 131, 1089, 1085, 3, 1087 + 132, 1261, 1106, 1280, 27 + 133, 27, 8, 784, 785 + 134, 27, 1261, 8, 785 + 135, 667, 893, 1227, 1226 + 136, 27, 784, 28, 1104 + 137, 27, 1104, 1106, 1280 + 138, 23, 41, 40, 30 + 139, 12, 11, 31, 13 + 140, 12, 38, 23, 31 + 141, 31, 15, 13, 972 + 142, 40, 980, 662, 1242 + 143, 40, 976, 30, 980 + 144, 40, 1242, 973, 980 + 145, 36, 666, 1226, 667 + 146, 36, 1271, 894, 893 + 147, 34, 664, 1213, 789 + 148, 2, 1085, 1082, 1084 + 149, 36, 34, 35, 30 + 150, 22, 10, 779, 11 + 151, 22, 778, 11, 779 + 152, 22, 778, 24, 11 + 153, 22, 1210, 778, 783 + 154, 22, 660, 1210, 671 + 155, 48, 8, 45, 1 + 156, 45, 9, 674, 7 + 157, 12, 24, 660, 661 + 158, 9, 781, 10, 8 + 159, 9, 10, 781, 783 + 160, 22, 9, 783, 10 + 161, 9, 783, 1209, 671 + 162, 9, 673, 671, 1209 + 163, 18, 779, 10, 11 + 164, 17, 19, 38, 31 + 165, 21, 18, 11, 788 + 166, 18, 11, 788, 779 + 167, 18, 1282, 788, 1092 + 168, 18, 1092, 1102, 1282 + 169, 1, 1261, 1260, 48 + 170, 42, 14, 670, 960 + 171, 42, 47, 16, 961 + 172, 43, 651, 45, 652 + 173, 48, 676, 1260, 1 + 174, 45, 47, 44, 49 + 175, 7, 9, 8, 45 + 176, 48, 9, 45, 8 + 177, 20, 975, 976, 977 + 178, 17, 46, 31, 13 + 179, 25, 13, 10, 16 + 180, 10, 28, 781, 780 + 181, 18, 10, 780, 28 + 182, 48, 9, 782, 673 + 183, 22, 779, 10, 783 + 184, 9, 8, 782, 781 + 185, 25, 26, 1099, 46 + 186, 25, 18, 1102, 28 + 187, 18, 780, 10, 779 + 188, 7, 14, 13, 15 + 189, 27, 1104, 28, 1103 + 190, 27, 47, 8, 49 + 191, 28, 1102, 1104, 1281 + 192, 28, 8, 781, 784 + 193, 28, 780, 1102, 1281 + 194, 47, 1311, 963, 962 + 195, 47, 45, 8, 49 + 196, 16, 47, 45, 8 + 197, 27, 47, 49, 1105 + 198, 667, 897, 664, 1227 + 199, 26, 27, 8, 25 + 200, 47, 27, 8, 26 + 201, 47, 49, 1105, 1083 + 202, 47, 963, 1311, 1098 + 203, 23, 663, 39, 665 + 204, 39, 1212, 789, 791 + 205, 39, 38, 791, 790 + 206, 32, 34, 39, 30 + 207, 23, 41, 663, 662 + 208, 23, 40, 31, 30 + 209, 7, 22, 11, 10 + 210, 12, 22, 24, 11 + 211, 24, 11, 786, 38 + 212, 24, 791, 38, 786 + 213, 24, 778, 1211, 791 + 214, 24, 778, 660, 1211 + 215, 33, 32, 1093, 1094 + 216, 33, 29, 1094, 1096 + 217, 35, 29, 33, 1096 + 218, 35, 33, 1097, 1096 + 219, 35, 36, 30, 41 + 220, 35, 41, 30, 978 + 221, 4, 1276, 1231, 925 + 222, 5, 927, 926, 928 + 223, 3, 929, 859, 862 + 224, 1086, 2, 3, 1085 + 225, 1085, 3, 1087, 2 + 226, 3, 1086, 1089, 1303 + 227, 3, 1086, 1303, 862 + 228, 33, 1097, 1298, 892 + 229, 36, 39, 663, 664 + 230, 35, 41, 978, 982 + 231, 1226, 1271, 1244, 666 + 232, 35, 36, 41, 982 + 233, 40, 41, 980, 30 + 234, 12, 24, 38, 11 + 235, 41, 30, 978, 980 + 236, 23, 40, 662, 661 + 237, 41, 980, 979, 1243 + 238, 41, 1243, 979, 666 + 239, 861, 863, 3, 1 + 240, 969, 1310, 1309, 1084 + 241, 1081, 1086, 2, 1 + 242, 1, 1106, 49, 1081 + 243, 1, 44, 2, 49 + 244, 37, 3, 926, 929 + 245, 2, 49, 44, 969 + 246, 27, 1261, 1, 8 + 247, 37, 6, 1, 863 + 248, 49, 1084, 969, 2 + 249, 1086, 861, 1, 1081 + 250, 2, 49, 1082, 1081 + 251, 1, 49, 2, 1081 + 252, 676, 45, 1, 48 + 253, 6, 652, 45, 1 + 254, 43, 6, 45, 44 + 255, 1294, 1261, 1081, 1333 + 256, 2, 1082, 49, 1084 + 257, 1298, 892, 1097, 1299 + 258, 1, 6, 2, 44 + 259, 27, 49, 1106, 1105 + 260, 4, 37, 659, 6 + 261, 49, 1081, 1106, 1082 + 262, 17, 20, 1100, 1091 + 263, 33, 792, 32, 793 + 264, 653, 863, 1, 6 + 265, 34, 793, 32, 790 + 266, 32, 1284, 1095, 787 + 267, 32, 38, 790, 787 + 268, 32, 792, 1284, 787 + 269, 44, 961, 962, 47 + 270, 26, 46, 13, 958 + 271, 14, 13, 974, 959 + 272, 26, 958, 13, 959 + 273, 24, 786, 11, 778 + 274, 22, 779, 783, 778 + 275, 35, 29, 978, 30 + 276, 1097, 35, 892, 33 + 277, 19, 29, 20, 1090 + 278, 29, 978, 30, 975 + 279, 20, 29, 30, 975 + 280, 12, 13, 31, 15 + 281, 29, 1315, 975, 1090 + 282, 29, 1315, 1090, 1096 + 283, 37, 654, 6, 863 + 284, 4, 6, 659, 657 + 285, 1, 3, 2, 37 + 286, 7, 10, 16, 8 + 287, 48, 674, 673, 675 + 288, 27, 1105, 1106, 1104 + 289, 25, 1102, 18, 1101 + 290, 7, 672, 22, 671 + 291, 1261, 8, 48, 1 + 292, 965, 44, 961, 962 + 293, 25, 1102, 1101, 1103 + 294, 7, 672, 671, 674 + 295, 19, 30, 20, 29 + 296, 11, 779, 786, 788 + 297, 24, 786, 778, 791 + 298, 30, 32, 34, 33 + 299, 3, 926, 862, 1303 + 300, 7, 671, 22, 9 + 301, 7, 22, 10, 9 + 302, 14, 960, 1240, 670 + 303, 42, 966, 1239, 651 + 304, 27, 1106, 49, 1 + 305, 42, 44, 45, 47 + 306, 18, 1092, 1101, 1102 + 307, 25, 18, 28, 10 + 308, 47, 49, 1083, 962 + 309, 1303, 1293, 1334, 862 + 310, 37, 6, 654, 659 + 311, 37, 929, 925, 930 + 312, 1221, 1265, 930, 1329 + 313, 46, 977, 972, 958 + 314, 31, 973, 972, 976 + 315, 25, 1103, 28, 1102 + 316, 10, 779, 780, 781 + 317, 26, 1312, 963, 1098 + 318, 670, 42, 675, 674 + 319, 12, 668, 40, 661 + 320, 25, 1099, 26, 1103 + 321, 47, 962, 1083, 1311 + 322, 25, 46, 1099, 1101 + 323, 7, 15, 12, 672 + 324, 7, 672, 674, 669 + 325, 14, 959, 974, 960 + 326, 15, 973, 40, 668 + 327, 7, 674, 14, 669 + 328, 966, 43, 42, 965 + 329, 14, 15, 669, 974 + 330, 15, 669, 974, 1241 + 331, 20, 977, 46, 1100 + 332, 30, 976, 978, 980 + 333, 17, 21, 11, 38 + 334, 17, 20, 46, 1100 + 335, 1244, 979, 1243, 666 + 336, 1330, 1271, 1244, 1226 + 337, 40, 980, 973, 976 + 338, 1316, 1096, 895, 1299 + 339, 655, 6, 652, 656 + 340, 43, 42, 651, 966 + 341, 3, 861, 1, 1086 + 342, 43, 656, 6, 964 + 343, 969, 49, 44, 962 + 344, 26, 959, 963, 958 + 345, 19, 32, 33, 1094 + 346, 25, 10, 13, 11 + 347, 39, 665, 663, 664 + 348, 33, 1256, 896, 892 + 349, 46, 1313, 958, 1099 + 350, 964, 657, 970, 1237 + 351, 3, 1086, 1085, 1089 + 352, 30, 976, 975, 978 + 353, 40, 662, 661, 668 + 354, 41, 980, 1243, 662 + 355, 34, 793, 790, 789 + 356, 39, 24, 791, 38 + 357, 24, 665, 791, 1211 + 358, 778, 1211, 1210, 660 + 359, 21, 1095, 1091, 1092 + 360, 788, 1282, 1283, 1092 + 361, 21, 1095, 1092, 1283 + 362, 1213, 664, 897, 1227 + 363, 789, 1213, 34, 1255 + 364, 1212, 665, 1211, 791 + 365, 32, 1093, 1284, 792 + 366, 38, 786, 790, 787 + 367, 1271, 36, 666, 1226 + 368, 980, 1243, 662, 1242 + 369, 1315, 975, 1090, 1314 + 370, 969, 968, 1087, 1309 + 371, 962, 1084, 969, 49 + 372, 966, 1238, 1239, 651 + 373, 44, 965, 967, 969 + 374, 17, 1101, 1091, 1100 + 375, 7, 12, 15, 13 + 376, 33, 32, 792, 1093 + 377, 29, 975, 1315, 981 + 378, 20, 1091, 1090, 1100 + 379, 46, 1100, 1099, 1101 + 380, 39, 664, 1212, 665 + 381, 24, 661, 665, 660 + 382, 1284, 1095, 787, 1283 + 383, 19, 1091, 1095, 1094 + 384, 20, 1100, 1090, 1314 + 385, 29, 1315, 1096, 981 + 386, 40, 1242, 662, 668 + 387, 22, 660, 778, 1210 + 388, 24, 1211, 660, 665 + 389, 39, 789, 1212, 664 + 390, 12, 661, 660, 672 + 391, 40, 668, 973, 1242 + 392, 23, 38, 39, 30 + 393, 32, 792, 790, 793 + 394, 20, 1314, 977, 1100 + 395, 12, 661, 672, 668 + 396, 973, 1241, 1242, 668 + 397, 21, 1283, 787, 1095 + 398, 18, 780, 788, 1282 + 399, 31, 30, 40, 976 + 400, 17, 1092, 21, 1091 + 401, 31, 976, 972, 977 + 402, 17, 20, 31, 46 + 403, 23, 665, 661, 663 + 404, 7, 16, 14, 45 + 405, 7, 13, 10, 11 + 406, 18, 780, 1282, 1102 + 407, 16, 959, 961, 963 + 408, 33, 1097, 1096, 1094 + 409, 16, 961, 959, 960 + 410, 35, 981, 982, 978 + 411, 1086, 1294, 861, 1081 + 412, 27, 8, 28, 784 + 413, 28, 780, 784, 781 + 414, 23, 663, 36, 39 + 415, 25, 1103, 1101, 1099 + 416, 42, 45, 16, 47 + 417, 974, 1240, 1241, 669 + 418, 962, 969, 1084, 1310 + 419, 14, 669, 670, 1240 + 420, 17, 1091, 1101, 1092 + 421, 977, 1100, 1314, 1313 + 422, 22, 783, 671, 1210 + 423, 39, 24, 665, 791 + 424, 18, 780, 779, 788 + 425, 9, 1209, 783, 782 + 426, 19, 29, 33, 30 + 427, 23, 661, 662, 663 + 428, 17, 21, 19, 1091 + 429, 21, 788, 11, 786 + 430, 32, 1095, 1284, 1093 + 431, 21, 38, 32, 787 + 432, 963, 1311, 1098, 1312 + 433, 25, 26, 16, 8 + 434, 46, 1100, 1313, 1099 + 435, 15, 668, 669, 1241 + 436, 26, 1098, 1099, 1312 + 437, 27, 1103, 1105, 1104 + 438, 19, 21, 1095, 1091 + 439, 960, 1240, 670, 1239 + 440, 47, 1098, 1311, 1083 + 441, 26, 1103, 1098, 1105 + 442, 48, 676, 675, 673 + 443, 12, 15, 668, 672 + 444, 958, 1099, 1313, 1312 + 445, 34, 789, 39, 664 + 446, 7, 45, 8, 16 + 447, 9, 1209, 782, 673 + 448, 23, 41, 36, 663 + 449, 12, 23, 40, 31 + 450, 6, 964, 657, 970 + 451, 15, 668, 672, 669 + 452, 8, 26, 16, 47 + 453, 14, 670, 669, 674 + 454, 10, 783, 779, 781 + 455, 28, 1281, 1104, 784 + 456, 28, 780, 1281, 784 + 457, 22, 783, 9, 671 + 458, 19, 38, 31, 30 + 459, 47, 961, 962, 963 + 460, 26, 1103, 1099, 1098 + 461, 17, 13, 25, 46 + 462, 48, 45, 9, 674 + 463, 26, 47, 27, 1105 + 464, 25, 16, 10, 8 + 465, 18, 1102, 28, 780 + 466, 1106, 1261, 1, 27 + 467, 23, 24, 661, 665 + 468, 13, 972, 959, 958 + 469, 7, 669, 15, 672 + 470, 1209, 1208, 782, 673 + 471, 17, 25, 1101, 46 + 472, 32, 787, 790, 792 + 473, 6, 653, 654, 655 + 474, 25, 27, 28, 1103 + 475, 9, 674, 671, 673 + 476, 1281, 1280, 1104, 784 + 477, 7, 10, 13, 16 + 478, 11, 779, 778, 786 + 479, 48, 1260, 1208, 782 + 480, 48, 673, 1208, 676 + 481, 37, 5, 2, 6 + 482, 654, 1222, 653, 863 + 483, 965, 43, 42, 44 + 484, 28, 1103, 1104, 1102 + 485, 783, 671, 1210, 1209 + 486, 1081, 1333, 1261, 1106 + 487, 1294, 860, 861, 1081 + 488, 1, 860, 1081, 861 + 489, 29, 1090, 1094, 1096 + 490, 47, 1083, 1105, 1098 + 491, 47, 963, 16, 961 + 492, 27, 1104, 1280, 784 + 493, 48, 9, 673, 674 + 494, 19, 30, 32, 38 + 495, 47, 49, 962, 44 + 496, 1081, 1261, 1, 1106 + 497, 27, 1261, 784, 1280 + 498, 962, 1310, 1084, 1083 + 499, 42, 670, 675, 651 + 500, 653, 676, 1260, 1327 + 501, 6, 653, 655, 652 + 502, 23, 665, 39, 24 + 503, 19, 1090, 20, 1091 + 504, 652, 653, 1, 6 + 505, 45, 652, 651, 675 + 506, 15, 973, 974, 972 + 507, 49, 1084, 1082, 1083 + 508, 41, 1243, 666, 662 + 509, 35, 30, 33, 29 + 510, 36, 41, 666, 663 + 511, 40, 980, 41, 662 + 512, 962, 1310, 1083, 1311 + 513, 35, 1316, 1270, 981 + 514, 859, 37, 654, 930 + 515, 35, 29, 1096, 981 + 516, 895, 1097, 35, 892 + 517, 1316, 981, 1096, 1315 + 518, 1236, 4, 657, 970 + 519, 1276, 658, 1231, 1331 + 520, 4, 927, 1276, 925 + 521, 4, 1236, 1276, 970 + 522, 4, 37, 6, 5 + 523, 667, 1227, 893, 897 + 524, 42, 47, 961, 44 + 525, 43, 45, 6, 652 + 526, 43, 656, 964, 1238 + 527, 6, 964, 970, 967 + 528, 4, 1276, 971, 970 + 529, 36, 1271, 979, 982 + 530, 36, 663, 666, 667 + 531, 44, 965, 969, 962 + 532, 5, 967, 971, 968 + 533, 4, 970, 971, 967 + 534, 4, 5, 971, 927 + 535, 6, 655, 659, 657 + 536, 2, 969, 1087, 1084 + 537, 10, 28, 8, 781 + 538, 654, 1232, 1221, 930 + 539, 1276, 1275, 927, 971 + 540, 4, 1276, 927, 971 + 541, 928, 1275, 1088, 5 + 542, 968, 1308, 1275, 1088 + 543, 43, 966, 964, 965 + 544, 44, 965, 961, 42 + 545, 895, 35, 1316, 1270 + 546, 19, 1090, 1091, 1094 + 547, 1096, 1097, 895, 1299 + 548, 35, 1316, 981, 1096 + 549, 35, 36, 982, 894 + 550, 38, 786, 791, 790 + 551, 13, 974, 959, 972 + 552, 1, 652, 45, 676 + 553, 21, 786, 787, 788 + 554, 36, 893, 1226, 1271 + 555, 22, 660, 671, 672 + 556, 37, 929, 927, 925 + 557, 964, 1238, 656, 1237 + 558, 37, 1, 6, 2 + 559, 1086, 1293, 861, 1294 + 560, 2, 1085, 1081, 1082 + 561, 23, 30, 36, 41 + 562, 27, 785, 784, 1261 + 563, 45, 1, 49, 44 + 564, 6, 654, 659, 655 + 565, 48, 785, 1260, 782 + 566, 8, 785, 781, 784 + 567, 43, 1238, 651, 656 + 568, 19, 21, 38, 32 + 569, 33, 1256, 892, 1298 + 570, 36, 666, 41, 979 + 571, 1, 1261, 860, 1260 + 572, 12, 40, 23, 661 + 573, 35, 895, 1316, 1096 + 574, 33, 792, 1256, 1285 + 575, 33, 1298, 1093, 1285 + 576, 32, 1095, 1093, 1094 + 577, 1284, 792, 1093, 1285 + 578, 1213, 1212, 789, 664 + 579, 33, 1285, 1256, 1298 + 580, 33, 1094, 1093, 1097 + 581, 34, 33, 32, 793 + 582, 667, 897, 36, 664 + 583, 41, 979, 980, 978 + 584, 34, 33, 793, 896 + 585, 7, 12, 11, 22 + 586, 17, 18, 25, 11 + 587, 7, 45, 14, 674 + 588, 35, 33, 34, 896 + 589, 35, 1270, 895, 894 + 590, 667, 897, 893, 36 + 591, 43, 651, 652, 656 + 592, 6, 654, 653, 863 + 593, 653, 863, 1222, 864 + 594, 8, 782, 781, 785 + 595, 35, 1096, 1097, 895 + 596, 44, 965, 964, 967 + 597, 22, 660, 24, 778 + 598, 36, 667, 664, 663 + 599, 29, 975, 981, 978 + 600, 33, 1093, 1298, 1097 + 601, 1256, 793, 896, 1255 + 602, 23, 39, 38, 24 + 603, 790, 34, 39, 32 + 604, 34, 897, 1213, 664 + 605, 1213, 34, 1255, 897 + 606, 48, 675, 45, 674 + 607, 34, 36, 664, 39 + 608, 1097, 895, 1299, 892 + 609, 43, 44, 965, 964 + 610, 897, 36, 664, 34 + 611, 1, 1261, 1081, 860 + 612, 42, 965, 961, 966 + 613, 1083, 962, 49, 1084 + 614, 4, 927, 37, 5 + 615, 2, 44, 6, 967 + 616, 666, 1271, 1244, 979 + 617, 26, 963, 47, 1098 + 618, 979, 36, 666, 1271 + 619, 42, 651, 1239, 670 + 620, 969, 1309, 1087, 1084 + 621, 17, 19, 31, 20 + 622, 1, 863, 864, 860 + 623, 37, 859, 3, 929 + 624, 964, 1237, 656, 657 + 625, 20, 1090, 29, 975 + 626, 41, 662, 666, 663 + 627, 676, 1260, 1208, 48 + 628, 39, 665, 1212, 791 + 629, 49, 1105, 1083, 1082 + 630, 1, 676, 1260, 653 + 631, 17, 13, 31, 11 + 632, 33, 1256, 792, 793 + 633, 36, 896, 893, 894 + 634, 33, 792, 1285, 1093 + 635, 34, 793, 789, 1255 + 636, 1265, 929, 1266, 859 + 637, 3, 926, 929, 862 + 638, 929, 862, 926, 1266 + 639, 2, 1085, 1084, 1087 + 640, 37, 859, 929, 930 + 641, 42, 961, 16, 960 + 643, 1275, 1088, 5, 968 + 644, 1275, 1336, 928, 1088 + 645, 5, 927, 928, 1275 + 646, 2, 44, 967, 969 + 647, 4, 659, 658, 657 + 648, 45, 1, 8, 49 + 649, 41, 982, 979, 978 + 650, 36, 41, 982, 979 + 651, 35, 894, 895, 892 + 652, 1316, 1335, 895, 1270 + 653, 4, 967, 6, 970 + 654, 1086, 1303, 862, 1293 + 655, 927, 5, 926, 37 + 656, 1316, 1335, 1299, 895 + 657, 1303, 926, 1304, 1089 + 658, 4, 659, 1231, 658 + 659, 864, 653, 863, 1 + 660, 654, 1222, 863, 1221 + 661, 659, 925, 1231, 1232 + 662, 653, 864, 1327, 1260 + 663, 3, 861, 1086, 862 + 664, 2, 1087, 969, 968 + 666, 3, 862, 859, 861 + 667, 4, 6, 657, 970 + 668, 39, 790, 791, 789 + 669, 780, 1281, 1282, 1102 + 670, 21, 18, 788, 1092 + 671, 49, 1105, 1082, 1106 + 672, 5, 2, 968, 1087 + 673, 33, 896, 1256, 793 + 674, 928, 926, 5, 1089 + 675, 4, 659, 37, 925 + 676, 863, 37, 654, 859 + 677, 1221, 930, 1265, 859 + 678, 1265, 859, 930, 929 + 679, 1213, 1227, 897, 1328 + 680, 896, 34, 36, 897 + 681, 1332, 1256, 1285, 1298 + 682, 5, 2, 6, 967 + 683, 1213, 1328, 897, 1255 + 684, 1303, 926, 862, 1266 + 685, 1, 864, 1260, 860 + 686, 1303, 1334, 1266, 862 + 687, 36, 894, 1271, 982 + 688, 1333, 1261, 1106, 1280 + 689, 45, 652, 675, 676 + 690, 37, 659, 930, 925 + 691, 6, 655, 657, 656 + 692, 968, 1088, 1087, 1308 + 693, 4, 1231, 1276, 658 + 694, 659, 925, 1232, 930 + 695, 4, 659, 925, 1231 + 696, 2, 968, 969, 967 + 697, 45, 1, 44, 6 + 698, 1, 863, 860, 861 + 699, 1089, 1304, 1088, 928 + 700, 929, 1266, 859, 862 + 701, 1294, 860, 1081, 1261 + 702, 9, 783, 781, 782 + 703, 5, 928, 1089, 1088 + 704, 6, 967, 44, 964 + 705, 861, 859, 3, 863 + 706, 1086, 861, 1293, 862 + 707, 1, 653, 1260, 864 + 708, 1271, 1270, 982, 894 + 709, 894, 892, 35, 896 + 710, 36, 667, 1226, 893 + 711, 1275, 1088, 1308, 1336 + 712, 653, 1222, 1327, 864 + 713, 1276, 1331, 1236, 658 + 714, 970, 1237, 657, 1236 + 715, 1085, 1081, 1086, 2 + 716, 1304, 1336, 1088, 928 + 717, 1089, 928, 926, 1304 + 718, 5, 1088, 1089, 1087 + 719, 658, 1236, 4, 657 + 720, 6, 964, 656, 657 + 721, 785, 8, 48, 1261 + 722, 5, 1087, 968, 1088 +*ELSET, ELSET=EALL +beam, +*NSET, NSET=NALL +nset_csys0, +*SURFACE, NAME = load, TYPE = ELEMENT + 41, S3 + 255, S2 + 309, S3 + 312, S2 + 487, S1 + 559, S3 + 571, S3 + 593, S3 + 622, S3 + 636, S4 + 660, S3 + 662, S3 + 666, S3 + 677, S4 + 685, S3 + 686, S3 + 698, S3 + 700, S3 + 701, S2 + 705, S2 + 706, S3 + 712, S3 +*SURFACE, NAME = fixed, TYPE = ELEMENT + 45, S3 + 123, S3 + 135, S3 + 146, S3 + 257, S2 + 336, S2 + 348, S3 + 523, S3 + 554, S3 + 569, S3 + 589, S3 + 601, S4 + 608, S3 + 633, S3 + 651, S3 + 652, S3 + 656, S3 + 679, S3 + 681, S2 + 683, S3 + 708, S2 + 709, S2 +**HWNAME LOADCOL 1 HM_IC_Load_Cols_1 +**HWCOLOR LOADCOL 1 3 +*INITIAL CONDITIONS, TYPE = TEMPERATURE +NALL,2.000000E +*****