Сегерлинд Л. - Применение метода конечных элементов (1050674), страница 45
Текст из файла (страница 45)
дятся против часовой стрелки) Коэффициенты, которые встречаютса при вычислении матриц элемента Учетверенная площадь элемента Удвоенная площадь элемента Величпна, значение которой определнет адрес последней ячейки памяти, отводимой лля (Ф) в одномерном массиве А ЕР РН! Т1 ТЕЕ ГО ГР НСЕ Х1. Х2, ХЗ Х(1), Х(2), Х(З) ТЕ 72, УЗ У(1), У(2), Ч (3) В н'С писываетсн в виде [Щь=[ОБД!)[ОР). Подпрограмма М(ГЕТЕР необходима для рещення иествцяонариых задач. У«аэыг Юьяра л Величина. значение которой определяет ад- рес послелней ячейки памяти для (р) в од- номерном массиве А Велпчнна, значение которой определяет ал- рес послелней нчейкн памяти, отводимой для [Х) в одномерном массвве А ЛЕНЕ 18.4.
Программа ТОКЕ)О)т) Программа ТОЕБГОГ4 вычисляет сдвнговые напряжения в упруп>м стержне иехругоеого сечения. Теоретической основой пргь граммьг служит материал, изложенный в рэзл. 6.1. При составле. иии программы попользовалась положения, обсужденные в гл. 7. Вхолные данные предсщвлеиы титульной картой, картой с параметрамн программы, походными даннымн элементов, значении. мн уэловьж снл н граничными значениями искомой величины, считываемыми подпрограммой ВПУЧАЕ, а также дополнительным комплектом походных данных элементов.
Титульная карта считывается в 18-й строке. а карта параметров програьгмы — в Ю-й строке. Переменные НР, НЕ, НВГЧ н 1РСН была определены выще в списке обозначений. Другие переменные, содержащееся в карте параметров, слцхующнег О Модуль сдвига, Н)смз 81„Длина стержня. см ТОЩИЕ Првкладываеыый крутящий момент, Н-см РСТ Удвоенное число осей симмнгрви. РСТ равно 1.6. если осей симметрии не существует Исходные данные элемента считываются в строке 46.
Формат этих данных совпадает с форматом лдииых, перфорнруемых программой ОЕ!О. Подпрограмма ВОРЧА!„обравгение к кпкрой происходит в строке 82, также имеет операторы ввода данных. Этн операторы обсуждаются в равд. !8.2. Исходныс данные элементе, необходимые для вычисления сдвнпжых напряжений, счп. тывзются в строке 92. Этн данные для элементов должны быть точной копией данныгч нспольауемых прн составлении глобальной матрицы жесткости.
Величины, укааываюпгие где ньгенно внутри одномерного мас. сава А хранятся узловые значения нсноыой величины, узловые силы н глобальная матрица жесткости, находятся в игроках 2Гь 27 и 28 Эти величавы были определены в списке обозначений. Входные данные длв варианта ьадачн о вручении, покззанного на фнг. 18.9, прелставлеиы полноспю, эа исключением нсхолных данных элементов, в табл. 18.2. Карта уаловых спл, которая считывается прн обращении к программе ВОЧЧАГ состоит только из данжо числа — 1.
так как узловые силы отсутствуют в задаче в О Й ФА СФ Ф Ю Ю Ю и $ Фа вв ~$ З" 'Г Тп Ю юаснйо Гсазв Фнонт.позвю. вв.тйайьь*зййзв,таппара ййзвв,папе ° 1 ННЮ1 твф~ф «11$ ея Фаа,э!,еазв,ааз °,сасз,о 11в ЫВ ТН ЬВГ ФОФЬ11 ФТГЬТГ ВГФГ ФОФГФГ СЕГСФ«ФВ ВЮГВ ВГ ЮЮНФФНО Я.аз ФФЕГТНВ1ОН О Н Е Саазйо. НВ Е Еаз с с 1 ' ''*''' ' Т С Сс 'Ф О *ГО+ Г С 11$ 5 11515 ю А СГ т«е 111 с» В Ге Фнезе са с НЕФОФФН ВТ111Е а с нмйаа Ф $ сааон* 1О в ознтеаа Фйо знпзмзиззсп ОГ т е сса й Естсн а О .«$ иа-а.
ЮФТГНТ ФГ Т111Е Ф В О ° »ЕЬО1ГСЗ нзтезпз ва тзтье $51 заноса а в $1«1' «ГГвв 1 ФЭГФ зьнав н .СФНА 5 Газ,ма вннкнь о Г*--':-' В$51 ю. 1 с ОГ 1не Ф.ОВФФ. $11 Гаезс 11$1 В сзаааь еОнсе НФ Взе с знаот Фнн ес о ФВ1нз ВГ еаенйнт пл а Сп Г па=а Е са снььтзс Г «» 111 НФ мз Г еы н нна ию еаевнз вви мтазс 5 Еа 11 ФФ ФВФИ ВФФФ $$11 ОФЕФФГФФЬ с с ЮААФвс ОГ заев» ннюпзззз з та вн сапааа св неаз Вьнв В=111 4 11 .Н 11 11 11 $1 ВФФЬФ ФФФФФ О1 ВН ФФ! Т СВ11 Е НООТГФС По НВ НАО11ай ОГ ТНЕ $Ф$% Ф ЕВАНА В с ыс потам. Фйсз 11 ФФАА Фь ао.нан,мьь ваВаазвн ПГ 1 Е Ф\ЕНЕ«Г йЕЬО\а ВВН ОГ Т Е ЕЬЕ"ЕНГ ВАТА ФВЬЬ Н1 Е Веаавн 11 нв. «$ $1 Гй ввт *1 та Уггбпэ гмьтмгмг щ "б 1 Еш с 11Е 1»Е (Вчй СОНЕТ)( ешиезм Фг нк мюж чн щз )э 111 »,з ы пеуп гпп 1 ЗП' »1»спгпщг г ыгщеп гг,-еф пк ге з фп «гонг ям ггыпвг внщк О»1»пг- г с»11 «(-и гп» 11 11 Г »11 ГЩПШЕ»Щ Е 1 гм.аз»пег «в Фнггог аг гве твн егзгь понг ггпепг 111 пг»е ° 1.ФЕЕ 1 зг гп гп»йрыпгг,гггмгы.аммм» Фгыке Э1 гегяггп 1 Е Ф 1 в 1»пгбп) 1(11(1111 ) е з Ф з е Ф з ж 11 влив с пк еыге Егаеи тб»1»м течгг ггмэг ° гщ 11 ° ° 1 Е*1 )»1 Ф, 1 1 е ° 1(гщжг(п 1„1 1 гы ше ыэВФ»~ мнме гы~ й г(- 11 пгг ~ггз»м %юг йй Ещ о крученш» стержня.
Набор ааданных узловь1х эяаченнй везомой величины, слпьшаемьи подпрограмгюй ВОТ'»»А1 может оканчиваться лнбо нулем, лнбо целым стрнцвтельным чнслом. Точно тан же онанчнгаепя карта угловых снл. В строке 1бб программы вредусмотрена возможность перфорнРовавня результннтов элемента (сдвиговых напряжений) для дальнейшего непользования пх прн вычнсленнн согласованных Результавтов элемента, т.
е. прн вычнслгнян узловых значений напряжезнй согласно методике. наложенной в гл. б. Программа СО)(БТК выгнслення согласованных узловых значений Описана в следующем разделе Программа СОР(БТЕ связана с ревлнзацней теорнн силасованиых результантсв элемента прн вычяслсння узловмх зкаченнй ре.
аультантов элемента. Зга программа составлена длв сонме»гпого нспользсвання лнбо с ТОК51ОУ(, лнбо БТЙЕББ (программа реше. ння двумерных задач теорнн упругсстн). Программой преду»иотрен ввод следующнх псрфшшрт: тнтульной карты, карти параметров при раммы, стандартных результзнтсв элементов, негодных данных элементов. Стандартнымн результантамн элемента служат значения напряжений, отперфарнрсванные программой ТОЙ51ОИ нлн БТКЕББ. Исхслные данные элеыентов в точности совпада»от г асг ходнымн Ланнымн элементов, нспользуемымн программой ТОК51ОВ нлн БТКЕББ. Прнмер комплекта исходных данных цля программы СОВБТК не будет приведен шгесь, так как перфснарты с этимн дзннымн (за нсключеннем двух карт) перфорнруются другнмн прсграммамя. Зта программа может быль также использована для вычнслення узловых значеннй результантов элементов в задачах течення грунтовых вод я переноса тепла.
Однако программы для этих двух прнклвдных Областей дснжны быть моднфяцнроааны таа. чтобы они моглн перфорнровать выходные данные. Прнзедеы обозначения некоторых нанбслге вав»ных переменных. кошрые попользуются в СОЙБТК н которые не былн опре. деланы ранее: БТК Согласованные результанты элемента. Их чнсло должно быть не больше 7 ЛСБТК Велвчнна, которая определнет адрес пгклезней ячейка памятн, Опюднмой для хранения стандартных значснвй напрягзений в одномерном масснве А Площадь отдельного элемента Вектор нагрузки дгя элемента МатРица ленпфнровання алемента Оату тагыыю1111 тыююю Е.—..::-=-": - '- =-- а юю а ю гата на юаюыгюю сан' с аыа11 1 гы 1 ю ганю 1\1м ам.ан.юы а ююгютюаз 1 ю ° \югюыыгпм а юаю сапам ютаюаю а еа111 111111 11 юп =1.1 ° аю 1 1 са 1 111«1.
п.пюа «аюаююнню а югл 11аюп смсююютпн м пе гыюпя аю ююаюю ын юю паям юю «ааааа а ююю 1 юю 11 йыыю 1 1 1 1 П 11 юю-ю 11-11 тю-ж 11111, 1 1 з 11111.л пяю ыююггг. а нзы ы» м юым 1 юмюп1гжз гата ыю ща 3 застав ° 1 1 а 1аы 11з сага 1 юса м юсю 11-11 а ы 11 юююу 111 юаню а аю1 1 1 . юа смюгт 11 ыю.аыю сю ня 11111 сю папан 1 Р ар таиаь111Л 1,11юп1юав 1111,ю юныас\ \зм 1Й 18.6. ИЮМСН Программа Р1утй4СН может быть нгпсльзована для анаюнэз двумерных беэвнхреаыз теченнй идеальных жндкосгей, которые рассматрвнаютсн а задачах о грунтовых водвх. Программой предусмотрен ввод следующих перфокарт: тнтульнсй карты, карты параметров програыыы, карты с нсхсднымн даннымн элементов, узлгаых пюл н заданных узловых эняченнй искомой велнчщпю, счнтывземых пря обращенпн к ВПТ4ГА1 В спнске обозначений не естречэлнсь следующие переменные: РХХ, РЪТ Коэффкцненты проницаемостн е напрянленнн осей х н у соответственнО .
'т'НЛ, 'ТЕСТ Компоненты вектора скорости течения в ияцравленнн осей х я у ссютветственно ' В табл. 18.3 представлены вхолные ланные для задачн о теченнн грунтовых вод, обсужденной в равд. 9.2. Сеть рзэбяення нн эдементы показана ня фяг. 9.3. Значения уэлспых снл опрелеляются количеством есды, просочювпейся вдоль русла ренн и выкаченной в точках расызложення двух языков. Вдоль двух гряннц сбласгн зеданы узловые значения нгкомой веля шны. Карта с чнглом — 1 для номера элемента была ясепзьтовака для завершения работы программы (сюрок» 84), твк кзк реэультвнты элеыентон не представляли ннтереса Н псзтому не вычнсляляы.
Программа вычисляет компоненты скорссгн, ссютвстствующне течению грунтовых вод и прямолинейному потоку. Компоненты скорости произвольного потенциального течення нс могут бьнь вычнсаены с поыощью этой программы. 18.7. ТОНЕЛТ Нрограмма ТПНЕАТ вычнсляет ряспрелеленне температуры в двумерном теле, н» границе которого лнбо пронсходнт конвектненый теплообмен, лнбо задека температура окружающей среды. Входные денные аналогичны тем, которые ясгюльзуются н пре. дыдущнх программах. Вйодятсн слелующне перфокэртыю тнтульная карта, кэрта параметров программы, карты походных данных элементов.
узловые силы я звдэкные узловые эначення искомой нелнчнны, вызываемые подпрограммой ВЮутЛЕ, нарты походных данных злемеюмов нлн карта с чнслом — 1. Приведем сбоэннчення нанболее важных переменных, юторые еще не былн определены: КХХ. КУУ Коэффнпненты тенлопроводноста в направления осей х н у соответственно Н Коэффнцнент теплсобмена ливлвтвн алтптаювю!ввтввт!Юапаа Ею т!Тсвюаа лл вп юата а ОО,ЕО вп спп,Еапат аевсапан ст и сю«тян. еп«пва ню - нкя в юиеаа тсйнтт 11.
е - нл лт ст пп тннз Р т - нте еастптт 1 $ с т Оейест1ОТ нен - ювллЕюй. п - и гветс!Гт Ен веяептен. $ ЕГРОГ тн тнт $$1$е сап $«6 т е св т«ОЮ Рлн е1еяз йе 611 нтт11$ае 3 санапгпм с Гввнаттюта.11.п!е.$1 *" спп$*таст О Онлеез лл 1 11птп Ю16м О тне слннй меспВ $1 т $3П1$ се вл От тя тттае а«в сета ейвт к тпепс л 1пае Р«тт и $ Гс ттае$1 !ей Гп л и 1Г$ ° т Г$ ОГ1Я Г $$. апе В 1 н От т е с$.6ем ОЮППЗЮ н пюм ° в Р.ийа ГФсе $1$и« С с $ н а«6 ес«О Р«1«$ и есе Г т е 1В е т к'1=1 $Р ает*$ $ «т и 1$ г-*г'юа- ю тг- е Рнп 5 Ет $1 31-$ с в1П В та Рт С11$ Ст Главе с $«вей $6« О е1ент«1 РвюРлитез т та т .е ляв н естВя йт 1 ° Е О5 1 ЮГЛТГ,Е П Г СРВ ЛР О т=пга $»$1 л '.
Ссн4«п с ллт 1саттса В О ТОЕОЮТ ОГ Гтв зтз е $5 ев ТЕО 3 сла ВО ие1ЮО6 1 и с 1 «Р л лат ОО аа 13 1 не Рспп !$$' е ! $. ° Е.та.ес. г п.п Вв ВР Гвй+Юайа йт«фею% Гй«ЕГЯГ е т юфтп1тт с Рплтсе Га«3 «ет. са ~ япяп н. е тп Глюп еамяа свт 1 в епееп Ва ВЫЕ1 1 ° $ Г ГЬ «1*«ыг» зз гн««)« ° «гг« с с««снщмсн сг «м «т«гст синчн«нгз с и» гг-«з з«г«=«» «г з<з««ы«з сы«=«эй г сгг««) «5 с«з «)-и «нг ««З « г « «« ю-нг.г«-«з гг-«« тзг тан«гз«*)йю» нп> нн«п гиг ««ж««г««с.ам нг«.«ж«.«г«« Й сю мЬ,«*,«,1«г э',««,к«г,эг «ю Ззэ «- Л « О я $ й «н ТП«)Р 1БГРЕ Температура окружающей срепы в точках поверхности конвективного теплообмева ЗГИ«Е Величина, которая определяет адрес псслеаней ячейки памяти, отводимой в одномерном массиве А для средних по элементу значений тсмперат)ры Стороны элемента, па которых происходят конвектнаные потери тепла В втой программе предусмсчрева возможность перфорации (строка !41) средних по элементу значений температуры, с тмч чтобы можно было исполыовать этп значения для анализа температурных напряжений с помощью программы ЯТЙЕ88 (программа решения двумерных задач теории упругости).
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
