Полезные программы Python-3. Книга третья. Сергей Фёдорович Гаврилов

программы.

      # -*– coding: cp1251 -*-

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

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

      uu=" "

      print (uu)

      u=" Расчет по Радиусу и Прогибу "

      print (u)

      print (uu)

      u=" ...... ...... ...... ....."

      print (u)

      print (uu)

      u=" Введите Радиус "

      print (u)

      print (uu)

      R=0.00000000

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

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

      u=" Введите Прогиб "

      print (uu)

      print (u)

      print (uu)

      b=0.00000000

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

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

      k=R-b

      x=(R*R)-(k*k)

      a=(math.sqrt(x))*2

      x=(a/2)/k

      aur = (math.atan(x))*2

      au=aur*180/Pii

      yg = au

      # Вывод по Хорде Прогиб Радиус

      # Далее следует блок вывода результата ( смотри выше ).

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

      # Позволяет рассмотреть результаты расчета

      # ..... Конец листинга программы ..... ..

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

      Расчет по Прогибу и длине Дуги

      Листинг программы.

      # -*– coding: cp1251 -*-

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

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

      uu=" "

      print (uu)

      u=" Расчет по Прогибу и длине Дуги "

      print (u)

      print (uu)

      u=" ...... ...... ...... ...... "

      print (u)

      print (uu)

      u=" Введите стрелу Прогиба "

      print (u)

      print (uu)

      b=0.00000000

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

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

      print (uu)

      u=" Введите длину Дуги "

      print (u)

      print (uu)

      L=0.00000000

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

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

      print (uu)

      u=" Подождите – идет расчет "

      print (u)

      print (uu)

      bb=0.000000 # Текущий прогиб

      Rt=L/Pii # Начальный Текущий радиус

      xx=b/100000 # Приращение радиуса

      xb=b/1000 # Допустимая ошибка..

      db=10000000*b # Начальное значение ошибки

      ugr=0.0000000

      # Уточнение подбором

      while db >xb:

      # Далее Cдвиг – четыре пробела в начале каждой строки

      Rt=Rt+xx

      ygr=(2*Pii)* (L/((Rt+Rt)*Pii)) # Угол раствора хорды радиан

      x=math.sin(ygr/2)

      y=Rt*x

      a=y+y # Хорда

      t=((Rt*Rt)-(y*y))

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

      bb=Rt-z # Прогиб

      x=(b-bb)*(b-bb)

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

      # Далее Cдвига нет..

      db=db+xb

      xx=b/100000000 # Приращение радиуса

      xb=b/1000000

      # Уточнение подбором ( вторая ступень )

      while db >xb:

      # Далее Cдвиг – четыре пробела в начале каждой строки

      Rt=Rt+xx

      ygr=(2*Pii)* (L/((Rt+Rt)*Pii)) # Угол раствора хорды радиан

      x=math.sin(ygr/2)

      y=Rt*x

      a=y+y # Хорда

      t=((Rt*Rt)-(y*y))

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

      bb=Rt-z # Прогиб

      x=(b-bb)*(b-bb)

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

      # Далее Cдвига – нет …

      # Конец подбора ........

      # Gsf = 461030/2725231222…

      R=Rt

      yg=ygr*180/Pii # Угол раствора хорды град.

      # Далее Cдвига – нет….

      x=a/(R+R)

      k=math.asin(x)

      au=(k+k)*180/Pii # Угол на хорду в градусах

      # Вывод