DYNAMIC ANALYSIS PROGRAM FOR MACHINE FOUNDATION. 

 PROGRAM INFO / LIMITATIONS.
( Actual program is available in Super Civil CD )

The foundation should be rectangular.   
Foundation is considered as safe if frequency ratios are less than 0. 80 or 
more than 1. 20 .
The user should ensure that the eccentricity of foundation should not exceed 5 % of the least length in any horizontal direction.

 CALCULATION OF VERTICAL FREQUENCY 

A = L * B
M = MG / 9.81
MF = L * B * D * 2.5 / 9.81
W = MG + L * B * D * 2.5
AFNZ = (CZ * A / (4 * 3.142 ^ 2 * W / 9.81)) ^ .5
AFNZ = AFNZ * 60
RATIO = AFNZ / F

 CALCULATION OF VERTICAL AMPLITUDE 

WM = CIRCULAR FREQUENCY OF MACHINE IN RAD/SEC.
WM = F * 2 * 3.142 / 60
WNZ = CIRCULAR FREQUENCY IN VERTICAL DIRECTION.
WNZ = AFNZ * 2 * 3.141 / 60
AMPZ = VERTICAL AMPLITUDE IN MMS
AMPZ = PZ / ((W / 9.81) * (WNZ ^ 2 - WM ^ 2)) * 1000

 CALC. OF FREQUENCY FOR SLIDING & ROCKING MOTION IN X-Z PLANE (F0Y) 

IY = (B * L ^ 3) / 12
XB = (M * XI + MF * L * .5) / (W / 9.81)
YB = (M * YI + MF * B * .5) / (W / 9.81)
ZB = (M * ZI + MF * D * .5) / (W / 9.81)
MY1 = (MF * (L ^ 2 + D ^ 2)) / 12
MY21 = M * ((XB - XI) ^ 2 + (ZB - ZI) ^ 2)
MY22 = MF * ((XB - L * .5) ^ 2 + (ZB - .5 * D) ^ 2)
MY2 = MY21 + MY22
MY = MY1 + MY2
M0Y = MY + W / 9.81 * ZB ^ 2
F0Y = ((C0 * IY - W * ZB) / (4 * 3.142 ^ 2 * M0Y)) ^ .5
F0Y = F0Y * 60
RATIO = F0Y / F

 CALC. OF FREQUENCY FOR SLIDING & ROCKING MOTION IN Y-Z PLANE (F0X) 

IX = (L * B ^ 3) / 12
MX1 = (MF * (B ^ 2 + D ^ 2)) / 12
MX21 = M * ((YB - YI) ^ 2 + (ZB - ZI) ^ 2)
MX22 = MF * ((YB - B * .5) ^ 2 + (ZB - .5 * D) ^ 2)
MX2 = MX21 + MX22
MX = MX1 + MX2
M0X = MX + W / 9.81 * ZB ^ 2
F0X = ((C0 * IX - W * ZB) / (4 * 3.142 ^ 2 * M0X)) ^ .5
F0X = F0X * 60
RATIO = F0X / F

 CALC. OF FREQUENCY FOR TWISTING MOTION @ Z AXIS (TORSION-> AFNX ) 

AFNX = (CS * A / (4 * 3.142 ^ 2 * W / 9.81)) ^ .5
AFNX = AFNX * 60
RATIO = AFNX / F

 CALC. OF COUPLED NATURAL FREQUENCY IN X-Z PLANE AFN1 & AFN2 

WTHETAY = ((F0Y / 60) ^ 2 * 4 * 3.142 ^ 2) ^ .5
ALFAY = MY / M0Y
WX = (CT * A / (W / 9.81)) ^ .5
WN11 = WTHETAY ^ 2 + WX ^ 2
WN12 = ((WTHETAY ^ 2 + WX ^ 2) ^ 2 - 4 * ALFAY * WTHETAY ^ 2 * WX ^ 2) ^ .5
WN1 = ((.5 / ALFAY) * (WN11 + WN12)) ^ .5
WN2 = ((.5 / ALFAY) * (WN11 - WN12)) ^ .5
AFN1 = WN1 * 60 / (2 * 3.142)
AFN2 = WN2 * 60 / (2 * 3.142)
RATIO1 = AFN1 / F
RATIO2 = AFN2 / F

 CALC. OF COUPLED NATURAL FREQUENCY IN Y-Z PLANE FM1 & FM2 

WTHETAX = ((F0X / 60) ^ 2 * 4 * 3.142 ^ 2) ^ .5
ALFAX = MX / M0X
WY = (CT * A / (W / 9.81)) ^ .5
WM11 = WTHETAX ^ 2 + WY ^ 2
WM12 = ((WTHETAX ^ 2 + WY ^ 2) ^ 2 - 4 * ALFAX * WTHETAX ^ 2 * WY ^ 2) ^ .5
WM1 = ((.5 / ALFAX) * (WM11 + WM12)) ^ .5
WM2 = ((.5 / ALFAX) * (WM11 - WM12)) ^ .5
FM1 = WM1 * 60 / (2 * 3.142)
FM2 = WM2 * 60 / (2 * 3.142)
RATIO1 = FM1 / F
COUPLED NATURAL FREQUENCY(2) @ Y-Z PLANE = "; FM2; " IN RPM
RATIO2 = FM2 / F'

 CALC. OF AMPLITUDE FOR HORIZONTAL & ROCKING MOTION IN X-Z PLANE 

MIX = PZ * YE + PY * ZE + MXV
MIY = PZ * XE + PX * ZE + MYV
MIZ = PX * YE + PY * XE + MZV
AX1 = (C0 * IY - W * ZB + CT * A * ZB ^ 2 - MY * WM ^ 2) * PX
AX2 = CT * A * ZB * MIY
AX3 = (W / 9.81) * MY * (WN1 ^ 2 - WM ^ 2) * (WN2 ^ 2 - WM ^ 2)
AMPX = (AX1 + AX2) / AX3 * 1000

 CALC. OF AMPLITUDE FOR HORIZONTAL & ROCKING MOTION IN Y-Z PLANE 

MIX = PZ * YE + PY * ZE + MXV
MIY = PZ * XE + PX * ZE + MYV
MIZ = PX * YE + PY * XE + MZV
AY1 = (C0 * IX - W * ZB + CT * A * ZB ^ 2 - MX * WM ^ 2) * PY
AY2 = CT * A * ZB * MIX
AY3 = (W / 9.81) * MX * (WM1 ^ 2 - WM ^ 2) * (WM2 ^ 2 - WM ^ 2)
AMPY = (AY1 + AY2) / AY3 * 1000

 CALC. OF AMPLITUDE FOR TWISTING MOTION ABOUT Z AXIS IN RADIANS 

WSAI = AFNX * 2 * 3.14 / 60
MZ1 = (MF * (L ^ 2 + B ^ 2)) / 12
MZ21 = M * ((XB - XI) ^ 2 + (YB - YI) ^ 2)
MZ22 = MF * ((XB - L * .5) ^ 2 + (YB - .5 * B) ^ 2)
MZ2 = MZ21 + MZ22
MZ = MZ1 + MZ2
ASAI = 1 / (WSAI ^ 2 - WM ^ 2) * MIZ / MZ