3. Перфокарта кастройки иа задачу (715) Кол. 1-5* Числ) узловых точек 6-10* Количество элементов 11 - 15* Количество граничных точек 16-£0* Количества случаев нагружения 21-25* Клнчество степеней свободы = 2 Количество различных материалов О печать исходных данных
26-30" 31-35*
1 пропуск печати исходных данных
(NP)
(NE)
(NB)
(NLD)
(NDF)
(NMAT)
4. Перфокарты, содержащие информацию о свойствах материалов (110,2F10.2) (по одной на каждый .материал) Кол. 1 - 10* Номер материала 11-20 Модуль Юнга 21-30 Коэффициент Пуассона
(0RT(N,1)) (0RT(N,2))
5. Перфокарты, содержащие информацию о коордипатах (по одной на каждую узловую точку) (I10,2F10.0) Кол. 1 - 10* Номер узла 11-20 Координата X 21-30 Координата Y
(C0RD(N,1)) (C0RD(N,2))
Перфокарты, содержащие информацию об элементах (615) (по одной иа каждый элгмент)
Кол. 4-5 * Номер элемента (N)
6-10* i (N0P(N,1))
11 - 15* / \ связи между элементами (N0P(N,2))
16-20» т ) (N0P(N,3))
21-25 Не используются (N0P(N,4))
26-30 * Номер материала (IMAT(N))
7. Перфокарты, содержащие информацию р граничных условиях (215) (по одной иа каждое граиичЕгое условие)
Кол. 1-5* Номер граничного узла
6-10* 01 - закрепление в направлении Y
10 - закрепление в направлении X
11 - закрепление в обеих направлениях
8. Перфокарты, содержащие информацию о нагрузке (I10,2F10.2) (по одной иа каждую точку)
Кол. 1-10* Нпмрп v4.:Ta
11-20 21-30
Номер узла Нагрузка по оси X Нагрузка по оси Y
(NBCII)) (NFIX(I))
R(l)
R(2)
Примечание. Перфокарты, содержащие информацию о нагрузке, заканчиваются информацией о нагрузке в поаеднем ya.te незавиЪимо от того, задана в нем нагрузка или иет.
* Означает, что десятичная точка в числе не набивается, остальные числа обязательно должны содержать десятичные точки.
Образец печати исходных данных и результатов Колонки
Одна из основных проблем, возникающих при практическом использовании метода конечных элементов, связана с огромной информацией, получаемой в результате счета, и с большими затратами времени на обработку выводимых на печать результатов. В значительной степени эту задачу облегчает использование при воспроизведении результатов автоматического самописца.
Одним из самых удобных способов представления результатов является построение изостат. Они могут быть построены либо во всей области, либо только в некоторой ее части, представляющей особый интерес. В простейшем случае напряжения в элементе усредняются по узловым значениям. При этом изо-статы вычерчиваются в виде набора прямолинейных в каждом элементе отрезков, согласующихся с этими узловыми значениями. На фиг. 20.4 приведен пример [7]. Недостаточная гладкость кривых вызывает сомнение в правильности решения. Программы, позволяющие вычерчивать гладкие кривые иа основе, данных, полученных в отдельных точках, можно приобрести у предприятий, производящих эти самописцы.
Второй способ использования самописца состоит в вычерчивании для каждого элемента в заданном масштабе векторов главных напряжений в соответствующих направлениях. Такой способ представления напряжений показан на фиг. 20.5.
Большинство применяемых в настоящее время элементов дают скачкообразно меняющиеся от элемента к элементу напряжения, хотя напряжения в двух соседних элементах колеблются относительно истинного. Для сглаживания разрывов напряжений при их описании (и в некоторых случаях для улучшения точности) обычно используются два способа усреднения. Первый заключается в усреднении напряжений по двум смежным элементам. Например, из двух треугольников составляется четырехугольник и в качестве напряжения в некоторой точке четырехугольника принимается среднее по этим двум треугольникам значение. При втором способе усреднения суммируются напряжения всех элементов, соединяющихся в рассматриваемом узле, и сумма делится на число этих элементов. Этот способ обычно дает совершенно гладкую кривую напряжений, достаточно точно описывающую их распределение во всей области, за исключением граничных точек или областей с высоким градиентом напряжений.
Блок-схема. Пример программы. Приведенная ниже программа определения вектора напряжений представляет собой неза-
Сетка С
VuHmugaJlO
Сетка А
Фиг 204 Линии уровней моментов косого моста (из Ргос. Inst. Civ. Eng., 821,
Aug. 1967).
1 «тии!
Фиг. 20.5. Автоматическое вычерчивание главных напряжений подпрограммой вычерчивания векторов.
висимую подпрограмму считывания с перфокарт, содержащих напряжения, и вывода на печать в соответствующем масштабе векторов для каждого элемента.
Все такие подпрограммы являются частью набора Calcomp Fortran.Типичная блок-схема приведена на стр. 502.
Программа 20-8
С Подпрограм.ма вычерчивания ве:ктора для Calcomp DIMENSION IBUF(IOOO)
С С С С
Вывод номера барабана и определение места во внешней памяти
