CC BY
40
УДК 624.071.22:531.011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ГИБКИХ НИТЕЙ С УЧЕТОМ ФИЗИЧЕСКОЙ НЕЛИНЕЙНОСТИ
Д. А. Тарасов, Д. Ю. Семенов
MATHEMATICAL MODELING OF STRESS-STRAIN STATE OF FLEXIBLE YARNS WITH CONSIDERATION OF PHYSICAL NONLINEARITY
D. A. Tarasov, D. Yu. Semenov
Аннотация. Актуальность цели. Объектом исследования являются основные силовые элементы механических систем, выполненные из гибких нитей, работающие по восприятию поперечных нагрузок. Предмет исследования - метод математического моделирования напряженно-деформированного состояния, одновременно учитывающий геометрическую и физическую нелинейность гибких нитей, что дает возможность исследовать работу конструкций за пределами упругой работы материала. Материалы и методы. Предложенный метод основан на теории математического моделирования с использованием аппарата дифференциального и интегрального исчисления функций одной и нескольких переменных, а также положениях сопротивления материалов. Результаты. Предложен метод нелинейного расчета напряженно-деформированного состояния гибких нитей с учетом реальных свойств материалов. Выводы. Сравнение параметров напряженно-деформированного состояния гибкой нити, полученного с помощью вычислительного эксперимента, с результатами, полученными другим методом расчета, подтвердило адекватность предложенного метода.
Ключевые слова: напряжения, деформации, гибкая нить, физическая нелинейность, геометрическая нелинейность.
Abstract. Background. The object of the study are the main power components of mechanical systems made of flexible threads, working on the perception of lateral loads. Subject of research - the method of mathematical modeling of stress-strain state (VAT) simultaneously takes into account geometric and physical nonlinearity flexible filaments that gives the opportunity to explore work structures beyond elastic behavior of the material. Materials and methods. The proposed method is based on the theory of mathematical modeling using the apparatus of differential and integral calculus of functions of one and several variables and positions of strength of materials. Results. The proposed method of nonlinear stress-strain state of flexible filaments taking into account the real properties of materials. Conclusions. Comparison of parameters of stress-strain state of flexible strands, obtained through computational experiments and also with results obtained by another method of calculation has confirmed the adequacy of the proposed method.
Key words: stress, strain, flexible filament, physical nonlinearity, geometric nonlin-
earity.
Введение
В инженерной практике широко распространены системы, в которых основными силовыми элементами, обеспечивающими прочность конструкции при внешнем воздействии, являются стальные канаты. К таким системам относятся воздушные линии электрических проводов, контактные провода электрифицированных железных дорог и трамваев, канаты висячих мостов и покрытий большепролетных производственных и общественных зданий, тросы канатных дорог и кабель-кранов и т. п. [1, 2].
В настоящее время определение напряженно-деформированного состояния (НДС) отдельных стальных канатов ведется с использованием известных методик расчета пологих гибких нитей, вошедших в нормативную литературу [3].
Актуальность настоящей работы вытекает из того, что в указанном нормативном документе расчет рассматриваемых элементов ведется без учета неупругих деформаций материала. Однако в последнее время все чаще находят применение синтетические канаты, материал которых подчиняется нелинейному закону деформирования. В связи с этим определение НДС таких элементов необходимо вести в физически нелинейной постановке задачи, так как отсутствует линейная зависимость между напряжениями и деформациями.
Метод математического моделирования НДС гибких нитей
Рассмотрим гибкую нить, расчетная схема которой представлена на рис. 1. При обозначении вертикальных опорных реакций ЯА и ЯБ примем реакцию ЯА на более высокой опоре. Условимся все параметры, относящиеся к начальному очертанию, обозначать с индексом «0», а к конечной линии равновесия - с индексом «1» [4].
Рис. 1. Расчетная схема гибкой нити: начальное состояние линии равновесия; - конечное состояние
линии равновесия от внешнего воздействия
При моделировании НДС гибких нитей приняты следующие допущения: материал непрерывно заполняет весь объем (гипотеза сплошности); свойства материала одинаковы во всех точках (однородность) и не зависят от направления (изотропность); поперечные нормальные сечения, плоские до деформации, остаются плоскими и после деформации (гипотеза Бернулли); поперечное сечение мало по сравнению с длиной; изменение поперечного сечения при растяжении не учитывается; из внутренних сил возникает лишь осевая сила - натяжения, а жесткость при изгибе и кручении отсутствует [1].
По поводу принятых гипотез следует отметить, что они являются общепринятыми при расчете гибких нитей на внешнее воздействие, и справедливость их подтверждается многолетней практикой эксплуатации конструкций и проведенными натурными и модельными экспериментами [1].
При определении деформированного состояния гибкой нити используются правила построения эпюры изгибающих моментов для балки. Линия равновесия под действием вертикальной нагрузки совпадает с эпюрой изгибающих моментов шарнирно опертой балки тем же пролетом, находящейся под действием той же нагрузки; при этом ординаты эпюры моментов уменьшены делением на величину распора и отложены от хорды АВ, соединяющей точки крепления. Математически это записывается так [4, 5]:
где Мо(ЯоА, Яов, х), М1(Я1А, Я1В, х) - функции изгибающего момента в шарнирно опертой балке пролетом I от действия начальной нагрузки и совместного действия начальной и дополнительной нагрузки, Н • м; Но, Н1 - распор, вызванный начальной нагрузкой и совместным действием начальной и дополнительной нагрузки, Н; х - текущая абсцисса (0 < х < I), м; в - угол наклона хорды АВ, град.
Распор является горизонтальной составляющей опорных реакций в точках крепления и равен по величине горизонтальной составляющей продольных усилий во всех сечениях гибкой нити [6].
Из выражений (1) и (2) видно, что построение линии равновесия сводится к определению распора, так как определение функции изгибающего момента в шарнирно опертой однопролетной балке не вызывает трудностей.
При определении внутренних усилий в однопролетной балке пролетом I поперечная сила равна сумме сил, расположенных по одну сторону от рассматриваемого сечения; изгибающий момент равен сумме моментов всех сил, расположенных по одну сторону от рассматриваемого сечения относительно центра тяжести сечения [4]. Математически внутренние усилия при изгибе можно записать следующим образом:
- при действии нагрузки, вызывающей начальное очертание:
Уо(х)
У1(х)
М0(Я0 А, Я В, х) Но
М1(Я1А, Я1В, х)
Н
+ х§р,
+ xtgP;
(1)
(2)
а(ЯоА,Яов,х) = ЯоА ■ (х > 0) + Яов • (х > I) +
о
о
(3)
Мо (Ra, RoB , x) = RAx(x > 0) + Rob (x -1)(x > l) +
rows( F0) x x
+ X [F0i(x- xFo )(x > xFo)] + Jq0(xx)(x- xx)dxx +Jq(xx)(x-xx)dxx; (4)
i=1 о о
- при совместном действии начальной и дополнительной нагрузки:
Wws(F Fo)
Qx(Ria,RlB,x) = R1A(x > 0) + RlB(x > l) + X [F (x > x^)] +
i=1
rows( Fj) x x x
+ X [F(x > xF1 )J + Jq0(x)dx + Jq1(x)dx + Jq(x)dx; (5)
i=1 ooo
M1 (R1 a, R1B, x) = Rux(x > 0) + R1B (x -1)(x > l) +
rows( Fo) rows( F1)
+ X [F0i (x - xF0,)(x > xF0,)] + X [F1i (x - xfJ(x > xfJ] +
1=1 1=1
+J q0 (xx)(x - xx)dxx + J q1 (xx)(x - xx)dxx + J q(xx)(x - xx)dxx, (6)
000
где Roa, Rob, R\A, R\B - опорные реакции при действии начальной нагрузки и совместном действии начальной и дополнительной нагрузки, Н; F0i, F1i - сосредоточенные начальные и дополнительные нагрузки, Н; xF0i, xF1i - абсциссы точек приложения сосредоточенных начальных и дополнительных нагрузок, м; q0(x), qi(x) - функция распределенной начальной и дополнительной нагрузки, Н/м; q(x) - функция равномерно-распределенной нагрузки от собственного веса, Н/м.
Функция равномерно-распределенной нагрузки от собственного веса постоянна при всех значениях абсцисс и равна
q( x) = pA, (7)
где р - объемный вес материала, Н/м3; А - площадь поперечного сечения, м2.
Распределенная нагрузка задается массивом точек контура для последующей линейной интерполяции функции. Предполагается, что координаты функции нагрузки от деформации гибкой нити не зависят.
В уравнения для поперечных сил и изгибающих моментов включены слагаемые от опорных реакций. Поскольку они не известны до момента определения внутренних усилий, их следует включать в параметры функций.
Определяются опорные реакции из условия равновесия балки:
- при действии нагрузки, вызывающей начальное очертание:
rows (Fo) l l
Roa + Rob + X F0l + J qo(x)dx + J q(x)dx = 0; (8)
i=1 о 0
Mo(RoA , RoB , l) = 0; (9)
- при совместном действии начальной и дополнительной нагрузки:
rows( Fo) rows( F1) III
RA + R1B + X F01 + X F\t + J qo(x)dx + J q1(x)dx + J q(x)dx = 0; (10)
1=1 1=1 000
МД Я1А, Я1В,
I) = о,
(11)
где I - пролет, м.
Для определения распора в точках крепления воспользуемся уравнением неразрывности деформаций [1]. Поскольку распор не известен до момента его определения, а он входит в уравнения для расчета начальной и конечной длины гибкой нити, все последующие выражения запишем в виде функций от распора:
¿о+дед)+Д4 = А(нД
(12)
где ¿о - начальная длина, м; ДЬ(Н1) - упругая деформация, м; ДА - температурная деформация, м; Ь^Н) - конечная длина, м.
В работах [7-9] получены выражения по определению всех членов уравнения неразрывности деформаций, однако определить удлинение гибкой нити при учете физической нелинейности материала по формуле, основанной на законе Гука, уже нельзя:
Н — Н
ДЬ(Н) = ^
1 ЕА
1 +
а< Яа , Я,в , х) Н,
+tgв
ёх,
(13)
где Но, Н\ - распор, вызванный начальной нагрузкой и совместным действием начальной и дополнительной нагрузки, Н; Е - модуль упругости материала, Па; А - площадь поперечного сечения, м2; ^(Я^, Я\В, х) - функция поперечной силы в шарнирно опертой балке пролетом I от совместного действия начальной и дополнительной нагрузки, Н; х - текущая абсцисса (о < х <I), м; в - угол наклона хорды АВ, град.
Учет физической нелинейности материала при моделировании НДС гибких нитей
Для определения удлинения гибкой нити необходимо по диаграмме деформирования, зная напряжения, найти относительную деформацию, а по ней определить искомое удлинение.
В качестве расчетных диаграмм состояния материала, определяющих связь между напряжениями и относительными деформациями, могут выступать любые виды диаграмм, например, такие как криволинейные или кусочно-линейные, но самое главное, отвечающие истинному поведению материала. При этом должны быть указаны основные параметрические точки диаграммы, каковыми являются максимальные напряжения и соответствующие относительные деформации, граничные значения и т.п.
В качестве рабочей диаграммы состояния материала, определяющей связь между напряжениями и относительными деформациями, примем кусочно-линейную диаграмму деформирования, представленную на рис. 2, так как с ее помощью можно аппроксимировать любую кривую.
При такой диаграмме растягивающие напряжения в зависимости от относительных деформаций удлинения материала определяются по формуле
а =
Ее при е < е1;
а1 + Е1(е-е1) при е1 <е<е2;
(14)
а + Е (е-е ) при е <е<е
п п \ п' г п п
где Еп - модуль упрочнения материала на п-м сегменте диаграммы деформирования, Па.
Рис. 2. Диаграмма деформирования
Модуль упрочнения материала на п-м сегменте диаграммы деформирования можно определить из уравнения
Е =ап+1 -ап
еп+1 еп
(15)
Для нахождения напряжений умножим обе части выражения (13) на Е
величину
М(Нг)Е = (И - И0)Е
Ь0
ЕАЬ0
1 +
Ql(Rl а , Яв, х) И,
+ tgP
ёх .
После сокращения и преобразования получим:
а(Н1) =
1 АЬ
1 +
Q1(Я1А , Я1В , х) + ?£р ^
Н
ёх .
(16)
(17)
В параметры функции напряжений введен распор, вызванный совместным действием начальной и дополнительной нагрузки, так как он не известен на момент определения напряжений.
Зная напряжения и то, что относительная деформация равна
дед)
£ = -
Ьо
■, найдем из выражения (14) удлинение гибкой нити:
дед)=
Е
при ст(Н1) <01;
Ьо
о(Н1) — о
при о1 < о(Н1) < о2;
(18)
Ьо
о(Н1) — о
при °п < о(Н1) < °п
После определения удлинения, вызванного изменением внешнего воздействия, а также нахождения начальной и конечной длины гибкой нити по формулам, представленным в работах [1, 3, 4], появляется возможность, используя численные методы, решить нелинейное уравнение (12) и найти значение распора.
После нахождения распора появляется возможность определить продольное усилие, возникающее в поперечных сечениях по длине гибкой нити.
Функция продольного усилия равна [3, 4]
21( х) ^д/н^^сЯ^А^Я^^Н^ .
(19)
Различные, но статически эквивалентные локальные нагрузки приводят к возникновению в гибкой нити одного и того же напряженного состояния, если не учитывать местные напряжения вблизи точки приложения нагрузки. Этим состоянием является линейное (одноосное) напряженное состояние. Известно, что при простом растяжении (а гибкие нити могут воспринимать только растягивающие усилия) во всех сечениях, параллельных оси рассматриваемого элемента, нормальные и касательные напряжения равны нулю. Таким образом, при простом растяжении в каждой точке гибкой нити главные площадки перпендикулярны и параллельны ее оси, а главные напряжения на них соответственно равны [1]
Т{ х)
о1(х) = —-; о2(х) = о3(х) = о .
(2о)
Проверка адекватности предложенного метода
В качестве оценки адекватности разработанного метода математического моделирования проведем сравнительный анализ НДС гибкой нити, полученного разными методами расчета. Для этого рассмотрим предложенный метод и метод конечных элементов (МКЭ) на примере задачи, расчетная схема которой представлена на рис. 3.
У
Рис. 3. Расчетная схема гибкой нити:
- - начальное состояние линии равновесия; - конечное состояние
линии равновесия от внешнего воздействия
Гибкая нить с опорами, расположенными на одном уровне пролетом I = 6 м, сечением А = 0,00031416 м2 нагружена равномерно-распределенной нагрузкой по всему пролету. Первоначальная стрела провеса и нагрузка, вызывающая начальное очертание, отсутствуют. Учет физически нелинейного поведения материала моделируется кусочно-линейной зависимостью между напряжениями и деформациями. Параметры закона нелинейного деформирования, представленного на рис. 4, выражены в напряжениях а1 = 300 МПа, а2 = 500 МПа, а3 = 700 МПа и соответствующих относительных деформациях е1 = 0,00204, е2 = 0,02, е3 = 0,03.
Сохранить закт в файл о Текущий закон Все законы проекта
Подтвердить I Отменить | Справка
Рис. 4. Закон нелинейного деформирования материала
Примем за основные критерии сравнения при определении НДС гибкой нити максимальное продольное усилие, возникающее в сечениях, расположенных на опорах, и прогиб в середине пролета.
Для этого решим ряд задач, последовательно приравнивая равномерно-распределенную нагрузку #1 к значениям в 3, 9 и 2о кН/м для получения напряжений на каждом сегменте заданной диаграммы деформирования.
Расчет МКЭ выполнен в программном комплексе «ЛИРА-САПР» версия 2о15 релиз Я4. Конструкция моделировалась универсальными пространственными стержневыми конечными элементами с учетом физической и геометрической нелинейности (тип 41о), разбивка производилась на 1оо конечных элементов. Расчет велся шаговым нелинейным процессором, предназначенным для решения физически и геометрически нелинейных, а также контактных задач. На каждом шаге матрица жесткости формировалась в системе координат «нового положения» с учетом изменения касательного модуля деформации.
В качестве примера на рис. 5 и 6 представлены результаты моделирования НДС гибкой нити, загруженной равномерно-распределенной нагрузкой #1 = 3 кН/м.
85.1 85.2 85.3 85.4 85.4 85.5 85.5 85.6 85.7
Нелинейное загружение 1 Мозаика N
Единицы измерения - кН
Рис. 5. Эпюра продольных усилий
Рис. 6. Эпюра перемещений
В результате математического моделирования НДС предложенным методом для ряда значений равномерно-распределенной нагрузки, действующей на гибкую нить с заданными физическими и геометрическими параметрами, получены соответствующие значения максимального продольного усилия Тт1(о) = Т^/) и прогиба в середине пролета w(ll2).
В табл. 1 представлены данные, рассчитанные с помощью МКЭ и разработанного метода моделирования.
Таблица 1
Сопоставление результатов расчетов
Критерий оценки Метод расчета Равномерно-распределенная нагрузка, кН/м
3 9 20
Продольное усилие в точках подвеса 7\(0), Т1(1), кН Предложенный метод 85,8 134,3 187,9
МКЭ 85,7 136 185
Погрешность, % 0,12 1,25 1,54
Прогиб в середине пролета ^(//2), м Предложенный метод 0,158 0,308 0,507
МКЭ 0,162 0,307 0,517
Погрешность, % 2,47 0,32 1,93
Заключение
Анализ параметров НДС гибкой нити, полученных с помощью предложенного метода математического моделирования, засвидетельствовал достаточно высокую согласованность с результатами расчета в программном комплексе «ЛИРА-САПР», реализующего МКЭ.
Расхождения значений параметров НДС, взятых за основные критерии сравнения на всех участках принятой диаграммы деформирования, составили не более 5 %, что подтверждает адекватность разработанного метода и принятых допущений при моделировании.
Представленный метод определения напряженно-деформированного состояния гибких нитей с учетом физической нелинейности позволяет решать задачи с учетом реальных свойств материалов, не подчиняющихся закону Гука, что характерно для механических систем с основными силовыми элементами, выполненными из синтетических канатов.
Библиографический список
1. Тарасов, Д. А. Математическое моделирование напряженно-деформированного состояния стальных канатов при действии поперечных статических и ударных нагрузок : дис. ... канд. техн. наук : 05.13.18 / Тарасов Д. А. - Пенза, 2016. - 146 с.
2. Тарасов, Д. А. Конструкция защитного сооружения от удара для железнодорожных переездов / Д. А Тарасов, В. В. Коновалов, В. Ю. Зайцев // Региональная архитектура и строительство. - 2014. - № 1 (18). - С. 111-117.
3. Пособие по проектированию стальных конструкций (к СНиП 11-23-81* «Стальные конструкции») / ЦНИИСК им. Кучеренко Госстроя СССР. - М. : ЦИТП Госстроя СССР, 1989. - 148 с.
4. Тарасов, Д. А. Математическое моделирование напряженно-деформированного состояния стальных канатов / Д. А. Тарасов, В. В. Коновалов, В. Ю. Зайцев // Вестник Саратовского государственного технического университета. - 2013. -№ 4 (73). - С. 215-221.
5. Тарасов, Д. А. Комплекс программ моделирования напряженно-деформированного состояния стальных канатов / Д. А. Тарасов // XXI век: итоги прошлого и проблемы настоящего плюс. - 2013. - № 12 (16). - С. 114-120.
6. Тарасов, Д. А. Компьютерное моделирование определения реакций опор гибких барьеров I В. В. Коновалов, Д. А. Тарасов, В. Ю. Зайцев, Н. В. Байкин I/ Известия Самарской государственной сельскохозяйственной академии. - 2о12. - № 3. -С. 72-79.
7. Тарасов, Д. А. Комплекс программ моделирования работы стальных канатов при действии поперечной динамической нагрузки I Д. А. Тарасов II XXI век: итоги прошлого и проблемы настоящего плюс. - 2о15. - № 2 (24). - С. 48-55.
8. Тарасов, Д. А. Параметрическая оптимизация стальных канатов при действии поперечного удара I Д. А. Тарасов II XXI век: итоги прошлого и проблемы настоящего плюс. - 2о15. - № 2 (24). - С. 56-62.
9. Тарасов, Д. А. Математическое моделирование оптимизации параметров несущих элементов, выполненных из стальных канатов I Д. А. Тарасов, В. В. Коновалов, В. Ю. Зайцев II Интеграл. - 2о12. - № 6. - С. 118-12о.
Тарасов Денис Александрович
кандидат технических наук, начальник отдела инженерного анализа конструкций и защиты интеллектуальной собственности,
ЗАО «ЦеСИС НИКИРЭТ» (г. Пенза) E-mail: den517375@ya.ru
Семенов Дмитрий Юрьевич магистр, инженер-конструктор отдела приводной и специальной техники, ЗАО «ЦеСИС НИКИРЭТ» (г. Пенза) E-mail: kb@cesis.ru
Tarasov Denis Alexandrovich
candidate of technical sciences,
head of department of engineering analysis
of structures and protection
of intellectual property,
Closed Joint-stock company
«Center of special engineering constructions
of research and design Institute
of radio electronic technology»
Semenov Dmitry Yuryevich
master, engineer-designer of department
of driving and special equipment,
Closed Joint-stock company
«Center of special engineering constructions
of research and design Institute
of radio electronic technology»
УДК 624.о71.22:531.о11 Тарасов, Д. А.
Математическое моделирование напряженно-деформированного состояния гибких нитей с учетом физической нелинейности I Д. А. Тарасов, Д. Ю. Семенов II Модели, системы, сети в экономике, технике, природе и обществе. - 2о17. -№ 3 (23). - С. 175-185.
