Механика деформируемых тел. Упругость. Деформация

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

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

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

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

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

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

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

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

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

На основе численных расчетов при разных значениях радиусов и свойствах материалов изучено сопряжение форм профилей контактирующих цилиндров в плоской постановке c привлечением точного решения Мусхелишвили и формул, полученных T.T. Loo и N.Y. Troy. Выполнен анализ влияния кривизны на напряженно-деформированное состояние в зоне контакта цилиндров. Проведено сопоставление и показана разница между решениями Герца и T.T. Loo, N.Y. Troy.

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

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

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

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

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

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

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

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

Асимптотическим методом изучается краевая задача плоской микрополярной теории упругости в тонкой прямоугольной области. Построены внутренняя (одномерная) модель и погранслой. Изучается задача сращивания, при помощи которой краевые условия плоской задачи на кромках прямоугольника перераспределяются между внутренней (одномерной) и погранслойной задачами. Построенная внутренняя одномерная модель трактуется как прикладная модель изгиба микрополярных упругих балок. Показывается идентичность прикладных моделей микрополярных балок, построенных на основе асимптотического метода и метода гипотез.

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

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

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

Исследована трехмерная контактная задача для составного упругого шероховатого клина при скользящей заделке внешней грани. Клин составлен из двух клиновидных слоев разного угла раствора, соединенных скользящей заделкой. Предполагается, что один из слоев несжимаем. Для вывода интегрального уравнения контактной задачи используется метод интегральных преобразований Фурье и Конторовича–Лебедева. Ядро интегрального уравнения контактной задачи зависит от вспомогательной функции, удовлетворяющей интегральному уравнению Фредгольма второго рода. При помощи метода нелинейных граничных интегральных уравнений определена область контакта и основные механические характеристики. Ранее аналогичная задача рассматривалась для однородного клина.

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

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

Представлена методика определения точек бифуркации для различных видов оболочек – цилиндрической оболочки с двумя типами винтовой анизотропии при продольном ее сжатии и круглой выпуклой мембраны. На первый взгляд столь различные задачи объединяет два фактора: 1) тип прикладной теории, на основе которой выведены уравнения равновесия (прикладная теория Киргоффа–Лява для непологих оболочек); 2) методы интегрирования нелинейных и линеаризованных задач. При численном анализе второй задачи (задачи устойчивости сферического купола) показано, что экспериментально проявляющаяся чувствительность сферического купола к несовершенствам связана с большим количеством близко расположенных точек бифуркации по неосесимметричным модам. Показано, что целенаправленным внесением небольших технологических изменений в форму купола можно добиться устранения большинства из этих точек бифуркации для обеспечения работы оболочки в осесимметричном режиме.

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

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

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

Рассмотрена задача об устойчивости нелинейно-упругого цилиндра с внутренними напряжениями, источником которых является клиновая дисклинация. На основе модифицированного полуобратного представления нелинейной теории упругости получены уравнения нейтрального равновесия для двух моделей материалов: полулинейного и Блейтца и Ко общего вида. На основе бифуркационного подхода изучено существование решений линеаризованных краевых задач. Получены критические значения кратности удлинения/укорочения цилиндра от параметра дефекта при растяжении и сжатии.

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

Рассмотрена контактная задача о взаимодействии абсолютно жесткого шара (штампа) с внутренней поверхностью трехслойного сферического основания. Предполагается, что внешняя поверхность сферического основания закреплена, слои имеют различные упругие постоянные и между собой жестко соединены. Задача сведена к решению интегрального уравнения (ИУ) первого рода, трансформанта ядра которого построена в явном аналитическом виде. Решение ИУ построено с помощью методов: симптотического, прямых коллокаций и конечных элементов. Проведен расчет распределения контактных напряжений, размеров области контакта и перемещения штампа. Проведено сравнение результатов расчетов, полученных этими методами.

Получены интегральные уравнения (ИУ) трехмерных контактных задач для упругого полупространства, составленного из двух клиновидных слоев, соединенных скользящей заделкой. Клиновидный слой, примыкающий к слою, в который вдавливается штамп, несжимаем (резино-металлическое сочленение). Внешняя грань этого слоя свободна от напряжений либо подчинена условиям скользящей заделки. Для решения вспомогательных краевых задач о действии заданной нормальной нагрузки применен метод комплексных интегральных преобразований Фурье и Конторовича — Лебедева, позволивший свести их к системам ИУ Фредгольма второго рода, решения которых затем вошли в ядра ИУ контактных задач. Для решения контактных задач использован метод Галанова.

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

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

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

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

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

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

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

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

Для обеспечения проектирования головных обтекателей летательных аппаратов, рис.1, представлены возможности многопараметрического расчета и оптимизации напряженно-деформированного состояния (НДС) элементов узла клеевого соединения тонкостенной керамической оболочки вращения с металлическим кольцом (опорным шпангоутом), рис.2. Система основных уравнений, описывающих статическое взаимодействие двух соосных оболочек, соединенных внахлестку с помощью промежуточного упругого слоя, преобразована к виду одного обыкновенного дифференциального уравнения второго порядка относительно вектора искомых функций. Компонентами искомой вектор-функции являются меридиональные и окружные перемещения керамической оболочки и обечайки шпангоута. Граничные условия тоже записаны в перемещениях. Таким образом, обеспечена возможность выполнения расчетов НДС с применением алгоритмов матричной прогонки [1] к решению полученной краевой задачи. При этом оптимизация напряженно-деформированного состояния обеспечивается путем изменения и подбора переменной по длине соединения толщины оболочек и толщины промежуточного слоя, а также путем изменения характеристик жесткости шпангоута и клея.

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

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

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

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

Предложен вариант метода асимптотического осреднения композиционных материалов с термоупругими характеристиками, позволяющий вычислять эффективные коэффициенты линейного теплового расширения (КЛТР) композитов. Сформулированы локальные задачи термоупругости на ячейке периодичности (ЯП) композитов. Дана вариационная формулировка задач термоупругости на ЯП, для численного решения которых применен метод конечных элементов. Для программной реализации этого метода использован программный комплекс, разработанный в Научно-образовательном центре «Симплекс» и на кафедре «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Представлены примеры численного решения локальной задачи термоупругости для композиционного материала на основе керамических волокон и полимерной матрицы. Рассчитаны эффективные КЛТР композиционных материалов с пространственным расположением керамических волокон и полимерной матрицей при различных температурах. Показано, что наличие процессов термодеструкции полимерной матрицы при высоких температурах приводит к немонотонной зависимости КЛТР от температуры. Предложенный алгоритм позволяет вычислять КЛТР композитов практически с произвольными структурами армирования волокнами и матрицами, претерпевающими физико-химические превращения при высоких температурах. В отличие от большинства известных приближенных методов расчета КЛТР разработанный метод дает возможность находить точные в математическом смысле значения этих коэффициентов.

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

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

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