Technique solution of the problem of the pipeline precast structure frequency response search with application of the finite element method

Full Text

Представленная в рамках данной статьи работа проводится при финансовой поддержке Правительства Российской Федерации (Минобрнауки России) по комплексному проекту 2012-218-03-120 «Автоматизация и повышение эффективности процессов изготовления и подготовки производства изделий авиатехники нового поколения на базе Научно-производственной корпорации «Иркут» с научным сопровождением Иркутского государственного технического университета» согласно постановлению Правительства Российской Федерации от 9 апреля 2010 г. № 218. 16 Математика, механика, информатика has been done. Frequency responses of collapsible pipelines constructions hydraulic system pressure line of a modern airplane before and after modernization on the presented method are presented; decrease of structures dynamic stresses in the operating frequency range. This work continues research of dependence of frequencies of natural vibration and tensely deformed condition ofpipeline systems from features of its detail of construction. Keywords: the pipeline, frequency response, hydraulic system, resonance, damping, dynamic tensely deformed condition, contact problem, finite element method. Изделия авиационной техники сегодня представляют собой сложнейший комплекс систем и оборудования, интегрированный в планер. В полете все элементы современного самолёта должны работать как одно целое, обеспечивая соответствие заявленным при его проектировании характеристикам. Однако учёт совместной работы элементов систем и планера при различных вариантах совместного нагружения является сложнейшей задачей, которая решается с применением итерационных подходов, с неоднократными доработками и модификациями тех или иных элементов и согласованием целого ряда параметров. Трубопровод - сложнонагруженная конструкция, исследование работы которой - комплексная задача, существующая на стыке исследования экспериментальных данных натурных испытаний, результатов традиционных численных решений и высокопроизводительных численных решений (рис. 1). Традиционные методы, описанные в [1-5], являются громоздкими и сложноприменимыми в производстве. Аналитические же методы в большей степени позволяют рассматривать отдельные элементы сборных конструкций и представляют известную сложность, связанную с определением всех возмож ных типов нагрузок и действующих напряжений, что не позволяет рассматривать масштабные сборные конструкции, где велико взаимовлияние элементов друг на друга, особенно при динамическом нагружении. Экспериментальные же исследования представляются слишком дорогостоящими ввиду необходимости наличия соответствующего оборудования и программных средств для обработки потока данных эксперимента (например, при тензометрировании динамических испытаний) и используются в основном для подтверждения заявленных характеристик и при отработках конструкций. Численные же методы решения подобных задач представлены системами инженерного анализа, основанными на методе конечных элементов. Методы численного решения задачи определения динамического напряженно-деформированного состояния трубопроводных систем описаны в серии работ [6-15]. Основные тенденции развития методики [6] заключаются в более глубокой проработке и уточнении средств рассмотрения экспериментальных данных натурных испытаний, а также в более точном определении параметров надёжности системы через сравнение частот отклика и возбуждения сборных конструкций. Рис. 1. Модель трубопроводной системы 17 Вестник СибГАУ. № 2(54). 2014 Согласно методике [8], необходимо рассмотреть задачу определения частотного отклика конструкции трубопроводной системы для определения критических частот работы трубопровода. По результатам работы необходимо сравнить критические частоты с рабочими частотами системы и устранить совпадающие частоты путём изменения конфигурации модели трубопровода или точек его крепления к элементам каркаса планера. Однако обнаружено, что при анализе конечно-элементной модели сборной конструкции трубопровода по заданному диапазону частот (анализ Transient static) в системе инженерного анализа ANSYS возможно разрушение контактов модели, что обусловлено нелинейным увеличением энергии колебаний постоянной амплитуды с увеличением их частоты. Таким образом, необходимо определить функцию колебаний точек крепления трубопроводов, позволяющую исследовать вибрации трубопровода с постоянной энергией и переменной частотой. Выполнен поиск функции перемещений в виде f (t) = Â(t) x sinra(t); (1) где t - время, с; A (t) - закон изменения амплитуды колебаний по времени; œ(t) - закон изменения частоты действия внешней силы. При этом амплитуда функции (1) уменьшается со временем, а частота возрастает. Модель сборной конструкции, подвергаемая динамическому анализу, может быть рассчитана при условии сохранения целостности контактов во время действия любых возникающих в процессе приложения исходной нагрузки сил. Опытным путём установлено, что максимальные амплитуды колебаний системы не вызывают нарушения контактных связей при амплитуде колебаний, равной 2 мм. Эффекта разрушения контактов в модели сборной конструкции можно избежать путём применения закона постоянства энергии колебательной системы при автоколебаниях: m ха2 x A2 2 - = const, где m - масса системы, кг; œ - частота колебаний системы; A - амплитуда колебаний системы. В любой момент времени t m xffljj x A02 = m xro2(t) x A2(t) 2 = 2 . (2) Таким образом, амплитуда зависит от частотного диапазона исследования: Л xa0 A(t ) = ra(t ) (3) где œ0 - начальное значение частоты (условно принятое равным 1); А0 - условное начальное значение амплитуды колебаний конструкции, не приводящее к разрушению контактов в модели. Частота по условиям задачи поиска частотного отклика зависит от времени по линейному закону ю = k • x + b . (4) Раскрывая уравнение (3) с учетом (4), подставляя граничные условия, получаем K (5) где K - коэффициент пропорциональности, зависящий от принятых начальных значений амплитуды, частоты и шага по времени. Для рассматриваемого расчётного случая значение К = 0,109 соответствует оптимальным с точки зрения задания переменных нагрузок в конечно-элементной модели трубопровода значениям начального времени t0 и начального значения амплитуды œ(t0) (исходные параметры для этих значений: t0 = 0,1 с, A0 = 2 мм). Из условия известных значений начальной и конечной частот вращения авиадвигателя как источника возбуждения колебаний n(t1) ^ 0 и n(t2) = 700 об/мин, получено: ro(t) = 700 x t x 2 xn . (6) 60 Согласно приведенной методике, функция (1) с учетом вышеуказанных условий имеет вид (рис. 2) 0,109 x sin(73,303 x t) f (t ) =t (7) Рис. 2. График функции f(t) с постоянной энергией колебаний в системе при увеличении частоты с заданными параметрами частотно-временных характеристик Функция (7) позволяет задавать перемещения точек опор трубопровода, анализировать сборные конструкции при постоянном значении энергии систем без разрушения её контактных взаимодействий и применима для инженерных анализов в стандартных КЭ-пакетах. Частотный отклик, полученный в результате анализа модели по представленной методике, отображен на рис. 3. Рис. 3. Частотный отклик модели трубопровода до модернизации 18 Математика, механика, информатика Таким образом, однозначно определён первый резонансный пик всей конструкции при частоте 125,44 Гц, что находится в диапазоне рабочих частот. Следовательно, модернизацию необходимо выполнить таким образом, чтобы исключить возможность появления резонансных пиков в диапазоне рабочих частот. Выполненный анализ частотных откликов модели с модернизированными кронштейнами подтвердил способность конструкции эффективно демпфировать динамические колебания. Отсутствие ярко выраженных пиков динамических напряжений и снижение общего фона напряжений при колебаниях системы с частотой, при которой в исходной конструкции прослеживались ярко выраженные пики напряжений, наглядно представлены на рис. 4. Рис. 4. Частотный отклик модели после модернизации Таким образом, уточнена методика определения частотного отклика протяженных сборных конструкций. Применение описанной методики позволит избежать итеративного поиска зависимости возбуждающей силы от времени и впоследствии эффективнее определять параметры любых решателей в стандартных КЭ-пакетах, применяемых для динамического анализа протяженных сборных конструкций с применением нелинейной контактной задачи теории конечных элементов.

About the authors

Igor Olegovich Bobarika

Irkutsk State Technical University

Email: MegusMC@mail.ru
Candidate of Engineering Science, associate professor of the Department of Aircraft construction and operations

Mikhail Sergeevich Yakhnenko

Irkutsk aircraft factory - the subsidiary of JSC “Corporation “Irkut”

Email: holtfor@mail.ru
Candidate of Engineering sciences, design engineer of the Department of detail designing

References

Supplementary files