Напряженно-деформированное состояние рулонированных оболочек тема автореферата и диссертации по механике, 01.02.03 ВАК РФ
Корх, Алексей Иванович
АВТОР
|
||||
кандидата технических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Киев
МЕСТО ЗАЩИТЫ
|
||||
1984
ГОД ЗАЩИТЫ
|
|
01.02.03
КОД ВАК РФ
|
||
|
ВВЕДЕНИЕ.
1. РАЗРЕШАЮЩИЕ СООТНОШЕНИЯ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ЦИЛИНДРИЧЕСКОЙ РУЛОНИРОВАННОЙ ОБОЛОЧКИ.
1.1. Постановка задач
1.2. Построение разрешающей системы уравнений
1.3. Решение системы интегральных уравнений методом коллокации
1.4. Проскальзывание слоев
1.5. Начальная погибь
1.6. Предварительный натяг'витков оболочки.
1.7. Предварительная опрессовка оболочки
1.8. Оболочка с витком,ослабленным щелевым разрезом.
2. ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ СМЕШАННОЙ ЗАДАЧИ ПЛОСКОЙ
ТЕОРИИ УПРУГОСТИ ДЛЯ ИЗОТРОПНОГО ТЕЛА.•
2.1. Система интегральных соотношений,связывающих граничные значения напряжений и перемещений.
2.2. Особенности решений на концах участка интегрирования.
2.3. Алгебраизация системы интегральных уравнений для кусочно-прямолинейной границы
2 Л. Алгебраизация интегральных соотношений для кусочно-круговой границы
2.5. Формирование системы линейных алгебраических уравнений . ^
2.6. Определение компонентов напряженно-деформированного состояния в области
3. ЧИСЛЕННАЯ РЕАЛИЗАЦИЯ ПРОЧНОСТНЫХ ЗАДАЧ . Ш
3.1. Алгоритмизация расчета напряженно-деформированного состояния оболочки . . . ш
3.1.1.Математическое обеспечение . III
3.1.2.Апробация численно-аналитической методики расчета
3.2. Алгоритмизация смешанной задачи плоской теории упругости для изотропного тела
3.2.1.Программная реализация метода граничных интегральных уравнений.Вычислительный комплекс K0NTAK.T.
3.2.2.Решение тестовых задач
4. НАПРЯЖЕННО-ДЕФОРМИРОВАННЫЕ СОСТОЯНИЯ МНОГОСЛОЙНЫХ РУЛОНИРОВАННЫХ ОБОЛОЧЕК
4.1. Исследование влияния трения и толщины листового проката на напряженное состояние оболочки.
4.2. Влияние не целого числа витков на напряженно -деформированное состояние
4.3. Учет технологической погиби
4.4. Предварительный натяг витков рулонированной стенки.
4.3. Опрессовка оболочки
4.6. Рулонированная стенка,ослабленная щелевым разрезом.
4.7. Расчет ловушки
4.8. Определение коэффициентов интенсивности напряжений.
Одной из важнейших задач народного хозяйства в одиннадцатой пятилетке является снижение материалоемкости и повышение надежности и эффективности сооружений и конструкций.
В области естественных и технических наук сосредоточить усилие на решении следующих важнейших проблем: развитие математической теории, повышение эффективности ее использования в прикладных целях.". (Из материалов ХХУ1 съезда КПСС,стр.146).
Решение этих задач предусматривает использование современных конструктивных форм, учет реальных условий эксплуатации и физических свойств материала конструкций, что приводит к необходимости рассмотрения более сложных математических моделей, описывающих работу конструкций или их элементов.
К числу наиболее распространенных конструкций в различных отраслях техники относятся тонкостенные трубопроводы. Те из них, которые работают в условиях высоких внутренних давлений,как правило, изготовляют из многослойных труб,имеющих ряд преимуществ перед однослойными с монолитной стенкой ОМ, [1*01.
В Основных направлениях экономического и социального развития СССР на 1981-1985 г.г. и на период до 1990 года, принятых ХХУ1 съездом КПСС, перед производственниками и учеными поставлена задача по увеличению объема добычи и транспортировки природного газа. Решение этой задачи связано с созданием магистральных трубопроводов большого диаметра (1022 мм и более),способных работать в условиях высоких рабочих давлений.
Один из конструктивных вариантов таких труб предложен в Институте электросварки им.Е.О.Патона АН УССР. Это многосекционная конструкция,состоящая из спирально-свернутых четырех-,пяти-слойных оболочек с внутренним диаметром 1420 мм и общей толщиной стенки 16-20 мы, сваренных между собой кольцевыми швами СП . Разработаны варианты и однослойных труб с многослойными рулони-рованными вставками СРаботы по повышению надежности маги-стальных трубопроводов из рулонированных оболочек проводит и ВНИИСТ по строительству магистральных трубопроводов.Опытные партии таких труб уложены в Западной Сибири и на Украине [101 .
В химической, нефтехимической,нефтеперерабатывающей и ряде других отраслей промышленности в технологических процессах,связанных с применением оборудованияработающего при высоких давлениях и температурах,применяются сосуды высокого давления,изготовленные с применением рулонирования. Их внутренний диаметр достигает 5000 мм, а рабочее давление 500-700 МПа Г42.11.
Многослойные рулонированные конструкции имеют ряд преимуществ перед однослойными. Их отличает сравнительно простая технология, возможность изготовления изделий практически любых размеров, существенное сокращение трудозатрат и снижение материалоемкости. Наиболее важными с точки зрения практических приложений являются такие качества рулонированных конструкций как большая стойкость к хрупкому разрушению и развитию магистральных трещин.
Применение таких оболочечных конструкций в технике требует исследования их прочностных свойств и изучения влияния различных факторов на них.
Обширный круг идей и практических результатов по общей теории многослойных оболочек и пластин нашел отражение в работах С.А.Амбарцумяна Г?1 ,[81 , А.Я.Александрова £31 , В.В.Боло-тина Ш1, И.А.Буякова [2,11 , А.Т.Василенко И 2,2,1,С 381 , В.В.Васильева 1, А.С.Вольмира 12.71» И.И.Воровича [» А.Л.Гольденвейзера 1зо] , Э.И.Григолюка и П.П.Чулкова
Ъ?1 , А.М.Григоренко, А.Н.Гузя Гз91 , Л.Либреску , Н.А.Кильчевского, Г.М.Куликова С6Э1 , В.Н.Москаленко и
Ю.Н.Новичкова С?91, Ю.Н.Новичкова £861 , Ю.В.Немировского и А.Н.Андреева [91 , В.Г.Пискунова 14СИ1 , А.П.Прусакова 1103"1, А.О.Рассказова [4051, А.Ф.Рябова 11091,[Ц01 , А.С.Сахарова [1161,
A.Г.Терегулова 1>М81, В.Й.Феодосьева [12141, Л.П.Хорошуна
B.Е.Чепига ЫП , Ю.Н.Шевченко [12 81 ,[551, В.М.Шопа и др.
Непосредственное применение теории и методов расчета слоистых оболочек к расчету спирально-свитых, т.е. замена рулонирован-ной стенки слоистой не позволяет выявить характерные особенности деформирования рулонированной оболочки. Существенное отличие между слоистой и рулонированной оболочками состоит в том,что в многослойной стенке рулонированной конструкции из-за несимметричности слоев (их неконцентричности) и наличия концентраторов в местах крепления навивки при нагружении внутренним равномерно распределенным давлением,кроме нормальных сил контактного взаимодействия между слоями,возникают также и касательные. Между слоями действуют силы кулоновского трения [681, при этом вдоль линии контакта возможны участки как со сцеплением,так и с проскальзыванием, причем положение и размеры этих зон заранее неизвестны. При отсутствии трения в такой системе внутренний пакет слоев не напряжен и работает как заполнитель,что,естественно,не свойственно многослойным в традиционном смысле оболочкам.
Кроме того,за счет смятия микро-,макрошероховатостей поверхностей контакта,а также за счет опрессовки неплотно свитых участков,деформативность оболочки в процессе нагружения увеличивается, что сказывается и на ее напряженном состоянии. Зависимость контактной нормальной податливости от сжимающих нормальных напряжений, действующих между слоями,нелинейна и определяется,как правило, экспериментально [711^761.
Проявляется влияние и ряда других факторов,влияющих на напряженно-деформированное состояние оболочки: преднапряжения навивки, начальной погиби, щелевых дефектов,температуры.
В ряде известных теоретических исследованиях,посвященных расчету рулонированных конструкций,авторы принимали различные расчетные модели с привлечением упрощающих гипотез. Упрощения расчетной модели продиктованы сложностью взаимодействия между слоями и недостатками применяемого расчетного аппарата.
Из-за сложности замеров деформаций в промежуточных слоях не получили развития и экспериментальные работы. Так,в работах, описанных в Е№, [96] ,[12.5], замеры деформаций проводились по свободным торцам слоев и наружной поверхности внешнего слоя.
Ю.М.Сагидаевым [1141 проведены экспериментальные исследования эволюции напряженно-деформированного состояния модели руло-нированного сосуда вплоть до его разрушения,опирающееся на данные тензометрирования наружных слоев.
М.И.Бураком [2.01 напряженное состояние слоев оболочки оценивалось по данным замеров деформации на поверхности небольшого сквозного круглого отверстия в стенке оболочки,расположенного вдали от продольных сварных швов и закрытого накладкой со стороны внутренней поверхности оболочки.
В работе В.Я.Баш, В.А.Коваленко и В.В.Пискуна экспериментально исследовалась модель рулонированной оболочки в виде тонкой трубы, скрепленной фольгой, на поверхности которой наносилась сетка и впоследствии по ней измерялись деформации. Полученные авторами результаты сопоставлялись с результатами решений по методике расчета тел вращения Ю.Н.Шевченко, М.Е.Бабешко В.В.Пискуна £1281 без учета проскальзывания.
Исследованию напряженного состояния спирально-многослойных оболочек поляризационно-оптическим методом на моделях из оптически чувствительного материала ЭПСА посвящены работы В.И.Савченко, Т.Ю.Кепича, А.М.Нуркиянова и К.И.Мазура 150"] ,\Л11 , ,
З] в Получены экспериментальные зависимости для контактного сближения поверхностей слоев из этого материала за счет макрошероховатостей и величины коэффициентов трения [ИЗ]. Определены усредненные кольцевые напряжения в слоях навивки для трех-,пяти-слойной моделей при некоторых значениях коэффициента трения [511, [87]. Рассматривалась и концентрация напряжений возле технологических сквозных отверстий в трехслойной оболочке и коэффициенты интенсивности напряжений в устьях щелевых дефектов [87].
Результаты экспериментальных исследований рулонированных оболочек носят частный характер. Даже в совокупности они не дают возможности сформулировать какие-либо общие количественные закономерности»описывающие наиболее важные стороны процесса деформирования этих конструкций.
В ряде исследований к расчету рулонированных конструкций применен приближенный подход, не позволяющий учесть ряд факторов, существенно влияющих на напряженное состояние. Так,в работах А.Ф.Безвербного , Г.С.Писаренко и О.Е.Богинича [100], Ю.М.Сагидаева [ЦБ] , Б.В.Яблонского \да"3 , а также [89] ,[90] рулокированная оболочка заменяется многослойным длинным цилиндром. Контактное сближение поверхностей слоев за счет макрошероховатостей не учитывается. В такой модели силы трения между слоями отсутствуют, и поэтому радиальные и кольцевые напряжения расчитывались по формулам Ламе [ИЗ] , а в £43] применялась и теория замкнутых тонкостенных оболочек. Аналогичная модель руло-нированной оболочки исследовалась авторами Л.И.Гурьевой [401 , М.А.Нетребским, Ю.М.Сагидаевым и Г.В.Раевским [83], Д.Драгое \\Ъ7,,с. 167-180 ] , с учетом предварительного натяга навивки.
Г.С.Писаренко и А.Е.Бабенко [901 ,[99], П.А.Вислобицким [^ рассмотрена безмоментная плотно свитая бесконечная спиральная оболочка, работающая на растяжение. Касательные силы взаимодействия предполагаются равными кулоновским силам трения на границах контакта всех слоев,что возможно при малых коэффициентах трения, когда в контакте имеет место проскальзывание. В расчетной модели не учитывается контактная податливость слоев. Для получения разрешающих уравнений используются уравнения без-моментной теории цилиндрической оболочки, толщина слоев полагается пренебрежимо малой по сравнению с радиусом кривизны витков. П.Г.Пимштейн №1 также рассматривает расчетную модель сосуда высокого давления в виде безмоментной спиральной оболочки,в которой в пределах всей зоны контакта имеет место проскальзывание. Напряженно-деформированное состояние внутренней обечайки определяется как для замкнутого кольца по формулам Ю.А.Шиманского, а для расчета навивки используются уравнения равновесия элемента цилиндрической оболочки в форме В.В.Новожилова. В расчете учитывается контактная податливость микро-,макронеровностей поверхностей слоев. Аналогичный подход к расчету рулонированной оболочки имеет место в исследованиях П.П.Бородавкина и А.М.Синю-кова , с. 148-1561 .
Исследованию влияния величины нормальной контактной податливости микро-,макронеровностей поверхностей слоев на напряженное состояние многослойной оболочки в целом посвящены исследования П.Г.Пимштейна и В.Н.Жуковой [951 , Е.Г.Борсука [971. Зависимость контактного сближения поверхностей от нормальных усилий для листового проката толщиной 6 мм, полученная экспериментально, и эмпирическая зависимость контактной теплопроводности от контактного давления приведены в [951 . Напряженное состояние оболочки от действия внутреннего давления и температуры определено для многослойного цилиндра по формулам Ламе. В работе [972 приводятся экспериментальные данные о сближении поверхностей слоев и изменении обьема межслойного пространства на единичной площади в зависимости от контактного давления и величины силы натяга навивки.
Модель рулонированной оболочки в виде плотной многолистной цилиндрической поверхности с постоянным радиусом кривизны и пренебрежимо малой толщиной слоев рассмотрена в работах Л.А.Ильина и Н.АЛобковой [45] Д^бЗ , А.Р.0динца[881 в предположении упругой работы материала витков. Контактная податливость микро-, макрошероховатостей поверхностей контакта в работах не учитывается. Применены уравнения безмоментной теории,на основе которых задача сведена к системе обыкновенных дифференциальных уравнений. Л.А.Ильин и Н.А.Лобкова рассмотрели вопрос о влиянии величины коэффициента трения и угла перехлеста краев навивки на напряженное состояние пятислойной оболочки.
А.Р.Одинцом [88] в разрешающие уравнения введена кусочно-постоянная функция приведения,позволяющая ограничить величину касательных сил взаимодействия и рассмотреть влияние частичного проскальзывания в зоне контакта. Значение этой функции определяется из условий минимума потенциальной энергии деформации.
Рулонированная оболочка,наружные слои которой представлены тонкими трубами,а промежуточные - косослоем, рассмотрена В.Н.Гор-деевым, М.А.Микитаренком и В.В.Перельмутером в работе £321 . Напряженно-деформированное состояние оболочки определялось методом конечных элементов (ШСЭ) в упругой и упруго-пластической стадиях работы.с учетом одностороннего характера связей между витками и кулоновских сил трения. Как и в [881, введено предположение, что проскальзывание имеет место на тех участках,где касательные силы взаимодействия выше кулоновских сил трения.
Н.П.Мельниковым и Н.М.Потаповым в работе £771 пакет слоев представлен трансверсально изотропной средой с приведенными характеристиками. Решение построено по МКЭ.
Анализ литературных источников показал, что теория расчета многослойных рулонированных оболочек в настоящее время находится на этапе разработки математической модели, дающей возможность выявить наиболее важные особенности процесса деформирования этих существенно нелинейных объектов.
Из сказанного следует,что разработка методики расчета многослойных рулонированных цилиндрических оболочек,учитывающей особенности взаимодействия слоев, является актуальной с точки зрения теории и практики проблемой.
В настоящей диссертационной работе предложена,развита и использована для расчетов реальных объектов численно-аналитическая методика,основанная на исследованиях Г.Б.Ковнеристова [58] в области механики деформирования многослойных систем. При этом подходе деформации наружной поверхности слоя оболочки представлены деформацией упругого слоя толщиной,равной половине толщины витка оболочки, сложенной с деформацией срединной поверхности витка, в результате чего контактная задача оказывается сведенной к интегральному уравнению относительно неизвестных сил контактного взаимодействия между слоями. На основе такого подхода открывается возможность учесть контактную податливость макрошероховатостей поверхностей контакта,учесть проскальзывание и другие факторы, влияющие на напряженно-деформированное состояние рулокированной оболочки.
В настоящей работе принята модель многослойной рулонирован-ной бесконечно-длинной оболочки в виде гибкой растяжимой нити у упругой вдоль длины и свернутой в спираль, концы которой надлежащим образом закреплены,образовав замкнутые крайние витки.На вогнутой и выпуклой поверхностях витков спирали имеют место ми-кро-,макронеровности. Оценка общего напряженного состояния на основе введенных технических гипотез рассматривается в качестве первого этапа расчета оболочки. На втором этапе полученные результаты корректируются путем наложения данных расчета напряженного состояния в областях,содержащих концентраторы и щелевые дефекты. С этой целью развита,предложенная Г.Б.Ковнеристовым,методика построения решений для составных кусочно-однородных многосвязных областей с позиции смешанной плоской задачи теории упругости путем сведения задачи к интегральным уравнениям Г561.
Замена краевой задачи для системы уравнений в частных производных интегральным или интегро-дифференциальным уравнением является одним из наиболее эффективных приемов в математической физике. В рамках круга проблем теории упругости этот способ понижения размерности функциональных уравнений успешно применяется, начиная с работ Буссинеска и Герца.
Применительно к двумерным задачам теории упругости значительные успехи этого направления связаны с именами М.А.Садовского [85] , Г.В.Колосова, Н.И.Мусхелишвили [80] ,[&ll и большого числа их последователей,с высокой эффективностью использовавших аппарат теории функций комплексного переменного.
В пространственных задачах теории упругости и термоупругости вопросы получения и исследования интегральных уравнений получили освещение в трудах В.Д.Купрадзе [?0] , посвященных распространению на эту область математической физики идей теории потенциала.
Другой путь к получению интегральных уравнений теории упругости указали работы С.М.Белоносова [i4] , И.И.Воровича [£91 , Я.С.Уфлянда [122] и др.ученых,предложивших использовать различного рода интегральные преобразования.
В основе работ А.Я.Александрова [51 , Ю.В.Верюжского [25], Н.А.Кильчевского [54] , Г.Б.Ковнеристова [56] , Т.Круза Ш1 и других ученых, посвященных получению разрешающих интегральных уравнений,лежит идея использования для этой цели теоремы Бетти о взаимности работ. Этот подход применяется и в настоящей работе. К достоинствам этого подхода следует отнести ярко выраженный физический смысл уравнения и наглядность его получения.
Для решения интегральных уравнений задач теории упругости разработано большое количество методов. Каждый из них учитывает специфику рассматриваемой области, необходимую точность результата и,иногда, доступность применяемой методики.Сначала,в работах Герца , Н.М.Беляева [153 , И.Я.Штаермана [130] применялись и развивались аналитические идеи теории потенциала. Усложнение задач связано с применением аппарата теории аналитических функций, где основное интегральное уравнение для области,ограниченной простой замкнутой линией,ищется в виде комплексных рядов с привлечением свойств аналитических функций,как описано в монографиях Н.И.Мусхелишвили [81] , Г.И.Савина [111], А.С.Космода-мианского [66] , С.Г.Лехницкого [?з] , Э.И.Григолюка и Л.А.Фильш-тинского№] , с привлечением методов неопределенных коэффициентов, малого параметра, сведения к задаче линейного сопряжения с помощью метода Карлемана и других методов. С другой стороны,применялись численные методы решения смешанных задач,где неизвестные плотности представлялись кусочно-постоянными функциями [5] , [2.5] , [56] . Подход Б.Н.Жемочкина и А.П.Синицына[4£] , давал, в свое время, наглядное представление решения интегрального уравнения смешанной задачи в виде кусочно-постоянной функции. Кусочно-линейные представления плотностей имеются в работе [5"б] , естественное развитие численных методов привело к сплайновой аппроксимации плотностей
Представление плотностей в виде полиномов Чебышева характерно для монографий В.В.Панасюка, А.П.Дацышина и М.П.Саврука
С91] ДШ1, Н.И.Васильева, Ю.А.Клокова и А.Я.Штерстена , и др. Классический подход Фредгольма положен в основу численного решения задач в исследованиях В.З.Партона и П.И.Перлина .
Следует также упомянуть обзор по контактным задачам, решаемых методом теории функций комплексной переменной и теории потенциала, под общей редакцией Л.А.Галина .
В разработку аппарата решения краевых задач на основе численной реализации интегральных уравнений теории упругости значительный вклад внесла большая группа зарубежных ученых: П.Бенерджи, Р.Баттерфилд, С.А.Бреббия, Т.Крузе, Э.О.Оливейра, Ф.Риццо и др.
Из всех существующих численных методов решения интегральных уравнений в настоящей работе принят метод Мультона-Каландия [¿<9], где неизвестные и известные функции плотностей в интегральном уравнении, составленном для одного прямолинейного или кругового участка интегрирования, представляются в кусочно-полиномиальном виде при помощи тригонометрических интерполяционных полиномов Лагранжа и полиномов специального вида.по Чебышевским узлам интерполяции с возможными особенностями решений в узловых и угловых точках контура. Забегая вперед, следует отметить,что такой вычислительный алгоритм приводит к исключительно точным результатам.
Целью настоящей работы является разработка практической численно-аналитической методики расчета многослойных рулонированных цилиндрических оболочек на внутреннее равномерное давление, анализ на ее основе напряженно-деформированного состояния оболочки,а также областей содержащих концентраторы. На первом этапе расчета определяется напряженное состояние оболочки в целому использованием метода сведения к интегральным уравнениям контактной задачи теории оболочек; на втором - на основе интегральных соотношений для напряжений и перемещений на контуре исследуемой области уточняется напряженно-деформированное состояние вблизи концентраторов и щелевых дефектов оболочки в постановке смешанной плоской задачи теории упругости.
Разработка методики расчета предусматривает решение следующих вопросов:
- построение разрешающих уравнений задач статики многослойных рулонированных оболочек;
- разработку численных алгоритмов прочностных задач рулонированных оболочек;
- развитие метода граничных интегральных уравнений плоской теории упругости на основе кусочно-полиномиального представления функций плотностей с учетом возможных особенностей их поведения в узловых и угловых точках контура исследуемой области;
- разработку численно-аналитических алгоритмов расчета составных кусочно-однородных многосвязных областей в постановке смешанной плоской задачи теории упругости на основе интегральных уравнений;
- исследование достоверности численно-аналитической методики на решении ряда контрольно-тестовых задач;
- применение разработанной методики к решению задач,представляющих практический интерес.
Диссертация выполнена в соответствии с общим планом научных исследований,выполняемых на кафедре строительной механики и Проблемной научно-исследовательской лаборатории тонкостенных пространственных конструкций Киевского ордена Трудового Красного Знамени инженерно-строительного института под руководством профессора, доктора технических наук Г.В.Йсаханова. Значительная помощь в разработке теоретических вопросов была оказана доцентом,кандидатом технических наук Ковнеристовым Г.Б.
Диссертация состоит из введения, четырех глав,заключения, списка использованной литературы и приложения.
ЗАКЛЮЧЕНИЕ
В диссертационной работе разработана двухстадийная численно-аналитическая методика определения напряженно-деформированного состояния многослойных рулонированных цилиндрических оболочек . На первом этапе расчета на основе физически обоснованных упрощающих статических и кинематических гипотез в постановке плоской контактной задачи определяется двухосное напряженное состояние оболочки в зонах, удаленных от торцевых сечений ; на втором - в постановке плоской смешанной задачи теории упругости уточняются напряженно-деформированные состояния вблизи концентраторов и щелевых дефектов .
Получены разрешающие уравнения статики многослойных рулониро-ванных цилиндрических оболочек . Учтены податливость макрошероховатостей поверхностей слоев и проскальзывание в контакте . Учтена начальная погибь, предварительный натяг и опрессовка рулонированной стенки оболочки . Предусмотрена возможность расчета оболочки с частично разрушенной стенкой .
Разработаны численные алгоритмы расчета рулонированных оболочек, основанные на использовании метода коллокации. Для описания процесса деформирования слоев оболочки из физически нелинейного материала . применялся метод И.А.Биргера .
Развиты вычислительные средства метода граничных интегральных уравнений для решения плоской задачи теории упругости на основе кусочно-полиномиальной аппроксимации плотностей , ориентированные на исследование составных многосвязных областей с учетом заданной асимптотики решений в узловых и угловых точках контура .
Решением модельных задач подтверждена достоверность результатов, получаемых на основе разработанной численно-аналитической методики. Сопоставлениями установлено , что результаты расчетов как на этапе определения общего напряженно-деформированного состояния оболочки , так и на этапе исследования концентрации напряжений правильно отражают действительную работу рулонированных конструкций .
Решены новые задачи о напряженно-деформированном состоянии многослойных рулонированных оболочек ; выявлены особенности поля напряжений , вносимые предварительной опрессовкой рулонированной стенки , предварительным натягом витков спирали , ослаблением продольным щелевым разрезом или ловушкой , влиянием погиби ; определены параметры концентрации напряжений вблизи концентраторов, в частности , в зоне рассверленной вершины ловушки .
Разработанная методика и полученные на ее основе результаты могут быть использованы при проверочных расчетах на прочность рулонированных трубопроводов и сосудов высокого давления , а также для практических расчетов широкого класса объектов и элементов конструк- . ций .
1. A.c. 232924 (СССР). Способ изготовления корпусов сосудов и труб / Б.Е.Патон и др.-Опубл. в Б.И., 1969, № 2.
2. Айнбиндер А.Б., Камерштейн А.Г. Расчет магистральных трубопроводов на прочность и устойчивость. М.: Недра,1982.-341с.
3. Александров А.Я., Куршин JI.M. Многослойные пластинки и оболочки. В кн.: Труды УП Всесоюзн.конф. по теории оболочек и пластинок.Днепропетровск, 1969.-М.: Наука,1970,с.714-721.
4. Александров А.Я. Решение основных задач теории упругости путем численной реализации метода интегральных уравнений. -В кн.: Успехи механики деформируемых сред.- М.,1975', с.3-24.
5. Александров А.Я., Соколов Ю.М. Пространственные задачи теории упругости. М.: Наука, 1978. - 462 с.
6. Александров В.И. Осесимметричная контактная задача для упругого бесконечного цилиндра. Изв. АН СССР. ОТН механ. и машиностроение, №5, 1962.
7. Амбарцумян С.А. Некоторые вопросы развития теории анизотропных слоистых оболочек. Изв. АН Арм.ССР,серия физ.-мат.наук, 1964, т.17, №3, с.15-27.
8. Амбарцумян O.A. Общая теория анизотропных оболочек. М.: Наука, 1974. - 448 с.
9. Андреев А.Н., Нешровский Ю.В. К теории упругих многослойных анизотропных оболочек. Изв. АН СССР. МТТ,1977,№5, с.87-96.
10. Аникин Е.А. Технология укладки трубопроводов 0 1420 из многослойных труб. Строительство трубопроводов, 1984, № 6, c.I9r2I.
11. Бахвалов Н.С. Численные методы, т.1. М.: Наука,1975 - 631с.
12. Баш В.Я., Коваленко В.А., Пискун В.В. Исследование упруго-пластического деформирования многослойных рулонированных цилиндрических обечаек. Проблемы прочности,1983,№ 8, с.93-96.
13. Безвербный А.Ф. О равнопрочности витков многослойного руло-нированного цилиндра. Проблемы прочности, 1982,1 6,с.64-67.
14. Белоносов С.М. Основные плоские статические задачи для од-носвязных и двусвязных областей.- Новосибирск: Сиб.отд. АН СССР, 1962. -231 с.
15. Беляев Н.М. Местные напряжения при сжатии упругих тел.-Ленинград: Сб.инженерные сооружения и строительная механика, 1924.
16. Болотин В.В., Новичков Ю.Н. Механика многослойных конструкций. М.: Машиностроение,1980, -375 с.
17. Бородавкин П.П., Синюков А.М. Прочность магистральных трубопроводов. М.: Недра, 1984, -245 с.
18. Брич З.С., Капилевич Д.В., Котик С.Ю., Цагельский В.И. ФОРТРАН ЕС ЭВМ. -М.: Статистика,1978. -264 с.
19. Броек Д. Основы механики разрушения. Пер. с англ. В.И.Дорофеева. -М.: В.ш.,1980. -388 с.
20. Бурак М.й. Экспериментальное определение напряжений в слоях многослойной рулонированной оболочки.- Проблемы прочности, 1982, № 6, с.67-71.
21. Буяков И.А. Нелинейные уравнения типа Тимошенко многослойных анизотропных оболочек.- Механика композиционных материалов, 1979, №3, с.501-507.
22. Василенко А.Т. К расчету по уточненной модели ортотропныхи слоистых оболочек переменной толщины.- Прикладная механика, 1977, т.13, К*?, с. 28-36.
23. Васильев В.В. Теория слоистых оболочек. В кн.: Теория пластин и оболочек: Докл.ХШ Всесоюзной конф. по теории пластин и оболочек - Таллин, 1983,ч.1, с.163-168.
24. Васильев Н.Й., Клоков Ю.А., Шкерстена А.Я. Применение полиномов Чебышева в численном анализе.- Рига: Зинатне,1984.- 240 с.
25. Верюжский Ю.В. Численные методы потенциала в некоторых задачах прикладной механики. Киев: Вища школа,1978. - 184с.
26. Вислобицкий П.А. Распределение усилий в трехслойной рулони-рованной оболочке. Проблемы прочности,1981, №2, с.73-76.
27. Вольмир A.C. Гибкие пластинки и оболочки. -М.: Гостехиздат, 1956, 419 с.
28. Ворович И.И. Некоторые математические вопросы теории пластин и оболочек. В кн.: Механика твердого тела: Труды П Всесоюзн. съезда по теор. и прикл.механике,вып.3. -М.: Наука,1966,с'. II6-I36.
29. Ворович И.И., Александров В.М., Бабенко В.А. Неклассические смешанные задачи теории упругости. М.: Наука,1974. - 455с.
30. Гольденвейзер А.Л. Асимптотический метод построения теории оболочек. В кн.: Материалы I Всесоюзной школы по теории и численным методам расчета оболочек и пластин. -Тбилиси,1975.
31. Гольденвейзер А.Л. Теория упругих тонких оболочек.- М.:Наука, 1976. -512 с.
32. Гордеев В.Н., Микитаренко М.А., Перельмутер В.В. О расчете спирально-многослойного сосуда высокого давления в упруго-пластической стадии. Проблемы прочности,1979,№ 7,с.99-104.
33. Градштейн И.С., Рыжик И.M. Таблицы интегралов,сумм,рядов и произведений. М.: Наука, 1971. - 1100 с.
34. Григолюк Э.И., Грингауз М.Г., Фильштинский Л.А. К решению двумерных задач теории упругости для областей с кусочно-гладкой границей. Докл. АН СССР,1984, т.275, М, - с.831-834.
35. Григолюк Э.Й., Фильштинский Л.А. Перфорированные пластины и оболочки. М.: Наука, 1970. -556 с.
36. Григолюк Э,И., Чулков П.П. Нелинейные уравнения пологих многослойных оболочек регулярного строения. Инж.журн.МТТ,1967, № I, с.163-169.
37. Григоренко Я.М., Василенко А.Т. Об учете неоднородности и деформации поперечного сдвига по толщине в слоистых оболочках. Прикладная механика,1977,т.13,№ 10, с.36-42.
38. Гузь А.Н. и др. Цилиндрические оболочки ослабленные отверстиями." Киев: Наук.думка,1974. -271 с.
39. Гурьева Л.И. Напряженное состояние цилиндрической оболочки, скрепленной обмоткой. -В кн.: Эксперимент и техника высоких газовых и твердофазных давлений, М.,1978,с.38-48.
40. Демкин Н.Б. Фактическая площадь касания твердых поверхностей. М.: Из-во АН СССР,1962.
41. Жемочкин Б.Н., Синицын А.П. Практические методы расчета фундаментных балок и плит на упругом основании без гипотезы Винклера. М.: Гос.изд.литер, по строит.,архитектуре и строит, материалов,1962. -238 с.
42. Ильин Л.А., Лобкова H.A., Нехотящий В.А. Исследование напряженно-деформированного состояния многослойной рулонированной стенки. Хим. и нефт.машиностроение,1979,№9, с.10-12.
43. Ильин Л.А., Лобкова H.A., Нехотящий В.А., Стариков Н.П. Схематизация многослойной рулонированной стенки сосуда анизотропным цилиндром,с усредненными свойствами навивки.— Прикладная механика,1979,т.15,МО, с.58-63.
44. Ильин JI.А., Лобкова H.A. Упругая деформация и проскальзывание слоев рулонированной цилиндрической оболочки при нагру-жении внутренним давлением. Прикладная механика,1981,т.17, № 6, с.59-63.
45. Ильин Л.А., Лобкова H.A. Исследование напряженного состояния рулонированной цилиндрической оболочки с учетом трения. При> кладная механика,1982,т.18, № II, с.45-56.
46. Исаханов Г.В., Коснеристов Г.Б., Корх А.И. Решение задачи прочности многослойной рулонированной цилиндрической оболочки. Проблемы прочности,1984, М, с.49-52.
47. Каландия А.И. Математические методы двумерной упругости. -М.: Наука, 1973. 304с.
48. Кепич Т.Ю., Нуркиянов A.M. Исследование многослойных конструкций с концентраторами поляризационно-оптическим методом.— В кн.: Экспериментальные методы исследования деформаций и напряжений. Киев: ИЭС им.Е.О.Патона, 1983, с.107-117.
49. Кепич Т.Ю., Мазур К.И., Нуркиянов A.M., Савченко В.И. Исследование спирально-многослойных оболочек поляризационно-опти-ческим методом. В кн.: Тез.докл. Xlii Всесоюзн.конф.по теории пластин и оболочек. - Таллин,1983, ч.Ш, с.43-48.
50. Кильчевский H.A. Основы аналитической механики оболочек. -Киев: Из-во АН УССР,1963. 354 с.
51. Кильчевский H.A., Издебская Г.А., Киселевская Л.М. Лекции по аналитической механике оболочек. Киев: Вища школа,1974.
52. Коваленко В.А., Пискун В.В., Шевченко Ю.Н. Напряженное состояние многослойного с зазорами между слоями цилиндра при опрессовке. Тепловые напряжения в элементах конструкций, 1980, вып.20, с.15-19.
53. Ковнеристов Г'.Б. Интегральные уравнения контактной задачи теории упругости для заглубленных штампов. Киев: Сб.научных трудов Киев.инж.строит.ин-та,1962, с.20.
54. Ковнеристов Г.Б. Оеесимметричная контактная задача для цилиндрической оболочки, В кн.: Сопротивление материалов и теория сооружений.- Киев: Буд1вельник,1975,вып.27, с.119-124.
55. Ковнеристов Г.Б. Модель Герца-Тимошенко в контактных задачах стержней,плит и оболочек.- В кн.: Статика сооружений.-Киев, КИСИ: Будiвельник,1978, с.116-121.
56. Ковнеристов Г.Б., Корх А.И. Интерполяционные полиномы в смешанных задачах теории упругости.Задачи для областей с кусочно-круговой границей.- Киев,1984. 66 с. - Рукопись предст. Киев.инж.-строит.ин. Деп. в УкрНИИНТИ 29 мая 1984г.,№1112 Ук-84 Деп.
57. Ковнеристов Г.Б., Корх А.И. Напряженное состояние плоских анизотропных областей с полостями,усиленными накладками. -В кн.: Ш Научно-технич.конф. "Совершенствование эксплуатации и ремонта корпусов судов": Тез.докл. Калининград,1984, с.157-159.
58. Ковнеристов Г.Б., Корх А.И., Мельниченко Г.И., Гетун Г.В. Интерполяционные полиномы в смешанных задачах теории упругости. Задачи для областей с кусочно-прямолинейной границей.
59. Киев,1984. 74 с. - Рукопись предст. Киев.инж.-строит.ин. Деп. в УкрНИИНТЙ 29 мая 1984г.1166 Ук-Деп.
60. Ковнеристов Г.Б., Мельниченко Г.И., Гетун Г.В. Интерполяционные полиномы в смешанных задачах теории упругости. -Смешанные задачи механ.деф.тела. В кн.: П Всесоюзная научная конференция: Тез.докл. Днепропетровск,1981, с.86-87.
61. Ковнеристов Г.Б., Мельниченко Г.И., Гезун Г.В. Решение задач прочности для тонколистовых конструкций,усиленных накладками в зонах локальных ослаблений. В кн.: Тез.докл.
62. П Научно-технической конференции "Совершенствование эксплуатации и ремонта корпусов судов". Калининград, 1981,с.95.
63. Королев Е.М., Борсук Е.Г., Раевский Г.В., Ромашевский В.Б.
64. Экспериментальные исследования напряженного состояния рулонированных сосудов высокого давления. В кн.: Химическое машиностроение.Рулонные сосуды высокого давления. - М.: НИИ химмаш, 1973,вып.63, с.3-23.
65. Корх А.И. Учет предварительного натяга и начальной погиби в задаче прочности рулонированной оболочки. Киев,1984.15 с. Рукопись предст.Киев.инж.-строит.ин. Деп. в УкрНИИНТИ 12 января 1984г., № 532 Ук-Д84.
66. Космодамианский A.C. Плоская задача теории упругости для пластин с отверстиями,вырезами и выступами. Киев: Вища школа, 1975. - 227с.
67. Крагельский И.В., Виноградова И.Э. Коэффициенты трения.- М.: Машгиз,1962. 220 с.
68. Крагельский И.В., Добычин М.Н., Камбалов B.C. Основы расчета на трение и износ. М.: Машиностроение,1977. - 525 с.
69. Куликов Г.М. К теории многослойных пологих оболочек конечного прогиба. Изв. АН СССР.МТТ,1979, №3, с.188-191.
70. Купрадзе В.Д. Методы потенциала в теории упругости.- М.: Физматгиз, 1963. 472 с.
71. Левина З.М., Решетов Д.М. Контактная жесткость деталей машин. М.: Машиностроение,1971. - 264 с.
72. Лейбензон Л.С. Курс теории упругости.- М., Л.: Гостехиздат, 1947. 464 с.
73. Лехницкий С.Г. Теория упругости анизотропного тела. 2-е изд., перераб. и доп. - М.: Наука, 1977. - 416 с.
74. Либреску Л. К уточненной линейной теории упругих анизотропных многослойных оболочек. ч.1. Механика полимеров,1975, №6, с.1038-1050; ч.П. - Механика полимеров, 1976, №1, с.100--109.
75. Лобкова H.A. Исследование зависимости напряженно-деформированного состояния рулонированной цилиндрической оболочки от перехлеста краев навивки. Прикладная механика,1984, т.20, №6, с.59-63.
76. Максак В.Н. Предварительное смещение и жесткость механического контакта.- М.: Наука,1975, 60 с.
77. Мельников Н.П., Потапов Н.М. Об одной физической модели, учитывающей действительную работу многослойных сосудов высокого давления. Строит.механика и расчет сооружений,1980, Ш, с.10-15.
78. Метод граничных интегральных уравнений. Вычислительные аспекты и приложения в механике / Круз Т., Риццо Ф. -М.: Мир, 1978. 215 с.
79. Москаленко В.Н., Новичков Ю.Н. Изгиб толстых многослойных оболочек. Инж.журнал. МТТ,1968, Ш 3, с.149-153.
80. Мусхелишвили Н.И. Сингулярные интегральные уравнения. М.: Наука, 1962. - 511 с.
81. Мусхелишвили Н.й. Некоторые основные задачи математической теории упругости. М.: Наука, 1966. - 707 с.
82. Натансон И.П. Конструктивная теория функций. М.; Л.: Гос-техиздат, 1949. - 688 с.
83. Нетребский М.А., Сагидаев Ю.М., Раевский Г.В. Оптимальные предварительные напряжения в многослойных трубах. Прикладная механика, 1971, №9, с.116-120.
84. Нехотящий В.А., Нетребский М.А., Тертышник К.Г. Исследование рулонированных конструкций на многослойных цилиндрических оболочках. В кн.: Хим.машиностроение. Рулонированные сосуды высокого давления.- М.: НИИхиммаш,1973, вып.63, с.95-99.
85. Новацкий В. Теория упругости. М.: Мир, 1975. - 872 с.
86. Новичков Ю.Н. Осесимметричная деформация многослойной цилиндрической оболочки с учетом проскальзывания между слоями.-Сб.науч.трудов Моск.энер.ин. М.: МЭИ,1975,вып.227, с.109-118.
87. Нуркиянов А.М. Применение поляризационно-оптического метода к исследованию рулонированных конструкций с концентраторами и дефектами. Автореферат дис. канд.ф-мат.наук Киев, 1983. - 21 с.
88. Одинец А.Р. Напряженно-деформированное состояние спирально-многослойных цилиндрических оболочек при воздействии внутренним давлением. Прикладная механика,1982,т.18,№5, c.III-114.
89. ОСТ 26-01.221-70. Сосуды многослойные рулонированные стальные высокого давления.-М.:Минхиммаш,1970.
90. ОСТ 26-1046-74. Сосуды и аппараты высокого давления:Нормы и методы расчета на прочность.- М.: Минхиммаш, 1974.
91. Панасюк В.В., Саврук М.П., Дацышин А.П. Распределение напряжений около трещин в пластинках и оболочках.-Киев:Наукова думка, 1976. 443 с.
92. Партон В.З., Перлин П.И. Интегральные уравнения теории упругости. М.: Наука, 1977. -312 с.
93. Пашковский С. Вычислительные применения многочленов и рядов Чебышева. М.: Наука,1983. - 384с.
94. Пимштейн П.Г. К теоретическому расчету рулонированных сосудов высокого давления. В кн.: Химическое машиностроение. Рулонирован. сосуды высокого давления. - М.: НИИхиммаш,1973, вып.63, с.59-71.
95. Пимштейн П.Г., Жукова В.Н. Расчет напряжений в многослойном цилиндре с учетом особенностей контакта слоев. Проблемы прочности,1977, №5, с.71-77.
96. Пимштейн П.Г., Борсук Е.Г. Экспериментальные исследования напряженного состояния рулонированной оболочки под действием внутреннего давления. Проблемы прочности,1978,№9,с.71-77.
97. Пимштейн П.Г., Борсук Е.Г. Способ определения среднего натяга слоев при изготовлении многослойных сосудов. Хим. и нефт.машиностроение,1979, № 9, с.18-19.
98. Писаренко Г.С., Бабенко А.Е. Напряженно-деформированное состояние трехслойной рулонированной цилиндрической оболочки под внутренним давлением. Проблемы прочности, 1977,}f°3, с.3-5.
99. Писаренко Г.С., Бабенко А.Е. Напряженно-деформированное состояние многослойной рулонированной оболочки под действием внутреннего давления. Проблемы прочности,1977, №7,с.3-6.
100. Писаренко Г.С., Богинич O.E. Приближенный расчет напряжен-но-деформировэнного состояния рулонированной цилиндрической оболочки. Проблемы прочности,1976, №2, с.10-15.
101. Пискунов В.Г. Об одном варианте неклассической теории многослойных пологих оболочек и пластин. Прикладная механика,1979, т.15, № II, с.76-81.
102. Привалов И.И. Введение в теорию функций комплексного переменного. М.: Наука, 1977.
103. Прусаков А.П. Нелинейные уравнения изгиба пологих многослойных оболочек. Прикладная механика, 1971, т.7,№3.
104. Развитие теории контактных задач в СССР /Под общ.ред. Л.А.Галина и др. М.: Наука, 1976.
105. Рассказов А.О. К теории многослойных ортотропных пологих оболочек. Прикладная механика,1976, т.12, № II, с.50-56.
106. Юб.Рвачев В.Л., Проценко B.C. Контактные задачи теории упругости для неклассических областей. Киев: Наукова думка,1977. - 235 с.
107. Рыжов Э.В. Контактная жесткость деталей машин. М.: Машиностроение, 1966.- 195с.
108. Рыжов Э.В., Колесников Ю.В., Суслов А.Г. Контактирование твердых тел при статических и динамических нагрузках. -Киев: Наукова думка,1982. 172 с.
109. Рябов А.Ф. Основные уравнения теории многослойных оболочек с учетом деформации поперечного сдвига. В кн.: Сопротивление материалов и теория сооружений. - Киев:Буд1-вельник,1965,вып.3, с.17-27.
110. ПО.Рябов А.Ф. Розрахунок багатошарових оболонок. Киев: Буд1вельник,1968. - 96 с.
111. Ш.Савин Г.Н. Распределение напряжений около отверстий. -Киев: Наукова думка, 1968. 887 с.
112. П2.Саврук М.П. Двумерные задачи упругости для тел с трещинами. Киев: Наукова думка,1981. - 324 с.
113. Савченко В.И., Кепич Т.Ю., Нуркиянов A.M. Применение по-ляризационно-оптического метода к исследованию рулонированных оболочек. Вiсник КГУ.Математика i мехашка,1983, вып.25, с.122-128.
114. Сагдаев Ю.М. Експериментальне досл!дження м1цност1 сшрально-багатошарово! оболонки. В кн.: Onip MaTepia-л1в i reopifl споруд. - Крив: Буд!вельник,1966,вып.4, с.107-112.
115. Сагидаев Ю.М. Расчет спирально-многослойной оболочки. -Прикладная механика,1967, т.З, №5, с,61-68.
116. Стечкин С.Б., Субботин Ю.Н. Сплайны в вычислительной математике. М.: Наука, 1976. - 248 с.
117. Терегулов А.Г. К теории многослойных анизотропных оболочек. Сб.: Исследования по теории пластин и оболочек.вып.6,7. Казань: йз-во Казан.ун-та, 1970, с.762-767.
118. Тимошенко С.П., Гудьер Дж. Теория упругости. М.: Наука, 1975. - 575 с.
119. Труфяков В.И. Некоторые вопросы повышения несущей способности и долговечности сварных конструкций.- В кн.: Надежность и долговечность машин и сооружений. Киев,Наукова думка, 1983, вып.З, с.3-12.
120. Усенко В.Г. Состояние и перспективы совершенствования аппаратов и трубопроводов высокого давления. -. Хим. и нефт. машиностроение,1979, № 9, с.3-4.
121. Уфлянд Я.С. Интегральные преобразования в задачах теории упругости. М.; Л.: Из-во АН СССР,1963. - 367 с.
122. Фадеев Д.К., Фадеева В.Н. Вычислительные методы линейной алгебры. М.: Наука,1963. - 734 с.
123. Феодосьев В.И. Геометрически нелинейные задачи теории пластин и оболочек. В кн.: Труды У1 Всесоюзн.конф. по теории оболочек и пластинок. Баку,1966. - М. : Наука,1966, с.971-976.
124. Химическое машиностроение. Рулонированные сосуды высокого давления. Сб.трудов Моск.' НИИхиммаша, - M. : 1973,вып.63.
125. Хорошун Л.П. О построении уравнений слоистых пластин и оболочек. Прикладная механика, 1978, т.14,№ 10, с.3-21.
126. Чепига В.Е. К уточненной теории слоистых оболочек. -Прикладная механика, 1976, т.12, № II, с.45-49.
127. Шевченко Ю.Н., Бабешко М.Е., Пискун В.В. и др. Решение осе-симметричной задачи термопластичности для тонкостенных и толстостенных тел вращения на ЕС ЭВМ. Киев: Наукова думка, 1980. - 196с.
128. Шопа В.М., Полевой Б.Н., Величкович C.B. Расчет пакета цилиндрических оболочек с учетом трения. Докл. АН УССР, 1982, А, № 5, с.50-55.
129. Штаерман И.Я. Контактные задачи теории упругости. М.; Л.: Гостехиздат, 1949, - 270 с.
130. Яблонский Б.В. Напряженное состояние многослойной конструкции при навивке ленты на цилиндр. Прикладная механика, 1971, т.7, № 2, с.130-133.
131. IDR-AGOß 5). CioCLOV. íleciftejufce suß pïesiune:
132. A no-tiza. sic&£r de "teusiune defownaiie.- Bumfc&ci* ,Ь138335Ц p.
133. DfoiTetm P.R.A. Plane zl^ss analysis a деиелизиб inieß^cut method . X ASO E , Eng. Mec-Pi. $>iir., >f96S, Fe£>*., p. ?9-85.
134. Potó 1н-к)К , Atuje\ric Ánd'zo. Bounola/2¿) efemerub ïneikod in plane elasíicbtj. g-fc^a.-ires-fcu.,№1,0*8, M§ 5-6 (i-ъ).
135. Rt^BO M. An. búñ^jxi execution opp^otJi to boundary ш¿ue pto^tems of ctocsstcat eàxsios-èaMcs Q.ua/fc¿. Appt. Matt., 496?, 25,4.
136. Manááctircclüe G.F. Methods of -ike ih-ecmj о} anoc&jiic -functions tu -the Игго'гу of &tcLsticxt^ . Compte oc. A not : Metk., T^nds, and. Appt. Be'üiin., 4983 , 2.8Q-2L95\
137. Ko-Beujaskt &h.oichi, Appticaiions oj- е^дсиЬток method (-boLLndcLruj element method ).— "^bgailpe , 0. ¿oc . Maie* .got., 9ap.", 1983,3a, MS 363 , 4293-1303 .