Научная статья на тему 'Расчет прочности оболочечных конструкций АПК бункерного типа'

Расчет прочности оболочечных конструкций АПК бункерного типа Текст научной статьи по специальности «Механика и машиностроение»

CC BY
117
22
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
тонкостенная конструкция бункерного типа / конечный элемент / условия сочленения цилиндра с конусом / концентрация напряжений / the thin-walled design of the hopper type / the final element / the conditions of the joint of the cylinder with the cone / the stress concentration

Аннотация научной статьи по механике и машиностроению, автор научной работы — Клочков Юрий Васильевич, Николаев Анатолий Петрович, Фомин Сергей Денисович, Соболевская Татьяна Алексеевна, Андреев Александр Сергеевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Клочков Юрий Васильевич, Николаев Анатолий Петрович, Фомин Сергей Денисович, Соболевская Татьяна Алексеевна, Андреев Александр Сергеевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

CALCULATION OF THE STRENGTH OF SHELL STRUCTURES OF AGROINDUSTRIAL COMPLEX OF BUNKER TYPE

Thin-walled elements are widely used in engineering structures of agriculture (agribusiness). They are structural elements of grain cleaning units, bunker structures for various purposes, as well as pipeline systems. Therefore, the development of modern numerical algorithms for strength calculations of the bunker type structure is an urgent task. In this article, an algorithm for calculating the strength of the thinwalled engineering structure of the bunker type as a structure consisting of two types of shells: cylindrical and conical. To determine the stress-strain state of such an aggregate, a numerical finite element method is used with the adoption of a curvilinear finite element developed by the authors in the form of a quadrangular fragment of the middle surface of the shell with nodal unknowns in the form of displacements and their first derivatives. The stiffness matrix of the used finite element is developed on the basis of the theory of thin shells with the use of the direct normal hypothesis. The kinematic parameters of the cylindrical shell are taken as the main unknown in the nodes at the junction of the cylindrical and conical shells. To obtain the expressions of the nodal parameters of a conical shell using the basic nodal unknown of the cylindrical shell used the kinematic conditions and the equilibrium equations of statics. The obtained relations are implemented in the matrix form and used in the transformation of stiffness matrices and vectors of nodal loads of finite elements of the conical shell adjacent to the boundary of the joint with the cylindrical shell. On the example of calculating the strength of a thin-walled bunker-type structure consisting of a cylinder and a truncated cone, the efficiency of the developed algorithms is shown and the correctness of the obtained numerical values of strength indicators is confirmed. The analysis of the numerical results allowed us to conclude that the stress concentration coefficient in the zone of joint of the cylinder with the cone increases proportionally to the angle of inclination of the cone to its axis.

Текст научной работы на тему «Расчет прочности оболочечных конструкций АПК бункерного типа»

***** ИЗВЕСТИЯ *****

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА: № 2 2019

НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

УДК 539.3 DOI: 10.32786/2071-9485-2019-02-34

РАСЧЕТ ПРОЧНОСТИ ОБОЛОЧЕЧНЫХ КОНСТРУКЦИЙ АПК

БУНКЕРНОГО ТИПА

CALCULATION OF THE STRENGTH OF SHELL STRUCTURES OF AGROINDUSTRIAL COMPLEX OF BUNKER TYPE

Ю.В. Клочков1, доктор технических наук, профессор А.П. Николаев1, доктор технических наук, профессор С.Д. Фомин1, доктор технических наук, доцент Т.А. Соболевская1, кандидат технических наук А.С. Андреев1, аспирант, ассистент

Yu.V. Klochkov1, A.P. Nikolaev1, S.D. Fomin1, T.A. Sobolevskaya1, A.S. Andreev1

1Волгоградский государственный аграрный университет

1 Volgograd State Agrarian University

Дата поступления в редакцию 11.02.2019 Дата принятия к печати 11.05.2019

Received 11.02.2019 Submitted 11.05.2019

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

Thin-walled elements are widely used in engineering structures of agriculture (agribusiness). They are structural elements of grain cleaning units, bunker structures for various purposes, as well as pipeline systems. Therefore, the development of modern numerical algorithms for strength calculations of the bunker type structure is an urgent task. In this article, an algorithm for calculating the strength of the thin-walled engineering structure of the bunker type as a structure consisting of two types of shells: cylindrical and conical. To determine the stress-strain state of such an aggregate, a numerical finite element method is used with the adoption of a curvilinear finite element developed by the authors in the form of a quadrangular fragment of the middle surface of the shell with nodal unknowns in the form of displacements and their first derivatives. The stiffness matrix of the used finite element is developed on the basis of the theory of thin shells with the use of the direct normal hypothesis. The kinematic parameters of the cylindrical shell are taken as the main unknown in the nodes at the junction of the cylindrical and conical shells. To obtain

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

the expressions of the nodal parameters of a conical shell using the basic nodal unknown of the cylindrical shell used the kinematic conditions and the equilibrium equations of statics. The obtained relations are implemented in the matrix form and used in the transformation of stiffness matrices and vectors of nodal loads of finite elements of the conical shell adjacent to the boundary of the joint with the cylindrical shell. On the example of calculating the strength of a thin-walled bunker-type structure consisting of a cylinder and a truncated cone, the efficiency of the developed algorithms is shown and the correctness of the obtained numerical values of strength indicators is confirmed. The analysis of the numerical results allowed us to conclude that the stress concentration coefficient in the zone of joint of the cylinder with the cone increases proportionally to the angle of inclination of the cone to its axis.

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

Key words: the thin-walled design of the hopper type, the final element, the conditions of the joint of the cylinder with the cone, the stress concentration.

Исследование выполнено при финансовой поддержке РФФИ и Администрации Волгоградской области в рамках научного проекта № 18-41-343003 р мол а

Цитирование. Клочков Ю.В., Николаев А.П., Фомин С.Д., Соболевская Т.А., Андреев А.С. Расчет прочности оболочечных конструкций АПК бункерного типа. Известия НВ АУК. 2019. 2(54). 285-297. DOI: 10.32786/2071-9485-2019-02-34.

Citation. Klochkov J. V., Nikolaev A. P., Fomin S. D., Sobolevskaya T. A., Andreev A. S. Calculation of strength of shell structures of AIC bunker type. Proc. of the Lower Volga Agro-University Comp. 2019. 2(54). 285-297. (in Russian). DOI: 10.32786/2071-9485-2019-02-34.

Введение. Неотъемлемой частью многих инженерно-технических объектов агропромышленного комплекса (АПК) являются тонкостенные конструкции бункерного типа, представляющие собой цилиндрические оболочки, сочлененные с конусом, такие как зерноочистительные агрегаты типа ЗАВ - 40, зерносушительные комплексы, сило-сы для хранения комбикормов, зерновых и масличных культур и другие. Расчёты на прочность и устойчивость таких объектов аналитическими методами достаточно трудоёмки и возможны лишь для сравнительно узкого класса задач. Поэтому в настоящее время применяют в основном численные методы расчёта [20, 21, 6, 5, 8, 3, 13], подразумевающие применение современной вычислительной техники и программного обеспечения. Наибольшую популярность среди численных методов сейчас имеет метод конечных элементов (МКЭ) в различных формулировках [1, 2, 7, 4, 9-12, 14-19, 22]. Современному инженеру предлагается достаточно широкий выбор вычислительных комплексов, многие из которых имеют иностранное происхождение. Однако большинство таких комплексов являются весьма дорогостоящими, а в сопроводительной документации ответственность за те или иные инженерные решения, принятые на основе результатов использования данных комплексов, целиком перекладывается на пользователя. Поэтому актуальной остаётся задача создания отечественных компактных программных модулей по конечно-элементному расчёту тонкостенных конструкций АПК бункерного типа в виде цилиндрической оболочки, сочленённой с усечённым конусом.

Материалы и методы. Применялись численные методы расчёта с использованием современной вычислительной техники и программного обеспечения. Использовался наиболее популярный среди численных методов - метод конечных элементов (МКЭ).

Результаты и обсуждение. Геометрия тонкостенной конструкции бункерного типа. Тонкостенная конструкция АПК бункерного типа представляет собой цилиндр, сочленённый с усеченным конусом. Срединная поверхность бункера описывается радиус-векторами:

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

- для цилиндрическом части:

Д0 = хТ + Язтб] + Исозвк (1)

- для конической части:

Я0 = (Ьс + х)Т + г(х)зтв] + г(х)созвк, (2)

где R - радиус цилиндра, 9 - угол, отсчитываемый от оси Оz против хода часовой стрелки в плоскости, перпендикулярной оси Ох; Lc - длина образующей цилиндра.

Входящая в (2) функция г(х) представляет собой радиус вращения конуса:

г(х) = Я — хЬда, (3)

где а - угол наклона образующей конуса к оси Ох.

Вычисляя частные производные от (1) и (2) по х и 9, можно найти векторы локального базиса для цилиндра:

а° = Я° = Т; а%= Я0в = Ясозв] — Изтвк; а0 = а° х а°/|а° х а°| = этв] + созвк (4)

и конуса

а1 = Т + гхзтв] + гхсозвк; а'® = г(х)зтв] + гхсозвк; а'0 = (—г(х) • гх1 + г(х)зтв] + г(х)созвк)/{г(х)^1 + г£), (5)

где гх = —1да.

Соотношения (4), (5) могут быть представлены в матричной форме: {а0} = [С0] {Т} ; {Т} =[С°]-1(а0};

3X1 ЗхЗ 3X1 3X1 ЗхЗ 3X1

{а0'} = [к0] {Т}; {Т} = [к°]~ча0'1 (6)

3X1 3X3 3X1 3X1 3X3 3X1

где {а0}Т = {а°1а°2а0 }; {Т} = {Т/7?}; {а0'}т = {а°'а°'а0'}.

1X3 1X3 1X3

На основании (6) можно установить матричную зависимость между векторами локальных базисов цилиндра и конуса:

{а0} = [С0]^0]"1^0'} = Ы0]{а0'}. (7)

3X1 3X3 3X3 3X1 3X3 3X1

Дифференцированием (4) и (5) можно найти производные векторов локального

базиса:

= 0; а°,х = 0; а° = 0; (8)

а°д = 0; а°,0 = —Язтб] — Ясозвк; а0в = созв] — 51П0А:; я?'* = 0; а°'х = гхсозв] — гхзтвк;

= ((—+ гхзтв] + гхсозвк) / (г(х)+

+(—г(х)гхТ + г(х)зтв] + г(х)со5 вк) • I--I; (9)

, у(г(х))2^1+г,2 у

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

= гвсозв] — гвзт6к; а.20 = — г(х)зтв] — г(х)созвк; а°д' = (г(х)сояв/ — г(х)зтвк)/(г(х)^1 + г£). Соотношения (8), (9) могут быть представлены в матричном виде:

{й°} = [0]®; ВД = В]Ш; {а0;} = [к°]т; ю = и]®,

или с учётом (6):

{Й°в} = [С°][С°]-1(а0} = {М}{а0}; {а«'} = [^[П"1^0'} = {у}{а0'};

{а°'} = [К°][К°]-Чй°} = Ы{й0}, (10)

где {а%} = {а1,ха2,ха,х}; {а°в} = {а1,еа2,еа°е}; {а,х'} = {а1,ха2,ха,х'}; = {а1,ва2,ва°в}-

При деформировании конструкции бункерного типа точка её срединной поверхности А0 переместится в точку А. Вектор перемещения данной точки и его производные по х и 0 записываются в локальном базисе точки А0:

х> = иа® + ра.® + wd0; х>х = + + £ха0;

v|в = tla% + t$a% + tвa0, (11)

где и, V, w - меридиональная, кольцевая и нормальная компоненты вектора перемещения, ¿х , ¿е - компоненты производных векторов перемещений.

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

V + = V + <а0,

Откуда

V = V + <(а — а0), (12)

где V - вектор перемещения точки, расположенной во внутреннем слое бункерной конструкции на расстоянии £ от её срединной поверхности; а - орт нормали в точке срединной поверхности в деформированном состоянии.

Ковариантные компоненты тензора деформаций могут быть определены соотношением механики сплошной среды (Седов Л. И. Механика сплошной среды, 1976):

= {дар — АЗА (13)

где даР = да • др; д^р = • 9р. (14)

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

= = (Д0 + (й°)а = а°а+ (а°а; (15)

да = = + У)а = д о + = д ° + + Ь2аа°2 + + ((а — а0),а.

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

Конечный элемент оболочечной конструкции бункерного типа. Конечным элементом бункерной конструкции выбирается четырёхугольный фрагмент её срединной поверхности с узлами т, п, р, s, расположенными в углах фрагмента. С целью реализации процедуры численного интегрирования по площади четырёхугольного фрагмента по квадратуре Гаусса в конечном элементе, наряду с глобальной х, 0 используется локальная система координат п, изменяющаяся в пределах от -1 до 1.

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

{Щ}Т = \{илу}Т{рЛ}Т{™Л}Т\; (16)

1X36 11X12 1X12 1X12 )

{^Г=(17)

1X36 11X12 1X12 1X12 ) где {гЛ}Т = {г^г^г?г^г^гр^}; = {г^г^г^г^г^^в}. Здесь

1X12 1X12

под z понимается компонента вектора перемещения и, V или w.

Между (16) и (17) можно установить матричную зависимость:

{Щ} = [Г] (18)

36X1 36х3636х1

Глобальные координаты х и 0 точки внутренней области четырехугольного конечного элемента выражаются через глобальные координаты узлов посредством матричной зависимости:

х = [ф}т{ху}; 9 = {^]Т{0У}, (19)

1X4 4x1 1X4 4X1

где {хуУ = {хтх"хрх5}; {ву}Т = {0т0"0р05}; {^>}Т = {^Ч^г^з^ч} - матрица-строка, содер-

1X4 1X4 1X4

жащая билинейные функции локальных координат:

_ (1-Р (1-4). _ (1+0 (1-4).

(20)

Для получения элементов входящей в (18) матрицы [Т] необходимо найти производные (19) по и п:

вл = {ф,ц}т{ву}; = к/М. (21)

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

г = (22)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1X1212X1

Выполняя операцию дифференцирования (22) по х и 0, можно получить интерполяционные выражения для производных компонент вектора перемещения первого и второго порядка:

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

г,а = +

= ({^„тЛЧа^ + {^¿пУ^.аЛ.р + Л.а^.ц) + {^.т^Л.аЛ.Ц + + {^лГ^Я2?}. (23)

Для компоновки матрицы жесткости и столбца узловых усилий четырёхузлового конечного элемента используется функционал Лагранжа:

ф = / И/ {^к - / , (24)

V Р

где (^д] = {^11^222^12}; = {а11а22а12} - деформации и напряжения в точке произ-

вольного слоя оболочки; [и]т = {и V ш) - компоненты вектора перемещения точки срединной поверхности оболочечной конструкции бункерного типа; {Р)т = {рг р2 Рз) - столбец внешней поверхностной нагрузки.

Дальнейшая процедура получения матрицы жесткости и столбца узловых усилий осуществляется по стандартному для МКЭ варианту в формулировке метода перемещений [7, 4, 9-12, 22, 17-19, 14-16].

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

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

~х>с = рк или

иа% + ра°2 + и/а° = и'а?' + р'а°2' + и/'а0'. (25)

Принимая во внимание (7), равенство (25) можно преобразовать к виду

и (й^а?' + + ^За0') + и (й^а?' + + ^За0') +

(й^а?' + + = и'а?' + р'а02' + и?'а0'. (26)

Из (26) можно получить следующие три равенства:

и' = ий.^ + vd21 + wdз1; V' = ийЧ2 + г^22 + wdз2;

иг' = ийЧ з + + . (27)

Применяя к (25) операцию дифференцирования по 0, можно записать следующие выражения:

и,^ = ^ или ^а? + + Ьв а0 = + + ^ва0'. (28)

Учитывая (7) из (28), можно получить ещё три равенства:

Ьд' = + £^21 + Ьд dз1; Ьд' = + £^22 + Ьд dз2;

г'д = + + 1д dз0з. (29)

Из (29) можно выразить производные компонент векторов перемещений конуса по координате в в узлах на границе пересечения через узловые неизвестные цилиндра:

290

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

и,

1',в = Л(и, и,х, и,в, V, и,х,

= /2 (и, и,х, и,в, V, V,х, ъ,в, ^ ™,х, (30)

= /з (и, и,х, и,0, и, V*, ш, ш,х, Жд).

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

рха0 = г^а0' или

+ + а° )а0 = + + ^а0')а0', откуда

= ^. (31)

Из (31) можно получить выражение для производной нормальной компоненты вектора перемещения по х конуса через узловые неизвестные цилиндрической части:

= Л(и, и,х, и,в, V, и,х, ™х, w|в). (32)

Статическое условие сочленения заключается в равенстве суммы сил в узлах на окружности сочленения:

с = 0 или

N¿0% + N¿0% +ЫХ а0 = -Л^'а?' - N/(2°' - Л^а0'. (33)

Из (33) можно записать равенство = — .

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

= /5 (и, Чв,™,в). (34)

Компонента и'х в узлах на окружности сочленения остаётся свободно варьируемой.

На основании (27), (30), (32) и (34) компонуются матрицы преобразований [PR], на которые умножаются матрицы жесткости и столбцы узловых усилий четырёхузло-вых конечных элементов усечённого конуса, которые граничат с окружностью сочленения.

Результаты и обсуждение. Была рассчитана тонкостенная бункерная конструкция, состоящая из цилиндра и усечённого конуса, нагруженная внутренним давлением q = 5,0 МПа. Торец цилиндрической части был шарнирно закреплён, торец усечённого конуса был свободен от закреплений (Рис. 1). Были приняты следующие исходные данные: длина образующей цилиндрической части Lc = 1,00 м; радиус цилиндра R = 0,9 м; толщина стенок бункера t = 0,02м; длина образующей конуса Lc = 0,6 м; угол а = - рад.

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

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

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

Координата х, м Напряжения о, МПа Число конечных элементов Аналитическое решение о, МПа

24 48 96

0,0 м < 80,89 81,0 81,0

< 81,17 81,13 81,13

сред °к 81,07 81,066 81,066 81,066

сред °к 225,05 225,04 225,04 225,00

1,0 м < 637,06 638,03 637,51

< -473,39 -475,85 -475,38

113,41 113,69 113,54

4Н -219,73 -220,47 -220,33

1,42 м < -0,13 -0,016 -0,001 0,000

< 0,14 0,016 0,001 0,000

165,82 165,85 165,86

152,33 152,30 152,29

сред °к 159,08 159,08 159,07

Рисунок 1 - Расчётная схема тонкостенной конструкции бункерного типа

На рисунке 2 представлены эпюры меридиональных напряжений на внутренней и наружной поверхностях бункерной конструкции. Как видно из рисунка, концентрация напряжений возникает в зоне сочленения цилиндрической части с усеченным конусом, что необходимо учитывать при выполнении прочностных расчётов бункерной конструкции.

Несомненный интерес представляет анализ влияния угла а на коэффициент концентрации напряжений в зоне сочленения цилиндрической части с усеченным конусом. На рисунке 3 представлен график изменения коэффициента концентрации напряжений К, равного отношению о^ при х = 0,0 м к при х = в зависимости от угла а.

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА: НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

183.47 161.1

Рисунок 2 - Эпюры нормальных напряжений °хх, МПа

30 35 40 45 50 55 60

Рисунок 3 - Зависимость коэффициента концентрации напряжений К от угла а

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

Как видно из рисунка 3, коэффициент концентрации напряжений К увеличивается с увеличением значения угла а.

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

Библиографический список

1. Игнатьев, А.В. Расчет тонких пластин по методу конечных элементов в форме классического смешанного метода с исключением перемещений конечных элементов как жесткого целого [Текст] / А. В. Игнатьев, В. А. Игнатьев, Е.А. Гамзатова // Известия высших учебных заведений. Строительство. - 2018. - № 3 (711). - С. 5-13.

2. Игнатьев, А. В. Анализ изгибаемых пластинок с односторонними связями по методу конечных элементов в форме классического смешанного метода / А. В. Игнатьев, В. А. Игнатьев, Е. А. Гамзатова // Известия высших учебных заведений. Строительство. - 2018. - № 8 (716). - С. 5-14.

3. Каюмов, Р.А. Потеря устойчивости листовых косоугольных консолей [Текст] / Р.А. Каюмов, И.З. Мухамедова, Г.Ф. Хазиева // Известия Казанского государственного архитектурно-строительного университета. - 2018. - № 3 (45). - С. 129-135.

4. Каюмов, Р. А. Моделирование процесса деформирования и оценка несущей способности системы грунт - тонкостенная конструкция [Текст]/ Р. А. Каюмов, Ф. Р. Шакирзянов, С.С. Гаврю-шин // Известия высших учебных заведений. Машиностроение. - 2014. - № 6. - С. 20-24.

5. Ким, А.Ю. Сравнение экспериментального и численного исследования большепролетного пневматического линзообразного сооружения [Текст] / А.Ю. Ким, С.В. Полников // Научное обозрение. - 2016. - № 15. - С. 36-41.

6. Пятикрестовский, К.П. О программировании нелинейного метода расчета деревянных конструкций [Текст] / К.П. Пятикрестовский, В.И. Травуш// Academia. Архитектура и строительство. - 2015. - № 2. - С. 115-119.

7. Расчетное обоснование механической безопасности стадионов к чемпионату мира по футболу 2018 года [Текст] / А. М. Белостоцкий, П. А. Акимов, А. А. Аул, Д. С. Дмитриев, Ю. Н. Дядченко, А. И. Нагибович, К. И. Островский, А. С. Павлов // Academia. Архитектура и строительство. - 2018. - № 3 - С. 118-129.

8. Хайруллин, Ф.С. Расчет тонких оболочек с использованием аппроксимирующих функций различного порядка [Текст] / Ф.С. Хайруллин, Д.Д. Мингалиев // Вестник Технологического университета. - 2017. - Т. 20. - № 14. - С. 102-104.

9. Чепурненко, А.С. Расчет цилиндрических гофрированных конструкций при помощи полуаналитического метода конечных элементов [Текст]/ А.С. Чепурненко, Б.М. Языев, М.С. Турко // Строительство и техногенная безопасность. - 2018. - № 12 (64). - С. 49-56.

10. Belostotsky, А.М. Adaptive finite-element models in structural health monitoring systems [Текст]/ А.М. Belostotsky, P.A. Akimov, O.A. Negrozov, N.O. Petryashev, S.O. Petryashev, S.V. Sherbina, D.K. Kalichava, T.B. Kaytukov // Инженерно-строительный журнал. - 2018. - № 2. - С. 169-178.

11. Continuum mechanics based bilinear shear deformable shell element using absolute nodal coordinate formulation, Trans. ASME. [Tekst] / H. Yamashita, Antti I. Valkeapaa, P. Jayakumar, H. Syqiyama // J. Comput. and Nonlinear Dyn. - 2015. - Vol. 10. - № 5. - P. 051012/1-051012/9.

12. Finite element simulatio of metallic cylindrical sandwich shells with graded aluminum tabular cores subjected to internal blast loading, [Tekst] / Shiqiang Li, Guoxing Lu, Zhihua Wang, Longmao Zhao, Guiying Wu // Int. J. Mech. Sci. - 2015. - № 96-97. - P. 1-12.

13. Galishnikova, V.V. Analysis of frame buckling without sidesway classification [Tekst] / V.V. Galishnikova, P.Ja. Pahl // Строительная механика инженерных конструкций и сооружений. 2018. - Т. 14. - № 4. - С. 299-312.

14. Hanslo, P. Tangential differential calculus and the finite element modeling of a large deformation elastic membrane problem [Tekst] / P. Hanslo, G. Larson Mats, F. Larson // Comput. Mech. - 2015. - Vol. 56. - № 1. - P. 87-95.

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

15. He, Xiaocong. Finite element analysis of torsional free vibration of adhesively bonded single-lap joints [Tekst] / He, Xiaocong // Int. J. Adhes.and Adnes. - 2014. - № 48. - P. 59-66.

16. Levyakov, S.V. Application of triangular element invariants for geometrically nonlinear analysis of functionally graded shells [Tekst] / S.V. Levyakov, V.V. Kuznetsov // Computational Mechanics. - 2011. - Vol. 48. - №. 4. - P. 499-513.

17. Nguyen, Nhung. Nonlinear, finite deformation, finite element analysise, ZAMP. [Tekst] / Nhung Nguyen, Anthonym Waas // Z. Angew. math. and Phys. - 2016. - Vol. 67. - №. 9. - P. 35/1-35/24.

18. Sartorato, Murilo. A finite element formulation for smart piezollectric composite shells: Mathematical formulation, computational analysis and experimental evaluation [Tekst] / Murilo Sartorato, Ricardo de Medeiros, Volnei Tita // Comros. Struct. - 2015. - №. 127. - P. 185-198.

19. Solodovnikov, A.S. Numerical study of strength properties for a composite material with short reinforcing fibers [Tekst] / A.S. Solodovnikov, S.V.Sheshenin // Moscow University Mechanics Bulletin. - 2017. - T. 72. - № 4. - P. 94-100.

20. Storozhuk, E.A. Stress Distribution Near a Circular Hole in a Flexible Orthotropic Cylindrical Shell of Elliptical Cross-Section [Tekst] / E.A. Storozhuk, S.M. Komarchuk // International Applied Mechanics. - 2018. - № 54(6). - C. 687-694.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

21. Storozhuk, E.A. Stress-Strain State Near a Hole in a Shear-Compliant Composite Cylindrical Shell with Elliptical Cross-Section [Tekst] / E.A. Storozhuk, I.S. Chernyshenko, A.V. Yatsura // International Applied Mechanics. - 2018. - № 54(5). - C. 559-567.

22. Zheleznov, L.P. Nonlinear deformation and stability of discretely reinforced elliptical cylindrical shells under transverse bending and internal pressure [Tekst] / L.P. Zheleznov, V.V. Kabanov, D.V. Boiko // Russian Aeronautics. - 2014. - T. 57. - № 2. - P. 118-126.

References

1. Ignat'ev, A. V. Raschet tonkih plastin po metodu konechnyh jelementov v forme klas-sicheskogo smeshannogo metoda s isklyucheniem peremeschenij konechnyh jelementov kak zhestkogo celogo [Tekst] / A. V. Ignat'ev, V. A. Ignat'ev, E. A. Gamzatova // Izvestiya vysshih uchebnyh zavedenij. Stroitel'stvo. - 2018. - № 3 (711). - P. 5-13.

2. Ignat'ev, A. V. Analiz izgibaemyh plastinok s odnostoronnimi svyazyami po metodu konechnyh ]lementov v forme klassicheskogo smeshannogo metoda / A. V. Ignat'ev, V. A. Ignat'ev, E. A. Gamzatova // Izvestiya vysshih uchebnyh zavedenij. Stroitel'stvo. - 2018. - № 8 (716). - P. 5-14.

3. Kayumov, R. A. Poterya ustojchivosti listovyh kosougol'nyh konsolej [Tekst] / R. A. Kayumov, I. Z. Muhamedova, G. F. Hazieva // Izvestiya Kazanskogo gosudarstvennogo arhitekturno-stroitel'nogo universiteta. - 2018. - № 3 (45). - P. 129-135.

4. Kayumov, R. A. Modelirovanie processa deformirovaniya i ocenka nesuschej sposobnosti sistemy grunt - tonkostennaya konstrukciya [Tekst]/ R. A. Kayumov, F. R. Shakirzyanov, S. S. Gavryushin // Izvestiya vysshih uchebnyh zavedenij. Mashinostroenie. - 2014. - № 6. - P. 20-24.

5. Kim, A. Yu. Sravnenie jeksperimental'nogo i chislennogo issledovaniya bol'sheproletnogo pnevmaticheskogo linzoobraznogo sooruzheniya [Tekst] / A. Yu. Kim, S. V. Polnikov // Nauchnoe obozrenie. - 2016. - № 15. - P. 36-41.

6. Pyatikrestovskij, K. P. O programmirovanii nelinejnogo metoda rascheta derevyannyh kon-strukcij [Tekst] / K. P. Pyatikrestovskij, V. I. Travush// Academia. Arhitektura i stroitel'stvo. - 2015. -№ 2.- P. 115-119.

7. Raschetnoe obosnovanie mehanicheskoj bezopasnosti stadionov k chempionatu mira po futbolu 2018 goda [Tekst] / A. M. Belostockij, P. A. Akimov, A. A. Aul, D. S. Dmitriev, Yu. N. Dy-adchenko, A. I. Nagibovich, K. I. Ostrovskij, A. S. Pavlov // Academia. Arhitektura i stroitel'stvo. -2018. - № 3 - P. 118-129.

8. Hajrullin, F. S. Raschet tonkih obolochek s ispol'zovaniem approksimiruyuschih funkcij razlichnogo poryadka [Tekst] / F. S. Hajrullin, D. D. Mingaliev // Vestnik Tehnologicheskogo univer-siteta. - 2017. - T. 20. - № 14. - P. 102-104.

9. Chepurnenko, A. S. Raschet cilindricheskih gofrirovannyh konstrukcij pri pomoschi polu-analiticheskogo metoda konechnyh jelementov [Tekst]/ A. S. Chepurnenko, B. M. Yazyev, M. S. Turko // Ctroitel'stvo i tehnogennaya bezopasnost'. - 2018. - № 12 (64). - P. 49-56.

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

10. Belostotsky, А.М. Adaptive finite-element models in structural health monitoring systems [Текст]/ А.М. Belostotsky, P.A. Akimov, O.A. Negrozov, N.O. Petryashev, S.O. Petryashev, S.V. Sherbina, D.K. Kalichava, T.B. Kaytukov // Инженерно-строительный журнал. - 2018. - № 2. - P. 169-178.

11. Continuum mechanics based bilinear shear deformable shell element using absolute nodal coordinate formulation, Trans. ASME. [Tekst] / H. Yamashita, Antti I. Valkeapaa, P. Jayakumar, H. Syqiyama // J. Comput. and Nonlinear Dyn. - 2015. - Vol. 10. - № 5. - P. 051012/1-051012/9.

12. Finite element simulatio of metallic cylindrical sandwich shells with graded aluminum tabular cores subjected to internal blast loading, [Tekst] / Shiqiang Li, Guoxing Lu, Zhihua Wang, Longmao Zhao, Guiying Wu // Int. J. Mech. Sci. - 2015. - № 96-97. - P. 1-12.

13. Galishnikova, V.V. Analysis of frame buckling without sidesway classification [Tekst] / V.V. Galishnikova, P.Ja. Pahl // Строительная механика инженерных конструкций и сооружений. 2018. - Т. 14. - № 4. - P. 299-312.

14. Hanslo, P. Tangential differential calculus and the finite element modeling of a large deformation elastic membrane problem [Tekst] / P. Hanslo, G. Larson Mats, F. Larson // Comput. Mech. -2015. - Vol. 56. - № 1. - P. 87-95.

15. He, Xiaocong. Finite element analysis of torsional free vibration of adhesively bonded single-lap joints [Tekst] / He, Xiaocong // Int. J. Adhes.and Adnes. - 2014. - № 48. - P. 59-66.

16. Levyakov, S.V. Application of triangular element invariants for geometrically nonlinear analysis of functionally graded shells [Tekst] / S.V. Levyakov, V.V. Kuznetsov // Computational Mechanics. - 2011. - Vol. 48. - №. 4. - P. 499-513.

17. Nguyen, Nhung. Nonlinear, finite deformation, finite element analysise, ZAMP. [Tekst] / Nhung Nguyen, Anthonym Waas // Z. Angew. math. and Phys. - 2016. - Vol. 67. - №. 9. - P. 35/1-35/24.

18. Sartorato, Murilo. A finite element formulation for smart piezollectric composite shells: Mathematical formulation, computational analysis and experimental evaluation [Tekst] / Murilo Sartorato, Ricardo de Medeiros, Volnei Tita // Comros. Struct. - 2015. - №. 127. - P. 185-198.

19. Solodovnikov, A.S. Numerical study of strength properties for a composite material with short reinforcing fibers [Tekst] / A.S. Solodovnikov, S.V.Sheshenin // Moscow University Mechanics Bulletin. - 2017. - Т. 72. - № 4. - P. 94-100.

20. Storozhuk, E.A. Stress Distribution Near a Circular Hole in a Flexible Orthotropic Cylindrical Shell of Elliptical Cross-Section [Tekst] / E.A. Storozhuk, S.M. Komarchuk // International Applied Mechanics. - 2018. - № 54(6). - P. 687-694.

21. Storozhuk, E.A. Stress-Strain State Near a Hole in a Shear-Compliant Composite Cylindrical Shell with Elliptical Cross-Section [Tekst] / E.A. Storozhuk, I.S. Chernyshenko, A.V. Yatsura // International Applied Mechanics. - 2018. - № 54(5). - P. 559-567.

22. Zheleznov, L.P. Nonlinear deformation and stability of discretely reinforced elliptical cylindrical shells under transverse bending and internal pressure [Tekst] / L.P. Zheleznov, V.V. Kabanov, D.V. Boiko // Russian Aeronautics. - 2014. - Т. 57. - № 2. - P. 118-126.

Информация об авторах Клочков Юрий Васильевич, заведующий кафедры «Высшая математика», Волгоградского государственного аграрного университета (РФ, 400002, г. Волгоград, Университетский проспект, 26), доктор технических наук, профессор.

ORCID: https://orcid.org/0000-0002-1027-1811. E-mail: klotchkov@bk.ru.

Николаев Анатолий Петрович, профессор кафедры «Прикладная геодезия, природообустрой-ство и водопользование», Волгоградского государственного аграрного университета (РФ, 400002, г. Волгоград, Университетский проспект, 26), доктор технических наук, профессор. ORCID: https://orcid.org/0000-0002-7098-5998. E-mail:anpetr40@yandex.ru.

Фомин Сергей Денисович, доцент кафедры «Механика», заведующий Центром наукометрического анализа и международных систем индексирования Волгоградского государственного аграрного университета (РФ, 400002, г. Волгоград, Университетский проспект, 26), доктор технических наук, доцент.

ORCID: https://orcid.org/0000-0001-7910-9284. E-mail: fsd_58@mail.ru.

НИЖНЕВОЛЖСКОГО АГРОУНИВЕРСИТЕТСКОГО КОМПЛЕКСА НАУКА И ВЫСШЕЕ ПРОФЕССИОНАЛЬНОЕ ОБРАЗОВАНИЕ

Соболевская Татьяна Алексеевна, доцент кафедры «Высшая математика», Волгоградского государственного аграрного университета (РФ, 400002, г. Волгоград, Университетский проспект, 26), кандидат технических наук.

ORCID: https://orcid.org/0000-0002-9167-075X. E-mail: moonway13@rambler.ru.

Андреев Александр Сергеевич, ассистент кафедры «Высшая математика», Волгоградского государственного аграрного университета (РФ, 400002, г. Волгоград, Университетский проспект, 26). ORCID: https://orcid.org/0000-0002-3763-0394. E-mail: aandreev.07.1988@gmail.com

Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов. Conflict of interest. The authors declare no conflict of interest.

УДК 519.245.53.08 DOI: 10.32786/2071-9485-2019-02-35

СИЛОВЫЕ И КИНЕМАТИЧЕСКИЕ ЗАВИСИМОСТИ РУЛЕВОГО УСИЛИТЕЛЯ

POWER AND KINEMATIC DEPENDENCE OF THE STEERING AMPLIFIER

Г.Д. Кокорев, доктор технических наук, доцент И.А.Успенский, доктор технических наук, профессор И.А. Юхин, доктор технических наук, доцент И.А. Афиногенов, аспирант

G.D. Kokorev, I.A. Uspensky, I.A. Yukhin, I.A. Afinogenov

ФГБОУ ВО «Рязанский государственный агротехнологический университет

имени П.А. Костычева»

Federal State Budget Educational Institution of Higher Education Ryazan State Agrarian University named after P.A. Kostycheva

Дата поступления в редакцию 04.03.2019 Дата принятия к печати 12.05.2019

Received 00.00.2019 Submitted 12.05.2019

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

The article discusses the method of assessing the power and kinematic dependencies of the steering amplifier using the mathematical models proposed by the author, using the example of practical calculation of the hydraulic steering amplifier of a Ural 4320 vehicle. ) and force tracking action (by the force of resistance to rotation). A theoretical and practical analysis of the proposed method and

i Надоели баннеры? Вы всегда можете отключить рекламу.