h=8#
F=100*10**3#
t=0.707*h#
A=4*60*t#
T=F/A#
print "T is %0.1f MPa "%(T)#
FOS=3#
Ta=95/FOS#
P=350*10**3#
h=12.5#
t=0.707*h#
l=P/(2*t*Ta)#
print "l is %0.0f mm "%(l)#
h=12#
t=0.707*h#
l=60#
Ta=80#
P=2*l*t*Ta#
P=P*10**-3#
print "P is %0.3f kN "%(P)#
from math import sin, pi, cos, sqrt
P=6*10**3#
e=150+(100/2)#
T=P*e#
A=200#
Td=P/A#
r=sqrt(2*50**2)#
Ixx=2*(100*50**2)#
Iyy=2*100**3/12#
IG=Ixx+Iyy#
Ts=r*T/IG#
Tmax=sqrt((Ts*sin(45*pi/180))**2+(Td+(Ts*cos(45*pi/180)))**2)#
Ta=80#
t=Tmax/Ta#
h=sqrt(2)*t#
h=3#
print "h is %0.0f mm "%(h)#
from math import sqrt
h=10#
t=10/sqrt(2)#
Ta=80#
x=((50*25)+(50*0))/(50+50)#
y=x#
ra=sqrt(x**2+37.5**2)#
Ixx=(7.07*50**3/12)+(50*7.07*(12.5**2))+(50*7.07*12.5**2)#
IG=2*Ixx#
e=100+(50-12.5)#
Tr=16.09*10**-3#
P=Ta/Tr#
P=P*10**-3#
print "P is %0.3f KN "%(P)#
from math import sqrt,pi
P=16*10**3#
l=300#
r=50#
M=P*l#
A=2*pi*r#
Ixx=pi*r**3#
sigb=M*r/Ixx#
Td=P/A#
Tmax=sqrt((sigb/2)**2+(Td**2))#
Ta=90#
t=Tmax/Ta#
h=sqrt(2)*t#
h=5#
print "h is %0.0f mm "%(h)#
sigut=415#
sige=sigut/3#
Ka=0.5#
Kb=0.85#
Kc=0.897#
SCF=1.5#
Kd=1/SCF#
FOS=2#
sige1=sige*Ka*Kb*Kc*Kd/FOS#
Pa=50*10**3#
h=10#
t=0.707*h#
l=Pa/(2*sige1*t)#
print "l is %0.0f mm "%(l)#
from math import sqrt,pi
l=300#
P=30*10**3#
T=P/(2*l)#
Ta=124#
t1=T/Ta#
h1=sqrt(2)*t1#
M=P*l#
Ixx=2*100*110**2#
sigb=M/Ixx*110#
#Let the allowable bending stress is Tab
Tab=200#
t2=sigb/Tab#
h2=t2/0.707#
h2=3#
print "h is %0.0f mm "%(h2)#
Ta=60#
l1=60#
l2=40#
P1=Ta*0.707*l1#
P2=Ta*0.707*l2#
P=80*10**3#
h=P/(P1+P2)#
h=20#
a=(P2*100)/(P1+P2)#
print " h is %0.0f mm "%(h)#
print "\n a is %0.0f mm "%(a)#
from math import sqrt
P=300*10**3#
l=500#
A=2*l#
Td=P/A#
T=(350-250)*P#
IG=(l**3*2/12)+(l*2*5**2)#
r=sqrt(250**2+5**2)#
Ts=T*r/IG#
Ts=Ts+Td#
Ta=110#
t=Ts/Ta#
h=t/0.707#
h=9#
print "h is %0.0f mm "%(h)#
t=30#
sigut=417#
sige=sigut/2#
Ka=0.5#
Kb=0.85#
Kc=0.897#
SCF=1.2#
Kd=1/SCF#
FOS=1.5#
sige1=sige*Ka*Kb*Kc*Kd/FOS#
Pa=60*10**3#
l=Pa/(sige1*t)#
print "l is %0.1f mm "%(l)#