save$bdn = {{True, "g"}, {False, "NAB"}} save$bodypos = {{xg, yg, zg, pitchg, yawg, rollg}, {xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB}} save$FileName = "Sub\\NAB" save$initpos = {{0, 0, 0}, {0, -0.757, 0}} save$mass = {{0, {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}}, {36.66, {{1.92, 0, 0}, {0, 2.6, 0}, {0, 0, 4.44}}}} save$mat = {{"Maraging Steel", 8., 195, 0.3, 1., -3, "Virgo Data"}, {"C-70 Steel", 7.8, 200, 0.3, 3., -4, "LIGO-G050087-00"}, {"Tungsten", 19.3, 411, 0.28, 1., -4, "Wikipedia"}, {"Copper Beryllium", 8.4, 134, 0.3, 5., -6, "www.alb-copperalloys.com"}, {"Sapphire", 3.98, 345, 0.3, 2., -7, "www.crystalsystems.com"}, {"Bolfur", 7.6, 157, 0.3, 1., -3, "www.alb-copperalloys.com"}} save$shape = {{"Doughnut[y]", {1.5, 0.8, 0.2, 0.2, 0.2}}, {"Cylinder[z]", {0.8, 0.14, 0.2, 0.2, 0.2}}} save$vspr = {{1, 1, 1750, 89.90865, True, -60, {0, 1, 0}, "NAB1", 10.}, {2, 1, 1750, 89.90865, True, -60, {0, 1, 0}, "NAB2", 10.}, {3, 1, 1750, 89.90865, True, -60, {0, 1, 0}, "NAB3", 10.}, {4, 1, 1750, 89.90865, True, -60, {0, 1, 0}, "NAB4", 10.}} save$wire = {{1, {-0.4189, 0., -0.05}, 2, {-0.4189, 0.0056, -0.05}, 3, 0.7455759544393915, 0.0007, 0, 0.0007, 89.90865, {0, -1, 0}, {0, 1, 0}, "g-NAB-1", 212025.91974076693, 0.005072885775047644}, {1, {-0.4189, 0., 0.05}, 2, {-0.4189, 0.0056, 0.05}, 3, 0.7455759544393915, 0.0007, 0, 0.0007, 89.90865, {0, -1, 0}, {0, 1, 0}, "g-NAB-2", 212025.91974076693, 0.005072885775047644}, {1, {0.4189, 0., 0.05}, 2, {0.4189, 0.0056, 0.05}, 3, 0.7455759544393915, 0.0007, 0, 0.0007, 89.90865, {0, -1, 0}, {0, 1, 0}, "g-NAB-3", 212025.91974076693, 0.005072885775047644}, {1, {0.4189, 0., -0.05}, 2, {0.4189, 0.0056, -0.05}, 3, 0.7455759544393915, 0.0007, 0, 0.0007, 89.90865, {0, -1, 0}, {0, 1, 0}, "g-NAB-4", 212025.91974076693, 0.005072885775047644}} save$wire2 = {{4, 1, 2, {0., 0.0056, 0.4189, 0.05, 0., 0.}, 3, {0.746, 0.0007, 0.02, 0.001, False, True}, 89.90865, {212025.91974076693, 0.005072885775047644}}} calc$combif = False calc$combif2 = False calc$combin = {xNAB + yNAB, xNAB - yNAB} calc$combout = {xNAB + yNAB, xNAB - yNAB} calc$eigenDET = 10 calc$eigenIS = 150 calc$eigenP = {True, True, True, True, True, True} calc$eigenPA = 0.5 calc$eigenTA = 0.2 calc$legend = True calc$lgndpos = {0.35, 0.3} calc$maxx = 2.5 calc$maxz = 2.5 calc$psdline = {} calc$psdname = {} calc$psdplot = {} calc$psdsel = {} calc$psdshow = {} calc$sflag = 0 calc$sni = {False, False, False, False, False, False} calc$sno = {False, False, False, False, False, False} calc$tfFi = {False, False, False, False, False, False} calc$tfFo = {False, False, False, False, False, False} calc$tfGi = {False, False, False, False, False, False} calc$tfGo = {False, False, False, False, False, False} calc$tfk = 1 calc$tfline = {} calc$tfname = {} calc$tfplot = {} calc$tfsel = {} calc$tfshow = {} calc$tnb = {False, False, False, False, False, False} calc$tnT = 300 calc$xrange = Automatic calc$ybottom = -2.5 calc$yrange = Automatic calc$ytop = 1. matKff = {{213775.91974076693, 0, 0, 0}, {0, 213775.91974076693, 0, 0}, {0, 0, 213775.91974076693, 0}, {0, 0, 0, 213775.91974076693}} matKffc = {{213775.91974076693 + 196.20259197407668*I, 0, 0, 0}, {0, 213775.91974076693 + 196.20259197407668*I, 0, 0}, {0, 0, 213775.91974076693 + 196.20259197407668*I, 0}, {0, 0, 0, 213775.91974076693 + 196.20259197407668*I}} matKfp = {{0, 0, 0, 0}, {212025.91974076693, 212025.91974076693, 212025.91974076693, 212025.91974076693}, {0, 0, 0, 0}, {10601.295987038347, -10601.295987038347, -10601.295987038347, 10601.295987038347}, {0, 0, 0, 0}, {-88817.65777940728, -88817.65777940728, 88817.65777940728, 88817.65777940728}} matKfpc = {{0, 0, 0, 0}, {212025.91974076693 + 21.202591974076693*I, 212025.91974076693 + 21.202591974076693*I, 212025.91974076693 + 21.202591974076693*I, 212025.91974076693 + 21.202591974076693*I}, {0, 0, 0, 0}, {10601.295987038347 + 1.0601295987038346*I, -10601.295987038347 - 1.0601295987038346*I, -10601.295987038347 - 1.0601295987038346*I, 10601.295987038347 + 1.0601295987038346*I}, {0, 0, 0, 0}, {-88817.65777940728 - 8.881765777940727*I, -88817.65777940728 - 8.881765777940727*I, 88817.65777940728 + 8.881765777940727*I, 88817.65777940728 + 8.881765777940727*I}} matKfv = {{0, 0, 0, 0}, {-212025.91974076693, -212025.91974076693, -212025.91974076693, -212025.91974076693}, {0, 0, 0, 0}, {-10601.295987038346, 10601.295987038346, 10601.295987038346, -10601.295987038346}, {0, 0, 0, 0}, {88817.65777940727, 88817.65777940727, -88817.65777940727, -88817.65777940727}} matKfvc = {{0, 0, 0, 0}, {-212025.91974076693 - 21.202591974076693*I, -212025.91974076693 - 21.202591974076693*I, -212025.91974076693 - 21.202591974076693*I, -212025.91974076693 - 21.202591974076693*I}, {0, 0, 0, 0}, {-10601.295987038346 - 1.0601295987038346*I, 10601.295987038346 + 1.0601295987038346*I, 10601.295987038346 + 1.0601295987038346*I, -10601.295987038346 - 1.0601295987038346*I}, {0, 0, 0, 0}, {88817.65777940727 + 8.881765777940727*I, 88817.65777940727 + 8.881765777940727*I, -88817.65777940727 - 8.881765777940727*I, -88817.65777940727 - 8.881765777940727*I}} matKpv = {{-491.76104232678995, 0, 0, 0, 0, -3.609568787884941}, {0, -6942.697007151903, 0, 0, 0, 0}, {0, 0, -491.76104232678995, 3.609568787884941, 0, 0}, {0, 0, -6.363430624915029, -17.310034384069695, 0, 0}, {0, 0, 0, 0, -87.54254984311416, 0}, {6.363430624915029, 0, 0, 0, 0, -1218.238392556581}} matKpvc = {{-491.76104232678995 - 0.0004841315987180795*I, 0, 0, 0, 0, -3.609568787884941 - 3.5535680087491487*^-6*I}, {0, -6943.259566389024 - 688.5914775489528*I, 0, 0, 0, 0}, {0, 0, -491.76104232678995 - 0.0004841315987180795*I, 3.609568787884941 + 3.5535680087491487*^-6*I, 0, 0}, {0, 0, -6.363430624915029 - 6.264704961570458*^-6*I, -17.311440782162663 - 1.7214786478889033*I, 0, 0}, {0, 0, 0, 0, -87.54254984311416 - 0.00008819354552270241*I, 0}, {6.363430624915029 + 6.264704961570458*^-6*I, 0, 0, 0, 0, -1218.3371088819695 - 120.83211126408442*I}} matKvv = {{491.76104232678995, 0, 0, 0, 0, -6.363430624915029}, {0, 6942.697007151903, 0, 0, 0, 0}, {0, 0, 491.76104232678995, 6.363430624915029, 0, 0}, {0, 0, 6.363430624915029, 22.092788841756374, 0, 0}, {0, 0, 0, 0, 87.54254984311416, 0}, {-6.363430624915029, 0, 0, 0, 0, 1223.021147014253}} matKvvc = {{491.76104232678995 + 0.0004841315987180795*I, 0, 0, 0, 0, -6.363430624915029 - 6.264704961570458*^-6*I}, {0, 6943.259566389024 + 688.5914775489528*I, 0, 0, 0, 0}, {0, 0, 491.76104232678995 + 0.0004841315987180795*I, 6.363430624915029 + 6.264704961570458*^-6*I, 0, 0}, {0, 0, 6.363430624915029 + 6.264704961570458*^-6*I, 22.09419523984934 + 1.721944145236065*I, 0, 0}, {0, 0, 0, 0, 87.54254984311416 + 0.00008819354552270241*I, 0}, {-6.363430624915029 - 6.264704961570458*^-6*I, 0, 0, 0, 0, 1223.1198633396707 + 120.83257676143157*I}} matGpv = {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}} matGvv = {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}} matMpv = {{0, 0, 0, 0, 0, 0}, {0, 0.03665999999999999, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}} matMvv = {{36.66, 0, 0, 0, 0, 0}, {0, 36.66, 0, 0, 0, 0}, {0, 0, 36.66, 0, 0, 0}, {0, 0, 0, 1.92, 0, 0}, {0, 0, 0, 0, 2.6, 0}, {0, 0, 0, 0, 0, 4.44}} allvars = {xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB} allparams = {xg, yg, zg, pitchg, yawg, rollg} allfloats = {hNAB1, hNAB2, hNAB3, hNAB4} subdefparams = {xg -> 0, yg -> 0, zg -> 0, pitchg -> 0, yawg -> 0, rollg -> 0} subminvars = {xNAB -> -3.0061690648916496*^-32, yNAB -> -0.7515999999999998, zNAB -> -1.4067265991291304*^-32, pitchNAB -> 1.0871075359068578*^-30, yawNAB -> 0., rollNAB -> -2.323144416745187*^-30, hNAB1 -> 1.2992781453899397*^-16, hNAB2 -> 1.2992781453899387*^-16, hNAB3 -> 1.2992781453899193*^-16, hNAB4 -> 1.2992781453899205*^-16} eigenf = {0.5336279607743916, 0.5828886417732024, 0.5886353204244023, 0.923513086053878, 2.190221346674201, 2.6414742305477867} eigenv = {{0., 0., -0.07965289705326095, 0.996822660251573, 0., 0.}, {-0.9999850433289613, 0., 0., 0., 0., -0.005469288653501728}, {0., 0., 0.5481759027063491, 0.8363630669105847, 0., 0.}, {0., 0., 0., 0., 1., 0.}, {0., 1., 0., 0., 0., 0.}, {0.0006624113339237442, 0., 0., 0., 0., -0.9999997806055883}} bodylist = {{{xg, yg, zg, pitchg, yawg, rollg}, 0, {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}, "Doughnut[y]", {1.5, 0.8, 0.2, 0.2, 0.2}}, {{xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB}, 36.66, {{1.92, 0, 0}, {0, 2.6, 0}, {0, 0, 4.44}}, "Cylinder[z]", {0.8, 0.14, 0.2, 0.2, 0.2}}, {{xg - 0.4189*Cos[rollg]*Cos[yawg] - 0.05*(Sin[pitchg]*Sin[rollg] + Cos[pitchg]*Cos[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB1)* (-(Cos[pitchg]*Sin[rollg]) + Cos[rollg]*Sin[pitchg]*Sin[yawg]), yg - 0.4189*Cos[yawg]*Sin[rollg] - 0.05*(-(Cos[rollg]*Sin[pitchg]) + Cos[pitchg]*Sin[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB1)* (Cos[pitchg]*Cos[rollg] + Sin[pitchg]*Sin[rollg]*Sin[yawg]), zg - 0.05*Cos[pitchg]*Cos[yawg] + (-0.007340086906449859 + hNAB1)* Cos[yawg]*Sin[pitchg] + 0.4189*Sin[yawg], pitchg, yawg, rollg}, 0, {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}, "Cylinder[y]", {0.1, 0.03, 0., 0, 0}}, {{xg - 0.4189*Cos[rollg]*Cos[yawg] + 0.05*(Sin[pitchg]*Sin[rollg] + Cos[pitchg]*Cos[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB2)* (-(Cos[pitchg]*Sin[rollg]) + Cos[rollg]*Sin[pitchg]*Sin[yawg]), yg - 0.4189*Cos[yawg]*Sin[rollg] + 0.05*(-(Cos[rollg]*Sin[pitchg]) + Cos[pitchg]*Sin[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB2)* (Cos[pitchg]*Cos[rollg] + Sin[pitchg]*Sin[rollg]*Sin[yawg]), zg + 0.05*Cos[pitchg]*Cos[yawg] + (-0.007340086906449859 + hNAB2)* Cos[yawg]*Sin[pitchg] + 0.4189*Sin[yawg], pitchg, yawg, rollg}, 0, {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}, "Cylinder[y]", {0.1, 0.03, 0., 0, 0}}, {{xg + 0.4189*Cos[rollg]*Cos[yawg] + 0.05*(Sin[pitchg]*Sin[rollg] + Cos[pitchg]*Cos[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB3)* (-(Cos[pitchg]*Sin[rollg]) + Cos[rollg]*Sin[pitchg]*Sin[yawg]), yg + 0.4189*Cos[yawg]*Sin[rollg] + 0.05*(-(Cos[rollg]*Sin[pitchg]) + Cos[pitchg]*Sin[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB3)* (Cos[pitchg]*Cos[rollg] + Sin[pitchg]*Sin[rollg]*Sin[yawg]), zg + 0.05*Cos[pitchg]*Cos[yawg] + (-0.007340086906449859 + hNAB3)* Cos[yawg]*Sin[pitchg] - 0.4189*Sin[yawg], pitchg, yawg, rollg}, 0, {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}, "Cylinder[y]", {0.1, 0.03, 0., 0, 0}}, {{xg + 0.4189*Cos[rollg]*Cos[yawg] - 0.05*(Sin[pitchg]*Sin[rollg] + Cos[pitchg]*Cos[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB4)* (-(Cos[pitchg]*Sin[rollg]) + Cos[rollg]*Sin[pitchg]*Sin[yawg]), yg + 0.4189*Cos[yawg]*Sin[rollg] - 0.05*(-(Cos[rollg]*Sin[pitchg]) + Cos[pitchg]*Sin[rollg]*Sin[yawg]) + (-0.007340086906449859 + hNAB4)* (Cos[pitchg]*Cos[rollg] + Sin[pitchg]*Sin[rollg]*Sin[yawg]), zg - 0.05*Cos[pitchg]*Cos[yawg] + (-0.007340086906449859 + hNAB4)* Cos[yawg]*Sin[pitchg] - 0.4189*Sin[yawg], pitchg, yawg, rollg}, 0, {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}, "Cylinder[y]", {0.1, 0.03, 0., 0, 0}}} bodylist /: bodylist::usage = "bodylist is a list of structures defining \ bodies. Each structure containing the following items:\n\tcoordinate list \ defining the body,\n\tmass,\n\tmoment of inertia \ tensor,\n\tshape,\n\tdimension" wirelist = {{{xg, yg, zg, pitchg, yawg, rollg}, {-0.4189, -0.007340086906449859 + hNAB1, -0.05}, {xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB}, {-0.4189, 0.01294008690644986, -0.05}, 0.7308957806264917, 212025.91974076693, 0.005072885775047644, 0.0001, 9.84485465598065*^-7, 0.0001}, {{xg, yg, zg, pitchg, yawg, rollg}, {-0.4189, -0.007340086906449859 + hNAB2, 0.05}, {xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB}, {-0.4189, 0.01294008690644986, 0.05}, 0.7308957806264917, 212025.91974076693, 0.005072885775047644, 0.0001, 9.84485465598065*^-7, 0.0001}, {{xg, yg, zg, pitchg, yawg, rollg}, {0.4189, -0.007340086906449859 + hNAB3, 0.05}, {xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB}, {0.4189, 0.01294008690644986, 0.05}, 0.7308957806264917, 212025.91974076693, 0.005072885775047644, 0.0001, 9.84485465598065*^-7, 0.0001}, {{xg, yg, zg, pitchg, yawg, rollg}, {0.4189, -0.007340086906449859 + hNAB4, -0.05}, {xNAB, yNAB, zNAB, pitchNAB, yawNAB, rollNAB}, {0.4189, 0.01294008690644986, -0.05}, 0.7308957806264917, 212025.91974076693, 0.005072885775047644, 0.0001, 9.84485465598065*^-7, 0.0001}} wirelist /: wirelist::usage = "wirelist is a list of structures defining \ wires. Each structure containing the following items:\n\tcoordinate list \ defining first mass,\n\tattachment point for first mass (local \ coordinates),\n\tcoordinate list defining second mass,\n\tattachment point \ for second mass (local coordinates),\n\tunstretched length,\n\tstretching \ elasticity,\n\ttorsional elasticity,\n\tloss angle of the stretching \ mode,\n\tloss angle of the bending mode,\n\tloss angle of the torsion mode," calcdone = True