% FILENAME: MC2003 SEPT 25.m % DESCRIPTION: Mode Cleaner Mirror Triple Pendulum % VERSION NO: MC file version 6 % NOTE: filename format RM.m % version number to be incremented plus date changed %%********************************************************************** %%COMMENTS%% %% PARAMETERS FOR THE MC SUSPENSION %% NEW DESIGN %% VERSION HAS NO LEVER_ARM INPUT %% Includes overide of the calculation of the lower blade frequency %% NUMBERS IN SLIDER GAIN MATCH THOSES CALCULATED FROM GEO_SR SUSPENSION %%********************************************************************** %**************************************************** %% Written by Calum Iain Torrie 7/98 %% modified by K.A. Strain 6/99 %% modified BY CALUM TORRIE DEC 2001 version 1 %% Norna/CIT - summer 2002 version 2 %% CIT - 01/Sept/2002 version 3 %% CIT - 29/Apr/2003 version 4 %% CIT - 01/May2003 (included masses and d's for a particular suspended suspension version 5 %% CIT - 25/SEPT/2003 (updated m1 and the d's) default CP version 6 %**************************************************** %% coordinates x = longitudinal = u_LIGO roll about this axis %% y = transverse = v_LIGO pitch about this axis %% z = vertical = w_LIGO yaw about this axis g =9.81; %%************************************************************************ %%%%%%%%%%%%%%%%%%%%%%%%%%% UPPER MASS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%*****1st ATTEMPT***** %% top mass is REPRESENTED by a rectangular BLOCK %% in reality it will be larger and less dense. %ux = 0.10; %dimensions of UPPER MASS (square) %uy = 0.18; %uz = 0.10; %den1 = 7000; %density (steel with holes) %m1 = den1* uy* uz* ux; %mass %I1x = m1*( uy^2+ uz^2)/12; %moment of inertia (transverse roll) %I1y = m1*( uz^2+ ux^2)/12; %moment of inertia (longitudinal pitch) %I1z = m1*( uy^2+ ux^2)/12; %moment of inertia (yaw) %m1_parameters = 'represented by a regtangular block'; %material1 = 'steel' %%%%%%%%%%%%%%%%%%%%%%% %%*****2nd ATTEMPT***** %% T-shaped calculated from mofi2.m % m1 = 3.1590; %T-shaped calculated from mofi2.m % I1x = 0.0207; % I1y = 0.0024; % I1z = 0.0206; %steel T piece, dimensions (x,y,z): % m1_parameters = 'T-shaped calculated from mofi2.m'; % material1 = 'steel' %%%%%%%%%%%%%%%%%%%%%%% %%*****3rd ATTEMPT***** %% FROM SOLIDWORKS ASSEMBLY OR ACTUAL MASS m1 = 3.130; %Actual shape calculated from SWorks 25 SEPT 2003 (assuming denisty of steel as 7800) %m1 = 3.110; %Updated Weight from lab-scales 24 SEPT 2003 CIT %??????????????????????????????????????? I1x = 0.0238; I1y = 0.0024; I1z = 0.0238; m1_parameters = 'Calculated'; material1 = 'combination steel+alum'; %%%%%%%%%%%%%%%%%%%%%%% pend.m1_parameters = m1_parameters; pend.material1 = material1; pend.m1 = m1; pend.I1x = I1x; pend.I1y = I1y; pend.I1z = I1z; %******************************************************************** %%%%%%%%%%%%%%%%%%%%INTERMEDIATE MASS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 1st Attempt - NOISE PROTOTYPE % Silica cylindrical mass ix = 0.075; %dimension of INTERMEDIATE MASS (cylinder) ir = 0.075; %den2 = 2202; %density (fused silica) %m2 = den2*pi* ir^2* ix ;%intermediate mass %I2x = m2*(ir^2/2); %moment of inertia (transverse roll) %I2y = m2*(ir^2/4+ix^2/12); %moment of inertia (longitudinal pitch) %I2z = m2*(ir^2/4+ix^2/12); %moment of inertia (yaw) %m2_parameters = 'Noise P-type: Silica Mass without flats and ears'; %material2 = 'silica'; % 2nd Attempt - CONTROLS PROTOTYPE % Aluminium "dummy" intermediate mass m2 = 2.979; %Actual shape calculated from SWorks 01/09/02 I2x = 0.00865; I2y = 0.00564; I2z = 0.00574; m2_parameters = 'Controls P-type: Calculated'; material2 = 'alum with holes + s/steel clamps'; pend.m2_parameters = m2_parameters; pend.material2 = material2; pend.ix = ix; pend.ir = ir; pend.m2 = m2; pend.I2x = I2x; pend.I2y = I2y; pend.I2z = I2z; %********************************************************************* %%%%%%%%%%%%%%%%%%%%%%%%%%%% TEST MASS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 1st Attempt - NOISE PROTOTYPE % Silica cylindrical test mass tx = 0.075; %dimensions of TEST MASS (cylinder) tr = 0.075; % den3 = 2202; %density (fused silica) % m3 = den3*pi*tr^2*tx; %test mass % I3x = m3*(tr^2/2); %moment of inertia (transverse roll) % I3y = m3*(tr^2/4+tx^2/12); %moment of inertia (longitudinal pitch) % I3z = m3*(tr^2/4+tx^2/12); %moment of inertia (yaw) % m3_parameters = 'Noise P-type: Silica Mass without flats and ears'; % material3 = 'silica'; % 2nd Attempt - CONTROLS PROTOTYPE % Aluminium "dummy" test mass m3 = 3.040; %Actual shape calculated from SWorks 01/09/02 I3x = 0.00876; I3y = 0.00589; I3z = 0.00565; m3_parameters = 'Controls P-tpye: Calculated'; material3 = 'alum with holes + s/steel clamps'; pend.m3_parameters = m3_parameters; pend.material3 = material3; pend.tx = tx; pend.tr = tr; pend.m3 = m3; pend.I3x = I3x; pend.I3y = I3y; pend.I3z = I3z; %******************************************************************** l1 = 0.445-0.15; %upper wire length l2 = 0.287-0.12; %intermediate wire length l3 = 0.28-0.06; %lower wire length pend.l1 = l1; pend.l2 = l2; pend.l3 = l3; %****************************************************************************************** nw1 = 2; % number of wires (= number of cantilevers if fitted) per stage (2 or 4) nw2 = 4; nw3 = 4; pend.nw1 = nw1; pend.nw2 = nw2; pend.nw3 = nw3; %*********************************************************************************** r1 = 180e-6; % radius of upper wire r2 = 100e-6; % radius of intermediate wire %%CONTROLS PROTOTYPE r3_parameters = 'Spring Steel Wires'; r3 = 75e-6; % radius of lower wire %%NOISE PROTOTYPE % r3_parameters = 'Fused-Silica'; % r3 = 75e-6; % radius of lower wire pend.r1 = r1; pend.r2 = r2; pend.r3_parameters = r3_parameters; pend.r3 = r3; %************************************************************************************ Y1 = 2.2e11; % Changed to YM for California Steel Wire by CIT 01/09/02, previously 1.65e11 Y2 = 2.2e11; % Changed to YM for California Steel Wire by CIT 01/09/02 %%CONTROLS PROTOTYPE Y3_parameters = 'Spring Steel Wires'; Y3 = 2.2e11; % Changed to YM for California Steel Wire by CIT 01/09/02 %%NOISE PROTOTYPE % Y3_parameters = 'Fused-Silica Fibres'; % Y3 = 7e10; % Youngs Modulus of lower wire (fused silica) pend.Y1 = Y1; pend.Y2 = Y2; pend.Y3_parameters = Y3_parameters; pend.Y3 = Y3; %********************************************************************** %********************************************************************** %blade design - upper blades mntb = (m1 +m2 +m3)/2;%total per blade mnb = m1/2;%uncoupled mass %opt(uncoupled mass, total mass per blade, max stress, max length, max width) %[uf,lnb,anb,hnb,stn] = opt(mnb,mntb,8e8,0.24,0.05); %ufc1 = uf; %***************************** % overide of the calculation of the UPPER blade frequency by CIT 01/09/02 ufc1 = 2.29; pend.ufc1 = ufc1; %***************************** % pend.l1b = lnb; % pend.a1b = anb; % pend.h1b = hnb; % pend.ufc1 = uf; % pend.st1 = stn; % pend.intmode_1 = 55*hnb*0.37^2/(0.002*lnb^2); %scaled from GEO blade %blade design - lower blades mntb = (m2 +m3)/4;%total per blade mnb = m2/4;%uncoupled mass %opt(uncoupled mass, total mass per blade, max stress, max length, max width) %[uf,lnb,anb,hnb,stn] = opt(mnb,mntb,8e8,0.12,0.018); %ufc2 = uf; %***************************** % overide of the calculation of the lower blade frequency by CIT 01/09/02 ufc2 = 3.22; pend.ufc2 = ufc2; %***************************** % pend.l2b = lnb; % pend.a2b = anb; % pend.h2b = hnb; % pend.ufc2 = uf; %***************************** % overide of the calculation of the lower blade frequency %pend.ufc2 = 2.8; %***************************** % pend.st2 = stn; % pend.intmode_2 = 55*hnb*0.37^2/(0.002*lnb^2); %scaled from GEO blade %****************************************************************************************** %********************************************************************** %d0 = 0.001; %height of upper wire break-off (above c.of m. upper mass) d0 =0.005; %changed to avoid close coupling of first pitch and long. modes 2/4/02 (actually 4.8mm SW, 25 SEPT 2003 CIT) %d1 = 0.001; %d1 = 0.004; %CHANGED Jan03 height of intermediate wire break-off (below c.of m. upper mass) d1 = 0.002; %MEASURED FROM SW 25 SEPT 2003 CIT d2 = 0.001; %height of intermediate wire break-off (above c.of m. of int. mass) d3 = 0.001; %height of lower wire break-off (below c.of m. intermediate mass) %d4 = 0.001; %d4 = 0.0025; % CHANGED Jan03 height of lower wire break-off (above c.of m.test mass) d4 = 0.001; % MEASURED SW 25 SEPT 2003 CIT pend.d0 = d0; pend.d1 = d1; pend.d2 = d2; pend.d3 = d3; pend.d4 = d4; %****************************************************************************************** % X direction separation su = 0.00; % 1/2 separation of upper wires si = 0.0285; % 1/2 separation of intermediate wires %si = 0.026; % sl = 0.005; % 1/2 separation of lower wires %sl = 0.004; pend.su = su; pend.si = si; pend.sl = sl; %****************************************************************************************** % Y direction separation %n0 = 0.03; % 1/2 separation of upper wires at suspension point % edited to fit in with Janeen's cage structure %n0 = 0.077; n0 = 0.0773; %MEASURED IN SW 25 SEPT 2003 CIT n1 = 0.1; % 1/2 separation of upper wires at upper mass %n2 = 0.04; % 1/2 separation of intermediate wires at upper mass n2 = 0.039; % MEASURED IN SW 25 SEPT 2003 CIT n3 = ir-0.0035+0.005; % 1/2 separation of intermedite wires at intermediate mass n4 = tr-0.0035+0.005; % 1/2 separation of lower wires at intermediate mass n5 = tr-0.0035+0.005; % 1/2 separation of lower wires at test mass pend.n0 = n0; pend.n1 = n1; pend.n2 = n2; pend.n3 = n3; pend.n4 = n4; pend.n5 = n5; %*********************************************************************************** %LENGTH CALCULATIONS tl1 = sqrt(l1^2 - (n0-n1)^2); tl2 = sqrt(l2^2 - (n2-n3)^2); tl3 = sqrt(l3^2 - (n4-n5)^2); pend.tl1 = tl1; pend.tl2 = tl2; pend.tl3 = tl3; %Overall length of pendulum from upper most wire break off to bottom edge of test mass l_total = tl1 + tl2 + tl3 + d0 + d1 + d2 + d3 + d4 + tr; pend.l_total = l_total; %Length of pendulum from upper most wire break off to centre of mass of test mass l_com = tl1 + tl2 + tl3 + d0 + d1 + d2 + d3 + d4; pend.l_com = l_com; %*********************************************************************************** % represents small loss bd = 0.01; % makes phases of open loop plots look nicer % CIT NAR DEC/JAN 2001/2002 % LEVER ARM NOT CURRENTLY CALLED FROM PARAMETER FILE THIS IS INCLUDED IN THE GAIN TRIANGLE of pendn.m %IMPORTANT FOR REFERENCE ONLY!!! % lever_pitch = 0.03; % lever_yaw = 0.08; % lever_roll = 0.06; %pend.lever_pitch=lever_pitch; %pend.lever_yaw=lever_yaw; %pend.lever_roll=lever_roll; %removd the "%" comment added by CIT Nov 2001 %pend %IMPORTANT FOR REFERENCE ONLY!!! %Gain triangles in pendn.m %NEW NUMBERS FROM NORNA AUGUST 2002 %gain = 0.06; %Gain triangle = (leverarm)^2 * (no. of coils) * gain %***Note that gains are linked in different directions due to common coils as below***** %***If gains in one channel only need to be changed then leverarm can be used*********** %gainzrtrl =gain; % vertical, z, pitch, rt, roll rl (coils on top of upper mass) %gaint = gain.*2; % transverse, t (coil on one end of upper mass) %gainlrz = gain; % longitudinal, l, yaw, rz (coils on long rear side of upper mass) %long = (1)^2 * 2 * gainlrz = 0.12 %pitch = (0.03)^2 * 2 * gainzrtrl = 1.08e-4 %vert = (1)^2 * 3 * gainzrtrl = 0.18 %yaw = (0.08)^2 * 2 * gainlrz = 7.68e-4 %trans = (1)^2 * 1 * gaint = 0.12 %roll = (0.06)^2 * 3 * gainzrtrl = 6.48e-4 %************************************************************************************** %VIEWER PREFERENCES (under edit on LTI viewer) %UNITS %Frequency in Hz using Log scale %Magnitude in absolute using log scale %Phase in degrees %TIME VECTOR %[0:0.01:50] %FREQUENCY VECTOR %logspace(-0.5,1.5,500) %**************************************** % EXTRACTING DATA FROM TRANSFER FUNCTIONS and IMPULSE RESPONSE %TRANSFER FUNCTION % File export - export to workspace %sys = tf(pendn_1(1,1)); [or equivalent i.e.sys = tf(pendn_1(6,10));] %damp(sys); %IMPULSE RESPONSE % right click in screen and select properties - characteristics - select settling time %set settling time 36.7% (1/e)