Все выпуски

На основе кусочно-квадратичной интерполяции получены полуаналитические аппроксимации потенциала двойного слоя вблизи и на границе двумерной области. Для вычисления интегралов, образующихся после интерполяции функции плотности, используется точное интегрирование по переменной $\rho=\left(r^2-d^2\right)^{1/2}$, где $d$ и $r$ — расстояния от наблюдаемой точки до границы области и до граничной точки интегрирования соответственно. Доказана устойчивая сходимость таких аппроксимаций с кубической скоростью равномерно вблизи границы класса $C^5$, а также на самой границе. Также доказано, что использование для вычисления интегралов стандартных квадратурных формул не нарушает равномерной кубической сходимости аппроксимаций прямого значения потенциала на границе класса $C^6$. При некоторых упрощениях доказано, что использование для вычисления интегралов стандартных квадратурных формул влечет отсутствие равномерной сходимости аппроксимаций потенциала внутри области вблизи любой граничной точки. Теоретические выводы подтверждены результатами численного решения задачи Дирихле для уравнения Лапласа в круговой области.

На основе кусочно-квадратичной интерполяции получены полуаналитические аппроксимации нормальной производной потенциала простого слоя вблизи и на границе двумерной области. Для вычисления интегралов, образующихся после интерполяции функции плотности, используется точное интегрирование по переменной $\rho =(r^{2} -d^{2} )^{1/2} $, где $d$ и $r$ — расстояния от наблюдаемой точки до границы области и до граничной точки интегрирования соответственно. Доказана устойчивая сходимость таких аппроксимаций с кубической скоростью равномерно вблизи границы класса $C^{5}$, а также на самой границе. Также доказано, что на границе аппроксимации по аналогии с точной функцией терпят разрыв, величина которого пропорциональна значениям интерполированной функции плотности, но могут быть доопределены на границе до функций, непрерывных или на замкнутой внутренней, или на замкнутой внешней приграничной области. Теоретические выводы о равномерной сходимости подтверждены результатами вычисления нормальной производной вблизи границы единичного круга.

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

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

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

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

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

Получено альтернативное представление для невязки, представляющее собой положительно определенную квадратичную форму от новых конечных разностей, заданных на границе. Элементы матрицы формы выражаются через многочлены Чебышёва, матрица обратима и такова, что обратная матрица имеет трехдиагональный вид. Эта особенность позволяет получить для спектра матрицы верхние и нижние оценки, не зависящие от размерности N. Данное обстоятельство позволяет провести исследование на качество аппроксимации для разных размерностей N и весовых коэффициентов ω∈[-1,1]. Показано, что наилучшее приближение дает параметр ω=0, а невязка стремится к нулю с ростом N.

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

Получено альтернативное представление для невязки, представляющее собой сумму двух положительно определенных квадратичных форм от новых конечных разностей, заданных на границе. Матрица первой формы имеет второй порядок и очевидный спектр. Элементы второй матрицы порядка N + 1 выражаются через многочлены Чебышева, матрица обратима и такова, что обратная матрица имеет трехдиагональный вид. Эта особенность позволяет получить для спектра матрицы верхние и нижние оценки, не зависящие от размерности N. Данное обстоятельство позволяет провести исследование на качество аппроксимации для разных размерностей N и весовых коэффициентов ω ∈ [−1, 1]. Показано, что наилучшее приближение дает параметр ω = 0, а невязка стремится к нулю с ростом N.

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