Где найти полный тестовый пример расчета для программы "Авторяд ЕС"?

Нубий-IV · 26.08.2025, 07:07

В книге Проектирование и расчет многоэтажных гражданских зданий и их элементов/Под ред. П. Ф. Дроздова— М.: Стройиздат, 1986 есть исходники старой программы для расчета многоэтажных зданий - "Авторяд ЕС".

Написана она на стареньком Фортране-IV, но собирается современными компиляторами, и вполне себе работает. Только неизвестно - правильно ли. В книге точно есть опечатки, типа пропавших переменных, скобок, перепутанных +/-, +/*, 0/O, I/1, плюс что-то добавил я при наборе с экрана. Часть глюков поймал компилятор, часть выловилась на несимметричных ответах в симметричных тестовых схемах.

Сейчас на примитивных тестовых схемах "один столб", "два параллельно", "крест" и т.п. получаются правильные ответы для сил, моментов и перемещений. А вот усилия в закладных/перемычках выглядят подозрительно. И чтобы их проверить, надо сначала вручную посчитать более сложную многосвязную схему по формулам, размазанным на десятках страниц. В книге приведен пример расчета, но злобный редактор заставил сэкономить одну страницу, и половина таблицы связей заменена многоточием. Разумеется, ответы со связями, добавленными "по среднему", не сходятся с книжными. Плюс описание ввода данных явно не соответствует версии исходников программы - не те входные массивы; возможно, ответы тут и не должны совпадать.

Кто-нибудь видел в книгах/журналах хоть один комплектный расчет, чтобы с ним свериться без ручной возни?

Программа с исходниками и примером из книги - во вложении. Есть мелкие отличия от книжной версии:

Поправлены найденные опечатки
Убраны неиспользованные метки/переменные/параметры - чтобы не ругался компилятор
Старые функции чтения/записи файлов закомментированы и заменены на современные - чтобы оно вообще компилировалось
Процедура MINV обращения матриц, которая, видимо, была библиотечной в фортране-IV, честно украдена с одного из сайтов по старым компьютерам
В паре мест добавлены проверки деления на ноль, которые в оригинале почему-то есть для одной оси, но пропущены для другой - чтобы проходили крестообразные тесты
Слегка поправлено форматирование исходников и результатов расчета

Программа консольная, читает с клавиатуры и выводит на экран. Разумеется, лучше сделать перенаправление:

Код:

[Выделить все]

AUTORDES < IN.TXT > OUT.TXT

Код:

[Выделить все]

&PAR ! Начало блока данных

H  = 6   ! Высота здания, м
N  = 2   ! Число столбов
MM = 1   ! Число рядов связей сдвига
KK = 5   ! Число расчетных сечений по высоте
MTK= 2   ! Признак печати матриц

! Массив модулей упругости столбов, кН/м2
E  = 2*0.1E8

! Массив толщин столбов, м
D  = 2*0.2

! Массив минимальных  координат Y столбов, м
Y1 = -4.5 +1.5

! Массив максимальных координат Y столбов, м
Y2 = -1.5 +4.5

! Массив минимальных  координат Z столбов, м
Z1 = 0 0

! Массив максимальных координат Z столбов, м
Z2 = 0 0

! Массив вертикально распределенных нагрузок, кН/м
P0 = 10 -10

! Массив горизонтальных нагрузок
!       My   Qy(my)     qy   ay   z0    Mz   Qz(mz)     qz   az   y0   MK
!      кНм       кН   кН/м    -    м   кНм       кН   кН/м    -    м   кН
UVK =    0        0      0    1    0     0        0      0    1    0    0

! Массив связей сдвига
!     N1  N2  S[м/кН]  Ymax  Ymin  Zmax  Zmin
SVE =  1   2    1e-10   1.5  -1.5     0     0

&END ! Конец блока данных

Код:

[Выделить все]

*************************************************** ИСХОДНЫЕ ДАННЫЕ ****************************************************

  ВЫСОТА ЗДАНИЯ                              H   = 6.000    
  ЧИСЛО СТОЛБОВ                              N   = 2
  ЧИСЛО СВЯЗЕЙ СДВИГА                        MM  = 1
  ЧИСЛО РАСЧЕТНЫХ СЕЧЕНИЙ ПО ВЫСОТЕ ЗДАНИЯ   KK  = 5
  ПРИЗНАК ПЕЧАТИ МАТРИЦ ЖЕСТКОСТИ И НАГРУЗОК MTK = 2

          E  - НАЧАЛЬНЫЕ МОДУЛИ УПРУГОСТЕЙ СТОЛБОВ
          D  - ТОЛЩИНЫ СТОЛБОВ
          Y1 - МАКСИМАЛЬНЫЕ КООРДИНАТЫ СТОЛБОВ ПО ОСИ OY
          Y2 - МИНИМАЛЬНЫЕ  КООРДИНАТЫ СТОЛБОВ ПО ОСИ OY
          Z1 - МАКСИМАЛЬНЫЕ КООРДИНАТЫ СТОЛБОВ ПО ОСИ 0Z
          Z2 - МИНИМАЛЬНЫЕ  КООРДИНАТЫ СТОЛБОВ ПО ОСИ 0Z
          P0 - ВЕРТИКАЛЬНЫЕ РАСПРЕДЕЛЕННЫЕ НАГРУЗКИ ПРИЛОЖЕННЫЕ К СТОЛБАМ
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА СТОЛБОВ
------------------------------------
 |         |      1    |      2    |
------------------------------------
 |   E     | 0.1000E+08| 0.1000E+08|
------------------------------------
 |   D     | 0.2000    | 0.2000    |
------------------------------------
 |  Y1     | -4.500    |  1.500    |
------------------------------------
 |  Y2     | -1.500    |  4.500    |
------------------------------------
 |  Z1     |  0.000    |  0.000    |
------------------------------------
 |  Z2     |  0.000    |  0.000    |
------------------------------------
 |  P0     |  10.00    | -10.00    |
------------------------------------

          UVK-МАССИВ ВНЕШНЕЙ ГОРИЗОНТАЛЬНОЙ НАГРУЗКИ
------------------------------------------------------------------------------------------------------------------------
 |         |      1    |      2    |      3    |      4    |      5    |      6    |      7    |      8    |      9    |
------------------------------------------------------------------------------------------------------------------------
 |     1   |  0.000    |  0.000    |  0.000    |  1.000    |  0.000    |  0.000    |  0.000    |  0.000    |  1.000    |
------------------------------------------------------------------------------------------------------------------------
------------------------------------
 |         |     10    |     11    |
------------------------------------
 |     1   |  0.000    |  0.000    |
------------------------------------

          SVE-МАССИВ ХАРАКТЕРИСТИК СВЯЗЕЙ СДВИГА
------------------------
 |         |      1    |
------------------------
 |   I     |  1.000    |
------------------------
 |   J     |  2.000    |
------------------------
 |   S     | 0.1000E-09|
------------------------
 |  YIJ    |  1.500    |
------------------------
 |  YJI    | -1.500    |
------------------------
 |  ZIJ    |  0.000    |
------------------------
 |  ZJI    |  0.000    |
------------------------

************************************************* РЕЗУЛЬТАТЫ РАСЧЕТОВ **************************************************

  КООРДИНАТА ЦЕНТРА ЖЕСТКОСТЕЙ ПО ОСИ OY  EY=  0.000    
  КООРДИНАТА ЦЕНТРА ЖЕСТКОСТЕЙ ПО ОСИ OZ  EZ=  0.000    

          A - МАТРИЦА ЖЕСТКОСТЕЙ
                    ПОСЛЕДНИЙ СТОЛБЕЦ ОТНОСИТСЯ К БИМОМЕНТУ
------------------------------------------------
 |         |      1    |      2    |      3    |
------------------------------------------------
 |     1   | -3611.    |  3611.    |  0.000    |
 |     2   |  3611.    | -3611.    |  0.000    |
 |     3   |  0.000    |  0.000    |  0.000    |
------------------------------------------------

          BO-МАТРИЦА ВНЕШНЕЙ НАГРУЗКИ
------------------------------------------------------------
 |         |      1    |      2    |      3    |      4    |
------------------------------------------------------------
 |     1   |  0.000    |  0.000    |  5556.    |  0.000    |
 |     2   |  0.000    |  0.000    | -5556.    |  0.000    |
 |     3   |  0.000    |  0.000    |  0.000    |  0.000    |
------------------------------------------------------------

          N - НОРМАЛЬНЫЕ УСИЛИЯ В СТОЛБАХ ОТ ДЕЙСТВИЯ ГОРИЗОНТАЛЬНЫХ А ТАКЖЕ УДЕЛЬНО-НЕРАВНЫХ ВЕРТИКАЛЬНЫХ НАГРУЗОК
          T - БИМОМЕНТ(ПОСЛЕДНИЙ СТОЛБЕЦ)
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА СТОЛБОВ, КРОМЕ ПОСЛЕДНЕГО СТОЛБЦА
                    ПО ВЕРТИКАЛИ ОТЛОЖЕНЫ НОМЕРА РАСЧЕТНЫХ СЕЧЕНИЙ ОТ ВЕРШИНЫ ЗДАНИЯ ДО ОСНОВАНИЯ
                    ЗНАК - СООТВЕТСТВУЕТ СЖАТИЮ, ЗНАК + РАСТЯЖЕНИЮ
------------------------------------------------
 |         |      1    |      2    |      3    |
------------------------------------------------
 |     1   |  0.000    |  0.000    |  0.000    |
 |     2   |  1.153    | -1.153    |  0.000    |
 |     3   |  2.307    | -2.307    |  0.000    |
 |     4   |  3.461    | -3.461    |  0.000    |
 |     5   |  4.615    | -4.615    |  0.000    |
------------------------------------------------

          Q(I,J) - ПЕРЕРЕЗЫВАЮЩИЕ УСИЛИЯ В СВЯЗЯХ СДВИГА
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА РАСЧЕТНЫХ СЕЧЕНИЙ ОТ ВЕРШИНЫ ЗДАНИЯ ДО ОСНОВАНИЯ
------------------------------------------------------------------------
 |         |      1    |      2    |      3    |      4    |      5    |
------------------------------------------------------------------------
 СВЯЗЬ 1- 2|  29.03    |  26.18    |  17.25    |  7.147    |  0.000    |
------------------------------------------------------------------------

          MY - МОМЕНТЫ В СТОЛБАХ ОТ ИЗГИБА В ПЛОСКОСТИ OXY
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА СТОЛБОВ
                    ПО ВЕРТИКАЛИ ОТЛОЖЕНЫ НОМЕРА РАСЧЕТНЫХ СЕЧЕНИЙ ОТ ВЕРШИНЫ ЗДАНИЯ ДО ОСНОВАНИЯ
------------------------------------
 |         |      1    |      2    |
------------------------------------
 |     1   |  0.000    |  0.000    |
 |     2   |  3.460    |  3.460    |
 |     3   |  6.920    |  6.920    |
 |     4   |  10.38    |  10.38    |
 |     5   |  13.84    |  13.84    |
------------------------------------

          MZ - МОМЕНТЫ В СТОЛБАХ ОТ ИЗГИБА В ПЛОСКОСТИ OXZ
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА СТОЛБОВ
                    ПО ВЕРТИКАЛИ ОТЛОЖЕНЫ НОМЕРА РАСЧЕТНЫХ СЕЧЕНИЙ ОТ ВЕРШИНЫ ДО ОСНОВАНИЯ
------------------------------------
 |         |      1    |      2    |
------------------------------------
 |     1   |  0.000    |  0.000    |
 |     2   |  0.000    |  0.000    |
 |     3   |  0.000    |  0.000    |
 |     4   |  0.000    |  0.000    |
 |     5   |  0.000    |  0.000    |
------------------------------------

          VY   - ПРОГИБЫ ЗДАНИЯ В ПЛОСКОСТИ OXY
          VZ   - ПРОГИБЫ ЗДАНИЯ В ПЛОСКОСТИ OXZ
          TETA - УГОЛ ЗАКРУЧИВАНИЯ ЗДАНИЯ В РАДИАНАХ
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА РАСЧЕТНЫХ СЕЧЕНИЙ ОТ ВЕРШИНЫ ДО ОСНОВАНИЯ
------------------------------------
 |         |      1    |      2    |
------------------------------------
 |  VY     !-0.3691E-04|  0.000    |
------------------------------------
 |  VZ     !  0.000    |  0.000    |
------------------------------------
 |  TETA   !  0.000    |  0.000    |
------------------------------------

          ПРОДОЛЬНЫЕ УСИЛИЯ В СТОЛБАХ ОТ ДЕЙСТВИЯ ВЕРТИКАЛЬНЫХ И ГОРИЗОНТАЛЬНЫХ НАГРУЗОК
                    ПО ГОРИЗОНТАЛИ ОТЛОЖЕНЫ НОМЕРА СТОЛБОВ
                    ПО ВЕРТИКАЛИ ОТЛОЖЕНЫ НОМЕРА РАСЧЕТНЫХ СЕЧЕНИЙ ОТ ВЕРШИНЫ ЗДАНИЯ ДО ОСНОВАНИЯ
------------------------------------
 |         |      1    |      2    |
------------------------------------
 |     1   |  0.000    | -0.000    |
 |     2   |  13.85    | -13.85    |
 |     3   |  27.69    | -27.69    |
 |     4   |  41.54    | -41.54    |
 |     5   |  55.39    | -55.39    |
------------------------------------

************************************************************************************************************************

Так что заполняем по образцу исходный файл, запускаем программу, смотрим ответы. Подробности в книге.

Программа собирается компилятором gfortran из комплекта GCC.

Просто компилятор можно взять тут:
https://github.com/skeeto/w64devkit
Кто захочет поотлаживаться (точки останова, просмотр переменных, пошаговое выполнение) - можно поставить Code::Blocks, у него gfortran идет в комплекте:
https://www.codeblocks.org/

Если удастся выловить всех глюков - тема может переехать в готовые программы, на память.

Нубий-IV · 26.08.2025, 09:49

Если кто вживую видел ЭВМ ЕС (например, не я) - можно поностальгировать.

Есть эмулятор старого железа Hercules:
https://github.com/SDL-Hercules-390/hyperion
Есть эмулятор старых терминалов:
https://x3270.bgp.nu/
Был сайт истории компьютеров it.museum-ic.ru. Сайт помер, но еще доступен в интернет-архиве. В его разделе загрузок была ссылка на образы софта под ЕС, она пока жива:
https://disk.yandex.ru/d/_vBDcqgNFjWTkw

Как запустить:

Качаем всех троих, распаковываем.
Из софта для ЕС выбираем версию по вкусу - например, cbm401.zip, - распаковываем.
В папке CBM-R401 в файле конфигурации CBM2165.cnf подправляем пути к образам дисков/лент.
Запускаем эмулятор с файлом конфигурации:
Код:
[Выделить все]
```
c:\Programs\Hercules\hercules.exe -f c:\Programs\Hercules.OS\CBM-R401\CBM2165.cnf
```
Запускаем первый терминал wx3270.exe - он изображает оператора ВЦ.
Подключаемся к localhost:3270.
В консоли эмулятора даем команду на загрузку с диска 240:
В консоли оператора отвечаем на вопросы системы при загрузке:
Запускаем второй терминал wx3270.exe - он изображает пользователя.
Подключаемся к localhost:3270.Логинимся как SYSPRG/SYSPRG
Что дальше - должен помнить пользователь ЭВМ ЕС.
Для остальных - руководства есть в интернете.

Еще есть готовая сборка, но с оригинальной буржуйской системой MVS:
https://www.prince-webdesign.nl/inde...3-8j-turnkey-5

Распаковать
Запустить mvs.bat, дождаться загрузки до логотипа
Запустить wx3270.exe, подключиться к localhost:3270
Залогиниться как HERC01/CUL8TR
Дальше см. канал moshix.
Например, Fortran compilers on IBM MVS 3.8 (по нему я запускал старый фортран-IV, для сравнения с современным):

Похожие темы
Тема	Автор	Раздел	Ответов	Последнее сообщение
Какой язык перспективен для инженера-конструктора с условием	The_Mercy_Seat	Программирование	705	17.03.2021 14:19
Помогите найти пример расчета заанкеренной шпунтовой стенки	Агиделька	Поиск литературы, чертежей, моделей и прочих материалов	1	26.03.2013 16:30
Необходим пример расчета прочность и устойчивость водоводов и др.	snipi	Поиск литературы, чертежей, моделей и прочих материалов	0	21.07.2011 08:51
где можно найти пример расчета сопряжения подкрановой металлической балки с жб колонной???	нага	Поиск литературы, чертежей, моделей и прочих материалов	0	03.06.2011 15:14
Помогите пожалуйста найти пример расчета	light	Конструкции зданий и сооружений	2	22.05.2011 10:02