Рассматривается статический расчет рамы в МКЭ-Ф. Также приводятся результаты расчета в студенческой версии ANSYS. Следует отметить, что в реализации МКЭ-Ф балочный элемент отличается от Beam188. Вероятно поэтому очень сильного сходства результатов для рамы добиться не удалось (разница в 5%). Необходимо также проверить программу на наличие багов и проанализировать задачу о консольной балке, хоть там и получились замечательные результаты в сравнении с формулами сопротивления материалов (рассмотрена тут).
Задача: Расчет рамы в МКЭ-Ф


Постановка во входном файле МКЭ-Ф
nodes 3 0,0,0 0.5,0,0 0.5,0.5,0 elems_113 2 113,1,2,0.0001,2e11,7800,0.333e-8 113,2,3,0.0001,2e11,7800,0.333e-8 bcfix 1 1,1,0,0,0 bcforce_stat 1 10,3,100,0,0
Результат


Увеличение числа узлов
Сделаем по одному промежуточному узлу.
nodes 5 0,0,0 0.25,0,0 0.5,0,0 0.5,0.25,0 0.5,0.5,0 elems_113 4 113,1,2,0.0001,2e11,7800,0.333e-8 113,2,3,0.0001,2e11,7800,0.333e-8 113,3,4,0.0001,2e11,7800,0.333e-8 113,4,5,0.0001,2e11,7800,0.333e-8 bcfix 1 1,1,0,0,0 bcforce_stat 1 10,5,100,0,0
Будет так:

Результат работы МКЭ-Ф почти не поменялся, только увеличилось число выводимых степеней свободы. Думаю это создаст трудности с выводом результатов при расчетах трехмерных рам. Не стоит выводить все подряд по всем узлам загромождает экран быстро.
