circular_plate.py

Created by tlake
Created on September 11, 2023
1.05 KB
Test program from Victor 4900 manual
from math import *

def f(t,r,y,d,p):
    return (.476*t/r**2)*sqrt(y/(d*(1-p**2)))

print("THIS PROGRAM USES THE FORMULA\n")
print("F=(.476T÷R*2)SQRT((Y÷D(1-P*2))\n")
print("TO FIND THE FUNDAMENTAL FREQ-")
print("UENCY, F, OF A CIRCULAR CLAMP-")
print("ED PLATE,\n")
print("WHERE P = POISSON'S RATIO")
print("      D = DENSITY")
print("      Y = YOUNG'S MODULUS")
print("      T = THICKNESS")
print("      R = RADIUS\n")
print("THE PROGRAM CALCULATES F FOR A")
print("RANGE OF RADII (R) FROM .1 TO")
print("1 IN STEPS OF .1, FROM 1 TO 10")
print("IN STEPS OF 1, AND FROM 10 TO")
print("100 IN STEPS OF 10.\n\n")

print("ENTER P")
p=float(input())
print("ENTER D")
d=float(input())
print("ENTER Y")
y=float(input())
print("ENTER T")
t=float(input())

print("\nRADIUS  FUNDAMENTAL FREQUENCY\n")
for r in range(1,10):
    print("{:5.1f} {:20.3f}".format(r/10.0,f(t,r/10.0,y,d,p)))
    
for r in range(1,10):
    print("{:5.1f} {:20.3f}".format(r,f(t,r,y,d,p)))
    
for r in range(1,11):
    print("{:5.1f} {:20.3f}".format(r*10.0,f(t,r*10.0,y,d,p)))