Раздел: Документация
0 ... 42 43 44 45 46 47 48 ... 87 etable.mziil 1 ,smisc, 12 etable,mziil2,smisc, 18 etable,mziil3,smisc,24 etable,mziil4,smisc,30 etable,mziil5,smisc,36 etable,mziil6,smisc,42 etable,mziil7,smisc,48 etable,mziil8,smisc,S4 etable,mziil9,smisc,60 etable,mzj,smisc,66 plls,mzi,mzj,3 etable,szi,ls,2 etable,szj,ls,32 plls,szi,szj,3 finish /exit ! построение эпюры изгибающих напряжений !.на верхи, поверхности балки ! окончание работы программы В результате работы программы в файле beam.grph будут содержаться четыре рисун. ка: деформированная и недеформированная формы балки (рис. 3.5 в), эпюра поперечных сил (рис. 3.5 г), эпюра изгибающих моментов (рис. 3.5 д) и эпюра изгибающих напряжений (рнс. 3.5 е). ш* за ш twt т. г Рис. 3.5, в АШ 30 2002 5TM5SS 8«в -I HSS MZJ Kttt -«20* -2Й» 21?» Ш lmz В FLOT BO. 5TEr*»J *ШК*1 ж .-.»«•* ssi ..я*»" SB ..и-*" СЗ ...я** -.n«»J ..IS"*? Puc. 3.5, e Как уже отмечалось выше, описанный метод построения эпюр имеет существе) недостаток: эпюры изгибающих моментов на участках с постоянной распределенной 1 „0й оказываются линейными, хотя известно, что оии должны быть представлены поэтичными параболами. Это объясняется тем, что графическое построение эпюр в ой программе опирается иа команду [plls] постпроцессора /POST1. Эта команда дей-Дтельно предназиачеиа для графического показа заданных характеристик в виде эпюр, ь"гоаДики этих характеристик по команде [plls] лииейиы в пределах одного конечного мента и строятся по значениям характеристик иа концах элементов. Кроме того, пользование опции KEYOPT(9) = 9 для элементов ВЕАМЗ при построении эпюр по коман-folls] ие Я3*1 Уточнения формы эпюры, т. к. для команды [plls] иет доступа к данным в вяти промежуточных точках, устанавливаемых по опции KEYOPT(9) = 9. 3.1.3. Определение реакций опор и усилий в стержнях плоской фермы Требуется определить опорные реакции и усилия в стержнях фермы, показанной иа рис. 3.6 а. Исходные данные: а — 2 м ; Pt = 10 кН; Рг = 20 кН. а)
6) -г пб ЮкН л 0 1 п2 е2 пЗ Рис. 3.6 Отнесем ферму к декартовой системе координат Оху (рис. 3.6. б) и рассчитаем координаты ее узлов nl, .... п7 (в м): nl— (0,0), п2— (2,0), пЗ — (4;0), п4— (6;0), п5 — {4;2/S), пб- {2;4/л[з), п7- {о;б/л[3~). Как известно, ферма — это конструкция, состоящая из стержней, соединенных между собой шарнирами (узлами фермы). Поскольку силы, действующие иа ферму, прилагаются в шарнирах, то все стержни фермы испытывают только осевые усилия растяжения (сжатия). Следовательно, подходящими конечными элементами для стержней ферм являются стержневые элементы LINK1 (в случае плоских ферм) и LINK8 (в случае простраиствеи-вых ферм). Тогда узлы конечно-элементной модели будут совпадать с узлами фермы, а каждый стержень фермы будет отдельным элементом модели. Данную задачу можно решать методами теоретической механики и тогда ие существенны никакие определяющие параметры стержней, кроме их длин. Однако для МКЭ элементы LINK1 (или LINKS) требуют задать свойство материала (модуль Юнга ЕХ) и константу элемента (площадь поперечного сечения AREA). Примем произвольно, что для *** стержней фермы ЕХ = 210П (Н/м2), AREA = 110~4 (м2). Однако отметим, что 1 Рассматриваемой статически определимой задаче значения этих параметров ие будут "Иять иа итоговые искомые величины. Ниже приводится листинг программы, предназначенной для решения описанной за-*чи в пакетном режиме /BATCH. /batch /com, ! файл plate girder /com, ! определение реакций опор и усилий в стержнях плоской фермы ! вход в препроцессор /show,girder,grph ! вывод графики в специальный файл girder.grph /title.plate girder ! задание заголовка plate girder е>, 1 Jink 1! выбор типа КЭ (стержневой КЭ linkl) и присвоение ему номера 1 г> 1,1е-4! произвольно назначенная площадь поперечного сечеиия ! элементов 1 (area =1-10 ) mp,ex,l,2el 1 ! произв. иазиачеииый модуль Юига материалов ! элементов / (ЕХ = 2е11) s3 = sqrt(3)! определяем узлы по их координатам иа плоскости n, 1,0,0$n,2,2,0$n,3,4,0$n,4,6,0$n,5,4,2/s3$n,6,2,4/s3$n,7,0,6/s3 С*** Здесь и далее: $ — знак, соединяющий несколько С*** отдельных операторов е,1,2$е,3,4$е,4,5$е,5,6$е,6,7$е,1,7$е,1,6$е,6,2$е,2,5$е,5,3 ! связываем узлы по их ! номерам элементами linkl текущего ! типа (type = I) с текущими наборами констант (real = 1) ! и свойств материалов (mat = У) finish! построение сетки КЭ (модели фермы) закончено /solu! вход в решатель /solution antype.static ! статический анализ конструкции ! ***** задание граничных условий: ******* ! определяем шарнирные опоры в узлах nl, п7 (задаем условия !закрепления) d,l,uy,0! смещение иу = 0 в узле 1 d,l,ux,0„7,6! смещение те = 0 в узлах nl и я7 ! задаем силы в узлах пЗ и пб f,3,fy,-2e4! fy = -20 кН в узле 3 f,6,fx,le4! fx= 10 кН в узле 6 solve! решаем систему линейных алгебраических уравнений finish! решение задачи закончено I ************************ /postl ! обработка результатов и их графическое представление Inapt ! подавление ненужного вывода set, first pldisp.l ! показ деформированной формы фермы n el =11 ! количество элементов /output,girder,res ! направляем вывод в файл girder.res ! ****** 1 вариант вывода результатов ******* prrsol! печатаем опорные реакции presol,smisc,l ! печатаем усилия в стержнях /output ! ****** 2 вариант вывода результатов ******** /pnum.node.O ! ие нумеровать узлы при графическом показе /pbc.u,, 1! показывать условия закрепления /vscale,,2, 1! установка масштаба показа векторов /pbc,f„l! показывать силы /pnum,elem, 1 ! нумеровать элементы eplot! нарисовать элементы etable,fi,smisc,l ! составление таблицы усилий в стержнях etable,fj,smisc,l plls,fi,fj! показывать эпюру усилий в стержнях *do,i, 1 ,n el! вход в цикл по элементам esel,s,elem,,i,i ! выбор элемента с номером / etable,forax,smisc,l ! выбор в таблице усилия в стержне i *get,n,eIem,i,etab,forax ! определяем п по указателю /gopr! активизация вывода *vwrite,i,n! вывод иа печать номера стержня i и усилия л (5x,,i = ,f4.0,5x,,n = ,610.4) /nopr *enddo ! конец цикла по элементам 0 ... 42 43 44 45 46 47 48 ... 87
|