Расчет параметров численным интегрированием

Данная программа разработана для расчета момента инерции сложного сечения. Для расчета сечение было разбито по горизонтали на одинаковые по высоте участки ( шаги ). Длинны участков были замеряны и внесены в таблицу. Затем запускаем программу и последовательно вводим из таблицы значения длинн участков. Программа расчитывает момент инерции участка и площадь участка – затем суммирует эти величины. При вводе нуля на запрос длинны – программа показывает результат работы. Программу легко переделать на расчеты интеграции по другим законам.

# -*– coding: cp1251 -*-

# Алгоритм расчета момента инерции и площади интегрированием

import math # Подключили математич модуль

u=" Алгоритм временного расчета Jx интегрированием,"

uu=" "

uuu=" ...... ...... ....... ....... ...... ....... ....... "

print (uu)

print (uu)

print (u)

print (uuu)

print (uu)

u=" Ввести Шаг вдоль оси Y-Y "

print(u)

t=0.00

t=input( ) # Вводим число

t=float(t) # Принудительно в вещественное число

print (uu)

J=0.0000

J=float(J)

Jx=0.0000

Jx=float(Jx)

S=0.0000

S=float(S)

Ss=0.0000

Ss=float(Ss)

d=0.00000001

d=float(d)

L=0.0

L=float(L)

n=1

n=float(n)

while d>0:

# Далее строки начинаются с четырех пробелов.

L=L+t # Суммируем смещение по оси

ss=str(n) # Преобразуем число в строку

Lt=str(L) # Преобразуем число в строку

u=" Ввести высоту по оси Y-Y Ввод № "

u=u+ss+" Смещение по оси = "+Lt

print(u)

d=0.00

d=input( ) # Вводим число

d=float(d) # Принудительно в вещественное число

print (uu)

J=t*d*d*d/12 # Находим момент инерции элемента

Jx=Jx+J # Суммируем J

S=t*d # Нашли площадь элемента

Ss=Ss+S # Суммируем S

n=n+1 # Счетчик числа вводов

# Далее строки начинаются без пробелов.

Pii=math.pi # Вытащили число " Пи "

# Rk=math.sqrt(a) # Квадратный корень из " а "

# sur=math.asin(x) # Угол А в радианах

# su=sur*180/Pii # Угол А в градусах

print (uu)

uuuu=" ,,,, ,,,, ,,,, ,,,,, ,,,,, ,,,, ,,,,, ,,,,, ,,,,, "

print (uuuu)

print (uu)

u=" Площадь сечения = "

ss=str(Ss) # Преобразуем число в строку

u=u+ss

print (u)

print (uu)

u=" Момент инерции по оси Х-Х = "

ss=str(Jx) # Преобразуем число в строку

u=u+ss

print (u)

print (uu)

print (uu)

u=" ...... ...... ...... Конец рассчета ...... ...... ...... "

print (u)

print (uu)

print (uu)

input( ) # Ожидание нажима Ентер

print (uu)

# .............

# Конец программы.

Загрузка...