Механика \ Основы численной технологической механики

Определение напряженно-деформированного состояния простой конструкции при статическом действии нагрузок, страница 2

 3   .8783E-05   .9999E-01   .6884E-06  -.9998E-01   .5817E-06  -.3333E-01  -.1057E-05  -.2197E-08  -.6316E-08   1

 4   .8782E-05   .9999E-01   .6757E-06  -.9998E-01   .5568E-06  -.3333E-01   .1070E-05  -.3607E-08   .9071E-10   1

 5   .8782E-05   .9999E-01   .5909E-06  -.9998E-01   .7098E-06  -.3333E-01   .2335E-08   .1031E-05  -.6067E-08   1

 6   .8783E-05   .9999E-01   .6136E-06  -.9998E-01   .7150E-06  -.3333E-01   .1005E-07  -.1031E-05   .3105E-08   1

 7   .8783E-05   .9999E-01  -.1180E-08  -.9999E-01   .1281E-07  -.3333E-01   .3816E-08   .2942E-08  -.4119E-08   1

Как видно из результатов решения, нормальные напряжения по оси Y=R равны значению распределенной нагрузки q=0.1=.9998E-01 в двенадцатом узле:

(R****CИГMA 13 .9998E-01)  и других: 14, 15, 10 и так далее по столбцам (см. ПAPAMETPЫ OБЬEMHOГO HДC).


3.1.2. Решение задачи


Приведем описание решения задачи. Расчетная схема одной четвертой симметричной части амортизатора приведена на рис. 3.2.

Жирными точками  показана нумерация узлов, для элементов — это обведенные окружностью числа. Справа внизу показана нумерация степеней свободы по направлениям осей X, Y, Z. На верхнюю грань (элементы 1-3) действует равномерно – распределённая нагрузка q=P/Sверх , где Р - сила сжатия, а Sверх - площадь верхней поверхности амортизатора. Поскольку выделена четвертая часть амортизатора, то в плоскостях симметрии необходимо задать ограничения на свободу перемещений. Так, в плоскости X,Z перемещения всех узлов, принадлежащих ей, не могут происходить по оси X. Следовательно, все степени свободы этих узлов необходимо ограничить. Аналогично, в плоскости симметрии X, Y перемещения узлов, принадлежащих ей, не могут происходить по оси Z. Кроме того нижние узлы необходимо закрепить в пространстве плоскости X, Z задав ограничение на свободу движения узлов по оси Y.  Такими cтепенями свободы узлов будут следующие:

  1   22   34   91  112  124  125  128  131  134  137  140  143  145  157  181

193  194  197  200  203  205  207  210  213  216  219  222  225  226  228  231

234  237  238  240  243  246  249  252  255  258  261  264  267  270  273  276

279  282  285  288  291  294  297  300  303  306  309  312  315  318  321  324

327  328  329  330  332  333  335  336  338  339  341  342  344  345  347  348 .

Последовательность расчета

1.  Задаются общие параметры решения задачи

NO-кол. эл-тов; KYZ-кол. узлов для чтения координат; MAXY-количество узлов в  выделенном фрагменте конструкции; MXM­- количество точек интегрирования в ней; KUV- количество степеней закрепления; NMAT-количество материалов в объекте;  NELH-параметр нелинейности; MAT()-номера материалов по элементам конструкции (расчетной схемы).

READ(2,2)IWRITE,NO,KYZ,MAXY,MXM,KUV,NMAT,NELH,(MAT(I),I=1,NO)

2.  Чтение массива NBIB()-типов элементов и чтение массива MI(J,I)- матрицы индексов READ (2,28) ((MI(J,I),J=1,20),I=1,NO) ! 20 – количество степеней свободы элемента.

3.  Чтение-печать координат угловых узлов элементов по два узла сразу с 1-й строки (N1=KYZ/2). В цикле: DO 23 I=1,N1  READ(2,39)N,X(N),Y(N),Z(N),M,X(M),Y(M),Z(M) ! ЧТЕНИЕ WRITE(1,39) N, X(N), Y(N), Z(N), M, X(M), Y(M), Z(M)  ! ПЕЧАТЬ

4.  Чтение неподвижных степеней свободы, количество которых равно KUV): 199 READ(2,3)(NU3(I),I=1,KUV)! NU3-нoмepa cтeпеней свободы с нулевым значением перемещений  WRITE(1,51) KUV, (NU3(I), I=1,KUV).

5.  Задание свойств материала в точках интегрирования при решении упругой задачи, если решается нелинейная задача, то чтение этих данных пропускается, а NELH равно 0.    DO 41 M=1,NMAT       READ(2,40)Epro,Pyas,Plot,ALt  В этом операторе: Epro-продольный модуль; Pyas-коэффициент Пуассона; Plot–плотность; Alt–коэффициент линейного теплового расширения.

Скачать файл