Численно-аналитические методы исследования концентрации напряжений в элементах конструкций при пространственном напряженном состоянии тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

ВВЕДЕНИЕ. Основное содержание работы и обзор литературы.

1. ИССЛЕДОВАНИЕ КОНЦЕНТРАЦИИ НАПРЯЖЕНИЙ В ИЗОТРОПНЫХ

ТЕЛАХ, СОДЕРЖАЩИХ СИСТЕМУ НЕПЕРЕСЕКАЮЩИХСЯ СФЕРИЧЕСКИХ КОНЦЕНТРАТОРОВ

1.1. Введение.

1.2. Представление напряжений и перемещений пространственного напряженного состояния через аналитические функции.

1.3. Периодическая задача для упругой среды с бесконечной системой сферических полостей.

1.4. Решение основных задач теории упругости для пространства с двумя полостями. /

1.5. Численная реализация задач.

1.5. Анализ напряженного состояния.

2. ИССЛЕДОВАНИЕ КОНЦЕНТРАЦИИ НАПРЯЖЕНИЙ В ТРАНС-ВЕРСАЛЬНО-ИЗОТРОПНЫХ ТЕЛАХ, СОДЕРЖАЩИХ СИСТЕМУ НЕПЕРЕСЕКАЮЩИХСЯ ЭЛЛИПСОИДАЛЬНЫХ ПОЛОСТЕЙ ИЛИ ВКЛЮЧЕНИЙ.

2.1.Представления перемещений и напряжений через аналитические функции.

2.2. Построение аналитических функций в случае двух концентраторов

2.3. Формула переноса.

2.4. Решение первой основной задачи теории упругости для пространства с двумя эллипсоидальными полостями.

2.5. Решение второй основной задачи для упругой среды с двумя эллипсоидальными полостями.

2.6. Определение напряжений.

2.7. Периодическая система полостей.

2:8. Исследование напряженного состояния.

2.9. Некоторые особенности численной реализации исследуемого класса задач.

3. РЕШЕНИЕ ОСЕСИММЕТРИЧНЫХ ЗАДАЧ ТЕОРИИ УПРУГОСТИ ДЛЯ ИЗОТРОПНЫХ ЦИЛИНДРИЧЕСКИХ ТЕЛ.

3.1. Осесимметричная задача теории упругости для бесконечного кругового цилиндра.

3.2. Осесимметричная задача для упругой среды с бесконечной цилиндрической полостью.

3.3. Определение напряженного состояния в сплошном цилиндре и в бесконечной среде с цилиндрической полостью.;.

3.3.1. Напряженное состояние в сплошном цидиндре.

3.3.2. Напряженное состояние в бесконечной упругой среде с цилиндрической полостью.

3.4. Осесимметричная задача для бесконечного полого цилиндра

3.5. Смешанные задачи для бесконечного цилиндра и упругой среды с цилиндрической полостью.

3.6. Смешанная задача для бесконечного полого цилиндра.

3:7. Особенности численной реализации задачи.

4. СТАТИЧЕСКИЕ И СТАЦИОНАРНЫЕ ДИНАМИЧЕСКИЕ ЗА

ДАЧИ ОПРЕДЕЛЕНИЯ НАПРЯЖЕННОГО СОСТОЯНИЯ В ИЗОТРОПНЫХ И ТРАНСВЕРСАЛЬНО-ИЗОТРОПНЫХ ЦИЛИНД РИЧЕСКИХ ТЕЛАХ

4.1. Осесимметричная задача для кругового цилиндра и бесконечной упругой среды с цилиндрической полостью.

4.2. Бесконечный сплошной цилиндр, находящийся под действием внешней нагрузки.

4.3. Осесимметричная задача для бесконечного полого цилиндра.

4.4. Напряженное состояние в полом двухслойном цилиндре.

4.5. Смешанная задача для сплошного цилиндра.

4. 6. Смешанная задача для двухслойного трубопровода.

5. ИССЛЕДОВАНИЕ ЛОКАЛЬНОГО ТРЕХМЕРНОГО НАПРЯЖЕН НОГО СОСТОЯНИЯ В МНОГОСЛОЙНЫХ ПЛАСТИНАХ С ОТВЕРСТИЯМИ

5.1. Постановка задачи.

5.2. Определение эффективных характеристик многослойной пластины.

5.3. Решение плоской задачи теории упругости.

5.4. Определение межслойных напряжений методом конечных элементов.

5.5. Анализ полученных результатов.

Основное содержание работы и обзор литературы

В современном машиностроении важнейшими критериями при проектировании и анализе работы конструкций в сложных эксплуатационных условиях являются прочностные, усталостные и жесткостные свойства отдельных ее элементов. Одним из основных факторов, влияющих на эти свойства, является концентрация напряжений, представляющая собой локальное увеличение напряжений, возникающих как в процессе эксплуатации конструкции, так и при ее изготовлении.

Основные виды источников концентрации напряжений (концентраторов) можно условно разделить на три группы: а) геометрические разрывы в материале, образованные полостями; б) неоднородность материала; в) разрывный характер нагрузки.

Каждый из этих видов связан с определенным классом задач учитывающих специфику моделируемого объекта и требующих индивидуальных методов их решения.

Наличие полостей может быть связано как с присутствием технологических дефектов в элементе конструкции, так и с особенностями структуры материала.

Неоднородность материала имеет множественный характер. В диссертационной работе рассматриваются два типа неоднородностей. К первому относятся жесткие включения, которые, как и в случае полостей могут моделировать либо технологические дефекты, либо армирующие элементы композиционного материала. Второй тип неоднородности связан с многослойными композиционными материалами, образованными слоями с различными механическими характеристиками.

Разрывный характер нагрузок может возникать как при эксплуатации конструкции в результате контактного взаимодействия отдельных ее элементов на ограниченной области поверхности, так и в результате конструктивных особенностей этих элементов. Как отмечалось выше, правильный учет влияния перечисленных видов концентраторов на напряженное состояние в элементах конструкции играет важнейшую роль при оценке прочности и долговечности конструкции в целом.

В настоящей работе предлагаются численно-аналитические методы исследования концентрации напряжений, инициируемых перечисленными выше типами концентраторов в трехмерных телах с позиций линейной теории упругости. В основе предлагаемых подходов к исследованию индивидуального локального характера концентрации напряжений, вызванной данным типом концентраторов, лежит использование решений задач плоской и пространственной теории упругости для бесконечных тел, что существенно упрощает решение поставленной задачи. Следует отметить, что решение перечисленных типов задач связано со значительно большими сложностями, чем решение аналогичных плоских задач. Это связано как с отсутствием такого же широкого круга аналитических методов, которые разработаны для решения плоских задач, так и с существенно большим объемом вычислений при использовании приближенных численных методов.

В I, II главах источниками концентрации напряжений являются системы сферических, эллипсоидальных полостей или жестких включений в бесконечных изотропных и трансверсально-изотропных телах. В главах III и IV концентратором напряжений является разрывный характер внешней нагрузки, приложенной к цилиндрическим поверхностям. В главе V концентрация напряжений обусловлена разнородностью материалов, объединенных в единый массив.

В первой главе исследуется концентрация напряжений, вызванная системой непересекающихся и несоприкасающихся сферических концентраторов в виде полостей или жестких включений с центрами на одной прямой в бесконечном упругом теле. Основная задача заключается в построении аналитического решения и создании на базе этого решения устойчивого численного алгоритма, позволяющего исследовать взаимное влияние сферических концентраторов на величину концентрации напряжений. Получены решения первой и второй основных задач пространственной теории упругости для бесконечных тел, ограниченных конечным или бесконечным числом замкнутых непересекающихся поверхностей с центрами на одной прямой.

Общие вопросы связанные с решением трехмерных задач для многосвязных тел с достаточной полнотой были исследованы в работах Т.В. Бурчуладзе [31], Перлина П.И. [117,118], Шермана Д.И. [155]. В основе этих работ лежат результаты, полученные для аналогичных задач теории гармонических функций Н.М. Гюнтером [61], Н.И. Мусхелешвили [93], Д.И. Шерманом [156] и В.Д. Купрадзе [84].

Первая основная задача теории упругости для многосвязных тел решалась Д.И. Шерманом и Т.В. Бурчуладзе. Различными методами эти авторы свели задачу к решению сингулярных интегральных уравнений, показав их разрешимость. Бурчуладзе Т.В. получил аналогичный результат так же и для второй основной задачи. Перлин П.И. предложил метод решения задач теории упругости для случая, когда область ограничена двумя поверхностями. С помощью этого метода задача сводится к решению двух сингулярных уравнений. В работе [141], C.J1. Соболев предложил использовать альтернирующий метод Шварца при решении задач теории упругости для тел ограниченных системой замкнутых поверхностей. Описание алгоритма этого метода дано в монографии J1.B. Канторовича и В.И. Крылова [65]. По своей сути близким к этому методу является метод фиктивных областей, изложенный в работах J1.H. .Ясницкого [158,159], который так же может быть использован при решении подобного класса задач. Несмотря на теоретическую полноту исследования основных пространственных задач теории упругости для тел ограниченных несколькими поверхностями, фактическая численная реализация интегральных уравнений, полученных в работах перечисленных выше авторов, сопряжена с существенными трудностями вычислительного характера, которые значительно возрастают при сближении поверхностей. Этим можно объяснить отсутствие работ по исследованию концентрации напряжений, использующих методы, предлагаемые перечисленными выше авторами. Исключение составляют работы Перлина П.И. [119,120], где автор на основе предложенного им метода решает конкретные осесимметричные задачи для области, ограниченной эллипсоидом и сферой.

Различные формы решения задачи теории упругости для одного сферического концентратора рассматривались широким кругом авторов. Решение задачи в замкнутом виде для сферической полости при одноосном напряженном состояния на бесконечности принадлежит Саусвеллу Р.В. и приводится в монографии Ю.Н. Работнова [137]. В дальнейшем, общие решения задач для пространства, содержащего одну сферическую полость или одно инородное включение, были получены в работах А.Я. Александрова и Б.Л. Абрамяна[2], А.Я. Александрова и Ю.И. Соловьева [18], B.C. Вольперта [47], А.И. Лурье [90,91], Г.Я. Попова [130], Э. Штернберга [157]. В двух последних работах дан анализ и обширный обзор результатов, посвященных данной проблеме.

Первые работы по решению задач теории упругости для тел, ограниченных системой (конечной или бесконечной) сферических полостей с центрами на одной прямой в осесимметричной постановке принадлежат K.Mikata [182,183], Н. Miyamoto [179-181], Е. Tsuchida [201], М.А, Sadovsky Е. и Е.Sternberg [200]. Последние два автора получили решение задачи теории упругости для бесконечного тела, содержащего две сферические полости одинакового радиуса при осесимметричном однородном напряженном состоянии на бесконечности. Решение строилось в биполярных координатах и приводилось к бесконечной системе линейных алгебраических уравнений. Аналогичная задача была решена Н. Miyamoto [179,180], в предположении, что внешняя нагрузка возникает за счет центробежного поля инерционных сил. В дальнейшем этот же автор исследовал, приближенно, влияние более чем двух полостей на напряженное состояние в теле [181]. В работе [201] решена задача для бесконечного тела с несколькими полостями. Решение разыскивалось в форме Папковича и с помощью преобразования Лежандра сводилось к бесконечной системе линейных алгебраических уравнений. Для случая двух полостей проведен численный анализ. Задача для пространства с двумя полостями в осесимметричной постановке решалась так же в работе M.Hamada [170] с помощью итерационного метода, использующего решение задачи с одной полостью. В этой работе было проведено сравнение величины концентрации напряжений, полученной предлагаемым методом и методом конечных элементов. Эта же задача, с помощью метода конечных элементов рассматривалась в работе [134]. Осесимметричная задача для случая бесконечной периодической системы полостей рассматривалась в работах [171,181-183]. Следует отметить, что с формальной точки зрения решение задачи о пространстве с двумя сферическими концентраторами и решение задачи об упругой сфере с неконцентрической полостью или включением могут считаться эквивалентными.

Р.Н. Кауфман [79] рассматривала первую основную задачу теории упругости о всестороннем сжатии упругого шара с неконцентрической полостью. Решение уравнения равновесия в смещениях разыскивалось в форме рядов и с помощью формул переноса задача сводилось к бесконечной системе линейных алгебраических уравнений. Значительно позже аналогичная задачи рассматривалась в работах A.A. Капшивого [68,70] с помощью обобщенных аналитических функций Г.Н. Положего [124-126] и так же сводилась к решению беско-несчных систем линейных алгебраических уравнений. Этим же автором и тем же методом решалась и осесимметричная задача о пространстве с двумя полостями [69]. Как и предыдущие, эта задача была сведена к решению аналогичной бесконечной системы линейных алгебраических уравнений, что подтверждает формальную эквивалентность решения обоих типов задач. Приближенное решение этой же задачи на основе асимптотического метода возмущения формы границы получено в работе B.C. Колесова и др. [42]. Следует отметить работу J1.B. Баева и др. [28], в которой решение задачи для пространства, содержащего полые сферические включения с центрами на одной прямой, строится с помощью представлений Н. Miyamoto [179] в виде рядов. На базе полученного решения для случая взаимодействия полого включения и полости авторами проводится количественный и качественный анализ влияния межцен-тренного расстояния и отношения радиусов концентраторов на величину концентрации напряжений при осесимметричном нагружении. Численно-аналитическое решение для неоднороднй среды со сферической полостью или включением, основанное на методе возмущений, предлагется в работе A.B. Василенко и Н.Д. Панкратова [39].

Исследованию взаимного влияния системы сферических полостей на напряженное состояние в теле в неосесимметричной постановке посвящено значительно меньше работ. Общие вопросы решения задач статики упругого тела, ограниченного сферическими поверхностями, исследовал Головчан В.Т. [57,58]. В работе Е. Tsuchida [202] решается задача о пространстве, содержащем три непересекающиеся полости с центрами на одной прямой, подверженном равномерному растяжению на бесконечности в направлении перпендикулярном к оси симметрии. Решение сводится к двум бесконечным системам линейных алгебраических уравнений, первая из которых соответствует осесимметричной задаче, а вторая учитывает асимметрию напряженного состояния. Следует отметить работу G.Rodin [198], в которой на базе решения Эшелби для пространства с одним включением строится приближенное решение для пространства с N включениями. В качестве примера рассматривается задача одноосного растяжения пространства с двумя одинаковыми полостями. Проводится сравнение величины концентрации напряжений в полярных точках полостей, полученной предлагаемым методом и с помощью МКЭ. Удовлетворительное совпадение сравниваемых величин отмечается при отношении радиуса полости R к величине межцентренного расстояния /большем 0.4. Тем не менее, величины потенциальной энергии определенные с помощью этих двух методов отличаются не более чем на 5% при значениях R/f, близких к 0.5.

Большинство цитируемых выше работ можно разделить на две группы. В первой решение задач тем или иным способом сводится к бесконечным системам линейных алгебраических уравнений (СЛАУ) и проводится исследование этих систем, связанное с существованием и единственностью их решения. При этом, пути их численной реализации не указываются и вопросы, связанные с концентрацией напряжений не исследуются. Ко второй группе относятся работы, не затрагивающие вопросы существования и единственности полученных решений, но рассматривающие конкретные задачи по определению концентрации напряжений и использованию полученных решений для определения эффективных модулей композиционных материалов, моделируемых упругой средой со сферическими включениями. Однако, в большинстве работ исследования концентрации напряжений относятся к осесимметричным напряженным состояниям и ограничиваются диапазоном изменения параметра оо >Я//>0.4 при равных радиусах концентраторов). Исключение составляют работы [28,197], где этот диапазон шире.

Автором диссертации получены решения пространственных задач теории упругости для бесконечной среды содержащей: периодическую систему сферических полостей, две полости, два включения, полость и включение с центрами на одной прямой. На базе метода суперпозиций, использующего комплексные потенциалы в форме Колосова-Мусхелишвили, задача сводится к ряду бесконечных систем линейных алгебраических уравнений. Показано, что полученные системы вполне непрерывны, а число этих систем зависит от вида граничных условий. Предложен и применительно к конкретным задачам реализован эффективный алгоритм решения этих систем, позволяющий определять концентрацию напряжений при любой толщине перемычки между концентраторами. Основные результаты этой главы изложены в работах [13,14,15,50,104,106,109].

Во второй главе получены решения первой и второй основных задач теории упругости, позволяющие проводить исследование концентрации напряжений в трансверсально-изотропных телах, содержащих систему непересекающихся эллипсоидальных концентраторов в виде полостей или жестких включений. С помощью эллипсоидальных полостей или включений, в зависимости от расстояния между их центрами и геометрическими параметрами, можно моделировать широкий круг концентраторов, являющихся технологическими дефектами или армирующими элементами. Задачи подобного типа, как и задачи главы I, возникают при оценке прочностных, усталостных и жесткостных характеристик трехмерных элементов конструкций, выполненных из композиционного материала, моделируемого трансверсально-изотропной упругой средой.

Наряду с этим, построенные решения позволяют с высокой степенью точности исследовать концентрацию напряжений и в изотропном материале содержащем такие же концентраторы. Полученные решения могут быть использованы при оценке влияния разрыва волокон на прочность композиционного материала, моделируемого изотропной матрицей с абсолютно жесткими волокнами в виде сильно вытянутых эллипсоидальных включений.

Исследованию концентрации напряжений вблизи эллипсоидального концентратора в изотропных и трансверсально-изотропных телах посвящено значительное число работ, в которых получены как точные, так и приближенные решения. Среди них, прежде всего, следует отметить работы А.Я. Александрова [6,7,17,18], Ю.И. Соловьева [142,144], B.C. Вольперта [482,51], И.А. Кунина [81-83], Ю.Н. Подильчука [121-123], К.В. Соляник-Краса [146,147]. В этих работах различными методами получены общие решения первой и второй осно-вых задач теории упругости, с помощью которых можно исследовать концентрацию напряжений вблизи единичного эллипсоидального концентратора, представляющего собой полость или инородное включение. В работах [82,83] концентрация напряжений определяется для предельных случаев (игольчатая трещина или игольчатое включение) при общем виде анизотропии материала.

Вопросы, связанные с исследованием взаимного влияния на концентрацию напряжений нескольких эллипсоидальных концентраторов, особенно при их близком расположении, изучены значительно меньше. В работе В.И. Куща [85] получено решение, позволяющее определять эффективные упругие модули изотропной упругой среды, армированной периодически расположенными эллипсоидальными включениями. Влияние сфероидального включения или полости, заключенных в упругом изотропном шаре, на концентрацию напряжений рассматривалось в работах П.И. Перлина [119,120] и А.Г. Николаева [103]. Построение решения этих задач, с формальной точки зрения, аналогично построению решения задач о взаимодействии сферического и эллипсоидального кон-центркторов и может быть использовано для решения последних.

В работах N.Noda [187,188] задача о концентрации напряжений в окрестности двух вытянутых эллипсоидальных полостей или включений в изотропной среде с фокусами на одной оси, сведена к системе сингулярных интегральных уравнений. При этом, полученное решение иллюстрирует частный случай одноосного растяжения на бесконечности под углом к общей оси симметрии концентраторов.

Автором диссертации было получено общее решение первой и второй основной задачи теории упругости для трансверсально-изотропной среды, содержащей периодическую систему эллипсоидальных непересекающихся и несоприкасающихся полостей, а так же двух непересекающихся и несоприкасающихся эллипсоидальных полостей или жестких включений, при некоторых ограничениях на их взаимное расположение. Использование подхода, основанного на методе суперпозиций с использованием комплексных потенциалов С.Г. Лехницкого [88,89], позволило получить решение имеющее единую форму как для сферических, так и для эллипсоидальных концентраторов. Как и в главе I решение рассматриваемых задач сводится к бесконечным системам линейных алгебраических уравнений, численная реализация которых осуществлялась с помощью того же алгоритма, что и в первой главе. Основные результаты главы II изложены в работах [97,102,105,109,113,116,191 ].

В третьей и четвертой главах получены решения задач исследования осе-симметричного напряженного состояния в изотропных и трансверсально- изотропных бесконечных телах, ограниченных цилиндрическими полостями. Осе-симметричные задачи теории упругости для изотропных тел, ограниченных цилиндрическими поверхностями в различных постановках рассматривались в работах В.М. Александрова [19,20], A.A. Баблояна [26], A.B. Белоконя [27], Г.М. Валова [37], И.И. Воровича [44], Колтунова [74], А.И. Лурье [90,91] , Г.Н. Положего [126,127], Г.Я. Попова [130], H.A. Ростовцева [135] и ряда других авторов [34,45,62,71,73,92,147,148,152]. В перечисленных работах строятся общие аналитические решения основных задач теории упругости. Как правило, решения первой и второй основных задач получено с помощью функций напряжений Лява, а решение смешанной (контактной) задачи с помощью найденных решений первой и второй основных задач сводится к решению интегральных уравнений специального вида. М.Я. Беленький [34] получил формулы, связывающие компоненты осесимметричного напряженного состояния с двумя аналитическими в плоскости осевого сечения тела вращения функциями, которые можно использовать при решении краевых задач для сплошного цилиндра.

Задачи исследования напряженного состояния в цилиндрических телах от действия осесимметричной ступенчатой нагрузки решались М.А. Колтуновым [74] и S.Masao [150]. В этих работах определение напряжений и смещений с помощью функций напряжений Лява и преобразования Фурье [148] сводилось к вычислению несобственных интегралов. Полученные в этих работах результаты можно использовать для решения широкого класса осесимметричных задач, при любом виде нагрузки, аппроксимируемой кусочно непрерывными функциями. Наиболее сложным представляется решение смешанной задачи. Для сплошного цилиндра, заключенного в жесткую втулку конечной длины при тугой посадке, и пространства с цилиндрической полостью с жестким вкладышем большего диаметра, в работах [19,20,21,22,37,133,174] задача сводится к интегральным уравнениям с ядрами специального вида. В этих работах предлагаются методы приближенных решений полученных уравнений, однако, исследование концентрации напряжений не проводится. С практической точки зрения значительно больший интерес представляет задача о взаимодействии полого бесконечного цилиндра, находящегося под действием внутреннего давления, с жесткой втулкой. Задача в подобной постановке может быть рассмотрена как модельная для элемента трубопровода. Сведение этой задачи к интегральным уравнениям не представляет особенных затруднений. Но ядра полученных уравнений имеют достаточно сложный вид, что не позволяет получить решения интегральных уравнений теми же методами, что и в случае сплошного цилиндра или пространства с цилиндрической полостью. Исследование напряженного состояния, вызванного взаимодействием жесткой втулки и полого цилиндра, моделирующих реальные элементы трубопроводов, представляется целесообразным проводить с помощью метода граничных элементов и его модификаций [32,33,63,75,78,150].

Работ посвященных решению задач теории упругости для трансверсаль-но-изотропных тел значительно меньше. Из них отметим работы Ф.А. Вайн-штейна [35,36], Р.Я. Сунчелеева [143], А.Ф. Хрусталева [151], H.Kasano [175], N.Noda [186]. Три последние работы посвящены решению контактных задач для полых трансверсально-изотропных цилиндров. Практический интерес представляет работа [175], в которой задача сведена к сингулярному интегральному уравнению, решение которого ищется в виде рядов по полиномам Чебы-шева. Полученное решение позволило авторам провести анализ напряженного состояния для полого эпоксидного цилиндра, со стекло или графитовыми волокнами, при заданных постоянных перемещениях на кольцевом участке внешней поверхности. Этот же автор решал аналогичную задачу для полого изотропного цилиндра, которую привел к системе четырех парных интегральных уравнений [174]. Различные формы решения задач с однотипными граничными условиями свидетельствуют об отсутствии единого подхода к решению задач для изотропных и трансверсально-изотропных полых цилиндров.

В третьей и четвертой главах диссертации, на основе метода суперпозиций, разработан единый подход к решению задач теории упругости для изотропных и трансверсально-изотропных тел, ограниченных цилиндрическими поверхностями. С помощью полученных решений можно исследовать осесим-метричные напряженные состояния в сплошном бесконечном цилиндре, пространстве с цилиндрической полостью и в полом однослойном или многослойном цилиндрах как при заданных на поверхности усилиях или смещениях, так и при смешанных граничных условиях. Использование комплексных потенциалов при решении подобных задач методом суперпозиции сводится к отысканию аналитических функций в бесконечных полосах. Кроме того, этот же подход применяется в главе IV и к решению стационарно динамических задач термоупругости, когда внешние силовые возмущения и температурные поля распространяются вдоль оси симметрии цилиндрического тела с постоянной скоростью. Задачу о распространении возмущений с постоянной скоростью вдоль границы изотропной полуплоскости методами теории функций комплексного переменного рассматривали J1.A. Галин [53] и J. Radoc [197]. Задача о движущейся с постоянной скоростью нагузке вдоль полосы рассматривалась В.М.Александровым и Е.В.Коваленко [21]. Плоская задача о движущемся штампе с учетом тепловыделения рассматривалась в работе И.И. Воровича и др. [52]. Для трансверсально-изотропной полуплоскости решение задачи о движущемся с постоянной скоростью вдоль границы штампе было получено B.C. Вольпертом [49]. Используя решения плоской задачи, в диссертационной работе получено решение задачи об осесимметричном напряженном состоянии в изотропных и трансверсально-изотропных сплошных и полых цилиндрах вдоль оси симметрии которых распространяются с постоянной скоростью возмущения. Следует отметить, что впервые связь между решением волнового уравнения в плоском и осесимметричном случаях была установлена, повидимому, В.И. Смирновым и C.J1. Соболевым [138]. Формально, решения осесимметрич-ных стационарных динамических задач имеют ту же форму, что и решения осе-симметричных статических задач для трансверсально-изотропных тел и сводятся к отысканию аналитических функций в полосе. Это позволяет при решении задач обоих классов использовать одни и те же вычислительные алгоритмы. Основные результаты автора диссертации, относящиеся к 3 и 4 главам опубликованы в работах [10,11,101,107,108,109,115].

Пятая глава диссертации посвящена разработке метода исследования межслойных напряжений в многослойных пластинах вблизи ненагруженных кромок круговых и эллиптических технологических отверстий.

Одно из наиболее эффективных направлений оценки прочности многослойных пластин с отверстиями основано на понятии характеристического размера и использовалось в работах Б.Д. Аннина и В.Н. Максименко [24,25]. Вместе с тем, на несущую способность многослойной пластины с отверстиями оказывает влияние и трехмерное напряженное состояние вблизи ненагруженных кромок отверстий, учет которого выходит за рамки данного подхода. Достаточно полный обзор основных методов исследования напряженно- деформированного состояния пространственных неоднородных конструкций, к которым относятся многослойные пластины с вырезами, дан в монографии В.О. Каледина [66].

Наличие межслойных напряжений характерно для многослойных пластин состоящих из однонаправленно армированных слоев при их различной взаимной ориентации. В большинстве работ по исследованию межслойных напряжений рассматриваются прямолинейные свободные кромки. Впервые полный трехмерный анализ в многослойной полосе при одноосном растяжении с использованием основных уравнений теории упругости проводился N.J. Pagano [193] методом конечных разностей. Было выявлено наличие нормальных и касательных напряжений, действующих на границе разделов слоев в узкой зоне, прилегающей к свободной кромке. Ширина зоны (пограничного слоя) примерно равнялась толщине пластины, В дальнейшем задача о межслойных напряжениях вдоль прямолинейной кромки пластины, в силу независимости напряженного состояния от осевой координаты направленной вдоль кромки, была сформулирована как двумерная и подробно исследовалась методом конечных элементов в работах К.Гераковича [54,172] и S.S. Wang [203,204]. Уменьшение размерности задачи позволило более детально исследовать напряженное состояние в слое. Полученные этими авторами результаты дают хорошее качественное согласование с результатами экспериментальных исследований описанных в работах [54,136,168,169].

Влияние свободных криволинейных кромок на напряженное состояние в пограничном слое менее изучено. Решение этой проблемы развивается в двух направлениях.

К первому относятся работы основанные на построении приближенного решения трехмерной задачи для бесконечной многослойной области, состоящей из однородных ортотропных слоев, с круговым отверстием. Прежде всего, это работы Zanq К. [206] и Avrashi J. [161], в которых заранее задан характер поведения межслойных напряжений вблизи отверстия. Несмотря на то, что полученные решения удовлетворяют уравнениям равновесия, полнота этих решений не исследуется и о их достоверности можно судить только путем сравнения с результатами решения аналогичных задач другими методами. В работах A.S. Wang [203,204] с помощью комплексных потенциалов С.Г. Лехницкого задача о межслойных напряжениях вблизи прямолинейной кромки сводится к решению системы дифференциальных уравнений, решение которых позволяет исследовать характер сингулярности в угловых точках контуров слоев. К этому же направлению можно отнести работы L. Chien-Chung [166,167], но для пластин с круговыми вырезами.

Второе направление исследования межслойных напряжений вблизи криволинейных отверстий полностью связано с применением метода конечных элементов. В большинстве работ [160,163,166,167,176,199] исследуются меж-слойные напряжения для пластин с круговыми отверстиями при одноосном растяжении. Схемы армирования пластин выбираются такими, что эффективные модули упругости пластины соответствуют ортотропному материалу, а направление растяжения совпадает с одной из осей ортотропии. Основная задача заключается при этом в детальном исследовании пограничного слоя, что требует значительного числа используемых конечных элементов. Одним из способов приводящих к сокращению объема вычислений является использование суперэлементов, как это сделано в работе Jones R. [173]. Несмотря на значительные усилия, направленные на исследование межслойных напряжений, общей теории, позволяющей прогнозировать характер поведения этих напряжений до сих пор нет и основное направление исследований связано с построением максимально эффективных алгоритмов по расчету межслойных напряжений на базе метода конечных элементов. Основные результаты, полученные автором диссертации по решению этой проблемы, содержатся в работах [1,29,98-100,110112,114,190,192].

Актуальность проблемы. Широкий круг проблем конструирования, эксплуатации, оценки долговечности и надежности реальных элементов конструкций, в значительном ряде случаев, связан с необходимостью максимально адекватно учитывать трехмерность напряженного состояния этих элементов. Одной из основных задач в этом направлении является расчет концентрации напряжений, возникающей как при создании этих элементов (технологические дефекты, неоднородность материала, и т.д.), так и в процессе эксплуатации (внешние нагрузки, контактное взаимодействие). В силу ярко выраженного индивидуального характера источников концентрации напряжений не существует единого метода к их определению. Несмотря на то, что в ряде случаев найдены простые замкнутые (в виде формул) решения, в целом разработка методов расчета концентрации напряжений остается актуальной. Мощным аппаратом при решении данной проблемы являются методы конечных и граничных элементов, позволяющие с высокой точностью определять концентрацию напряжений. Однако, при необходимости исследования концентрации напряжений для заданного класса концентраторов, зависящего от геометрических параметров, изменяющихся в широком диапазоне, применение этих методов становится мало эффективным. Это связано со значительным ростом объема подготовки исходных данных, учитывающих особенности изменения геометрических параметров. Особенно, это относится к задачам исследования концентрации напряжений вызванной взаимодействием ряда концентраторов, в зависимости от расстояния между ними и их формы. Поэтому, наряду с МКЭ и МГЭ представляется перспективным развивать аналитические методы и численно-аналитические методы, позволяющие получать устойчивые решения данной проблемы для широкого диапазона изменения геометрических параметров концентраторов.

Цель работы заключается в разработке аналитических и численно-аналитических методов исследования концентрации напряжений различной природы в элементах конструкций, находящихся в условиях трехмерного напряженного состояния, и разработке на основе этих методов эффективных алгоритмов расчета как напряженного состояния в целом, так и концентрации напряжений.

Научная новизна работы.

• Проведено обобщение метода суперпозиций на решение пространственных задач теории упругости для многосвязных тел следующего вида: изотропная упругая среда, содержащая систему непересекающихся и несоприкасающихся сферических полостей или включений с центрами на одной прямой; трансверсально-изотропная упругая среда, содержащая систему непересекающихся и несоприкасающихся эллипсоидальных полостей или включений симметричных относительно общей оси.

• На основе полученных решений в виде бесконечных систем линейных алгебраических уравнений разработаны устойчивые, эффективные численные алгоритмы расчета концентрации напряжений, в окрестности полостей или включений. Определены области применения полученных решений в зависимости от геометрических и жесткостных параметров. Проведено исследование концентрации напряжений в широком диапазоне изменения этих параметров.

• Разработан общий подход к решению ряда основных осесимметричных задач теории упругости для изотропных и трансверсально-изотропных тел, ограниченных цилиндрическими поверхностями. Получены решения статических задач для бесконечного сплошного цилиндра, бесконечной цилиндрической полости в упругой среде, полого бесконечного многослойного цилиндра. Для этого же класса цилиндрических тел получены осесимметричные решения задач стационарной динамической термоупругости, в предположении, что внешние нагрузки и температурные поля распространяются в направлении оси цилиндра с постоянной скоростью.

• Разработан численно-аналитический метод расчета межслойных напряжений вблизи ненагруженных кромок круговых и эллиптических отверстий в многослойных пластинах из композиционных материалов. Исследованы границы применимости метода для пластин с круговыми и эллиптическими отверстиями.

Методы исследований основаны на: использовании теории функций комплексного переменного, позволяющей представить искомое в заданной области решение в виде аналитических функций заранее известной формы, соответствующих этой области; численных итерационных алгоритмах; применении метода конечных элементов.

Достоверность научных положений, результатов и выводов, содержащихся в диссертационной работе, основывается на:

• корректном использовании основных уравнений теории упругости, аппарата теории функций комплексного переменного и метода конечных элементов;

• исследовании сходимости и устойчивости разработанных численных алгоритмов;

• совпадении с известными численными результатами в частных случаях.

Практическая значимость и реализация результатов исследований заключается:

• в разработке численных алгоритмов расчета трехмерного напряженного состояния и концентрации напряжений, порожденной различного типа концентраторами;

• в численных результатах позволяющих проводить анализ влияния геометрических параметров концентраторов напряжений и характера внешней нагрузки на концентрацию напряжений;

• во внедрении результатов и пакетов программ в расчетную практику заинтересованных организаций.

Работа проводилась по договорам с ФГУП СибНИА им С.А. Чаплыгина и другими предприятиями Минавиапрома в соответствии с правительственными научно-техническими программами "Икарус-МАП", программой Минвуза

РСФСР "Полет", а так же выполнялась в рамках федеральной целевой программы "Государственная поддержка интеграции высшего образования и фундаментальной науки на 1997-2000 годы".

На защиту выносятся:

• разработанный метод определения концентрации напряжений в окрестности системы сферических или эллипсоидальных концентраторов с общей осью симметрии в изотропных и трансверсально-изотропнных средах, а так же алгоритмы и результаты расчета полученные на основе этого метода;

• единый метод определения осесимметричного напряженного состояния в изотропных и трансверсально-изотропных телах, ограниченных цилиндрическими поверхностями в статической и квазистатической постановке;

• численно-аналитический метод расчета межслойных напряжений в окрестности круговых и эллиптических вырезов в многослойных пластинах.

Апробация работы. Основное содержание работы и отдельных ее разделов докладывалось в "Летней научной школе по вопросам прочности конструкций из стекла и керамики "(Обнинск, 1976 г.); на Всесоюзной конференции по теории упругости (Ереван, 1979 г.); на научно-технической конференции по статической прочности авиационных конструкций в ЦАГИ (Жуковский 1980 г.); на Республиканской научно-технической конференции "Концентрация напряжений" (Донецк, 1983г.); на X и XI Всесоюзных конференциях "Конструкции и технология получения изделий из неметаллических материалов" (Обнинск, 1986 г., 1988г.) ); на Межвузовских научно-технических конференциях "Композиционные материалы в конструкциях глубоководных технических средств" (Николаев, 1989 г., 1991г.); на XI Дальневосточной научно-технической конференции "Повреждения и эксплуатационная надежность судовых конструкций" (Владивосток, 1990); на IX Бубновских чтениях по эксплуатационной конструктивной прочности судовых конструкций (Нижний Новгород, 1991); на международной научно-технической конференции "Расчетные методы механики деформируемого твердого тела" (Новосибирск, 1995 г.); на IY Всероссийской научно-технической конференции "Проблемы прочности и усталостной долговечности материалов и конструкций" (Новосибирск, 1997); на международной конференции "Всесибирские чтения по математике и механике" (Томск 1996 г.), на I, II, III Международных российско-корейских научно-технических конференциях CORUS "Научные основы высоких технологий" (Ульсан, Корея, 1997 г., Томск 1998 г., Новосибирск 1999 г.); на межвузовской научной конференции "Численно - аналитические методы решения краевых задач" (Новокузнецк, 1998 г.); на Всероссийской научной конференции "Математическое моделирование процессов в синергетических системах" (Улан-Удэ, 1999 г.); на III школе-семинаре СО РАН "Математические проблемы механики сплошных сред" (Новосибирск, (1999 г.), на международной научно-технической конференции RUSCO "Научные основы высоких технологий" (Новосибирск 2001 г.);

Публикации. По теме диссертации опубликовано 38 печатных работ. Результаты исследований автора, выполненные по заказам промышленности и отраслевых КБ, отражены более чем в 30 научно-технических отчетах.

ЗАКЛЮЧЕНИЕ

1. Разработан метод решения задач исследования концентрации напряжений в зоне взаимного влияния сферических и эллипсоидальных концентраторов в виде полостей или жестких включений с общей осью симметрии в изотропных и трансверсально-изотропных телах.

2. Разработаны и апробированы на примерах решения конкретных задач устойчивые численные алгоритмы определения полей напряжений в окрестности рассматриваемых типов концентраторов напряжений.

3. Исследованы особенности влияния геометрических параметров концентраторов, степени их близости и механических характеристик упругой среды, содержащей концентраторы, на величину концентрации напряжений, при заданных однородных напряженных состояниях на бесконечности.

4. Разработан общий подход к решению широкого класса задач осесимметричной теории упругости для тел ограниченных цилиндрическими поверхностями. С помощью методов теории функций комплексного переменного получены общие решения статических и стационарно динамических задач термоупругости для сплошных и полых изотропных и трансверсально-изотропных цилиндров.

5. Исследовано влияние скорости распространения осесимметричной внешней нагрузки вдоль поверхности сплошного изотропного цилиндра на напряженное состояние.

6. С помощью приближенного решения исследовано напряженное состояние в стенках двухслойного трубопровода вызванное взаимодействием с жесткой подкрепляющей втулкой конечной ширины.

7. Предложен численно-аналитический метод определения межслойных напряжений в многослойных пластинах из однонаправленно армированных слоев вблизи ненагруженных кромок круговых и эллиптических отверстий. Метод позволяет минимизировать число конечных элементов при анализе трехмерного напряженного состояния вблизи кромок отверстий. Показано влияние схемы . армирования многослойной пластины на величину межслойных напряжений.

1. Абраменко Л.В., Олегин ИЛ., Пустовой Н.В. Определение межслойных напряжений в многослойных пластинах с эллиптическими и круговыми отверстиями // Док. СОАН Вышей школы.-2000.-№1.-с.85-91.

2. Абрамян Б.Л., Александров А.Я. Осесимметричные задачи теории упругости //Механика деформируемого твердого тела. -Тр. II Всесоюзн. съезда по тео-рет. и прикладной механике.-М.:Наука.-1966. -С.7-37.

3. Александров А.Я. Некоторые зависимости между решениями плоской и осе-симметричной задач теории упругости и решение осесимметричных задач при помощи аналитических функций // Докл. АН СССР. 1959.- т. 129.- № 4." с.754-757.

4. Александров А.Я. решение осесимметричных задач теории упругости при помощи зависимостей между осесимметричными и плоскими состояниями// ПММ- 1961.- т.25.- вып.5.- с.912-920.

5. Александров А.Я. О решении осеисмметричных задач теории упругости методами теории функций комплексного переменного // Сб. «Некоторые проблемы математики и механики». (К 60-летию акад. Лаврентьева). Изд-во СО АН СССР.- 1962,- № 2,- с.15-24.

6. Александров А.Я. Решение пространственных задач теории упругости при помощи аналитических функций и контурных интегралов для некоторых неосесимметричных тел// Докл. АН СССР.- 1970.-т. 191.-№ 5.- с. 1-11.

7. Александров А.Я., Волъперт B.C. О применении одного метода решения осе-симметричных задач теории упругости к задаче о шаре и о пространстве с шаровой полостью // Изв. АН СССР. ОТН / Механика и машиностроение.-1961,-№6.- с.106-109.

8. Александров А.Я., Олегин И.П. Численный и аналитический методы исследования концентрации напряжений в задачах теории упругости. // Труды на-учно-техн.конф. по статической прочности авиационных конструкций. Из-дат.отдел ЦАГИ.- 1981.- с. 67-75.

9. Александров А.Я., Олегин И.П. Исследование пространственного напряженного состояния в окрестности системы сферических концентраторов в упругих телах. //Тезисы Республиканской научно-техн.конф. /Концентрация напряжений.- Донецк.- 1983.- с.19.

10. Александров А.Я., Олегин И.П. Исследование концентрации напряжений в окрестности системы сферических концентраций // Теор. и прикл. мех.- Ки-ев.-1985.- №16.-с. 3-9.

11. Александров А.Я., Олегин И.П. Пространственное напряженное состояние упругой среды, содержащей осесимметричную систему сферических концентраторов //Изв. АН СССР. МТТ.- № 2,- 1985.-е. 85-92.

12. Александров А.Я., Волъперт B.C., Соловьев Ю.И. Об одной задаче для упруго трехслойного эллипсоида // Актуальные проблемы механики деформируемых сред / Сб. науч. трудов. Днепропетровск.- ДГУ.- 1979.- с. 13-20.

13. Александров А.Я., Соловьев Ю.И. //Пространственные задачи теории упругости.- М.: Наука.- 1975. 462 с.

14. Александров В.М., Ромалис Б.Л. //Контактные задачи в машиностроении. М.: Машиностроение. -1986.- 173 с.

15. Аннин Б.Д., Максименко В.Н. Оценка разрушения пластин из композитных материалов с отверстиями //Механика композитных материалов. -1989. -№2. -С.284-290.

16. Аннин Б.Д., Максименко В.Н. Прогнозирование прочности слоистых композитов с отверстиями и повреждениями типа трещин //VII Всесоюз. съезд по теор. и прикл. механ.: Аннот. докл. -1991. -С. 19.

17. Аржаных И. С., Бондаренко Б.А. Приложение теории функций комплексного переменного к трехмерным задачам математической теории упругости // Труды ин-та математики им. Романовского А.Н. Уз.ССР. 1961. - вып.23. -С.15-19.

18. Баблоян A.A., Мелконян А.П. Осесимметричная задача для полого бесконечного цилиндра с периодически насаженными на него дисками //Изв. АН

19. АрмССР. Механика. -1969. - 22. - № 3. -С.4-9.

20. Белоконь A.B. Смешанные задачи теории упругости для цилиндрических тел // Контактные задачи и их инженерные приложения. М.: изд. НИИмаш. -1969.-С.

21. Бенерджи П., Баттерфилд Р. Методы граничных элементов в прикладных науках /М.: Мир. 1984.- 496с.

22. Бреббия К., Уокер С. Применение метода граничных элементов в технике. //М.: Мир. 1982.- 248с.

23. Беленький М.Я. Некоторые осесимметричные задачи теории упругости. ПММ.- I960. т.24. - вып. 3.-С.34-40.

24. Вайнштейн Ф.А., Хрусталев А.Ф. О напряженном состоянии трансверсаль-но-изотропного цилиндра //Изв. вузов. Математика. - 1961. - №1. С. 1116.

25. Вайнштейн Ф.А., Хрусталев А.Ф. Об одной смешанной задачи теории упругости для трансверсально-изотропного полого цилиндра //Изв. вузов. Математика. - 1963. - № 3. -С57-63.

26. Валов Г.М. Контактная задача об упругой и термоупругой осесимметричной деформации бесконечного сплошного цилиндра //Изв. АН СССР. Механика.- 1965.- № 5. С.60-67.

27. Ванин Г.Ф. Микромеханика композиционных материалов // Киев. -Наукова думка.-1985.-304 с.

28. Василенко A.T., Панкратова Н.Д. Численно-аналитическое решение задачи теории упругости для неоднородной среды со сферическими полостями или включениями //Изв. РАН. МТТ. -1993. -№1. -С.82-88.

29. Васильев В.В. Механика конструкций из композиционных материалов //М: Машиностроение.-1988.-270с.41 .Ватсон Г.Н. Теория бесселевых функций //М.: ИЛ. -1949. -457с.

30. Власов Н.М., Колесов B.C., Тисовский Л.О., Шацкий И.П. Концентрация напряжений в упругом шаре с неконцентрической сферической полостью //Мат. методы и физ.-мех. поля. (Киев). - 1989. -№ 30. - с.37-40.

31. Василенко А.Т., Панкратова Н.Д. Численно-аналитическое решение задачи теории упругости для неоднородной среды со сферической полостью или включением //Изв. РАН. МТТ. 1993. -№ 1. - с.82-88.

32. Ворович И.И., Александров В.М., Бабешко В.А. Неклассические смешанные задачи теории упругости /¡М.: Наука. 1974. - 456 с.

33. Ворогин Т. А. Контактные напряжения, возникающие при тугой посадке жесткой втулки на бесконечный цилиндр //Изв. АН СССР. ОТН. 1957. - №8. -С.51-58.

34. Верюжский Ю.В. Численные методы потенциала в некоторых задачах прикладной механики//Киев. Вища школа. - 1978,- 184с.

35. Волъперт B.C. Решение основных задач осесимметричной теории упругости для полой сферы //Строительная механика. Тр. Новосиб. ин-та инж. ж.д. тр-та,- 1967,- вып. 62.- С.96-102.

36. Волъперт B.C. Пространственная задача теории упругости для эллипсоида вращения и эллипсоидальной полости // МТТ. 1967. - № 3. -С.118-124.

37. Волъперт B.C. Решение некоторых плоских стационарных динамических задач для анизотропного тела при помощи функций комплексного переменного //Тр. Новосиб. ин-та инж. ж.д. тр-та. 1970. - вып.96. - С.

38. Волъперт B.C., Олегин И.П. Осесимметричное напряженное состояние пространства, содержащего систему сферических полостей или включений //Деп. ВИНИТИ № 3266-77. -15с.

39. Волъперт B.C. Концентрация напряжений около эллипсоидальной полости или включения в трансверсально-изотропном теле //Механика деформируемого тела и расчет сооружений. Тр. Новосиб. ин-та инж. ж.д. тр-та. -1972,- вып.137. - С.56-78.

40. Ворович И.И., Пожарский Д. А., Чебаков М.И. Задача о движущемся штампе при учете тепловыделения от трения // Прикл. мат. и мех. -1994. -т.58. -№3. -С.161-166.

41. Галин Л.А. //Контактные задачи теории упругости. М.: ГИТЛ. - 1953. -264с.

42. Геракович К. Кромочные эффекты в слоистых композитах //В сб. статей "Прикладная механика композитов".- М:Мир. 1989. - с.295-341.

43. Гобсон Е.В. Теория сферических и эллипсоидальных функций // М.: ИЛ. -1952. -476с.

44. Головчан В.Т. К решению граничных задач статики упругого тела, ограниченного сферическими поверхностями //Докл. АН УССР. 1974. - А. № 1. - С.61-64.

45. Головчан В.Т. Анизотропия физико-механических свойств материалов //Киев.-1987.-303с.

46. Головчан В.Т. Полные системы двояко- и троякопериодических частных решений уравнений теории упругости в сферических координатах //Прикл. мех. -Киев. -1995. -С.29-34.

47. Градштейн И.С., Рыжик И.М. Таблицы интегралов, суммы рядов и произведений // М.6 Физматгиз. 1962. -1109 с.

48. Григолюк Э.И., Филъшъинский JI.A. Перфорированные пластины и оболочки //М.: Наука.-1970.-554 с.6\.Гюнтер Н.М. Теория потенциала и ее применения к основным задачам математической физики //М. -JI: Гостехтеориздат. 1953. -411с.

49. Зайдес С.А., Дружинина Т.Я. Напряженное состояние при осесимметричном нагружении ограниченной части цилиндра //Сб. "Исслед. по мех. деформир. сред".- Иркут. политехи, ин-т. Иркутск. - 1991. -С.22-27.

50. Зиновьев Б.М., Карманова Т.Ф. К решению пространственных задач теории упругости анизотропного тела методом А.Я.Александрова // Межвуз. сб. науч. тр. Новосибирск. - 1984. -С.20-24.

51. Зиновьев Б.М. Один приближенный метод расчета тел с разрезами //Тр. Но-восиб. ин-та инж. ж.-д. трансп. -1972. -вып.137. -С.105-125.

52. Канторович Л.В., Крылов В.И. Приближенные методы высшего анализа //М.: Физматгиз. 1962. - 707с.

53. Капшивый A.A., Копыстра Н.П., Ломонос JJ.H. Осесимметричное напряженное состояние шара с неконцентрической шаровой полостью //Докл. АН УССР. 1980. - А. - № 9. - С.50-55.

54. Капшивый A.A., Ломонос Л.Н., Стоян H.H. Вторая основная задача об осе-симметричном напряженном состоянии пространства с двумя сферическими полостями //Докл.АН УССР. 1984. - А. - № 6. - С.42-46.

55. Капшивый A.A., Ломонос Л.Н., Стоян H.H. Первая основная задача об осе-симметричном напряженном состоянии шара с неконцентрической шаровой полостью //Докл.АН УССР. 1986. - А. -№ 3. -С.37-41.

56. Квитка А.Л., Ворошко П.П., Бобрицкая С.Д. Напряженно деформированные состояния тел вращения // Киев. -Наукова думка. -1977. -208 с.

57. Копейкин Ю.Д. Интегральные уравнения пространственных задач статики упругого тела// Прикладная механика. Киев. -1965. - т. 1. -вып. 5.-С.29-35.

58. Коган Б.И., Хрусталев А.Ф. Об одной осесимметричной задаче теории упругости для полого цилиндра //ПММ. 1958. - т.22. - вып.5. -С.84-92.

59. Колтунов М.А., Васильев Ю.Н., Черных В.А. Упругость и прочность цилиндрических тел //М.: Высшая школа. -1975. -525 с.

60. Крауч С., Сгарфилд А. Методы граничных элементов в механике твердого тела //М.: Мир. 1987. - 328с.

61. Кристенсен Р. Введение в механику композитов //М.: Мир. -1982. -334с.

62. Кауфман Р.Н. Сжатие упругого шара с неконцентрической шаровой полостью //ПММ.- 1964. т.28. - вып. 4. - с.787-790.

63. Князев П.Н. Интегральные преобразования //Минск.-изд. Вышейшая школа. -1969.-197 с.

64. Кунин И.А., Соснина Э.Г. Эллипсоидальная неоднородность в упругой среде //ДАН СССР. -1971. -т.199. -№3. -С.571-574.

65. Кунин И.А., Соснина Э.Г. Концентрация напряжений на эллипсоидальной неоднородности в анизотропной упругой среде //ПММ. -1973. -т.37. -вып.2. -С.306-315.

66. Кунин И.А., Миренкова Г.Н., Соснина Э.Г. Эллипсоидальная трещина и игла в анизотропной упругой среде // ПММ. -1973. -т.37. -вып.З. -С.531-542.

67. Купрадзе В.Д. Методы потенциала в теории упругости //М.: -Физматгиз. -1963.-472 с.

68. Кущ В.И. напряженное состояние и эффективные упругие модули среды, армированной периодически расположенными сфероидальными включениями //Прикл. мех. -Киев-1995. -31. -№ 3. с.32-39.

69. Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного //М.: -Наука. -1965. -716с.

70. Лебедев H.H. Специальные функции и их приложения //Изд.2. -М.:ФМЛ. -1963.-358с.

71. Лехницкий С.Г. Анизотропные пластинки //М.: ГИТТЛ. -1957. 391с.

72. Лехницкий С.Г. Теория упругости анизотропного тела //Изд. 2. -М.: «Наука». 1977. -416с.

73. Лурье А.И. Пространственные задачи теории упругости //М.-.ГИТЛ. 1955. -492с.91 .Лурье А.И. Теория упругости //М.:Наука. -1970. -940с.

74. Маховиков В.М. О решении пространственной задачи теории упругости для цилиндра //Инженерный журнал. АНСССР. -1964. -т.4. -вып.2. -С.268-275.

75. Максименко В.Н., Олегин И.П. Прогнозирование прочности композитных конструкций. // Учебное пособие. -НГТУ. Новосибирск. - 1994. -89с.

76. Нигирич Ю.Б., Олегин И.П. Осесимметричное напряженное состояние упругой среды в окрестности эллипсоидальных включений //Тез. док. X Всесоюз.конф. "Конструкции и технология получения изделий из неметаллических материалов'. Обнинск. - 1986. -С.35.

77. Николаев А.Г., Проценко B.C. Первая и вторая основные задачи для двухсвязных областей, ограниченных поверхностями сферы и сфероида //Прикл. мат. и мех. 1990. -т.54. - № 1- С. 65-74.

78. Олегин И.П. Пространственное напряженное состояние бесконечного тела, содержащего однопериодическую систему жестких включений //В сб. Механика деформируемого тела и расчет транспортных сооружений. Тр. Новосиб ин-та инж. ж.д. тр-та. - 1980.- С.67-74.

79. Олегин И.П. Напряженное состояние упругой среды в окрестности системы армирующих элементов. // Тез. док. XI Всесоюз.конф. "Конструкции и технология получения изделий из неметаллических материалов".- Изд-во НПО "Технология".- Обнинск. 1988. -С.43.

80. Олегин И.П. Пространственное напряженное состояние упругого тела в окрестности периодической системы сферических полостей //Динамика и прочность авиационных конструкций. Новосибирск. - НЭТИ. - 1989. -С.

81. Олегин И.П., Перепелкин И.Г. Напряженное состояние в полом бесконечном цилиндре, вызванное движущейся нагрузкой //Динамика и прочность элементов авиационных конструкций. Новосибирск. - НЭТИ. - 1990. -Новосибирск. - НЭТИ. -С.

82. Олегин И.П. Осесимметричное напряженное состояние в слоистом транс-версально-изотропном цилиндре // Межвузовский сборник научных трудов "Динамика и прочность конструкций". Новосибирск. -НГТУ. - 1994. -С.20-26.

83. Олегин И.П. Решение некоторых классов задач теории упругости для тел вращения методом суперпозиций // Научный вестник НГТУ. 1996. - № 2. -С. 57-69.

84. Олегин И.П., Расторгуев Г.И. Исследование кромочных эффектов в слоистых композитах с отверстиями // Избранные доклады международной конференции "Всесибирские чтения по математике и механике".- Том И. -Механика. Томск. - 1997. - С.34-38.

85. Олегин И.П. Решение пространственной задачи теории упругости для трансверсально-изотропного тела, содержащего периодическую систему эллипсоидальных полостей // Сибирский журнал индустриальной математики.- 1999.- том 2. №1. - С.117-122.

86. Олегин И.П. Определение напряженного состояния в трансверсально-изотропных цилиндрическихтелах // Научный вестник НГТУ. 2001. - № . -С.

87. Олегин И.П. Формулы переноса для функций Лежандра второго рода в задачах теории упругости // Док. СОАН Вышей школы.-2001.-№2.-С.67-74.

88. Партон В.З., Перлин П.И. Методы математической теории упругости //М.: Наука.- 1981.- 683с.

89. Перлин П.И. К обобщению на пространственный случай одного способа решения основных плоских задач теории потенциала и теории упругости //Докл. АН СССР. 1963.-т. 153.- №5.- С.1033-1036.

90. Перлин П.И. Решение первой основной осесимметричной задачи теории упругости для области, ограниченной эллипсоидом и сферой //Инженерный ж. 1964.- т.5. - №2,- С.275-280.

91. Перлин П.И., Ступак С.Ф. Решение задачи теории упругости для шара, ослабленного сфероидальной полостью //Изв. АН СССР. МТТ. - 1974 - № 5.- С. 165-170.

92. Подилъчук Ю.Н. Деформация упругого сфероида //Прикл. мех. -Киев. -1967.-т. 3.- вып. 12. -С.34-42.

93. Подилъчук Ю.Н., Соколовский Я.И. Напряженное состояние трансвер-сально-изотропной среды с анизотропным сфероидальным включением // Прикл. мех. -Киев. 1991.-т. 3. - вып. 27. -С.22-30.

94. Подилъчук Ю.Н. Обзор по точным решениям пространственных задач теории упругости для трансверсально-изитропных тел канонической формы //Прикладная механика. Киев. -1997-т. 33. - №. 10. -С. -.

95. Положий Г.М. Про одне штегральне перетворення узагальнених ана-л1тичных функцш //В1сник Кшв. ун-ту. сер. астроном., матем., мех. - 1959. - вып. 1,- №2.-С. 19-29.

96. Положий Г.Н. О (p,q) аналитических функциях комплексного переменного и некоторых их применениях //В сб. Исследования по современным проблемам теории функций комплексного переменного. - М.: Физматгиз. -1960. -С.483-515.

97. Положий Г.Н. О краевых задачах осесимметричной теории упругости. Метод р- аналитических функций комплексного переменного //Укр. Матем. Журнал/- 1963.-т. 15.- № 1.-С.25-45.

98. Положий Г.М., Капшивий О.О. Про розв 'язания осесиметричних задач теорп . пружносп для сюнченного цилшдра //Прикл. мехашка. 1961.- т. 7.- вып. 6.-С.616-626.

99. Положий Г.М., Чемерис B.C. До питания про застосувення р-аналггичних функщй в осесиметрш Teopii пружносп //Доповда АН УРСР/- 1958/- № 12. -С.1284-1287.

100. Положий Г.Н., Чемерис B.C. Об интегральных уравнениях осесиммет-ричной теории упругости //В сб. Исследования по современным проблемам теории функций комплексного переменного. М.: Физматгиз. - 1961. -С.399-412.

101. Попов Г.Я. К решению контактных (смешанных) задач теории упругости для бесконечного длинного кругового цилиндра //Изв. АН АрмССр. Сер. физ.-мат. наук. - 1964. - 17. - N 4. -С.

102. Рабинович А.Л. Об упругих постоянных и прочности анизотропных материалов //Труды ЦАГИ. 1946. - № 582. -С.

103. Работное Ю.Н. Механика деформируемого твердрго тела //М.Наука. -1979.-712 с.

104. Развитие теории контактных задач в СССР//М.: Наука. 1976.-496 с.

105. Ройтфарб ИЗ., Ворона Ю.В., Дао К.Т. Применение МГЭ к исследованию напряженно-деформированного состояния массивов, ослабленных полостями //Деп. в УкрНИИТИ (Киев)-1989. -N 2383-УК89. 17с.

106. Ростовцев H.A. Комплексные функции напряжений в осесимметричной контактной задаче теории упругости /ЛТММ. -1953. т. 17. - вып. 5. —С.611— 614.

107. Рыска Л. Влияние кромочных эффектов на напряженное состояние и прочность углеродных композитов при растяжениии //Мех. композит, материалов» (Рига). 1985. - №2. - С. 263-270

108. Свирский И.В. Методы типа Бубнова-Галеркина и последовательных приближений //М: Наука. -1968. -199с.

109. Смирнов В.К, Соболев С.Л. О применении нового метода к изучению упругих колебаний в пространстве с осевой симметрией //Тр. Сейсмолог. Инта АН СССР. 1933.- №29.- 48с.

110. Смирнов Л.Г., Федик И.И. Упругие напряжения в сфере с инородным эксцентрическим сферическим включением //Мат. методы и физ.-мех. поля. -1990. -№31. -С.79-83.

111. Снеддон КН., БерриД.С. Классическая теория упругости //М.:Физматгиз. 1961.-220с.

112. Соболев С.Л. Алгоритм Шварца в теории упругости //Докл. АН СССР. -1936. -4(13). -№6. -С.235-238.

113. Соловьев Ю.К Решение пространственной осесимметричной задачи теории упругости для многосвязных тел вращения при помощи обобщенных аналитических функций//ДАН СССР. 1966.- т. 169.- № 2.-С.310-313.

114. Соловьев Ю.К. Основная смешанная задача осесимметричной теории упругости //ПММ. 1966. - т. 30. - вып. 5. -С.956-962.

115. Соловьев Ю.К. Некоторые вопросы, связанные с решением пространственной осесимметричной задачи теории упругости при помощи обобщенных аналитических функций //СтроителЬная механика. Тр. Новосиб. ин-та инж. ж. д. тр-та. - 1967- вып. 62. -С.5-35.

116. Соляник-Красса КВ. Решение осесимметричной задачи в сферических координатах //Труды Ленинградского политех, ин-та. 1955. - № 178. -С. 17-31.

117. Соляник-Красса КВ. Решение осесимметричной задачи в эллиптических координатах //Труды Ленинградского политех, ин-та. -1958. -№ 195. -С.5-25.

118. Соляник-Красса КВ. Осесимметричная задача теории упругости // М: Стройиздат. -1987. -329с.

119. Уфлянд Я.С. Интегральные преобразования в задачах теории упругости. //Изд. 2-е. -Л.:Наука. -1967. -402с.

120. Сунчелеев Р.Я. Некоторые краевые задачи теории упругости для транс-версально изотропного цилиндра//Изв.АНУССР. -1966. -№5. -С.41-47.

121. Угодчиков А.Г., Хуторянский Н.М. Метод граничных элементов в механике деформируемого твердого тела // Изд-во Казанского ун-та. Казань. -1986.-259с.

122. Хрусталев А.'Ф., Байнштейн Ф.А. /Об одной смешанной задаче упругости для трансверсально-изотропного полого цилиндра //Изв.ВУЗ-ов. -Математика. -1963. -№3. -С. 162-165.

123. Чемерис B.C. Про метод наближеного розв'язання друпл основно осеси-метрично1 задач1 теорп пружност1 //Прикл. мехашка. 1963. - т. 9. - вып. 1. -С.68-76.

124. Чернопинский Д.И. Некоторые особенности расчета деформируемых цилиндров в модифицированных функциях Бесселя //Прикладная механика. -1978. -т. 14. -№6. -С.45-51.

125. Шардаков И.Н., Трояновский И.Е., Тру фанов И.H. Метод геометрического погружения для решения краевых задач теории упругости //Препринт. ИМСС УНЦ АН СССР. -Свердловск. -1984. -66с.

126. Шерман Д.И. Метод интегральных уравнения в плоских и пространственных задачах статической теории упругости //Трукды Всесоюзного съезда по теоретической и прикладной механике. 1960. - изд. АН М., Л. - 1962. -С.10-12.

127. Шерман Д.И. Пространственная статическая задача теории упругости с заданными смещениями на границе //ПММ. 1943- т.7. вып.5. -С.341-360.

128. Штернберг Э. Трехмерная концентрация напряжений в теории упругости //«Механика», Период, сб. переводов иностр. статей. 1958. - вып. 6(52), С.73-80.

129. Ясницкий JJ.H. Метод фиктивных канонических областей в механике сплошных сред //М: Наука. -1992. -125.

130. Ясницкий Л.Н. Композиция расчетной области в методе фиктивных канонических областей // МТТ. -.-№6. -С. 168-172.

131. Akasaka Т., Kenmochi К., Asamo К., Sato Т. Stress distribution around hole of angle-ply laminates under uni-axial tension //Progr. Sci and Eng.Compos. Proc. hlnt.Cont ICCM-IV. Tokyo. Oct 25-28. - 1982. - Vol 1. -P.397-404.

132. Bar-Yoseph Pinchas, Avrashi Jacob. Interlaminar stress analysis for laminated plates containing a curvilinear hole //Comput. and Struct. 1985.-vol. 21. -№ 5. -P.917-932.

133. Carlsson Leif. Interlaminar stress at a hole in a composite member subject to in-plani loading //J.Compos.Mater. 1983. -17-№3. -P.238-249.

134. Chen W.T. Stress concentration around spheroidal inclusion and cavities in a transversely isotropic material under pure shear //Trasection of the ASME . Ser. E. -J-l of Appl. Mech. -1970. -vol.37. -№1. -P.88-96.

135. Chen H.S., Acrivos A. The solution of the equations of linear elasticity for an infinite region containing two spherical inclusions // //Int. J. Solids and Structures. -1978. -v.14. -№2. -P.331-339.

136. Chien-Chung L., Chu-Cheng K. Method for calculating the interlaminar stresses in symmetric laminates containing a circular hole //AIAA Journal. -1992. -30. —№1. -P.197-204.

137. Chien-Chung L., Chu-Cheng K. A new method for interlaminar stresses i symmetrica. laminated plates containing a circular hole // Proc. Nat. Sci. Counc., Rep. China. A. -1994. -18. -№5. -C.485-494.

138. Czarnek R., Post D., Herakovich C.T. Edge effects in composites by moire interferometry //Exp. Techn. -1983. -№7. -P. 18-27.

139. Detcher A.K., Webber J.P.H. Edge effects in uniaxial compression testing of cross-ply carbon-fiber laminates //J. Compos. Mater. 1982. -16. -P.

140. Hamada M., Kodama J. Axisymmetricical tension of an infinite body containing two" spherical cavities //Bull. JSME. 1985. - 28. -№ 237. -P.408-413.

141. Hamada M., Kodama J., Satake K. Axisymmetricical tension of an infinite body containing a series spherical cavities 11 Bull. JSME. -1986. 29. -№ 247. -P.8-13.

142. Herakovich C .7!, Nagarkar A.P. , O'Brien D.A. Failure analysis of composite laminates with three edges //Modern Developments in Composite Materials and Structures. -ASME. -1979. -P. 53-62.

143. Jones R., Callinan R., Teh K.K.,Brown K.C. Analysis of multi-layer laminates using three-dimensional superelements //Int.J. Number. Meth.Eng. 1984. -vol.20.-№3.-P. 583-587.

144. Kasano H., Ogino K., Matsumoto H, Nakahara I. II Three dimensional contact problem of a rigid ring bonded to an elastic hollow cylinder. -Trans. Jap. Soc. Mech. Eng. -1982. -A 48.-№433. -p. 1160-1167.

145. Kasano H., Hara H., Matsumoto H., Nakahara I. Three dimensional stress in a transversely isotropic hollow cylinder with a rigid hoop //Trans. Japan Soc. Mech. Eng. 1986. - A52. -№ 477. - P.355-360.

146. Lucking W.M., Hoa S.V., Sankar T.S. The effect of geometry of interlaminarstress of o° /90° ^composite laminares with circular holes. Journal of Composite Materials. Vol. 17- March 1984, p. 47-56.

147. Masao S., Juhachi O. On the problems of stick and hollow cylinders under the axially symmetrical deformations //Mem. Fac. Technol. Kanazawa Univ. 1968. -v.5-№l. -P.47-54.

148. Mar James W, Lin Kuen Y. Fracture of Boron/Aluminum Composites with Discontinuities //J. Composite Materials. 1977, Vol. 11. -P.405-419.

149. Miyamoto H. On the problem of elastisity theory for an infinite region containing two spherical cavities same size. Proc. Fith Japan Nat. Congr. Appl. Mech. 1955, p. 125-129.

150. Miyamoto H. On the problem of elastisity theory for an infinite region containing three spherical cavities. Proc. Sixs Japan Nat. Congr. Appl. Mech. 1956, p.98-101.

151. Miyamoto H. On the problem of elasticity theory for an infinite region containing more then two spherical cavities same size. Trans. Japan Soc. Mech. Eng., 23, N131, 1957, p. 431-443.

152. Mikata Kenndji. //Res. Repts. Hachinohe Techn Coll. 1972. - №7. -P. 18-23.

153. Mikata Kenndji. IIRes. Repts. Hachinohe Techn. -Coll. 1973. -№ 8. -P. 26-31.

154. Molinar A., El Mouden M. The problem of elastic inclusions at finite concentration //Int. J. Solid and Struct. -1996. -vol. 33. -№ 20-22. -P. 1331-1350.

155. Nishioka T., Atluri S.H. Stress analysis of holes in angle-ply laminates // AIAA/ASME/ASCE/AHS 21st Struct, struct. Dyn.and Mater. Conf., Seatle, Wash. 1980. Collect, teen. pap.Part INew York, N.Y.,1980, P.295-301.

156. Noda N., Ashida F. Transient thermoelastic contact problem of a transversely isotropic cylinder// Ind.-Arch. -1988. -58. -№4. -P.241-251.

157. Noda N., Ogasawara N., Matsua T. Interaction effect between ellipsoidal cavities in an infinite body under uniaxial tension //Trans. Jup. Soc. Mech. Eng. A. -1996. -v.62. -№596. -P.1051-1058.

158. Noda N., Ogasawara N., Matsua T. Interaction effect between ellipsoidal inclusions is an infinite body under asymmetrical uniaxial tension // Trans. Japan Soc. Mech. Eng.A. -1998. -v/-64. -№622. -P. 1577-1582.

159. Olegin I. P. Researh of The Interlaminar Stress in the Composite Plates. //Abstracts The 1st Korea-Russia International Symposium on Science and Technology "KORUS'97". University of Ulsan Republic of Korea, 1997. P. 21.

160. Pagano N.J., Pipes R.B. Interlaminar stresses in composite laminates under uniform axial extension //J. Composite Materials. -1970. -v.4. -P.204-214.

161. Parhi H., Itou S., Atsumi A. Stress in transversely isotropic medium containing an infinite row of spherical cavities //Lett. Appl. and Eng. Sci. -1975. -3. -№6. -P.453-467.

162. Parnes R. Applied tractions on the surface of an infinite cylindrical bore //Int. J Solids and Structure. -1983. -19. -№2. -P.165-177.

163. Popov G.Ya. Problems of stress concentration in the neibourhood of spherical defect //Успехи мех. -1992. -15.-№ 1-2. -P.71-110.

164. Radoc J.R.M., On the solution of problems of dynamic plan elasticity //Quarterly of applied mathematics. 1956. - vol. 14. -№3. -P.44-51.

165. Rodin Gregory, Hwang Yuh-Long. On the problem of linear elasticity for an infinite region containing a finite number of non-interestiyg spherical ingomogen-ius //Int. J. Solids and Struct. -1991. -vol. 27. -№ 2. -P.145-159.

166. Rybicki Edmund F., Schmueser David W. Effect of Stacking Sequence and Lay-Up Angle on Free Edge Stresses Around a Hole in a Laminated Plate Under Tension //J. Composite Materials. -1978. Vol. 12. -P.300-310.

167. РОССИЙСКАЯ ГввУДА=С7й£НК4Ябиблиотека