Solution of a mixed nonaxisymmetric problem of the theory of elasticity for anisotropic bodies of revolution

Abstract


The paper developed a technique for solving mixed nonaxisymmetric problems of the theory of elasticity for bounded bodies of revolution made of a transversely isotropic material under the action of surface forces specified according to a cyclic law. The technique involves the development of the energy method of boundary states, which is based on the concepts of spaces of internal and boundary states, conjugated by isomorphism, which makes it possible to establish a one-to-one correspondence between the elements of these spaces. The internal state includes the components of the tensor of stresses, deformations, and the displacement vector. The boundary state includes efforts and displacements at the boundary of the body. The isomorphism of the state spaces is proved, which allows finding the internal state to be reduced to the study of the boundary state isomorphic to it. The basis is formed on the basis of the general solution of the boundary value problem of elastostatics for a transversely isotropic body of revolution. Orthogonalization of state spaces is carried out, where the internal energy of elastic deformation is used as scalar products in the space of internal states; in the space of boundary states, the work of external forces is used. Finally, finding the desired state is reduced to solving an infinite system of algebraic equations for the Fourier coefficients. The solution of the problem with mixed boundary conditions for a circular in plan cylinder of transversely isotropic coarse dark gray siltstone with anisotropy axis coinciding with the geometric axis of symmetry is presented. The solution is analytical and the characteristics of the stress-strain state have a polynomial form. Explicit and indirect signs of convergence of problem solutions and graphical visualization of the results are presented.

Full Text

Детали машин и механизмов, изготовленные из современных материалов, таких как эластомеры, поликристаллические металлы, керамика, а также горные породы и обладающие значительной анизотропией свойств, пребывают в условиях сложного кинематического взаимодействия с другими телами. Естественно, что силы взаимодействия этих тел не распределяются симметрично относительно, например, оси вращения цилиндрического тела или какой-либо плоскости. Исследование напряженно-деформированного состояния, возникающего в теле и которое носит несимметричный характер, является актуальной научной задачей для анизотропных тел. Решению краевых задач для трансверсально-изотропных тел посвящено множество работ. На сегодняшний момент исследуются частные аспекты данных задач для усложненных по геометрии и структуре материалов, например многосвязные, слоистые и др. Например, в работе [1] получены точные аналитические решения задач о равновесии полых и составных транстропных сфер, находящихся под действием внешнего или внутреннего давления. В работе [2] рассмотрена задача о деформации трансверсально-изотропного цилиндрического слоя под действием нормального давления. Полученные асимптотические формулы позволяют описать поведение слоя с разными жесткостями в трансверсальном и тангенциальном направлениях. В работе [3] с помощью преобразования Фурье решена смешанная краевая задача Дирихле - Неймана для уравнения Пуассона в области, ограниченной двумя параллельными гиперплоскостями. Решение записывается через построенную функцию Грина оператора Лапласа. В работе [4] предлагается подход к определению трехмерного напряженно-деформированного состояния (НДС) многослойного транстропного полупространства в случае воздействия на него нормальной нагрузки. Работа [5] посвящена решению контактной задачи для транстропного полупространства с неизвестной областью контакта. Задача сведена к интегральному уравнению относительно давления в зоне контакта, для решения которого применяется численный метод Галанова. Ряд работ посвящен расчету оболочек. Например, в работе [6] рассмотрены осесимметричные краевые задачи теории упругости для неоднородной трансверсально-изотропной конической оболочки. Решение строится методом асимптотического интегрирования уравнений теории упругости. В работе [7], применяя прямую формулировку метода граничных состояний, решены краевые задачи трехмерной анизотропной теории упругости. Приведены численные примеры для нескольких видов анизотропии. В работе [8] исследовалось предельное нагружение конструкций из трансверсально-изотропных материалов в условиях кусочно-линейной текучести. В работе [9] приводятся доказательства теорем существования и единственности решения упругопластической краевой задачи, основанной на теории пластического течения трансверсально-изотропных тел. Смешанные задачи в теории упругости рассматривались реже, чем задачи с однотипными граничными условиями, однако их исследование проводилось в приложении к различным направлениям механики. Например, в работе [10] рассматривался алгоритм численного решения смешанной задачи теории упругости для тела, имеющего одностороннее контактное взаимодействие с упругим полупространством. В исследовании [11] предложен аналитический метод решения смешанной плоской задачи теории упругости для двухслойной кольцевой области. В статье [12] осуществлены математический и численный анализ асимптотических решений трехмерных статических задач теории упругости со смешанными граничными условиями. В работе [13] Для решения смешанных задач использовались конечно-элементные технологии, основанные на смешанной формулировке, построенной с помощью функционала Рейсснера. В исследовании [14] метод граничных состояний применен для решения смешанных задач теории изотропной упругости, а в статье [15] он получил развитие на класс смешанных задач для неограниченной упругой среды с полостями. В исследовании [16] для односвязного изотропного тела построены параметрические решения, содержащие константы среды в неявном виде. В работе [17] метод граничных состояний применен для решения задачи о движении идеальной жидкости в осесимметричной трубе при смешанных граничных условиях. На сегодняшний момент для трансверсально-изотропных тел вращения средствами метода граничных состояний решена первая основная задача теории упругости при участии массовых сил [18; 19]. Решению второй основной задача теории упругости при одновременном действии на тело массовых сил посвящена работа [20]. По идентичной методике, что и во второй основной задаче, решены основная смешанная [21; 22] и контактная [23] задачи. Особенность решения данных задач заключатся в том, что полученное упругое поле удовлетворяет одновременно заданным граничным условиям и заданным массовым силам. Целью данной работы является решение неосесимметричных смешанных задач теории упругости для трансверсально-изотропных тел вращения. 1. Постановка задачи Рассматривается упругое равновесие трансверсально-изотропного конечного и односвязного тела вращения с осью анизотропии, совпадающей с геометрической осью симметрии. Требуется восстановить упругое поле в области V по заданным поверхностным усилиям на части границе и перемещениям на части границе (рис. 1). Естественно, что . Рис. 1. Трансверсально-изотропное тело вращения Fig. 1. The transversely isotropic body of revolution Статические поверхностные усилия и ненулевые перемещения точек границы изменяются по циклическому закону (синуса или косинуса). Объемные силы отсутствуют. 2. Определяющие соотношения для среды Для однородной трансверсально-изотропной среды в цилиндрических координатах , , имеют место следующие соотношения. Дифференциальные уравнения равновесия при отсутствии объемных сил [24; 25]: ; ; (1) , Соотношения Коши [24]: ; ; ; (2) ; ; . Уравнения совместности деформаций [26]: ; ; ; ; (3) ; . Обобщенный закон Гука [24]: ; ; (4) ; ; ; . Здесь: u, v, w - компоненты вектора перемещений u вдоль оси r, θ, z соответственно; , , , , , - компоненты тензора деформаций; , , , , , - компоненты тензора напряжений; и - модули упругости соответственно в направлении оси z и в плоскости изотропии; - коэффициент Пуассона, характеризующий сжатие вдоль оси r при растяжении вдоль оси z; - коэффициент Пуассона, характеризующий поперечное сжатие в плоскости изотропии при растяжении в этой же плоскости; и - модули сдвига в плоскости изотропии и перпендикулярной к ней. 3. Метод решения Для решения поставленной задачи используется метод граничных состояний (МГС) [27]. Основу метода составляют пространства внутренних и граничных Г состояний: ; . (5) Внутреннее состояние определяется наборами компонент вектора перемещений, теноров деформаций и напряжений (индекс k в правой части помещен наверх): . Воспользуемся при построении решения основных задач механики уравнением Клапейрона при отсутствии объемных сил [25; 28]: . (6) На основе равенства (6) можно назначить скалярные произведения в пространствах состояний. Скалярное произведение в пространстве внутренних состояний выражает внутреннюю энергию упругого деформирования (например, для 1-го и 2-го внутреннего состояний): , (7) причем в силу коммутативности состояний среды: . Граничное состояние определяется наборами компонент вектора перемещения точек границы , поверхностными усилиями : , , (8) где - компонента нормали к границе. В пространстве граничных состояний Г согласно (6) скалярное произведение выражает работу внешних сил по поверхности тела S (например, для 1-го и 2-го внутреннего состояний): , причем в силу тождества Бетти и соотношения Клапейрона: . В случае гладкой границы и в силу (6) оба пространства состояний являются гильбертовыми и сопряжены изоморфизмом. По определению каждому элементу соответствует единственный элемент , причем это соответствие взаимно-однозначное: . Это позволяет отыскание внутреннего состояния свести к построению изоморфного ему граничного состояния. После построения базисов пространств внутренних и граничных состояний, их необходимо проортонормировать. Ортонормирование базиса пространства осуществляется по разработанному рекурсивно-матричному алгоритму ортогонализации [29], где в качестве перекрестных скалярных произведений принимается (7). Алгоритм основан на процессе Грама - Шмидта, переписанном в форме, использующей лишь перекрестные скалярные произведения элементов исходного базиса, которые сведены в матрицу Грама. Если в процессе ортогонализации на k-м шаге встречается некоторый элемент базиса внутренних состояний , алгоритм на этом шаге выдаст «0» (нулевое), если этот элемент является линейной комбинацией элементов , …, . Для сохранения ортогональности выходных элементов и для предотвращения деления на ноль при ортонормировании алгоритм делает проверку на нулевые элементы и исключает их. На их место идут следующие элементы исходного базиса внутренних состояний, и процесс повторяется. Ортонормированный базис Г редуцируется из ортонормированного базиса внутренних состояний, используя выражения (8) и (1). Окончательно проблема сводится к разрешающей системе уравнений относительно коэффициентов Фурье, разложения искомых внутреннего и граничного состояний в ряд по элементам ортонормированного базиса: ; , или в развернутом виде: ; ; ; . (9) Ортонормированный базис позволяет для элементов базиса граничных состояний записать следующие выражения: ; (10) Представим слагаемые из (10) в следующем виде (нижний индекс v в развернутых выражениях опущен): ; и, подставляя их в (10), получим: . Группируя слагаемые и обозначая: (11) легко убедиться что . Преобразуем следующим образом: базисные упругие характеристики , заменяем заданными, и перебор будем осуществлять по индексу j, образуя тем самым матрицы коэффициентов: ; ; . Следует отметить, что матрица B является кососимметричной (, ). Матрица-столбец коэффициентов Фурье рассчитывается так: , (12) где N - число используемых элементов базиса. Окончательно решение имеет вид (9). Тестирование коэффициентов Фурье осуществляется подстановкой одного из базисных элементов с соответствующими ГУ в качестве заданного, при этом должны выполняться условия , где n - номер тестируемого элемента, остальные коэффициенты Фурье должны равняться нулю. 4. Формирование базиса Основную сложность формирования решения в МГС является конструирование базиса внутренних состояний, который опирается на общее или фундаментальное решение для среды; также возможно использование каких-либо частных или специальных решений. В работе [24] с помощью метода интегральных наложений установлена зависимость между пространственным напряженным и деформированным состоянием упругого трансверсально-изотропного конечного тела без полостей и определенными вспомогательными двумерными состояниями, компоненты которого зависят от двух координат z и y (переменных). При установлении зависимости используется следующий прием. Упругое тело, напряженное состояние которого требуется изучить, рассматривается как часть некоторого бесконечного цилиндра с осью , параллельной образующей цилиндра. С телом связана система координат . Меридианное сечение тела совпадает с плоскостью поперечного сечения бесконечного цилиндра с осями координат zy (направление плоскости zy, ось z общая для тела и цилиндра). Предполагается, что цилиндр находится в некотором двумерном напряженном состоянии, не меняющемся вдоль образующей. Компоненты этого состояния , , , , , , определяют плоскую деформацию с депланацией сечения цилиндра в плоскости zy. Эти же компоненты определяют напряженное состояние заданного упругого тела, так как оно является частью цилиндра. Для определения пространственного напряженного состояния тела рассматривается ряд цилиндров, отличающихся направлением образующей, или углом поворота относительно оси z. Представляя последовательно тело вырезанным из каждого такого цилиндра при (m - число цилиндров), редуцируется ряд напряженных состояний, суперпозиция которых и дает суммарное трехмерное состояние. В общем случае напряженное состояние каждого цилиндра в процессе поворота изменяется, и трехмерное состояние тела будет неосесимметричным (зависящим от угловой координаты ). Например, для компоненты вектора перемещения пространственного состояния имеет место выражение: Перейдя к пределу при , сумма заменяется интегралом. При последующей замене переменой интегрирования окончательно связь между перемещениями неосесимметричной деформации и плоской деформации соответствующего цилиндра имеет вид [24]: ; ; (13) ; , где a и b - пределы суммирования. Соответственно для компонент тензора напряжений имеют место следующие выражения: ; ; ; ; (14) ; ; ; ; ; ; ; ; ; ; ; ; ; . Здесь компоненты с верхним индексом pl соответствуют компонентам плоского вспомогательного состояния. В качестве плоских вспомогательных состояний используется плоская деформация и депланация, возникающие в сечениях цилиндра, имеющего в каждой точке плоскость упругой симметрии, параллельную плоскости zy (направление плоскости zy) [24]: ; ; ; ; ; ; ; , где константы , и определены упругими параметрами материала; ; и - комплексные корни характеристического уравнения [24]; функции - аналитические по своим переменным. Придавая функциям поочередно значения: , (15) можно сгенерировать набор плоских вспомогательных состояний (15) и по зависимостям (13) и (14) построить совокупность пространственных состояний, образуя базис пространства внутренних состояний (5). 5. Решение задач Исследуется упругое равновесие трансверсально-изотропного кругового в плане цилиндра из горной породы алевролита крупного темно-серого [28]. После процедуры обезразмеривания параметров задачи, аналогия которой приведена в работе [30], упругие характеристики материала составили: ; ; ; ; . Цилиндр занимает область . При построении базиса внутренних состояний необходимо стремиться к наибольшей простоте вида функций, описывающих компоненты упруго поля. Рассмотрим сначала базис, формируемый из левых частей выражений (13) и пределами суммирования и (аналогично и для напряжений): ; (16) В данном случае задача будет разрешима, если на поверхностях тела компоненты заданных сил , , содержат тригонометрические функции , , соответственно, например: ; (17) В противном случае коэффициенты Фурье (12) будут равны нулю. Если формировать базис из правых частей выражений (13): ; ; , (18) то решение существует, если силы , , содержат тригонометрические функции , , соответственно. Если в выражениях (16) и (18) использовать пределы суммирования и , то для заданных поверхностных сил вида (17) решения не существует, в этом случае приближенное решение задачи ищется для заданных функций вида или . В случае, когда заданные силы имеют вид , (19) уже необходимо использовать выражения (13) в полной мере с пределами суммирования и . При этом возможно получение не только приближенных, но и строгих решений. Если при sin и cos в выражении (19) разные коэффициенты, например , , то решения не существует. Это связано с одинаковыми коэффициентами (единица) при соответствующих функциях в базисных выражениях (13), (14). В этом случае целесообразно воспользоваться принципом независимости действия сил и решить две отдельные задачи, в каждой из которых заданы и , а полученные упругие поля сложить. В случае, когда заданные силы зависят от или , , то в выражениях (13), (16), (18) необходимо использовать пределы суммирования . Для последнего случая и приведем пример решения задачи. Цилиндр жестко закреплен по верхнему торцу. На нижнем торце действуют усилия (рис. 2). Рис. 2. Граничные условия к задаче для цилиндра Fig. 2. Boundary conditions to the problem for a cylinder После процедуры ортонормирования и исключения линейно-зависимых элементов базисный набор для компонент вектора перемещения представлен в табл. 1 (показано 6 элементов). В табл. 1 и 2 истинное значение показанной величины равно показанному значению, умноженному на Усеченная до матрица коэффициентов (11) представлена в табл. 2 (i - строка, j - столбец). При решении использовался базис из 71 элемента. Приведем значения для восьми коэффициентов Фурье: ; ; ; ; ; ; ; . График, иллюстрирующий «насыщение» суммы Бесселя (левая часть неравенства Бесселя), представлен на рис. 3. Это является косвенным признаком сходимости решения. Искомые характеристики упругого поля определяются по зависимостям (9). Проверка результата и оценка точности осуществляются сопоставлением заданных ГУ с восстановленными граничными условиями в результате решения (рис. 4). Здесь и далее на графиках заданные () и восстановленные () ГУ изображены в масштабе. Например, истинное значение на левом Таблица 1/Table 1 Перемещения ортонормированного базиса () Displacements of an orthonormal basis () U v w 0 Таблица 2/Table 2 Матрица коэффициентов () Coefficients matrix () 1 1 1 1 1 1 1 1 Рис. 3. Сумма Бесселя Fig. 3. Bessel sum графике рис. 4 равно значению на графике, умноженному на коэффициент . С целью выявления максимальной погрешности, угол q выбран таким образом, чтобы тригонометрические функции имели максимальные значения. На рис. 5 показана верификация граничных условий в зависимости от угловой координаты q для ненулевых компонент заданной распределенной силы на поверхности . Характеристики полученного напряженно-деформированного состояния, имеющие полиномиальный вид, представлены в виде изолиний на рис. 6 (в явном виде необозримы). Истинное значение показанной величины равно соответствующему значению на изолиниях, умноженному на коэффициент . На изолиниях рис. 6, а, показана область . На рис. 6, b, показаны изолинии некоторых компонент тензора напряжений в зависимости от угловой координаты q. Показано сечение . Полученные компоненты упругого поля удовлетворяют всем уравнениям теории упругости для трансверсально-изотропного тела (1)-(4). Таким образом, в работе сформулирована методика решения неосесимметричных краевых задач смешанного типа для трансверсально-изотропных тел вращения в случае, когда ненулевые усилия или перемещения точек границы заданы по циклическому закону. Методика представляет собой развитие метода граничных состояний, а именно представлен новый способ формирования Рис. 4. Верификация граничных условий на участках поверхности Fig. 4. Verification of boundary conditions on surface areas Рис. 5. Верификация граничных условий на участке границы S2 в зависимости от угловой координаты Fig. 5. Verification of the boundary conditions at the border section S2 depending on the angular coordinate , , , , , , , , , , , , , , , , , а , , , , b Рис. 6. Изолинии компонент полученного упругого поля Fig. 6. Isolines of the components of the obtained elastic field базиса внутренних состояний, основанный на общем решении задачи о плоской деформации. Исследован прием декомпозиции базиса под условия конкретной задачи. Предложенная методика, однако, не является общей для любого класса рассматриваемых областей (односвязных и многосвязных). Скорость сходимости рядов зависит от граничных условий и условий внутри области, а так же от геометрии тела. Решение является аналитическим и имеет полиномиальный вид, что позволяет легко проводить анализ полученных характеристик напряженно-деформированного состояния.

About the authors

D. A Ivanychev

Lipetsk State Technical University

References

  1. Фукалов А.А., Кутергин А.В. Точные аналитические решения задач о равновесии упругих анизотропных тяжелых тел с центральной и осевой симметрией и их приложения // Вестник Нижегородского университета им. Н.И. Лобачевского. - 2011. - № 4 (4). - С. 25-26.
  2. Стружанов В.В. Сагдуллаева Д.А. Осесимметричные деформации трансверсально-изотропного цилиндрического слоя под действием нормального давления // Вестник СПбГУ. - 2015. - Серия 1. - Т. 2 (60). - Вып. 3. - С. 426-430.
  3. Алгазин О.Д., Копаев А.В. Решение смешанной краевой задачи для уравнения Лапласа в многомерном бесконечном слое // Вестник МГТУ им. Н.Э. Баумана. Серия: Естественные науки. - 2015. - № 1. - С. 3-13. doi: 10.18698/1812-3368-2015-1-3-13.
  4. Круподеров А.В. Функции Грина для трансверсально-изотропных оснований // Вестник БНТУ. - 2011. - № 5. - С. 54-60.
  5. Пожарский Д.А., Давтян Д.Б. Трёхмерная контактная задача для трансверсально изотропного тела // Вестник ДГТУ. - 2013. - № 7/8 (75). - С. 22-26. doi: 10.12737/2016.
  6. Ахмедов Н.К., Мехтиев М.Ф., Шахвердиева Г.Н. Анализ осесимметричной задачи теории упругости для неоднородной трансверсально-изотропной конической оболочки // Известия вузов. Северо-кавказский регион. Естественные науки. - 2015. - № 2. - С. 5-11.
  7. Игумнов Л.А., Марков И.П., Пазин В.П. Гранично-элементное решение краевых задач трехмерной анизотропной теории упругости // Вестник Нижегородского университета им. Н.И. Лобачевского. - 2013. - № 1 (3). С. 115-119.
  8. Семыкина Т.Д., Цуканова Л.П. Расчет предельных нагрузок для конструкций из трансверсально-изотропных материалов // Вестник Воронежского государственного технического университета. - 2011. - Т. 7, № 4. - С. 233-236.
  9. Кодиров А.У. Решение задач для упругопластических трансверсально-изотропных тел // Бюллетень науки и практики. - 2015. - Т. 5, № 2. - С. 10-13. doi: 10.33619/2414-2948/39/01.
  10. Станкевич И.В. Численное решение смешанных задач теории упругости с односторонними связями // Математика и математической моделирование. - 2017. - № 5. - С. 40-53. doi: 10.24108/mathm.
  11. Божкова Л.В., Рябов В.Г., Норицина Г.И. Смешанная плоская задача теории упругости для двухслойной кольцевой области // Известия Московского государственного технического университета МАМИ. - 2011. - № 1 (11). - C. 217-221.
  12. Соболь Б.В. Об асимптотических решениях трехмерных статических задач теории упругости со смешанными граничными условиями // Вестник Нижегородского университета им. Н.И. Лобачевского. - 2011. - № 4 (4). - С. 1778-1780.
  13. Станкевич И.В. Математическое моделирование задач теории упругости с использованием МКЭ на основе функционала Рейсснера // Символ науки. - 2017. - № 4 (2). - С. 21-25.
  14. Пеньков В.Б., Пеньков В.В. Применение метода граничных состояний для решения основной смешанной задачи линейного континуума // Известия Тульского государственного университета. Серия: Математика. Механика. Информатика. - 2000. - Т. 6, № 2. - С. 124-127.
  15. Пеньков В.Б, Саталкина Л.В., Шульмин А.С. Основная смешанная задача для сферической полости в упругом пространстве // Известия Тульского государственного университета. Серия: Естественные науки. - 2014. - Вып. 1. - Ч. 1. - С. 207-215.
  16. Пеньков В.Б., Новикова О.С., Левина Л.В. Построение полнопараметрических аналитических решений в основной смешанной задаче эластостатики односвязного тела. Вестник Самарского государственного технического университета. Серия: Физико-математические науки. - 2018. - Т. 22, № 3. - С. 586-598. doi: 10.14498/vsgtu1603.
  17. Penkov V.B., Polikarpov M.V., Levina L.V. Efficient solutions of mixed-type axial symmetry problems for perfect fluids. Proceedings - 2020 2nd International Conference on Control Systems, Mathematical Modeling, Automation and Energy Efficiency, SUMMA 2020. - 2020. - P. 52-55. doi: 10.1109/SUMMA50634.2020.9280583.
  18. Иванычев Д.А. Метод граничных состояний в решении первой основной задачи теории анизотропной упругости с массовыми силами // Вестник Томского государственного университета. Математика и механика. - 2020. - № 66. - С. 96-111. doi: 10.17223/19988621/66/8.
  19. Ivanychev D.A. The solution of boundary value problems of various types with consideration of volume forces for anisotropic bodies of revolution // Вестник МГТУ им. Н.Э. Баумана. Серия: Естественные науки. - 2021. - № 4 (97). - С. 59-70. doi: 10.18698/1812-3368-2021-4-57-70.
  20. Иванычев Д.А. Метод граничных состояний в решении второй основной задачи теории анизотропной упругости с массовыми силами // Вестник Томского государственного университета. Математика и механика. - 2019. - № 61. - С. 45-60. doi: 10.17223/19988621/61/5.
  21. Иванычев Д.А. Метод граничных состояний при решении смешанной задачи теории анизотропной упругости с массовыми силами // Вестник Томского государственного университета. Математика и механика. - 2021. - № 71. - С. 63-77. doi: 10.17223/19988621/71/6.
  22. Solving the mixed problem of elasticity theory with mass forces for transversal-isotropic body. Proceedings - 2020 1st International Conference on Control Systems, Mathematical Modelling, Automation and Energy Efficiency, SUMMA 2020. - P. 56-61. doi: 10.1109/SUMMA50634.2020.9280697.
  23. Иванычев Д.А. Решение контактной задачи теории упругости для анизотропных тел вращения с массовыми силами // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2019. - № 2. - С. 49-62. doi: 10.15593/perm.mech/2019.2.05.
  24. Александров А.Я., Соловьев Ю.И. Пространственные задачи теории упругости (применение методов теории функций комплексного переменного). - М.: Наука, 1978. - 464 с.
  25. Новацкий В. Теория упругости. - М.: Мир, 1975. - 872 с.
  26. Лурье А.И. Теория упругости. - М.: Наука, 1970. - 940 с.
  27. Пеньков В.Б., Пеньков В.В. Метод граничных состояний для решения задач линейной механики // Дальневосточный математический журнал. - 2001. - Т. 2, № 2. - С. 115-137.
  28. Лехницкий С.Г. Теория упругости анизотропного тела. - 2-е изд. - М.: Наука, 1977. - 416 с.
  29. Саталкина Л.В. Наращивание базиса пространства состояний при жестких ограничениях к энергоемкости вычислений // Сборник тезисов докладов научной конференции студентов и аспирантов Липецкого государственного технического университета. - Липецк: ЛГТУ, 2007. - С. 130-131.
  30. Левина Л.В., Новикова О.С., Пеньков В.Б. Полнопараметрическое решение задачи теории упругости односвязного ограниченного тела // Вестник ЛГТУ. - 2016. - № 2 (28). - С. 16-24.

Statistics

Views

Abstract - 78

PDF (Russian) - 21

Cited-By


PlumX


Copyright (c) 2022 Ivanychev D.A.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies