Сегерлинд Л. Дж. - Применение метода конечных элементов (1051193), страница 47
Текст из файла (страница 47)
0091*1 3 ЕГ(<(-6( Яйа/6 ° ООВЗ=1.3 5 15н<!.з>*<В<!» ° в<а»с(<> с(з>»/айя ти5ейт!ои ог есемеит Рйойейт<15 тито тне ссоваа Вт!Ггиевв ндтатх сот<=1.3 ! 1 - "И 5 < 1 > оотвз*» исс З5-"ЗСГ>(З 1> МРЧ!1 15 А</5>=А(/5>ЧЕГ(<> ОО>ТЗЕ<,3 ЗЗ=М5<З> ЗЗ=ЗЗ ПЧ1 <Г ! зз> 17, <7, 16 16 /5*/65Н> (ЗЗ 1> ИР>11 Я</5> А</В>ЧЕ59<»,З> 17 СОМТ1НОЕ 7 соиттмое с ноо<Г<саттои Аио !о<оттен ог тме ВТВтен аг еооат!ОИВ с СЯ11 ВОТЧЯ<<А</65МЧ1> ° А<ЗСГ+1>,НР ИВН,НС1> СА<1 ОСНРВО<а<ЗС5Н ° 1> ИР НВМ> са<1 51чво(а</65Н>1>.а<зсг>1>.а<1>.нй,нвм,мсс.<о»> с ° ° ° ' ° ° ° ° с саасосаттои ог тне !семен! йевостаи<5 с !МРО< ог тне е<енеит одтд ОО 66 кк=1.иЕ йеяв<ти 1> ме< М5 х1,71 ° хг т2,х»,та 100 155 110 115 129 ЗХ,13<Е1Я«5' 130 ЕС ЕНЕИТ 190 ЯООЕ 15ГЬ ° 2 ТОЯВОЕ <5Г1 .с 'с с с с с с с с с с с с с ЯЕТйгЕЧАЕ ОГ ТНЕ ИООАС ЧАЕОЕ5 002011 3 11*ИЗ< 1! 25 РНТ<1! А(11! САССНСАТТОИ ОГ ТНЕ ЗНЕАИ 51йЕЗЗ ЧЬ<ОЕ5 АИО ТН1$ТТИО НОНЕИТ Гсй ОНЕ ОЕСЯЕЕ ОГ ТН15Т Яйг <хгетзекзет!ех( тг-кгет(-кзехг-х( ° НЗ> НОНЕМТ НОНЕИТЕЯЬ2 ° <РН(<1>ЕРН1<2>ЕРНТ<З>И/3 ° В(1>*Т2 Тэ В(2>*ЧЗ Т1 в<э> т( тг с(<! кз-кг С<2! 51 Хэ С<3>ех2 Х1 /г ЗСГ+ие1 ЭЗ=/ОТ<НЕ<НЕЙ я(/г>ев.в Я</3> 0,0 оо 211 <.3 Я(/г> ° Я(/2>-В<1! РНТ(1>/Яйг 21 А</3! А</3>еС(1> РН1(1>/АМ2 ВЬ СОИ<(ИОЕ МОМЕНТ НОНЕИТ РСГ ТнЕЯТЯ=ТОЯООЕ/НОМЕЙТ САЯСОЙАТТОЙ АМО ООТРОТ ОГ ТНЕ ТВОЕ НООАЙ ЧАЙОЕ5 Ий!ТЕ(!0,30> Т1<ЙЕ ЗВ ГОННАТ(1Й1////1Хе2044//1Хе12НИООАС ЧАСОЕЗ ! 00311=1 НР 3< Я<1> Я(1>е(НЕЯТЬ МЯ1ТЕ<10 32! (1 ° А(1> ° 1е1,ЙР> Эг ГОЯНАт<(Хе<Зег<!.Тезке(Э;Е<9 ° Чеэхе<З,ЕАЯ.Т<ЗХ,ТЗ,Е(Я,ве 1! ООТРОТ ОГ ТНЕ 5НЕАй ЗТЯЕ55 ЧВСОЕ5 А</г! ° ТАО2Те Я(ЭЗ! Й ТАОТХ ИЯ1ТЕ( 10 55>(1(СЕ ЬВ ГОЯМЯТ«х,/////<1,20ЯЯ//(хе(ЬНЕ<ЕЙЕИТ 51255525//схе 5ЯН ! <Аигк .
' <аогт ТАОНВХ > ООЬЬ| 1 ИЕ /г=/СГ4 /3 /С/ей( ° 1 Я</2! А</2> ТНЕАТА А(/З>=Я</3> ТНЕЯТЯ тнАЯ=$0йт<А(/2! 2>я(33>еег> 1Г<ТНАХ <Т ТЯОНАХ! ООТО90 ТЯОМАК=ТМЯХ 'КЕе1 90 Ий<ТЕ(ТО,05> Т.А<33>,4(32<еТИЯК ЬЬ ГОЯМЯТ (Тх 13,3<5< Е12,9! > 1Г<1РСН,ТО,О! СОТОЬЬ ИМ<ТЕ<<И.ВЬ> <.Я</З!.4<32>,ТНЯК 59 ГОМНАТ<13 2Х 6112 ° Ь! Ьв СОНТ1МОЕ ий(ТЕ<10 ЬТ> ТНЕАТА ТАОНАХ КЕ 67 ГОЯНЯТ( //// <Х,ЯОИЯИСЙЕ ОГ ТИ15Т ГОЬ ТНЕ ЯРРСТЕО ТО 2,ЬИ ОЕСЯЕЕ5/1Х,ЯЬНМЯХ1МОИ 5НЕАЯ 5<ЯЕ55 Гой ТМЕ ЯРРЯ<ЕО 25 ° </(кеигННАХ1МОН 5НЕАй 5ТЯЕ55 <5 <ОСАТЕО <и Е<ЕИЕНТ19> $ТОР ЕИО Учебные программы о кручении стержня.
Набор заданных узловых значений искомой величины, считываемых подпрограммой В1ЭУЪ'А1., может оканчиваться либо нулем, либо целым отрицательным числом. Точно так же оканчивается карта узловых сил. В строке 150 программы предусмотрена возможность перфорирования результантов элемента (сдвиговых напряжений) для дальнейшего использования их при вычислении согласованных результантов элемента, т.
е. при вычислении узловых значений напряжений согласно методике, изложенной в гл. б. Программа СОИЯТК вычисления согласованных узловых значений описана в следующем разделе. 18.5. СОНАТЕ Программа СОХЯТК связана с реализацией теории согласованных результантов элемента при вычислении узловых значений результантов элемента.
Эта программа составлена для совместного использования либо с ТОКЯО)ч), либо БТКЕЫ (программа решения двумерных задач теории упругости), Программой предусмотрен ввод следующих перфокарт: титульной карты, карты параметров программы, стандартных результантов элементов, исходных данных элементов. Стандартными результантами элемента служат значения напряжений, отперфорированные программой ТОКЯ1ОМ или БТКЕЬЯ.
Исходные данные элементов в точности совпадают с исходными данными элементов, используемыми программой ТОК81ОХ или ЬТКЕ88. Пример комплекта исходных данных для программы СОгчЬТК не будет приведен здесь, так как перфокарты с этими данными (за исключением двух карт) перфорируются другими программа'ми. Эта программа может быть также использована для вычисления узловых значений результантов элементов в задачах течения грунтовых вод и переноса тепла. Однако программы для этих 'двух прикладных областей должны быть модифицированы так, чтобы они могли перфорировать выходные данные. Приведем обозначения некоторых наиболее важных переменных, которые используются в СОИБТК и которые не были определены ранее: ЯТК Согласованные результанты элемента.
Их число должно быть не больше 7 ЛСБТК Величина, которая определяет адрес последней ячейки памяти, отводимой для хранения стандартных значений напряжений в одномерном массиве А АКЕА Площадь отдельного элемента ЕКМ Вектор нагрузки для элемента ЕСМ Матрица демпфирования элемента 15 15 20 55 Ой!РОТ! Е1ЕНЕИТ 1МРО( оГ т!715 амо Раяанетей само ЯЕАО(1М 1> Т1Т(Е 1 гойнат(гааа! яеао((м.г> нй,ме,йем.нс1 г Гоянат(6!3! са(со!а!(ом оГ н>!М(БЯБ Амо (н(т(аа(гат!Сн оГ тме с /СГ=ИР ИСС 355й 3СГ 2 4ЕИО=/СБН>ИР йВМ /С5(й ЗЕИОАИЕ НСА 00131 1,/С5Тй 15 А(!> О.О 1йРот ОГ тне соихемт!ОНА! еАенеит йе501(Аит5 0061=1 МЕ ЯЕАО(1М А> ИЕЕ (5ТЯ(Л 3*1 МС!> А ГОЯНАТ(13,2Х ВЕ(2 ° 6> 006К-"1 НС( 41 /Ено>(К-1! ИЕ>ИЕА Ь А(/1> 5ТЯ(К! С(!СО(АТ10й ОГ ТНЕ ((ЕНЕИТ >С> НАТВ!Х Айо ТНЕ 11ЕНЕ ООАОКК 1 НЕ аеао((И,Б> ме(.НБ.К(,т (.хг.тг,(3,73 3 Гойнат(613,67(а.а! го ай!а (КЕР>з>хз 71>х( тг-хг ° 71-(3 тг-х( ° 73>/г.
0051=1(3 ООВК=1 НС!. /!-"ЗЕМО>(К-1> МЕ>КК В еян(1 к>=АРБА а(/1>/3, ООБ/ !.3 СОЕГ=1 0 1Г(1 ЕО /> СОЕГ 2. 5 есм(! л-"соеГ аяеА/12 1ИБЕЯТ!ОМ ОГ Е(ЕМЕМТ РЯОРЕЯТ1Е5 !НТО ТНЕ СО(ОММ ТЕС 0071=1 3 11*й5(1! ОО>5/ 1 МСС Л /БГ>$3-1>РМР>!1 15 а(35> А(/5>>ЕЯН(! 3> 0017/=1 ° 3 //=ИБ(Л 43 33 11 >1 1Г(33>17 17 16 16 /5>35$Н>(33-1! МР+!1 а(35> А(/5>!есн(1 3> 17 СОНТ1НОЕ 7 СОНТ1ИОЕ АО СОИТ1ИОЕ $010(ТОН ОГ ТМЕ БТБТЕН ОГ ЕООА(10МБ СА(Е ОСНРВО(а(/СБЧ ° 1>,ЧР,ИВМ> Са(С 51ТВО(а(/СБМ>1>,А(/СГ>1! А(1> НР ИВН МСА 10 БТОР ЕМО 0!ОНИ Вестой а МТ (Я> ИАТй1Х ТОЯ а( > Ряосяан БОМБ(й((ймотмоотйот.тай!Во Тмйот.(АРБББ* о1немБ(он есн(з.3>.еян(3,7>'.МБ(3>.5(й($!',А(6505> СОННОМ/САЕ/Т1(АЕ(20! оата (м/60/,го/61/,101/В/ с» ° С МР йОИВЕЯ ОГ ИООЕБ ° НЕ НОНВЕЯ ОГ ЕСЕНЕИТ$ С ИВМ - ВАНО М1ОТН ИСА МОНВГй ОГ $ТЯЕ$5Е5 РЕ с 363 Учебные нрограмгеье 18.6.
Н 0МСН Программа Г1.1)МСН может быть использована для анализа двумерных безвихревых течений идеальных жидкостей, которые рассматриваются в задачах о грунтовых водах. Программой предусмотрен ввод следующих перфокарт: титульной карты, карты параметров программы, карты с исходными данными элементов, узловых снл и заданных узловых значений искомой величины, считываемых при обращении к В1гУЧА1.
В списке обозначений не встречались следующие переменные; РХХ, РТУ Коэффициенты проницаемости в направлении осей х и у соответственно 1ГЕ1,Х, ЧЕ1'г' Компоненты вектора скорости течения в направлении осей х и у соответственно В табл. 18.3 представлены входные данные для задачи о течении грунтовых вод, обсужденной в равд. 9.2. Сеть разбиения на элементы показана на фиг. 9.3. Значения узловых сил определяются количеством воды, просочившейся вдоль русла реки и выкачанной в точках расположения двух насосов. Вдоль двух границ области заданы узловые значения искомой величины.
Карта с числом — 1 для номера элемента была использована для завершения работы программы (строка 84), так как результанты элементов не представляли интереса и поэтому не вычислялись. Программа вычисляет компоненты скорости, соответствующие течению грунтовых вод и прямолинейному потоку. Компоненты скорости произвольного потенциального течения не могут быть вычислены с помощью этой программы. 18.7. ТРНЕАТ Программа ТОНЕАТ вычисляет распределение температуры в двумерном теле, на границе которого либо происходит конвективиый теплообмен, либо задана температура окружающей среды.
Входные данные аналогичны тем, которые используются в предыдущих программах. Вводятся следующие перфокарты: титульная карта, карта параметров программы, карты исходных данных элементов, узловые силы и заданные узловые значения искомой величины, вызываемые подпрограммой ВггУЧА1, карты исходных данных элементов или карта с числом — 1. -Приведем обозначения наиболее важных переменных, которые еще ие были определены: КХХ, КУУ Коэффициенты теплопроводности в направлении осей х и у соответственно Н Коэффициент теплообмена О о С О О М О О о О О Ъ О о О о 1 О а * о о о о О О О о о о о о о О о о о о О О а о О о о О о а о о о О о о ОО о о О О О О о с О о Ю « о О о О О О о со ° ° о ! О О $ ч $~ Ю Ф Э) Ъ Ф х с.~ Ф И ФЬ ) 3 О % И Ю; ае 65 га 75 ва вы>*тг-тз в<г>=чз-$1 В<з> ч<-чг с<ы =кз-кг с<г>=к<-кз с<з>=кг-к> ааа=<кг тз>кз ° <1>к< ° тг-кг ° $1-кз ° тг-к! Рз) г.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
