Formule del Roark

Gli script di programmazione raccontano storie che leggono dati e forniscono dati. Gli script bisogna interpretarli a livello di funzionamento e a livello di risultati. L'interprete di esecuzione degli script in questo caso è Python 3.2. I risultati dipendono invece dall'oggetto di studio. Fra un milione di scelte su cui lavorare un classico è l'analisi lineare di travi rettilinee. Il Roark fornisce una casistica ben strutturata da interpretare e implementare negli script. Sotto è il caso della trave a sbalzo

Tab 8 caso 1.0

Script in python

import math

#valori di input
l = 600. #[mm] lunghezza trave
W = 50. #[kN] carico concentrato
a = 0. #[mm] distanza dal vincolo A del carico
I = 416.67 #[Cm^4] momento d'inerzia
E = 210000. #[MPa] modulo elasticità acciaio

#calcolo
'''caso 8.1a pag. 100 (7a edizione in pdf-2542X_08 pag. 65'''
I *= 10000 #[mm^4]
W *= 1000 #[N]
RA = 0.
MA = 0.
yA = -W/(6*E*I) * (2*l**3 - 3*l**2*a + a**3)
ThetaA = (W * (l - a)**2) / (2 * E * I)
MB = -W * (l - a) / (1000*1000)
RB = W / 1000
yB = 0.
ThetaB = 0.
outputA = 'RA = {0:.1f}, MA = {1:.1f},yA = {2:.2f}, ThetaA = {3:.2f}'.format(RA, MA, yA, ThetaA)
outputB = 'RB = {0:.1f}, MB = {1:.1f},yB = {2:.2f}, ThetaB = {3:.2f}'.format(RB, MB, yB, ThetaB)
print('Reazioni in kN, momenti in kN*m, freccia in mm')
print(outputA)
print(outputB)