Бидерман В.Л. - Механика тонеостенных конструкций (1070924), страница 68
Текст из файла (страница 68)
481 МЮ вЂ” коэффициент Пуассона; Т [1: 12) = ) ~, д, Н, Мы Мг, Ты Тч, о~ы огз, опю ог, ~) — полный вектор состояния, выводимый на печать. В блоке обратной прогонки для определения начальных параметров в точке 5К использованы идентификаторы 6, 61, 62, которые представляют собой главный и вспомогательные определители системы краевых условий. Идентификатор 6! повторно использован при обратной прогонке для переноса граничного условия ь = 0 ~ри 5 = 5К; !пои! ( В!Π— 2', ...) — процедура ввода числовой информации; !пои! ('Р2 — 10', ...) — процедура вывода числовой информации; сог! ('гЬг', 1, пг); ...
— процедура, реализующая уход на метку 1, выт: сод ('ге!игп') ~ нолнение действий до метки и и возврат к оператору, следующему за процедурой соН ('гЬг', 1, т). ТЕКСТ ПРОГРАММЫ Ьей!и !п1едег Н, К, Р, Н1, 5, 5К, 51, К1, Е1, 1„И, !ОС, Р!.ОС, Д, ж! 7,.7,6; геа! ХО„ХК, Х, Х1, Н, НС, ЕР5„ЕРШОВ, 70, РО, ТЕТА, 5Т, СТ, Р, ТОЛ, ЕН, ЕНЗ, Ф, Я!4, М!О, МЮ1, 6, 61, 62, А; аггау Т [1: 12); тои! ('Л!Π— 2', Н, К, Ж1, 5К, 51, К1, 'Р2 — 10'„Л', К, Н1, 5К, 51, К1); ж2 Ьея!п аггау 2 [1: Н, 1: К), ОМЕГА [1: К, 1: К], !У [1: Н, 1: 5! ], У [1: К1, 1: 5К!, ВЕТА [1: Л'1, 0: 5К), С [1: Л'1, 1: Л71), В [1: Л'!1, Ц, У, У!, У2, Р, РО, Р1, Р2 [1: Н); сопппеп! до с, !5, ОЛ1ЕГА, 9, Ф', 765, У, 250, ВЕТА, 102, С, 4„В, 2, жЗ О, У, У1, У2, Р, РО, Р1, Р2, 5; !пои! ('Я!0 — 2', ХО, ХК, 1, ЕР5, Р!.ОС, ЕО, РО, 6, МЮ, ЕН, ЕНЗ, ж4 ТОЛ, Я, )Р2 — !О', ХО, ХК, 1., ЕР5, Р!.ОС, !.О, РО, 6, МЮ, ЕН, ЕНЗ, ТОЛ, г); МЮ1:= 1 — МЮ ": 2; Н:= (ХК вЂ” ХО)/5К; !.1:= О; 5:= 0; Х1:= ХО; ж5 М1: 1ог 7:== 1 а1ер 1 нпИ Н до 1ог 7:= 1 а[ер 1 ппИ К г(о %' [1, КХ5+ 1]:= с [1, У]; И 5 = 0 1!зеп ао 1о М 2; О:=- 1; 1ог 1: 1 а!ер 1 пп!!1 К вЂ” 1 бо ж7 1ог 7:= 1 а1ер 1 пп!!1 К до Ьей!п У И, 5):=- ОМЕГА [7, у]; О:= ()+ 1 епд; Ч 5 = 5К !Ьеп до!о ПРОГОНКА; М2: ХО:= Х1; Р:= 1; ж8 МЗ: !1 Р = К Ьйеп Д:=- 1 е!зе Д =- О; 1ог 1:= 1 з1ер 1 ппИ Н до У [!) '-= с.
[! Р); ж9 ИНТЕГРИРОВАНИЕ: НС:= Н7Р7ОС; ЕОС:= 0; ж10 М4: Х:= Х1; 1ог 1:= 1 з(ер 1 пи!11 Н бо О [!):= У [!); сей ('гЬг', Л45, Мб); 1ог 7:= 1 а[ер 1 ппШ Л' до Ьей!п РО [Г] ".=, Р [7]; Н [!):= У! [7]:= У [!) + НСХРО [7]!3 епб; Х:= Х! + НС(З; сой ('гЬг', ~Ч5, Мб); 1ог 7:=-. 1 з!ер 1 нп1!1 Л! до Ьей)п Р! [!]:=- Р [!); (! [!]:= У1 [!):= У [!] -'- НСХ РО [!)!6+ + НСХР! [!)76 епд; сод ('гЬг', Мб, Мб); 1ог 7:= 1 а[ер 1 нп1!! Лг до Ьей!и Р! [!):= Р [!]; (! [!):=- У1 [!]:= У [!) [- НСх, РО [!)!8-1 3 Х НСХ Р1 [! ) 78 епд; Х:= Х1+ НС)2; сод ('гЬг', гЧ5, Мб); 1ог 1:= ! а[ер 1 ппИ Н до Ъед(п Р2 [1]:=- Р [!); (! [7[:=- У! [!);= У [1] + НСХ ГО [!)!2— ЗХНСХР1 [!)(2+ 2ХНСХР2 [!] епд; Х;=- Х1+ НС; сой ('гбг', М5, Мб); !ог /:= 1 з!ер 1 ип1!! Н-до Ьеа(п Р1 [/]:= — Р [/]; У2 [/]:= У- [/]+ НСХРО [/]/б-г.
2ХНСХ ХР2 [/]/3+ НСХР1 [/]/6 епд; И;= 1; 1ог /:= 1 в1ер 1 ип!!! Н до Ьеи!и ВИНОЛИ:= АВЗ (0.2Х (Н Е = 0 1/ У! [/] = О !Ьеп У! [/]— У2 (/] е1зе 1 — У2 [/]/У1 [/])); '!! ЕЙНОН ЕРЯ !Ьеп Ьец1п НС:= — НС/2; Р1.0С:= 2ХРЕОС; ЕОС:= 2Х ЕОС; до!о М4 епд; !! ЕННОНХ64 ) ЕРЯ !Ьеп И:= О епд; Х1:= Л1+ НС; /.ОС:=- ЕОС+ 1; !ог /:=- 1 е!ер 1 ип!1! И бо У [/]:= У2 [1]; И ЛОС<' РЕОС !Ьеп Ьед(п !! И = 1 Л, ЕОС =- (НОС вЂ”:2)Х2 Л РЕОС» 1 !Ьеп Ъеи!и НС:=- 2ХНС; 10С:= ЕОС вЂ”:2; РЕОС:= РЕОС вЂ”;2 епд; до!о М4 епй еЬе йо!о М7; М5: И.Е! = 1 Йеп ио1о Мб; А .'= МЮХСТ/Н; Р [1]:= ОХ ( — АХ (/ [1] — ЯТХ (/ [2] + МЮ1Х СТ/~ЕН1'НХ Х(СТХ(/ [3]+ ДХЯТХФ)); Р [2]:= ОХ( — АХ(/ [2]+ !2ХМЮ!/ЕН3/НХ(/ [4]); Р [3 ]:= ОХ (ЕН/В Х (/ [1 ] + А Х (/ [3 ] + ДХ (МЮХ Я Т/Л Х Ф— ЛХОЛВ Р [4]:= ОХ(ЕНЗХ СТ 72/12/ИХ(/ [2] + ЯТХ (/ [3] -'- АХ (/ [4] — ДХСТХФ); Р [5]:= ОХ( — МЮХЯТ/ИХ(/ [1]+ СТХ(/ [2]+ МЮ1Х ЯТ/ЕН/КХ (СТХ (/ [3] + ДХ ЯТХ Ф)); Мб: сос1 ('тешит') М7: !пг 1:= 1 з!ер !'ип!!1 Л/ до с, [/, Р]:= У [/]„ !! Р( К !Ьеп Ъефп Х1: =ХО) Р:= Р+ 1; ио!о МЗ епд; Я:= Я+ 1; ОРТОГОНАЛИЗАЦИЯ: !ог (~:= 1 а!ер ! ип!!1 К до Ьефп (ог /;= 1 е!ер 1 ип!!1 О йо Ьефп ОМЕГА [/, Я:= 0; !ог /:= 1 в!ер 1 ип!!! Л' йо ОМЕГА [/, Ц:= ОМЕГА [/, О]+ Я [./, 1]Х хг [7, О] епд; И 9» 1 !Ьеп Ьеи!и !ог /:= 1 з1ер 1 ип1!1 Π— 1 йо Ьеи]п ОМЕГА [О, 9]:= ОМЕГА [О, О] — ОМЕГА [/, (с] ! 2; !ог 7:= 1 к!ер 1 ип!!1 /1/ !(о Е (/, !]] = 3 [,/, О] — Е [У, /]ХОМЕГА [/ О] епд опав; ОМЕГА [О, О]:= ЯЩТ (ОМЕГА Ф И); И !] = К !Ьеп ОМЕТА,Щ (]]:= 1; !ог,,/:= 1 в!ер 1 ип!!! /1/ до Я [У, !',1]:=- Х [,/;О,]/!ОМЕГА [(], 9] еп4 КОНЕЦ ОРТОГОНАЛИЗАЦИИ; до!о М1; ПРОГОНКА: 51:= 1; О:= 01:= 02:= ВЕТА [1, ЯК]:= О!/О; ВЕТА [2, ЯК]:= — 02/О; до!о М9; М8: Я:=Я вЂ” 1; Х:=-Х вЂ” Н; !с:=- 1; 1ог /;= 1 е!ер ! ипШ К вЂ” 1 до !ог /:= / ь!ер 1 ип!!1 К до э11 е12 *13 е!4 е15 е16 е18 483 Ье8!л ОМЕГА [1, У]:= [г [Ц, 5+ 1]; (4:= О+ 1 епд; [ог У:=- К вЂ” 1 з!ер — 1 ил61 1 до Ье8!и ВЕТА [У, 5]:= ВЕТА [У, 5+ 1] — ОМЕГА [У, К]; !! У< К вЂ” 1 !Ьеп ЬеИ!л [ог 1:= К вЂ” ! 8(ер — 1 ип!!! У+ ! до ВЕТА [У, 5]:=.- ВЕТА [У, 5] — ОМЕГА [У, 1]ХВЕТА [1, 5] епд; ВЕТА [У, 5): —.
ВЕТА [У, 5)!ОМЕГА [У, У] елд; !ог 1:= ! 8!ер 1 ип!1! М до [ог У:=.1 8!ер 1 ил!!1 К до 2 [1, У]:= [)У [1, КХ 5+ У]; М9: [ог 1:= 1 8(ер 1 ип(!! 5! до Ье8[л У [1]:= 2 [1, К]; !ог У:=-. 1 8!ер 1 ип!!1 Л!! до У [1]:= У [1] -!'- ВЕТЛ [У, 5]Х7 [1, У] елд; !! 5 = 5К (Ьеп 01:=- У [5); сод ('252', М5, М6); Т [1]:= ЛОХУ [1]) Т [2]:= У [2]1 Т [3):= Р01101ЕХ У [3]1 Т [4]:=- Р01ЕХ У [4]; Т [5]:= МЮХ Т !4] + РОХЕН3Х СТ1121ЙХТ [2); Т [6];=- СТХ Т [3]+ 5ТХ РОХ Ф11.01ЕП Т [7]:= МЮХ Т [6) + РОХЕН110 421ЕХТ [1); Т [8):= Т [6)1ТОЛН.О; Т [9):= Т [7)1ТОЛ110; Т [1О]:= 6Х Т [4]1(ТОЛХЕО) ] 2; Т [11):= 6Х Т [5]1(ТОЛХЕО) [2; Т [12]:= ЕОХ(У [5) — О!); !пои! ('Р2 — 10'„Х, Т); И 5 = 0'!Ьеп сод ('8(ор'); 80 !о М8 епд епд; е20 ф21 е23 е24 э25 484 П р и м е ч а н и я: Ф! — перечень целочисленных н дробных переменных; Э2 — ввод и вывод границ массивов: зз3 — команда распределения памяти (указанные длины массивов позволяют для сястемы 5-го порядка назначать да 50 узлов ортогонализацни); Ф4 — ввод и вывод числовой информации; цсб — начало прямой прогонки; зуб — запись матрицы артагональных решений Х в блочную матрицу %', Ф7 — запись верхнего треугольника матрицы ОМЕГА в столбец Я матрицы у; Эб — начала блока интегрирования; Ф9 — выбор начальных условий интегрирования; 910 — начало программы ннтегрвровзния методом Кутта — Мерсона; Ф1! конец программы Кутта — Мерсона! выход иа метку М7 после интегриревания на шаге ортогоналнзацин Еи за!2 на метке Мз псременйые параметры и нагрузочные члены описываются для данной оболочки как функции независимой переменной Х (см.
примеры расчета в 5 18); йз13 — вычисление производных компонентов вектора состояния; ф14 — запись частного решения в соответствующий столбец матрицы Х: если матрица Х сформирована не полностью, то возврат на метку МЗ; в противном случае выполниетси ортогонализация столбцов матрицы 7; ' м15 — описание блока ортогоиализацни; Э! б — нормирование частных решений однородной задачи; м17 — вычисление определителей системы грани ~ных условий; зтн процедуры конкретизируются в зависимости от способа закрепления оболочки (см.
примеры расчета в 4! 8); Ф18 — вычисление произвольных постоянных в точке ЯК (в конце интервала интегрирования); Ф !9 — запись столбца Я + 1 матрицы У в верхней треугольник матрицы ОМЕГА; 320 — программа решения системы линейных уравнений ()(з+Г) С(з) == ь(з+1) (матрица П вЂ” треугольная); 321 — вычисление настоянных ннтегрирования в точке Я; за22 — выбор матрацы ортоганальпых решений 2 (з)- из блочной матрицы )Гг; 323 — вычисление вектора полного решения краевой задачи; Ф24 — вычисление полного физического вектора состояния Т [1: 121, включающего размерные величины перемещений, внутренних снл н напряжений; вт25 — печать вектора Т: Список литературы 1.
Абрамов А, А. О переносе граничных условий для систем линейных обыкновенн дифференциальных уравнений (вариант метода прогонки). — «Журнал вычислительной математики и математической физики», 1961, т. 1, № 3, с. 542 — 545, 2. Аксельрад Э. Л. Изгиб и потеря устойчивости тоякостснпых труб при гидростатическом давлении. — «Изв.
АН СССР. Сер. Механика и машиностроение», 1962, № 1; с. 98 — 114. 3. Алексеев С. А. Основы теории мягких осесимметричных оболочек. — В кн. Расчет пространственных конструкций. Вып. Х. М., Гасстройнздат, !965. с. 5 — 37. 4. Алексеев С. А.
Задачи статики и динамики мягких оболочек. — Труды Ч! Всесоюзной конференции по теории оболочек и пластинок. М., «Наука», !966, с. 28 — 37. 5. Алумяэ Н. А. Теория упругих оболочек и пластинок. — В кп.: Механика в СССР за 50 лет, т. 1П, М., «Наука», 1972. с. 227 — 266. 6. Амбарцумян С. А.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
