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

На сегодняшний день существуют модели механики деформируемого твердого тела, в которых деформация среды описывается не только вектором перемещений, но также вектором поворота. В этих моделях, в отличие от классической (симметричной) теории, напряжённое состояние описывается несимметричным тензором напряжений. Как отмечено в работе [47], модель классической теории упругости хорошо совпадает с экспериментами, проводимыми с конструкционными материалами (сталь, алюминий, бетон) при напряжениях, остающихся в пределах упругости материала. Значительное различие между теорией и экспериментом возникает в тех случаях, когда существенными являются градиенты напряжения. Это имеет основное значение при концентрации напряжений вокруг отверстий и выточек. Расхождение между экспериментом и классической теорией упругости появляется также в задачах о колебаниях, при распространении волн и при вынужденных высокочастотных (ультразвуковых) колебаниях. Это происходит из-за того, что при высокочастотных колебаниях и достаточно малых длинах волн неизбежно сказывается влияние микроструктуры материала. Наконец, классическая теория упругости не описывает с необходимой точностью явления, происходящие в зернистых средах и при прохождении акустических волн через кристаллы, поликристаллические структуры и полимеры [47].

История развития учёта вращения в деформируемых телах рассматривается в монографии [15].

Первое упоминание об изучении вращения в трехмерном пространстве было сделано У. Гамильтоном в 1848 г. в его фундаментальной работе [94], посвящённой теории кватерниона.

О важности учёта моментных напряжений говорилось в работе В. Фойхта (1887г.) [126].

Э. и Ф. Коссера обобщили и развили работы Г. Кирхгофа, А. Клебша, П. Дюгема, В. Фойхта и в 1909 г. появилась теория, согласно которой каждая материальная точка континуума наделяется свойствами твёрдого тела путем учёта ротационных степеней свободы [72]. Согласно концепции братьев Коссера, при изучении напряжённого состояния твёрдого деформируемого континуума необходимо наряду с обычными напряжениями (сила на единицу площади) вводить в рассмотрение моментные напряжения (момент силы на единицу площади).

В 1911 году появилась работа Леру [108], в которой материальная точка наделялась способностью к микродеформации.

Отдельные работы по теории микрополярных сред публиковались в 20-е, 30-е и 40-е годы, но наибольшее развитие эти теории получили в конце 50-х — 60-х годов прошлого столетия в работах Дж. Эриксена и К. А. Трусделла [79], Э. Л. Аэро и Е. В. Кувшинского [1], В. А. Пальмова [49], В. Новацкого [116], [117], А. Г. Ерингена и Е. С. Сухуби [80].

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

1. Псевдосреда Коссера. Теория среды со «стеснённым вращением» .

Такую среду часто называют псевдоупругой средой Коссера или псевдосредой Коссера.

Описание основных положений данной теории можно найти в работах Э. Л. Аэро и Е. В. Ку вши некого [1], Р. Д. Миндлина (R.D.Mindlin [111]), Р. Д. Миндлина и Г. Ф. Тирстена (R.D.Mindlin and H.F.Tierstin [112]), Ю. Н. Немиша ([44], [45]), В. Т. Койтера (W.T.Koiter [100]), Н. Ф. Морозова ([42], [41]) Г. Н. Савина ([58], [55]), А. И. Каландии [20] и др.

В теории псевдосреды Коссера считается, что перемещения й точек этой среды и их жёсткие малые повороты со связаны зависимостью: — rotu. 2.

Таким образом, для псевдосреды Коссера имеется одна независимая кинематическая неизвестная — перемещения гГ, и в рассмотрение вводятся несимметричные тензора напряжений, а и моментных напряжений Д.

Этот вариант несимметричной теории понижает её полноту, так как число физических констант для изотропного упругого тела сокращается до четырёх. Например, часто используются Е — модуль Юнга, 7 — коэффициент Пуассона, I — постоянная, имеющая размерность длины и Вбезразмерная постоянная, называемая модулем изгиба ([111], [42], [41], [55], [20], [67]).

Кроме этого, получаемая структура уравнений такова [47], что если, в частности, на поверхности упругого тела заданы перемещения, то не удаётся произвольно задать нормальную составляющую вектора поворота.

Несмотря на эти недостатки, имеется достаточно большое число работ, посвящённых псевдосреде Коссера. Предложено несколько общих теорем, методов интегрирования и дано решение ряда задач. Так, Миндлин и Тирстен [112] обобщили представление Папковича-Нейбера для статических задач, а также получили фундаментальное решение в бесконечном упругом пространстве. Их теоретические выводы были проиллюстрированы несколькими примерами.

Боджи и Стернберг (D.В.Bogy and E. Sternberg [70]) рассматривали задачи для случая плоско-деформированного состояния. Были обобщены решения Файлона и задача о штампе на континуум псевдосреды Коссера. Особенно интересными являются следствия полученных в [70] результатов, касающихся сингулярных решений для плоско-деформированного состояния.

Можно отметить также решение задачи об изгибе кругового цилиндра [67], плоской граничной задачи о действии сосредоточенной силы на бесконечной плоскости с круговым отверстием [97], задачи для бесконечной упругой изотропной области, ослабленной конечным числом произвольно расположенных несоприкасающихся круговых отверстий [46], задачи о деформировании плоского кольца [82].

2. Теория среды Коссера. Эта теория была развита в 60−70-х годах прошлого столетия независимо несколькими исследователями: В. Новацким (W.Nowacki [47]), В. А. Пальмовым [49], H. Schaefer [122], Э. Л. Аэро и Е. В. Кувшинским [2] и др.

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

В этом варианте упругое поведение изотропной линейной среды характеризуется шестью упругими константами ([47], [49], [50]): две постоянные Ляме и четыре новые константы, характеризующие микроструктуру В случае квадратично-нелинейной среды количество новых констант увеличивается до девяти [15].

Во многих работах (напр. [50], [120], [115]) отмечается, что псевдосреда Коссера является следствием среды Коссера при условии стремления одной из новых упругих констант к бесконечности.

Известны точные аналитические решения ряда задач для среды Коссера, несмотря на значительные трудности при разрешении получающихся дифференциальных уравнений равновесия или движения. Например, найдена концентрация напряжений вблизи кругового отверстия [50], решены задачи о действии сосредоточенной силы и сосредоточенного момента в безграничном упругом пространстве ([120], [114]), о равновесии полупространства ([110], [76], [77]), изгиба пластины под действием сосредоточенной силы, представленной в виде распределенного на небольшой поверхности давления [8]. Рассмотрены неоднородная задача изгиба пластины из материала Коссера и материала с ограниченными локальными поворотами по модели Гриоли-Тупина [90] и задача для прямоугольного параллелепипеда [12].

Также следует отметить работу [6], в которой получено решение уравнений термоупругости, содержащих произвольные постоянные, работу [91], в которой рассмотрена обобщённая задача плоско-напряжённого состояния для упругого слоя, и построены решения для вектора перемещений и вектора бесконечно малого поворота в случае упругого слоя из материала Коссера и работу [59], где в случае тонкой пластинки для решения поставленной краевой задачи трёхмерной несимметричной теории упругости излагается асимптотический метод интегрирования определяющих уравнений трёхмерной несимметричной теории упругости.

В работе [103] предложен подход к построению точных аналитических решений некоторых одномерных и двумерных статических краевых задач в рамках несимметричной теории упругости для среды Коссера. С использованием данного подхода получены точные аналитические решения для следующего ряда задач: задача о сдвиге плоского бесконечного слоя (пластины), закреплённого по обоим краям, под действием силы тяжестизадача о растяжении бесконечной полосы (пластины), ослабленной центральным круговым отверстием (задача Кирша) — задача о кручении жёстко закреплённого по внешнему контуру кольца за счёт поворота внутреннего контура на фиксированный уголзадача о деформировании жёстко закреплённого по внешнему контуру кольца за счёт сдвига внутреннего контура на фиксированную величину. На основе анализа полученных решений введены экспериментально измеряемые макропараметры, откликающиеся на «моментное» поведение среды. Безразмерная форма записи полученных в [103] аналитических решений позволила наглядно установить принципиальное различие безразмерных моментных и классических решений. А именно, безразмерное моментное решение зависит от характерного геометрического размера, а классическое — нет.

3. Континуум Jlepy (градиентная модель). К понятию моментных напряжений приводит и учёт зависимости энергии деформаций от высших градиентов вектора перемещений. Впервые на целесообразность учета высших градиентов перемещений указал Jlepy [108]. Деформированное состояние при этом определяется двумя тензорами: тензором макродеформации второго порядка и градиентом микродисторсии третьего порядка. Градиент микродисторсии связан с вектором перемещений и не связан с вектором поворота. Следовательно, вращение частиц среды в этом случае является стеснённым [15]. Напряжённое состояние определяется объёмной плотностью внутренней энергии, через которую вычисляются тензор напряжений и тензор третьего порядка «двойных напряжений», антисимметричная часть которого является тензором моментных напряжений Д [57].

В случае физической нелинейности в определяющие соотношения этой модели входят, помимо констант Ляме, семь констант Ландау, определяющих нелинейность, и две новые константы, характеризующие микроструктуру [15]. Для линейной среды общее количество констант сокращается, как и в случае псевдосреды Коссера, до четырёх.

4. Микроморфная среда Миндлина-Эрингена. Данная теория развита Р. Д. Миндлиным в работе [113] и А. К. Эрингеном [80].

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

Меры деформаций микроморфной среды являются обобщением деформационных характеристик двух вышеназванных моделей — континуума Коссера и континуума Jlepy.

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

5. Прочие теории. Среди прочих можно выделить работу Э. Л. Аэро и Е. В. Кувшинского [2], где развита 45-константная теория. В работе [80] предлагается мультиполярная теория, количество физических констант которой определяется её степенью. Обобщение несимметричной теории на случай анизотропии было сделано H. Neuber [115], где упругое поведение новой анизотропной (в самом общем случае анизотропии) среды будет характеризоваться 171 упругой постоянной.

Различные аспекты моделей несимметричной среды можно найти также в работах [5], [9], [18], [37], [40], [54], [60], [92], [93], [100], [61], [124], [89], [86], [19], [119], [60], [65], [66], [68], [73], [75], [78], [81], [87], [88], [95], [98], [99], [121].

Известны обобщения несимметричной теории на случай термоупругости и больших деформаций. Также известны решения ряда динамических задач, например, систематическое изложение современной теории распространения и взаимодействия упругих волн в твёрдых телах с микроструктурой можно найти в работе В. И. Ерофеева [15].

Имеется всего несколько работ по реализации процедур идентификации моделей такого типа для нескольких конкретных материалов. Измерение констант упругости на основе статических экспериментов проводилось в работе [83]. Более точные динамические (в частности, ультразвуковые) эксперименты использовались в работах: [16], [57], [56] для идентификации моделей Леру и псевдосреды Коссера, [84], [14], [104] для идентификации линейной среды Коссера, [51] и [52] для нелинейной среды Коссера (для смеси).

Так в работах [57], [56] указывают на экспериментальную возможность оценить численное значение одной из новых констант моментной теории путем измерения скорости ультразвука в зависимости от частоты колебаний. Такие измерения были проделаны на лабораторной установке, использующей известный метод наложения ультразвуковых импульсов, модулированных высокой частотой. Поперечные ультразвуковые колебания (в виде радиоимпульсов) вводились нормально в прямоугольные образцы из поликристаллических металлов технической чистоты и принимались с противоположной стороны образца. С помощью электронно-счётного частотомера в опыте фиксировалась частота повторения радиоимпульсов, обратная величина которой была равна двойному времени пробега ультразвукового импульса в образце. Отсюда при известной базе прозвучивания (длине пути пробега луча) подсчитывались значения скорости звука с точностью порядка Ю-4.

С помощью динамических экспериментальных методов в [104] делается попытка определить механические характеристики несимметричной упругости. Используется так называемый метод размерного эффекта, основанный на зависимости от размера жёсткости при кручении и жёсткости при изгибе, которые имеют место при изгибе и кручении балок с круговыми поперечными сечениями. Экспериментальная аппаратура, используемая в эксперименте, содержит дисковидный постоянный магнит (для приложения крутящего момента к образцу), специальное кольцо Гельмгольца, которое производит однородное магнитное поле в результате протекания электрического тока. Крутящий момент образца пропорционален току в кольце Гельмгольца. Изгиб и кручение достигаются с помощью определённым образом ориентированного кольца Гельмгольца. Измерение угловых перемещений выполнялось с помощью интерференционного метода.

Что касается вопроса об экспериментальном обнаружении «момент-ного» эффекта, то в некоторых работах (напр. [111], [50]) сопоставление решений несимметричной теории с классическими решениями осуществляется на основе анализа показателя концентрации напряжений и его зависимости от характерного размера концентратора. Иллюстрируется факт возрастания (или убывания) показателя концентрации напряжений по сравнению с классической теорией при уменьшении характерного размера. Необходимо отметить, что данный факт представляется интересным, однако использование показателя концентрации в качестве экспериментально измеряемого параметра является весьма проблематичным. Так, например, попытка зафиксировать изменение показателя концентрации методом фотоупругости не увенчались успехом, так как разрешающая способность этого метода слишком низка на необходимом характерном размере концентратора [48]. Как установил Миндлин [111], порядок характерной постоянной материала /, появляющейся в теории псевдосреды Коссера и имеющей размерность длины, для идеальных кристаллов и аморфных материалов такой же, как длина трещиныдля поликристаллов и зернистых материалов она несколько больше размеров толщины. Несмотря на то, что фактор концентрации напряжений, измеренный методами фотоупругости, совпадает с расчётным значением, выполненным по методике классической теории упругости, максимальное напряжение из-за присутствия моментных напряжений в образце с надрезом должно быть меньше, чем-то, которое получается без учёта этих напряжений. Причина в том, что величина I для оптически активных материалов, применяющихся в фотоупругости, достаточно мала. Обычно это не учитывается.

Имеется ряд работ, в которых получены аналитические решения в рамках несимметричной теории упругости применительно к изгибу [100] и кручению [102] стержней различного поперечного сечения. Койтер [100], рассмотрев задачу изгиба бруса, имеющего прямоугольное сечение, с учётом влияния моментных напряжений, пришёл к выводу, что эффективным методом экспериментального определения новых характеристических постоянных упругого материала являются испытания на изгиб. Сопоставление моментного и классического решений осуществлялось на основе анализа изменения жёсткости на изгиб или кручение в зависимости от характерного размера. Да, жесткость с точки зрения экспериментальной реализации является хорошо измеряемым параметром. Но полагать, что экспериментальные измерения жёсткости при изгибе (кручении) позволят продемонстрировать «моментный» отклик среды, не приходится. Это связано с тем, что в этих задачах отсутствует одно из необходимых условий яркого проявления «моментного» поведения среды, а именно отсутствие высокого градиента напряжений. Эксперименты, описанные в [48], подтверждают это.

Так, в работе [123], учитывающей влияние толщины пластины на жёсткость, проведено несколько испытаний балки на изгиб. Исследовались пять балок разной толщины, изготовленные из дисперсноупрочнённого алюминиевого сплава. Полученные экспериментальные данные позволяют заключить, что влияние толщины на значение модуля Юнга почти незаметно, а влияние крутящего момента находится в пределах погрешности измерений.

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

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

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

Метод конечных элементов является одним из базовых численных методов, успешно применяемый в классической теории упругости. В настоящее время этот метод прочно вошел в расчётную практику. Наличие большого числа монографий, например [11], [17], [53], [118], а также обзоров по различным аспектам метода конечных элементов в рамках классической теории упругости позволяет в рамках данной работы не приводить обзора литературы в этой области с позиций анализа достоинств и недостатков этого метода. Поэтому остановимся на обзоре литературы по методу ко нечных элементов в рамках несимметричной теории упругости.

Идея распространить этот метод на решение задач моментной теории упругости появилась в 80-х годах прошлого столетия в работах R.S.Lakes, R.L.Benedict, N. Nakamura [105], [106]. В качестве конечных элементов выбирались плоские треугольные напряжённые элементы [105], а так* же постоянные по деформации трёхузловые треугольные элементы и че-тырёхузловые изопараметрические элементы [106]. Программа тестировалась на примере вычисления коэффициента концентрации напряжений в окрестности отверстия изотропного микрополярного материала, для которого существует точное аналитическое решение [105]. В [106] исследовалась задача о деформировании эллиптического контура в полосе упругого материала Коссера (микрополярного материала) под действием растягивающего усилия.

В последующих работах рассматривались вопросы применения различных типов конечных элементов для задач несимметричной теории упругости. В [109] применялись четырёхугольные конечные элементы для решения задач линейной теории упругости по определению показателей напряжённо-деформированного состояния в случае микрополярного упругого материала. Применение плоских конечных элементов с вращательными степенями свободы можно найти в работах [71], [13], [85], [96]. В [71] рассмотрены вопросы эффективности и точности применения треугольного конечного элемента с тремя узлами, в каждом из которых присутствуют две поступательные степени свободы и одна вращательная. В [13] приводится матрица жёсткости плоского треугольного элемента с вращательными степенями свободы в узлах, а численные примеры показывают достаточно высокую сходимость предлагаемого элемента. Линейно-упругий анализ в исследовании четырёхугольных смешанных конечных элементов проводился в [96] для моментной теории упругости в плоском случае. Новая сеточная конечно-элементная модель типа Вороного разработана в [85] для анализа микрополярного термоупругого напряжённого состояния в произвольно неоднородных материалах. Приведено несколько иллюстрирующих численных примеров для подтверждения эффективности данной модели в анализе полей температур и напряжений для микрополярных упругих материалов.

Ряд работ посвящен решению различных задач в рамках несимметричной теории упругости с применением метода конечных элементов. В [107] рассматривается локализация концевых нагрузок, действующих на плоскую полосу микрополярного материала. Конечно-элементный подход для реализации плоских и осесимметричных задач моментной теории упругости рассмотрен в [7]. При моделировании горного массива применялась модель континуума Коссера, а элементы Коссера представлены традиционными конечными элементами [74].

В ряде работ отмечается, что для определения наиболее «информативных» экспериментов с точки зрения идентификации материальных констант материала в различных разделах механики и физики оказывается эффективным аппарат анализа чувствительности.

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

Математические достоинства такого подхода заключаются в том, что нелинейная задача анализа чувствительности рассматривается при помощи методов, которые обладают преимуществами математических свойств линейных (при фиксированных переменных проектирования) операторов состояния [62].

В [3] изложены и проанализированы три метода расчёта чувствительности по проектным переменным: метод виртуальных нагрузок, пространства состояний, проектного пространства. На основании этого исследования авторы делают вывод о предпочтительном использовании метода пространства состояний.

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

Анализ чувствительности при наличии ограничений рассмотрен Э.Дж.Хогом и Я. Аророй в [4], обзор методов вычисления производных чувствительности применительно к конечно-элементным моделям конструкции приведён в работах [64], [125].

Анализ обзора литературы позволяет сделать следующие выводы.

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

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

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

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

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

3. Постановка и отработка экспериментов, позволяющих зафиксировать факт «моментного» поведения упругих материалов.

Научная новизна работы состоит в том, что:

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

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

3. Оценена «информативность» полученных численных решений с точки зрения проведения эксперимента для определения неизвестных констант в модели среды Коссера.

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

Достоверность полученных результатов подтверждается:

1. Сравнением отдельных численных результатов с известными аналитическими решениями.

2. Анализом реализованных экспериментов методами численного эксперимента.

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

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

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

3. Созданием новых экспериментов, направленных на установление фактов «моментного» поведения упругих материалов.

Апробация работы. Основные положения и результаты работы докладывались на: Всероссийской конференции молодых учёных «Математическое моделирование в естественных науках» (Пермь, 2000 г.), VIII Всероссийском съезде по теоретической и прикладной механике (Пермь, 2001 г.), XXX Summer School «Advanced Problems in Mechanics» (St. Petersburg, 2002), Всероссийской научной конференции «Современные проблемы математики, механики, информатики» (Тула, 2002 г.), конференции молодых учёных «Неравновесные процессы в сплошных средах» (Пермь, 2002 г.), Зимней школе по механике сплошных сред (Пермь, 2003 г.), III Всероссийской конференции по теории упругости (Ростов-на-Дону, 2003 г.), European Congress on Computational Methods in Applied Sciences and Engineering (Jyvaskyla, Finland, 2004 г.), XXI International Congress of Theoretical and Applied Mechanics (Warsaw, Poland, 2004r.).

Структура работы. Работа состоит из введения, четырёх глав, заключения и списка литературы.

В первой главе приведены основные положения несимметричной теории упругости: уравнения равновесия, геометрические и определяющие соотношения, вариационные уравнения.

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

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

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

Выводы по работе.

В данной работе получены следующие новые результаты:

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

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

3. Оценена «информативность» полученных численных решений с точки зрения проведения эксперимента для определения неизвестных констант в модели среды Коссера.

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

5. Выполнен анализ экспериментальных данных и приведены численные расчёты, определяющие программу экспериментальных работ по установлению фактов «моментного» поведения упругих материалов.