#!/usr/bin/bc

# implement the deflection equations given in the tslots.com catalog

scale = 3  # request 3 digits after the decimal

l = 60  # inches

# modulus of elasticity
e1 = 10000000	# tslots.com
e2 = 10200000	# 8020.net

# values for tlots.com products -- use with "e1"
i[1010] = .046
i[1020] = .321
i[1021] = .087  # 1020 on side
i[1030] = .986
i[1515] = .172   # 1515vl
i[1530] = 1.431  # 1530l
i[1531] = .408  # 1530l on side

# from 8020.net -- use with "e2"
i[2515] = .183   # 8020.net 1515ULS ( .936 lb/ft)
i[3515] = .255   # 8020.net 1515    (1.3428 lb/ft)
i[2530] = 1.382  # 8020.net 1530L   (2.011 lb/ft)
i[3530] = 1.233  # 8020.net 1530UL  (1.608 lb/ft)
i[4530] = 1.276  # 8020.net 1530ULS (1.590 lb/ft)


define deflections(f, l, p)
{

    fl3 = f * (l^3)
    if (p < 2000 ) {
	ei = e1 * i[p]
    } else {
	ei = e2 * i[p]
    }

    sd = (5 * fl3) / (384 * ei)
    sp = fl3 / (48 * ei)
    fd = fl3 / (384 * ei)
    fp = fl3 / (192 * ei)

    print " ", p, "	", sd, "\q	", sp, "\q		", fd, "\q	", fp, "\q\n"

}

define all_deflections(f, l)
{
    print f, "lbs\n"
#    r = deflections(f, l, 1010)
#    r = deflections(f, l, 1020)
#    r = deflections(f, l, 1021)
#    r = deflections(f, l, 1030)
#    r = deflections(f, l, 1515)
    r = deflections(f, l, 2515)
#    r = deflections(f, l, 3515)
#    r = deflections(f, l, 1530)
#    r = deflections(f, l, 2530)
#    r = deflections(f, l, 3530)
    r = deflections(f, l, 4530)
#    r = deflections(f, l, 1531)
    print "\n"
}

print "length is ", l, "\q\n"

    print "	  supported		   fixed\n"
    print "	 even	point		even	point\n"


r = all_deflections(400, l)
r = all_deflections(300, l)
r = all_deflections(200, l)
r = all_deflections(100, l)
r = all_deflections(67, l)
r = all_deflections(40, l)
r = all_deflections(30, l)
r = all_deflections(20, l)

quit