RU2265193C1 - Способ определения резонансной частоты и добротности колебаний объекта - Google Patents
Способ определения резонансной частоты и добротности колебаний объекта Download PDFInfo
- Publication number
- RU2265193C1 RU2265193C1 RU2004106169/28A RU2004106169A RU2265193C1 RU 2265193 C1 RU2265193 C1 RU 2265193C1 RU 2004106169/28 A RU2004106169/28 A RU 2004106169/28A RU 2004106169 A RU2004106169 A RU 2004106169A RU 2265193 C1 RU2265193 C1 RU 2265193C1
- Authority
- RU
- Russia
- Prior art keywords
- frequency
- parameter
- oscillations
- time
- function
- Prior art date
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Apparatuses For Generation Of Mechanical Vibrations (AREA)
- General Electrical Machinery Utilizing Piezoelectricity, Electrostriction Or Magnetostriction (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Изобретение относится к машиностроению, а именно к способам определения резонансной частоты и добротности колебаний объекта. Предлагаемый способ позволяет определять резонансную частоту и добротность колебаний объекта в процессе эксплуатации без влияния на режим работы объекта для проведения измерений, без остановок и задержек на резонансных частотах, что уменьшит вероятность поломок объекта в процессе эксплуатации. Предлагаемый способ также позволяет определять резонансную частоту и добротность колебаний объекта в процессе его запуска и остановки. Операции предлагаемого способа могут быть легко автоматизированы для проведения измерений, что во многих случаях позволит уменьшить время и квалификацию исследователя для определения необходимых параметров. Способ определения резонансной частоты и добротности колебаний объекта заключается в том, что для исследуемого объекта выявляют резонансную область частот (fн; fв), в процессе работы, когда объект находится под возбуждающим воздействием, измеряют и регистрируют параметр колебаний объекта в процессе работы как функцию времени x(t) при условии нахождения частоты возбуждающего воздействия в резонансной области частот [fн; fв] и по возможности близкому к постоянному уровню возбуждающем воздействии. По зарегистрированным значениям параметра колебаний объекта x(t) как функции времени t судят о резонансной частоте fр и добротности Q. 1 з.п. ф-лы, 1 ил.
Description
Изобретение относится к машиностроению, а именно к способам определения резонансной частоты и добротности колебаний объектов.
Известен способ [А.С. СССР №381020, G 01 H 15/00, Бюл. №21, 1973 г.] определения декремента затухания, по которому возбуждают резонансные колебания исследуемого объекта, затем возбуждают колебания на нерезонансной частоте при воздействии на объект увеличенной вынуждающей силы, обеспечивающей получение амплитуды колебаний, равной амплитуде резонансных колебаний, логарифмический декремент определяют по формуле
где ω - нерезонансная частота колебаний, Р - резонансная частота колебаний, В - постоянная, равная относительному повышению вынуждающей силы.
Недостатком известного способа является невозможность определения резонансной частоты и добротности в процессе эксплуатации без влияния на режим работы объекта для проведения измерений, без остановок и задержек на резонансных частотах, что увеличивает вероятность поломок объекта, а также большое время реализации способа. Недостатком известного способа является также необходимость проведения измерений на резонансных частотах в установившемся режиме, что может привести к разрушению объекта при превышении возбуждающего воздействия допустимого уровня.
Наиболее близким к предлагаемому способу является способ, приведенный патент РФ №2086943, G 01 M 7/02, Бюл. №22, 1997 г., в котором выявляют резонансную область частот (fн; fв), устанавливают возбуждающее воздействие фиксированного уровня на нижней границе частоты резонансной области fн, увеличивают частоту F возбуждающего воздействия со скоростью vf1, в процессе увеличения частоты F регистрируют амплитуды колебаний объекта как функцию частоты возбуждения f, определяют частоту fнв как частоту возбуждающего воздействия, выше которой наблюдается уменьшение амплитуды колебаний, прекращают увеличение частоты возбуждающего воздействия f; затем уменьшают частоту возбуждающего воздействия f со скоростью vf2, в процессе уменьшения частоты f регистрируют амплитуды колебаний объекта как функцию частоты возбуждения f и определяют частоту fвн как частоту возбуждающего воздействия, ниже которой наблюдается уменьшение амплитуды колебаний, прекращают уменьшение частоты f, повторяют операции 3-10 до нахождения резонансной частоты fp с требуемой точностью, регистрируют амплитуду резонансных колебаний Ар объекта на частоте fp, затем производят расстройку резонанса путем изменения частоты возбуждающего воздействия до некоторой случайно выбранной величины fp+Δf, регистрируют амплитуду колебаний Δд объекта на частоте fp+Δf, по совокупности значений Ар, Ад, fp, fp+Δf судят о добротности колебаний Q.
Недостатком известного способа является большая сложность нахождения резонансной частоты fp с требуемой точностью, другим недостатком является сложность поддержания постоянной амплитуды возбуждающего воздействия. При случайной расстройке частоты f в области малых Δf погрешность определения Δf существенно влияет на точность определения логарифмического декремента затухания.
Общим недостатком известных способов является то, что необходимость проведения измерений на резонансных частотах в установившемся режиме может привести к разрушению объекта при превышении возбуждающего воздействия допустимого уровня. Также существенным недостатком известных способов является то, что они не позволяют определять резонансную частоту и добротность колебаний объектов, для которых недопустимы резонансы на стационарных режимах. Известные способы не могут определять резонансную частоту и добротность колебаний объекта в динамических режимах.
Поставлена задача расширения возможностей использования способа, а именно определения резонансной частоты и добротности в процессе эксплуатации без влияния на режим работы объекта для проведения измерений, без остановок и задержек на резонансных частотах, что уменьшит вероятность поломок объекта, а также уменьшит время, необходимое для определения резонансной частоты и добротности колебаний.
Поставленная задача достигается за счет того, что в способе определения резонансной частоты и добротности колебаний выявляют резонансную область частот (fн; fв), в процессе работы, когда объект находится под возбуждающим воздействием, измеряют и регистрируют параметр колебаний объекта в процессе работы как функцию времени x(t) при условии нахождения частоты возбуждающего воздействия в резонансной области частот [fн; fв] (в качестве регистрируемого параметра колебаний объекта x(t) может использоваться функция колебательного перемещения объекта от времени, или функция скорости перемещения объекта от времени, или функция ускорения перемещения объекта от времени), по зарегистрированным значениям параметра колебаний объекта x(t) как функции времени t судят о резонансной частоте fр и добротности Q.
О величине резонансной частоты fр и добротности Q судят следующим образом: из зарегистрированных значений параметра колебаний объекта в процессе работы как функции времени x(t) выбирают интервал времени с монотонным изменением частоты возбуждающего воздействия внутри всей резонансной области частот [fн; fв], для выбранного интервала времени по зарегистрированным значениям параметра колебаний объекта x(t) определяют спектральную плотность параметра колебаний объекта x(f), по спектральной плотности параметра колебаний объекта x(f) судят о резонансной частоте fр и добротности Q колебательной системы объекта.
Сущность способа поясняется схемой, представленной на чертеже, где показаны в резонансной области частот зависимости спектральной плотности параметра колебаний объекта x(f) (сплошная линия) и параметра колебаний объекта x(t) как функция частоты f (ось времени пересчитана в соответствующие значения частоты возбуждающего воздействия F(t); зависимость параметра колебаний объекта x(t) от частоты f показана штриховой линией и для наглядности увеличена в 5 раз).
Определение резонансной частоты и добротности колебаний предлагаемым способом осуществляется следующим образом: выявляют резонансную область частот (fн; fв), в процессе работы, когда объект находится под возбуждающим воздействием, измеряют и регистрируют параметр колебаний объекта в процессе работы как функцию времени x(t) при условии нахождения частоты возбуждающего воздействия в резонансной области частот [fн; fв], по зарегистрированным значениям параметра колебаний объекта x(t) как функции времени t судят о резонансной частоте fр и добротности Q.
О величине резонансной частоты fр и добротности Q судят следующим образом: из зарегистрированных значений параметра колебаний объекта в процессе работы как функции времени x(t) выбирают интервал времени с монотонным изменением частоты возбуждающего воздействия внутри всей резонансной области частот [fн; fв], для выбранного интервала времени по зарегистрированным значениям параметра колебаний объекта x(t) определяют спектральную плотность параметра колебаний объекта x(f), о резонансной частоте fр колебательной системы объекта судят, например, по положению максимума спектральной плотности параметра колебаний объекта x(f), а о добротности Q колебательной системы объекта судят, например, по отношению резонансной частоты fр колебательной системы объекта к ширине зависимости спектральной плотности параметра колебаний объекта от частоты x(f) на общепринятом для определения полосы пропускания колебательных систем уровне .
Для пояснения способа определения резонансной частоты fр и добротности Q сделаем предположение, что нам известна частота возбуждающего воздействия F(t) как функция времени t. Тогда, по значениям частоты входного воздействия F(t) во всей резонансной области [fн; fв], можно моделировать воздействие на объект в виде какой либо функции той же частоты, например в виде гармонического колебания e(t)=e0·cos[2π·F(t)·t+φ].
Для модели воздействия e(t) и выделенного интервала времени [tн; tк] пребывания частоты F возбуждающего воздействия во всей резонансной области [fн; fв] находят спектральную плотность e(f) модели воздействия e(t) по формуле прямого преобразования Фурье
Для того же самого выделенного интервала времени [tн; tк] и зарегистрированных как функцию времени значений параметра колебаний объекта x(t) также по формуле прямого преобразования Фурье находят спектральную плотность x(f)
Реакция x(f) линейной системы в спектральной области связана с воздействием E(f) [Гоноровский И.С. Радиотехнические цепи и сигналы. Изд. 2-е. М.: Советское радио, 1971. - 672 с., ил.] соотношением
где G(f) - частотная характеристика колебательной системы объекта, E(f) - спектральная плотность возбуждающего воздействия.
Если модель e(t) возбуждающего воздействия E(t) выбрана так, что спектральная плотность e(f) модели для резонансной области частот [fн; fв] пропорциональна спектральной плотности возбуждающего воздействия E(f)
e(f)=K·E(f),
где К - коэффициент пропорциональности,
то выражение (3) будет иметь вид
x(f)=Gв(f)·E(f),
где Gв(f) - вспомогательная функция, связанная с G(f)- частотной характеристикой колебательной системы объекта соотношением
Cв(f)=G(f)/K.
Колебательные свойства объекта могут быть определены по его частотной характеристике G(f), которую в свою очередь можно определить из выражения (3)
Если резонансная область [fн; fв] выбрана так, что колебательная система объекта содержит только одну резонансную частоту, то амплитудно-частотную характеристику колебательной системы объекта в резонансной области [Гоноровский И.С. Радиотехнические цепи и сигналы. Изд. 2-е. М.: Советское радио, 1971. - 672 с., ил.] можно описать выражением вида
где а - обобщенная расстройка, характеризующая основные свойства колебательной системы объекта
здесь Δf=f-fp - абсолютная расстройка; fp - резонансная частота колебательной системы объекта.
Для определения параметров колебательной системы объекта: добротности Q и резонансной частоты fp, можно аппроксимировать модуль частотной характеристики |G(f)| модулем выражения (5), т.е.
Максимальное значение функции |Ga(fp)| достигается на частоте f=fp, поэтому частоту, на которой достигается наибольшее значение модуля частотной характеристики |G(f)|, можно считать равной резонансной частоте колебательной системы объекта.
Добротность Q колебательной системы объекта можно определить с помощью выражения (7). Для этого находят частоту f1, для которой значения |Ga(f1)| уменьшаются по сравнению с резонансной в K1 раз, и частоту f2, для которой значения |Ga(f2)| уменьшаются по сравнению с резонансной в К2 раз, т.е.
Из этих выражений можно определить добротность Q колебательной системы объекта по одной из следующих формул как
или
Если из-за помех или других причин положение максимума функции |G(f)| определить затруднительно, то резонансную частоту fp определяют из выражения
Влияние помех можно существенно уменьшить, если уровни K1 и К2 выбрать для таких частот f1 и f2, при которых функция |Ga(f)| имеет наибольшую крутизну. В этом случае частоты f1 и f2 определяют, приравнивая к нулю вторые производные K1 и К2 по частоте выражений (8). После подстановки найденных значений частот f1 и f2 в выражения (8) получают оптимальные значения K1 и К2 равными
В этом случае резонансную частоту fp и добротность Q определяют подстановкой значений K1 и К2 из выражения (12) в выражения (11) и (9).
Получаем, что
Если принять значения K1 и К2 на общепринятом для определения полосы пропускания колебательных систем уровне
то из выражений (11) и (9) следует
При изменении частоты возбуждающего воздействия F(t) как функции времени t по закону, близкому к линейному, и постоянном уровне гармонического возбуждающего воздействия спектральная плотность e(f) модели возбуждающего воздействия e(t) будет близка к функции постоянного уровня, т.е. будет мало зависеть от текущего значения частоты возбуждающего воздействия F(t), поэтому в качестве частотной характеристики Gв(f) для определения параметров колебательной системы объекта: добротности Q и резонансной частоты fp, можно использовать спектральную плотность x(f) параметра колебаний объекта x(t), которая по форме будет повторять ее.
Колебательная система с сосредоточенными параметрами и с одной степенью свободы, как известно, описываются линейным неоднородным дифференциальным уравнением (л.н.д.у.) второго порядка. Считая, что возбуждающее воздействие действует на колебательную систему с линейно изменяющейся во времени частотой ω=ωн+vft, свободный член дифференциального уравнения запишется в виде
С учетом сказанного, уравнение колебаний объекта можно будет записать в виде
где x(t) - уравнение (функция) колебания объекта от времени;
δ - коэффициент затухания колебательной системы;
ω0 - резонансное значение угловой частоты колебаний объекта;
А0 - значение амплитуды ускорения, вызванного действующим на объект возбуждающим воздействием;
ωн - начальное значение угловой частоты колебаний объекта;
vf - скорость изменения частоты действующего на объект возбуждающего воздействия;
t - время;
φ0 - начальная фаза в уравнении колебаний объекта.
Зависимости, показанные на чертеже, получены из исследования линейного неоднородного дифференциального уравнения (л.н.д.у.) второго порядка с параметрами: добротность Q=30, резонансная частота 2πfр=333 рад/с, скорость изменения частоты F возбуждающего воздействия vf=160 рад/с2, амплитуда ускорения, вызванная действующим на объект возбуждающим воздействием А0=100 м/с2.
Предлагаемый способ позволяет определять резонансную частоту и добротность колебаний объекта в процессе эксплуатации без влияния на режим работы объекта для проведения измерений, без остановок и задержек на резонансных частотах, что уменьшит вероятность поломок объекта в процессе эксплуатации. Предлагаемый способ также позволяет определять резонансную частоту и добротность колебаний объекта в процессе его запуска и остановки. Реализация способа не требует вмешательства в режим работы объекта, поэтому не требует соответствующих дополнительных затрат времени для выполнения необходимых замеров. Все необходимые для определения параметров колебательной системы объекта вычислительные преобразования могут быть выполнены на современных вычислительных средствах за малое по сравнению с продолжительностью работы объекта время, поэтому способ обладает существенно большим быстродействием по сравнению с известными способами. Операции предлагаемого способа могут быть легко автоматизированы для проведения измерений, что во многих случаях позволит уменьшить время и квалификацию исследователя для определения необходимых параметров.
Claims (2)
1. Способ определения резонансной частоты и добротности колебаний объекта, заключающийся в том, что выявляют резонансную область частот (fн; fв), в процессе работы, когда объект находится под возбуждающим воздействием, измеряют и регистрируют параметр колебания объекта, отличающийся тем, что измеряют и регистрируют параметр колебания объекта в процессе работы как функцию времени x(t) при условии нахождения частоты возбуждающего воздействия в резонансной области частот (fн; fв] и по возможности близкому к постоянному уровню возбуждающем воздействии, по зарегистрированному параметру колебания объекта x(t) как функции времени t судят о резонансной частоте fр и добротности Q.
2. Способ по п.1, отличающийся тем, что из зарегистрированных значений параметра колебаний объекта в процессе работы как функции времени x(t) выбирают интервал времени с монотонным изменением частоты возбуждающего воздействия внутри всей резонансной области частот [fн; fв], для выбранного интервала времени по зарегистрированным значениям параметра колебания объекта x(t) определяют спектральную плотность параметра колебания объекта x(f), по спектральной плотности параметра колебания объекта x(f) судят о резонансной частоте fр и добротности Q колебательной системы объекта.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2004106169/28A RU2265193C1 (ru) | 2004-03-02 | 2004-03-02 | Способ определения резонансной частоты и добротности колебаний объекта |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2004106169/28A RU2265193C1 (ru) | 2004-03-02 | 2004-03-02 | Способ определения резонансной частоты и добротности колебаний объекта |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2004106169A RU2004106169A (ru) | 2005-08-10 |
RU2265193C1 true RU2265193C1 (ru) | 2005-11-27 |
Family
ID=35844819
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2004106169/28A RU2265193C1 (ru) | 2004-03-02 | 2004-03-02 | Способ определения резонансной частоты и добротности колебаний объекта |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2265193C1 (ru) |
-
2004
- 2004-03-02 RU RU2004106169/28A patent/RU2265193C1/ru not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
RU2004106169A (ru) | 2005-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9897518B2 (en) | Method and apparatus for measuring damping in a workpiece | |
RU2362136C1 (ru) | Способ испытания конструкции на ударные воздействия | |
RU2265193C1 (ru) | Способ определения резонансной частоты и добротности колебаний объекта | |
RU2262671C1 (ru) | Способ определения резонансной частоты, добротности, амплитуды стационарных резонансных колебаний объекта и уровня возбуждающего воздействия | |
US6116088A (en) | Method of operating a machine for stress relieving workpieces | |
RU2264605C1 (ru) | Способ определения резонансной частоты, добротности, амплитуды стационарных резонансных колебаний объекта и уровня возбуждающего воздействия | |
RU2285247C2 (ru) | Способ определения резонансной частоты, добротности и амплитуды стационарных резонансных колебаний | |
RU2239165C1 (ru) | Способ определения механических параметров и параметров резонансных колебаний | |
RU2086943C1 (ru) | Способ определения логарифмического декремента колебаний | |
RU2642155C1 (ru) | Стенд для виброакустических испытаний моделей систем виброизоляции судовых энергетических установок машинного отделения судна | |
RU2292026C1 (ru) | Способ определения динамических характеристик механической системы | |
SU560168A1 (ru) | Способ определени декремента колебаний по ширине резонансного пика | |
US6530270B1 (en) | Method and apparatus for analyzing a drivetrain assembly | |
RU2797126C1 (ru) | Устройство для измерения прочности бетона | |
RU2624411C1 (ru) | Способ определения добротности механической колебательной системы | |
RU2466368C1 (ru) | Способ определения динамических характеристик тензометрического преобразователя давления (варианты) | |
RU2778631C1 (ru) | Способ диагностики технического состояния трубопровода путем анализа декремента модального затухания | |
SU619863A1 (ru) | Установка дл определени чувствительности пьезоаксеерометров к деформации объекта измерени | |
RU2619812C1 (ru) | Способ неразрушающего контроля скрытых дефектов в технически сложном элементе конструкции, к которому нет доступа, и устройство для его осуществления | |
SU1165937A1 (ru) | Фазовый способ определени характеристик рассе ни энергии колебаний | |
US20220146367A1 (en) | Resonant frequency vibrational test | |
SU1002896A1 (ru) | Способ определени характеристики внутреннего трени в гибком элементе | |
RU2108502C1 (ru) | Способ определения относительных коэффициентов демпфирования механических и электромеханических колебательных систем по ускорению | |
JP4376588B2 (ja) | 地盤探査装置 | |
SU798185A1 (ru) | Способ контрол процесса вибро-ОбРАбОТКи КОНСТРуКций и уСТРОйСТВОдл ЕгО ОСущЕСТВлЕНи |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20060303 |