ОТРАСЛЕВОЙ СТАНДАРТ
Распоряжением Министерства от 30 декабря 1980 г. № 087-16 срок введения установлен с 1 января 1982 г. Настоящий стандарт устанавливает методику выполнения измерений вектора линейного ускорения, перегрузок, действующих на пилотируемые летательные аппараты (ЛА) в процессе проведения их летных испытаний. 1. ОБЩИЕ ПОЛОЖЕНИЯ1.1. Вектор линейного ускорения , обусловленный действием на ЛА результирующей аэродинамической силы , силы тяги и силы отдачи при стрельбе из бортового оружия, определяется для характерной точки ЛА и выражается отношением результирующей к массе m ЛА, т.е. ____________________ * Все физические величины в формулах, приведенных в стандарте, выражены в единицах, установленных Международной системой единиц. 1.2. Под перегрузкой, действующей на ЛА, понимается отношение действующего значения вектора линейного ускорения к стандартному значению ускорения свободного падения gc т.е. (2) или с учетом формулы (1) (3) В условиях пробега ЛА по взлетно-посадочной полосе при взлете и приземлении результирующая сила в этом случае включает в себя силы реакции Земли. 1.3. Измерение вектора линейного ускорения, перегрузок производится измерительными преобразователями, представляющими собой низкочастотные линейные акселерометры (в дальнейшем изложении - акселерометры). Типы акселерометров, основные параметры и технические требования, предъявляемые к ним, устанавливаются в соответствии с ОСТ 1 03751-86. Допускается применение акселерометров, принцип действия которых основан на других методах преобразования, но удовлетворяющих требованиям настоящего стандарта. 1.4. Внешние условия применения акселерометров при выполнении измерений в процессе проведения летных испытаний ЛА характеризуются параметрами, установленными ГОСТ 4401-81, особенностями конструкции ЛА и условиями выполнения полета (высотой h, скоростью Vn, широтой места j, направлением полета ψ и др.). Термин «летные испытания» приведен в справочном приложении 1. Основные параметры внешних условий работы акселерометров приведены в рекомендуемом приложении 2. Примечание. Под внешними условиями понимаются такие условия, воздействие которых оказывает влияние на метод измерения и приводит к появлению методических погрешностей измерения. 1.5. Градуировка и поверка акселерометров, предназначенных для определения вектора линейного ускорения, перегрузок, действующих на ЛА, должны производиться в соответствии с требованиями ГОСТ 8.179-76 и ГОСТ 8.476-82, определяющих в государственной системе обеспечения единства измерений соответственно: государственный первичный эталон и общесоюзную поверочную схему для средств измерений постоянного линейного ускорения твердого тела в диапазоне 0,001 - 200 м/с2 и государственный специальный эталон и общесоюзную поверочную схему для средств измерений перемещения, скорости и ускорения при колебательном движении твердого тела в диапазоне частот 0,5 - 104 Гц. 1.6. Градуировка каждого акселерометра и определение метрологических характеристик при его эксплуатации выполняются предприятием-потребителем. Поверяемые метрологические характеристики должны выбираться из числа, устанавливаемых ГОСТ 8.009-84 и ОСТ 1 00181-75. 1.7. Требования, предъявляемые при установке акселерометров на ЛА, обусловлены, главным образом, необходимостью уменьшения влияния на результаты измерений непараллельности измерительных осей акселерометров относительно связанной системы осей координат ЛА, а также влияния деформации (за счет нагруженности, изменения температуры и др.) элементов конструкции ЛА, к которым крепятся акселерометры. 2. МЕТОДИКА ОПРЕДЕЛЕНИЯ ВЕКТОРА ЛИНЕЙНОГО УСКОРЕНИЯ, ПЕРЕГРУЗОК, ДЕЙСТВУЮЩИХ НА ЛА ПРИ ПРОВЕДЕНИИ ЛЕТНЫХ ИСПЫТАНИЙ2.1. Методика определения вектора линейного ускорения и перегрузок основана на измерении трех составляющих - проекций вектора на связанную систему осей координат ЛА, а именно: - вертикальной составляющей вектора линейного ускорения (перегрузки) ; - продольной составляющей вектора линейного ускорения (перегрузки) ; - поперечной (боковой) составляющей вектора линейного ускорения (перегрузки) 2.2. Знак линейного ускорения, действующего вдоль какой-либо оси координат, совпадает со знаком соответствующей силы. При этом отклонение массы чувствительного элемента вдоль этой оси противоположно направлению действия силы. Направление движения инертной массы чувствительных элементов при действии на ЛА положительно направленных составляющих векторов и показано на черт. 1. Черт. 1 2.3. Векторы линейного ускорения и перегрузок определяются по формулам: - модули векторов (4) (5) - направляющие косинусы для векторов (6) (7) При этом aj (nj) должны браться со своими знаками. 2.4. Компоновочная схема размещения акселерометров на ЛА представлена на черт. 2. Положение акселерометров относительно центра масс ЛА характеризуется координатами:
Схема установки акселерометра на ЛА Черт. 2 2.5. Значения составляющих векторов линейного ускорения, перегрузок в центре масс ЛА определяются по формулам: (8) (9) где - измеренные акселерометрами значения составляющих векторов линейного ускорения, перегрузок; - поправки на положение акселерометров относительно центра масс ЛА (за счет вращения ЛА относительно своего центра масс). Расчет поправок выполняется по формулам (10) (11) где ωx, ωy, ωz - составляющие вектора угловой скорости, направленные по осям связанной системы координат ЛА; - соответственно составляющие векторов углового ускорения ЛА. 2.6. Не параллельность измерительных осей акселерометров относительно осей связанной системы координат ЛА характеризуется установочными углами Daуст и Dgуст (черт. 3). Положительное значение Daуст принимается при наклоне измерительной оси акселерометра назад, отрицательное - при наклоне вперед. Положительное значение Dgуст принимается при наклоне измерительной оси акселерометра вправо (если смотреть с конца оси ОХ), отрицательное - при наклоне влево. При внесении поправки на установочный угол, получим: (12) Установочный угол акселерометра
Установочный угол акселерометра Установочный угол акселерометра
Черт. 3 (13) 2.7. Значения составляющих вектора линейного ускорения и перегрузок в скоростной системе осей координат (воздушной) определяются по формулам: (14) (15) где a - угол атаки ЛА; b - угол скольжения ЛА (ГОСТ 20058-80). При выполнении полета в плоскости О, Y, Ya, X, Xa (b = 0) значения составляющих вектора линейного ускорения, перегрузок в скоростной системе осей координат (воздушной) определяются по формулам: (16) (17) Методика измерения угла атаки при летных испытаниях приведена в рекомендуемом приложении 3. 2.8. При необходимости определения вектора линейного ускорения, перегрузок в других системах осей координат (в нормальной земной системе осей координат Оо, Xg, Yg, Zg; в скоростной системе осей координат О, Xn, Yn, Zn, связанной с пространственным углом атаки; в траекторией системе осей координат О, Xк, Yк, Zк и др.), пересчет вектора линейного ускорения, перегрузок должен производиться с использованием направляющих косинусов, составленных на основании углов, определяющих положение осей координат относительно друг друга (ГОСТ 20058-80). Взаимное положение связанной системы осей координат ЛА относительно земной системы осей координат, скоростной системы осей координат, траекторной и др., а также углы, характеризующие положение осей относительно друг друга, показаны на черт. 1 справочного приложения 4. 2.9. В расчетах при определении характеристик ЛА в ряде случаев необходимо знать действующее значение ускорения свободного падения ЛА. Ускорение свободного падения определяется по формуле где - коэффициент аппроксимации, равный 9,8235 м/с2; h - высота полета; rЗ - радиус Земли, равный 6360·103 м; ωv - угловая скорость ЛА относительно центра Земли за счет полета, рассчитанная по формуле
где ωЗ - угловая скорость вращения Земли, равная 0,728·10-4 рад/с; Vn - скорость полета ЛА; j - географическая широта полета; ψ - направление полета. Вывод формулы (18) приведен в обязательном приложении 5. 3. МЕТОДИКА ВЫБОРА АКСЕЛЕРОМЕТРОВ3.1. Методика выбора акселерометров основана на определении основных метрологических характеристик акселерометров, удовлетворяющих требованиям, предъявляемым к измерению составляющих вектора линейного ускорения, перегрузок при летных испытаниях ЛА. 3.2. Пригодность акселерометра для применения в зависимости от этапов и задач испытаний ЛА характеризуется допустимой суммарной погрешностью измерения, а также требуемым коэффициентом преобразования и диапазоном регистрируемых частот. 4. ПОВЕРКА И ГРАДУИРОВКА АКСЕЛЕРОМЕТРОВ4.1. При проведении поверки следует выполнять операции, указанные в табл. 1. Таблица 1
4.2. Перечень средств поверки и их нормативно-технические характеристики представлены в табл. 2 и рекомендуемом приложении 6. Таблица 2
4.3. При проведении поверки акселерометров должны соблюдаться нормальные условия измерений: - температура окружающей среды 293 К ± 5 К (20 °С ± 5 °С); - атмосферное давление, 101,3 КПа (760 мм рт. ст.); - относительная влажность (65 ± 15) %;, - уровень акустических полей, не более 60 дБ; - частота переменного тока сети питания (50 ± 0,5) Гц; - отклонение напряжения питания от номинального значения, указанного в технических условиях и эксплуатационной документации на акселерометр конкретного типа, не более ±2 %; - уровень внешних электрических и магнитных полей, а также действия вибрации должны быть в пределах норм, установленных в технических условиях и эксплуатационной документации на акселерометр конкретного типа и промежуточные преобразователи измерительной информации, входящие в комплект поверочных средств. 4.4. Подготовка и поверка образцовых поверочных и вспомогательных средств измерений, а также крепления поверяемых акселерометров на образцовых установках должны выполняться в соответствии методик поверки, применяемых при поверках. 4.5. Проведение поверки 4.5.1. Одновременно с внешним осмотром акселерометра производится также проверка его массы, габаритных и присоединительных размеров. 4.5.1.1. При внешнем осмотре должно быть установлено: - отсутствие электрических и механических повреждений; - отсутствие грязи или смазки на штепсельных разъемах; - наличие приспособления, которое препятствует вскрытию акселерометра без нарушения клейма или пломбы; - наличие пломбы; - наличие маркировки акселерометра; - наличие паспорта. 4.5.1.2. При взвешивании акселерометра погрешность оценки его массы должна быть не более ±5 % от измеряемой величины. 4.5.1.3. Габаритные и установочные размеры акселерометров должны соответствовать требованиям ОСТ 1 03751-86. 4.5.1.4. В случае несоответствия акселерометра хотя бы одному из указанных требований, он признается непригодным к применению (проводятся мероприятия по устранению недостатков). 4.6. Опробование акселерометра 4.6.1. Опробование акселерометра производится посредством поверочного лабораторного пульта, с помощью которого акселерометр подключается к источнику питания и осуществляется подача на акселерометр стимулирующих сигналов; выходной сигнал акселерометра подается на визуальный показывающий прибор. 4.6.2. При опробовании проверяется функционирование акселерометра при изменении его положения в гравитационном поле Земли и при подаче на него стимулирующих сигналов, определяется потребляемая сила тока и характеристика выходного сигнала, для чего, если в этом есть необходимость, производится легкое постукивание по акселерометру. 4.7. Определение сопротивления изоляции и проверка ее электрической прочности 4.7.1. Сопротивление изоляции проверяется между измерительными цепями и корпусом акселерометра (при этом все выходные контакты измерительной цепи должны быть соединены между собой накоротко), а также между измерительными цепями и цепями вспомогательных элементов акселерометра (цепями питания электромоторов, термостатов и др.). 4.7.2. Сопротивление изоляции должно быть не менее 20 МОм или значения, указанного в технических условиях на акселерометр. 4.7.3. Электрическая прочность изоляции проверяется между измерительными цепями и корпусом акселерометра, а также между измерительными цепями и цепями вспомогательных элементов акселерометра. 4.7.4. Проверка электрической прочности изоляции проводится на установке, позволяющей плавно повышать испытательное напряжение от нуля до заданного значения, указанного в технических условиях на акселерометр. Мощность установки (на стороне высокого напряжения) должна быть не менее 0,25 кВ·А, и измерение испытательного напряжения должно проводиться на стороне высокого напряжения. Форма кривой напряжения должна быть такой, чтобы отношение амплитудного значения напряжения к действующему находилось в пределах 1,34 - 1,48. Скорость изменения напряжения от нуля до заданного значения должна быть за время от 5 до 20 с. Уменьшение напряжения до нуля должно производиться с той же скоростью. 4.8. При проведении работ по поверке акселерометров должны соблюдаться «Правила технической эксплуатации электроустановок потребителей» и «Правила техники безопасности при эксплуатации электроустановок потребителей», утвержденные Госэнергонадзором, и требования, установленные ГОСТ 12.2.007.0-75. На каждый тип поверочной установки должны быть разработаны правила проведения работ на ней и правила техники безопасности. 4.9. Поверку акселерометров должны проводить лица, прошедшие курс обучения и имеющие право на проведение поверки. 4.10. Определение градуировочной характеристики акселерометров 4.10.1. Определение градуировочной характеристики акселерометра производится в процессе градуировки его, при этом определяются: - диапазон измерения (по выходу); - номинальная градуировочная характеристика, выражающая зависимость выходного сигнала либо в функции линейного ускорения (), либо в функции перегрузки : (19) (20) где - ускорения, воспроизводимые образцовыми установками; gc - стандартное ускорение свободного падения, равное 9,80665 м/с2; y - величина выходного сигнала. 4.10.2. При применении в процессе градуировки акселерометра поворотной установки (использование принципа поворота акселерометра в гравитационном поле Земли) должно быть известно местное значение ускорения свободного падения (gм), полученное с точностью до пятого знака. Значение линейного ускорения определяется в соответствии с наклоном поворотного стола пропорционально местному значению ускорения свободного падения, т.е. (21) Значение перегрузки в этом случае рассчитывается по формуле (22) где a - фиксированные значения угла наклона плоскости поворотного стола к уровню горизонта. В случае, если градуировка акселерометра производится на центрифуге, воспроизводимое постоянное ускорение, перегрузки определяются по формулам (23) (24) где ωст.ц - угловая скорость вращения стола центрифуги; Rст.ц - расстояние от оси вращения до центра масс чувствительного элемента акселерометра. Градуировка акселерометров должна производиться в соответствии с требованиями ОСТ 1 00108-73 и ОСТ 1 00181-75. Градуировочная характеристика средства измерений имеет вид степенного полинома и определяется оценкой коэффициентов a0; а1; a2, ..., aк: х = а0 + а1у + а2у2 + ... + aкук, (25) где x - входной сигнал средства измерения; y - выходной сигнал средства измерения; к - максимальная степень полинома, которым выражена градуировочная характеристика рассматриваемого средства измерения; aк - коэффициент степенного полинома. При обработке результатов измерений выбор степени полинома, аппроксимирующего градуировочную характеристику, осуществляется по критерию минимума СКО случайной составляющей погрешности. 4.11. Определение основной погрешности акселерометра 4.11.1. При нормировании метрологических характеристик индивидуально градуируемых средств основная погрешность включает: - погрешность образцовых средств измерений; - случайную составляющую погрешности; - нестабильность градуировочной характеристики во времени. 4.11.2. Случайная составляющая погрешности характеризуется рассеянием градуировочных точек относительно полученного аппроксимирующего полинома. СКО случайной составляющей погрешности определяется по формуле (26) где N - число измерений; к - степень полинома, аппроксимирующего градуировочную характеристику; xi - значение входной величины при градуировке для выходного сигнала; - значение входной величины, оцененное по полиному для того же выходного сигнала. Аппроксимация градуировочной зависимости осуществляется по методу наименьших квадратов. При обнаружении больших величин СКО случайной составляющей погрешности необходима тщательная проверка числового материала на наличие грубых погрешностей. Исключение грубых погрешностей измерений производится: - на предварительном этапе рассмотрения результатов измерения, если известно о некондиционной работе средства измерения; - выделение грубых погрешностей из ряда измерений с применением методов математической статистики. Для этого определяется принадлежность результатов градуировки к нормальному закону распределения с учетом количества результатов наблюдений (при N > 50, при N < 50). Грубые погрешности исключаются. При распределении случайной величины по закону, отличному от нормального (закон неизвестен), значение критерия для определения грубых, погрешностей определяется по формулам если xв < ; (27) если xв > . (28) Если любое из неравенств выполняется, то значение xв выделяется как грубая ошибка. Данный критерий использует аппроксимирующую систему распределения Пирсона. Аa/2 и А1-a/2 - квантили системы распределения Пирсона для уровня значимости a - вероятности ошибки первого рода (уровень значимости) - определяются в зависимости от третьего и четвертого центральных моментов ряда измерений. Оценки третьего (μ3) и четвертого (μ4) центральных моментов рассчитываются по формулам: (29) (30) Значения квантилей Аa/2 и А1-a/2 нормированных случайных величин системы распределения Пирсона в зависимости от и μ4 для a = 0,005 приведены в табл. 3 и 4. Таблица 3
Таблица 4
4.11.3. Нестабильность градуировочной характеристики оценивается как разница в значениях полиномов, определенных через определенный интервал времени. Необходимо сравнивать полиномы одной и той же степени. 4.12. Действия влияющих величин выражаются: - в сдвиге и повороте градуировочной характеристики, относительно градуировочной характеристики, полученной в нормальных условиях; - в изменении случайной составляющей погрешности. 4.12.2. Для каждого из значений влияющих величин производится аппроксимация градуировочной характеристики. В результате аппроксимации в общем случае получается совокупность градуировочных характеристик: yμ(x) = a0μ + a1μx + a2μ(x)2 + ... + aкμ(x)к (31) где μ = 1, 2, ..., с - номера значений влияющих величин, имеющих место при эксперименте. Функции влияния для каждого из коэффициентов aкμ определяются по методу наименьших квадратов в виде отклонений величины коэффициентов от значений, полученных в нормальных условиях. Обычно степень полинома, аппроксимирующего градуировочную характеристику, не превышает 2 - 3-й степени, а полином, определяющий функцию влияния, - 1 - 2-й степени. В результате аппроксимации получается функция влияния вида Daμy = b1μ(zμ - zμн) + b2μ(zμ - zμн)2, (32) где zμ и zμн - значение μ-й влияющей величины, измеренные и соответствующие нормальным условиям; b1, b2 - коэффициент функции влияния. Данная функция влияния получена при раздельном действии влияющих величин для каждой влияющей величины. 4.13. Определение динамических характеристик акселерометра Динамическая характеристика выражается передаточной функцией W(p), которая представляет собой отношение преобразования Лапласа y(p) выходного сигнала y(t) к преобразованию Лапласа x(р) входного сигнала x(t), т.е. (33) где р - комплексная переменная. Если р принимает чисто мнимое значение, т.е. р = jω, то передаточная функция выразится зависимостью W(jω) = B(ω) + jм(ω). (34) Функции B(ω) и М(ω) называются соответственно вещественной и мнимой частотными характеристиками. Выражение для частотной характеристики может быть представлено в виде W(jω) = A(ω)ejj(ω). (35) Вещественные функции A(ω) и j(ω) называются амплитудной и фазовой частотными характеристиками. 4.14. Оформление результатов поверки 4.14.1. Положительные результаты поверки акселерометров вносятся в паспорт акселерометра и специальный журнал прохождения периодических поверок данного акселерометра. 4.14.2. Результаты поверки и градуировки акселерометров, а также их метрологические характеристики заносятся в протокол, который приведен в рекомендуемом приложении 7. 4.14.3. Акселерометры, не соответствующие требованиям технических условий на них, признаются непригодными для эксплуатации; в их паспорта заносятся отметки о непригодности с указанием причин неисправности. 5. ТРЕБОВАНИЯ К РАЗМЕЩЕНИЮ И СПОСОБЫ КРЕПЛЕНИЯ АКСЕЛЕРОМЕТРОВ НА ЛЕТАТЕЛЬНОМ АППАРАТЕ5.1. Акселерометры на ЛА должны устанавливаться в его центре масс, и при этом направление осей акселерометров должно совпадать с направлением осей связанной системы координат ЛА. 5.2. Для снижения систематических ошибок измерения, обусловленных установкой акселерометра на ЛА, в общем случае допустимая непараллельность осей должна соответствовать требованиям технических условий на ценный тип акселерометра. В особых случаях, в частности, при определении аэродинамических и тяговых характеристик ЛА (Сx, Сy, Сp), допустимая не параллельность осей не должна превышать 5 - 10'. Если плоскость симметрии акселерометра совпадает с плоскостью симметрии ЛА или параллельна ей, то допуск на непараллельность продольных и вертикальных осей должен соответствовать требованиям технических условий на данный тип акселерометра. 5.3. Если акселерометры установлены вне центра масс ЛА, производится определение координат акселерометра относительно центра масс. Координаты акселерометров заносятся в кроссировочную таблицу измерений. Погрешность определения координат не должна превышать ± 20 мм. В связи с тем, что положение центра масс ЛА в полете изменяется, в ряде случаев при обработке результатов (при пересчете измеренных значений линейных ускорений, перегрузок) необходимо учитывать изменение положения центра масс относительно первоначального положения, установленного до полета. 5.4. Материал и жесткость конструкции кронштейнов выбираются такими, чтобы обеспечить заданную не параллельность осей акселерометра и ЛА в процессе полета при воздействии на конструкцию перегрузок и возможных перепадов температур. Кроме того, жесткость конструкции кронштейна не должна приводить к возникновению резонансных колебаний кронштейна с преобразователем. 5.5. Крепление кронштейнов с акселерометрами должно производиться к таким силовым элементам ЛА, деформации которых незначительны. 5.6. В целях обеспечения технологичности в эксплуатации целесообразно акселерометры закреплять на переходном кронштейне, устанавливаемом на нивелировочной площадке. Нивелировочная площадка на ЛА должна строго выставляться в плоскости симметрии продольной и поперечной осей связанной системы координат ЛА либо параллельно ей. Допустимая непараллельность установочной плоскости нивелировочной площадки относительно плоскости O, X, Z должна быть не более ± 15'. 5.7. Для установки нивелировочной площадки акселерометров производится выставление связанных осей ЛА посредством нивелирования в плоскости горизонта, соблюдая все правила нивелирования данного типа ЛА (не менее чем при четырех точках опоры). После чего устанавливается нивелировочная площадка акселерометров строго горизонтально посредством уровня или оптического угломера КО-10. 5.8. На нивелировочной площадке акселерометров наносятся риски, которые ориентируются относительно осей ОХ, OZ связанной системы координат ЛА. Кронштейн, на котором выставляются оси акселерометра относительно его базовой поверхности, также должен иметь риски, позволяющие выставить измерительную ось акселерометра относительно рисок нивелировочной площадки. 5.9. После установки акселерометров на ЛА производится регистрация установочных углов (Daуст и Dgуст) по отношению к связанной системе осей координат ЛА (последний должен находиться в плоскости горизонта). На нивелировочной площадке или в удобном месте на ЛА выбирается площадка для установки угломера, посредством которого при подготовке измерительных средств определяется угловое положение ЛА относительно горизонта и оценивается правильность работы акселерометра в гравитационном поле Земли. 5.10. Условия эксплуатации акселерометров на ЛА должны соответствовать требованиям ОСТ 1 03751-86. 5.11. Основные измерительные средства, с помощью которых производится установка и контроль установки акселерометров на ЛА: - квадрант оптический - КО-10 ГОСТ 14967-80. - нивелир высокоточный - НА-1 ГОСТ 10528-76; - теодолит - ТБ-1, ТТ-5 ГОСТ 10529-86; - рейка ГОСТ 11158-83; - рулетка - Р3-20 ГОСТ 7502-80. 6. АЛГОРИТМ ИЗМЕРЕНИЙ6.1. Общий алгоритм измерений должен соответствовать установленному согласно техническим условиям на средства измерений. ПРИЛОЖЕНИЕ 1Справочное ТЕРМИНЫ И ИХ ОПРЕДЕЛЕНИЯ
ПРИЛОЖЕНИЕ 2Рекомендуемое 1. ВНЕШНИЕ УСЛОВИЯ РАБОТЫ АКСЕЛЕРОМЕТРА1.1. Под внешними условиями понимаются такие условия, воздействие которых оказывает влияние на метод измерения и приводит к появлению методических погрешностей измерения. К ним относятся: - условна окружающей среды (воздушной среды, гравитационного поля Земли, вращение и форма Земли и др.); - направление полета и местонахождение ЛА; - угловые скорости и ускорения относительно центра масс, вызванные маневрированием ЛА; - угловые скорости и ускорения, вызванные упругими изгибными и крутильными деформациями конструкции ЛА; - статические, изгибные и крутильные деформации конструкции ЛА. 1.2. Воздушная среда земной атмосферы 1.2.1. Воздушная среда земной атмосферы характеризуется давлением, плотностью, температурой, влажностью, турбулентностью, скоростью течений воздушной среды, параметры которых меняются в зависимости от времени года и суток, координат места, метеорологических условий, солнечной активности и других факторов. 1.2.2. Вместо действительных характеристик атмосферы, являющихся случайными функциями времени и пространства, влиянием которых на работу акселерометров в ряде случаев можно пренебречь, а в ряде случаев влияние которых необходимо избежать, в качестве модели атмосферы принимают характеристики, соответствующие ГОСТ 4401-81. В ГОСТ 4401-81 вместо действительных характеристик приведены средние значения (математические ожидания) их при следующих граничных условиях: давление воздуха над уровнем моря (поверхность героида) на широте 45°, кПа (мм рт. ст.)........................................ 101,3 (760) температура, К (°С)........................................................... 288,15 (+15) плотность воздуха, кг/м3................................................... 1,225 скорость звука, м/с............................................................ 340,294 Действительная нестационарность параметров атмосферы приведет к некоторому разбросу результатов измерений. 1.3. Параметры Земли, влияющие на результаты измерений 1.3.1. Особенности геометрической формы Земли, строение поверхности, вращение и параметры гравитационного поля оказывают значительное влияние на работу акселерометров. 1.3.2. Точное математическое описание формы Земли с практической точки зрения неудобно, поэтому имеется ряд моделей, упрощающих эту задачу. Более точные модели предполагают, что Земля является эллипсоидом, симметричным полярному и экваториальному диаметрам и построенному по среднему уровню моря. 1.3.3. В практике для решения различных задач применяются следующие модели Земли: в США используется эллипсоид Кларка, в европейских странах - эллипсоид Бесселя, в Советском Союзе - эллипсоид Красовского со следующими размерами: - большая полуось (экватоиальный радиус Земли) a = 6378245 м; - малая полуось (полярный радиус Земли) b = 6356863 м; - эллиптичность (сжатие) 1.3.4. Для решения задач, не требующих высокой точности, Земля по форме принимается близкой к шару-сфероиду, описываемому простым математическим выражением. Клеро вывел зависимость между изменением силы тяжести от экватора к полюсам и сжатием сфероида. Его вывод был точен только для малых членов порядка величины сжатия, при такой точности сфероид и эллипсоид неразличимы. 1.3.4.1. За сфероид принимается эквипотенциальная поверхность гравитационного поля Земли, охватывающая всю ее массу. Начало координат 0 помещается в центре масс Земли, ось Z декартовых координат направлена по оси вращения, оси X и У расположены в плоскости экватора, ось X расположена в плоскости начального меридиана (черт. 1). 1.3.4.2. Расстояние r от центра Земли до точки М на поверхности сфероида (модуль геоцентрического радиуса-вектора ) находится из уравнения Клеро Используя сферические и полярные координаты, выполнив соответствующие преобразования в уравнении (1) и произведя упрощение, для определения r пользуются формулой r = a(1 - esin2j). (2) а - сферическая и полярная системы координат; б - направление вертикалей Черт. 1 В ГОСТ 4401-81 за условный радиус Земли принят радиус r3 = 6356766 м, при котором ускорение свободного падения и вертикальный градиент на среднем уровне моря наиболее близки к истинным на широте 45°32'33''. 1.3.4.3. Для решения практических задач в сфероидальной Земле используются следующие вертикали (см. черт. 1): - геоцентрическая вертикаль - линия, характеризующая положение точки и проходящая через центр Земли; - линия отвеса - линия, характеризующая направление вектора силы тяжести, являющегося суммарным вектором гравитационного ускорения и центростремительного ускорения, вызванного вращением Земли; - гравитационная вертикаль - линия, совпадающая с направлением поля тяготения Земли. 1.3.4.4. Угол между геоцентрической вертикалью и плоскостью экватора образует геоцентрическую широту j. Угол между линией отвеса и экваториальной плоскостью образует географическую широту j'. Угол между гравитационной вертикалью и экваториальной плоскостью является гравитационной широтой j''. Угол между географической и геоцентрической широтами определяется по приближенной формуле j' - j = esin2j. (3) На широте 45° разность j' - j достигает максимума, т.е. 11'30''. Угол между направлением силы тяжести и силы тяготения определяется формулой j' - j'' = 5,9'sin2j. (4) 1.3.5. Существенную роль в работе акселерометра играет вращение Земли. 1.3.5.1. Земля совершает сложное движение, включающее: - вращение вокруг своей оси с Запада на Восток с периодом в 23 ч 56 мин 4,091 с = 86164,091 с среднего солнечного времени, или 24 ч = 86400 с звездного времени; угловая скорость при этом соответственно равна
Вектор угловой скорости Земли направлен по оси вращения от южного полюса к северному в соответствии с правилами законов для правых систем координат; - годичное обращение вокруг Солнца со средней скоростью движения по орбите 29,893 км/с; - нутационное колебание земной оси с периодом около 18,6 года и амплитудой колебаний, не превосходящей 11''; - прецессионное движение относительно оси эклиптики с периодом 25800 лет; - движение вместе с солнечной системой относительно других звезд. Все перечисленные составляющие движения Земли, кроме суточного вращения, в измерениях и расчетах не учитываются, так как их влияние чрезвычайно мало. 1.3.6. Земля - огромное сферическое тело с массой М = 5,975·1024 кг создает вокруг себя центральное гравитационное поле тяготения. Это поле притягивает к себе любые материальные тела, обладающие массой. 1.3.6.1. Гравитационное поле, создавая силу притяжения, вызывает ускоренное движение свободных материальных тел. Если тело неподвижно на земной поверхности, то возникает сила тяжести. Для тел, неподвижных относительно поверхности Земли, распределение силы тяжести в пространстве рассматривается как поле и характеризуется его потенциалом или геопотенциалом, характеризующим потенциальную энергию тел, расположенных в данной точке. 1.3.6.2. Потенциал поля тяжести зависит от гравитационного поля Земли и центростремительного ускорения, вызванного вращением Земли: (5) При этом гравитационный потенциал
а кинетическая энергия вращения равна где r - радиус-вектор; G - универсальная гравитационная постоянная; М - масса Земли; ω - угловая скорость вращения Земли; J, к - имперические константы. Произведение гравитационной постоянной G на массу Земли М постоянно и для расчетов может быть принято равным 3,986004·1014 м3/с2. 1.3.7. Как уже отмечалось, эквипотенциальная поверхность поля тяжести образует сфероид Земли. Направления действия сил тяжести, в том числе и отвеса, перпендикулярны эквипотенциальной поверхности поля тяжести, и указаны на черт. 2. 1.3.8. Ускорение силы тяжести имеет радиальную и горизонтальную составляющие, лежащие в плоскости меридиана, т.е. (6) (7) 1.3.9. Кроме основного гравитационного поля имеются гравитационные аномалии, вызванные неоднородным распределением плотности и масс (девиация масс) в земном геоде (горы, впадины). В ряде случаев в практике очень важно знать величину g применительно к той лаборатории, в которой проводятся работы. Этот вопрос может быть решен только путем непосредственного измерения величины gм в данном месте. 1 - радиус-вектор точки М; 2 - гравитационная вертикаль; 3 - нормаль к эллипсоиду Земли Черт. 2 1.3.10. Для решения поставленных задач в качестве модели Земли можно принять сферу с равномерным распределением плотности, а гравитационное поле считается центральным, сила притяжения Земли характеризуется ускорением свободного падения. При этом условии зависимость ускорения свободного падения на среднем уровне моря от географической широты вычисляют по формуле gj = 9,80616(1 - 0,0026373cos2j + 0,0000059cos22j). (8) Ускорение свободного падения неподвижного тела на высоте определяют по формуле (9) Стандартное значение ускорения свободного падения принято равным gc = 9,80665 м/с2. ПРИЛОЖЕНИЕ 3Рекомендуемое МЕТОДИКА ИЗМЕРЕНИЯ УГЛА АТАКИ ЛЕТАТЕЛЬНОГО АППАРАТА1. Задача измерения угла атаки ЛА решается путем применения прямого метода с использованием измерительных преобразователей флюгерного типа и путем применения косвенного Метода, основанного на измерении угла тангажа и скоростного угла тангажа ЛА. 1.1. Особенностью работы преобразователей флюгерного типа является то, что на результаты их показаний оказывает значительнее влияние характер обтекания измерительной системы воздушным потоком. Наибольшее влияние обуславливается формированием скачков уплотнения на лопасти флюгеров при малых околозвуковых числах М и взаимодействия скачков уплотнения, образованных на лопасти флюгера, и элементов конструкции ЛА (носовой штанги) при больших околозвуковых и сверхзвуковых числах М. В качестве примера влияния условий обтекания измерительных преобразователей флюгерного типа, установленного на носовой штанге ЛА, может быть рассмотрен график, представленный на черт. 1. a = 22° a = 16° a = 8° a = 0° Черт. 1 Как видно из графика, влияние обтекания на показания преобразователей флюгерного типа имеет прямую зависимость от угла атаки и числа М. Изменение показаний преобразователя флюгерного типа практически характеризуется линейной зависимостью от изменения числа М и учитывается путем введения в результаты измерений аэродинамической поправки. 1.1.1. Угол атаки определяется по формуле a = aпр - Daаэр, (1) где aпр - показания измерительного преобразователя флюгерного типа; Daаэр - аэродинамическая поправка к преобразователю флюгерного типа, устанавливаемого на носовой штанге ЛА. 1.2. Косвенным способом угол атаки измеряется как разность между углом тангажа и скоростным углом тангажа , т.е. a = - . (2) 1.2.1. Угол тангажа может измеряться с помощью гировертикали, предварительно откорректированной на режиме горизонтального полета. Способ измерения угла тангажа иллюстрирован на черт. 2. Установочные углы:
Черт. 2 Изменение угла тангажа определяется по формуле
Поправка рассчитывается по формуле
где , - поправки направлений гировертикали соответственно за счет вращения и кривизны Земли; ψ - курс полета; j - широта полета; ω3 - угловая скорость вращения Земли; r3 - радиус Земли; DL, t - расстояние и время пролета от момента выключения коррекции. 1.2.2. Истинный угол тангажа: (3) где - гировертикали; - установочный угол гировертикали относительно средней аэродинамической хорды САХ (определяется при установке ПА в линию горизонта); - поправка к показаниям гировертикали, обусловленная вращением и кривизной Земли. 1.2.2.1. Введение поправки вызвано тем, что в гировертикали совмещение оси OY с местной вертикалью обеспечивается только при включенной коррекции в прямолинейном установившемся горизонтальном полете. При выключении системы коррекции гировертикаль сохраняет свое неизменное положение в пространстве. В полете с момента выключения коррекции местная вертикаль будет отклоняться от оси OY гировертикали на величину за счет вращения и кривизны Земли. Значение величины поправки за интервал времени Dt, прошедшего от момента выключения коррекции, определяется по формуле (4) Обычно изменение широты невелико и им можно пренебречь, тогда формула для вычисления будет иметь вид: (5) 1.2.3. Для измерения скоростного угла тангажа на участке прямолинейного полета, близкого к горизонтальному, определяется начальное значение на основании измерения вертикальной скорости Vg и воздушной скорости ЛА V. При отсутствии специального преобразователя вертикальная скорость Vy может быть получена на основании измеренных значений h и вертикальной составляющей ускорения, перегрузки. В дальнейшем определяется как сумма начального значения и приращения , вычисленного на основе измерения вертикальной составляющей ускорения и воздушной скорости численным интегрированием с интервалом 0,5 с, т.е. (6) При условии, что (7) где Dt = 0,5 с; V и - скорость полета и скоростной угол тангажа в начале интервала; - среднее значение вертикальной составляющей перегрузки на интервале, рассчитанное по формуле (8) 1.2.4. Определение представляет собой типичную задачу, решаемую методом последовательных приближений. Количество потребных приближений будет определяться заданной разницей Da значения угла атаки, полученного по результатам предыдущего и последующего приближений. 1.2.4.1. При измерении центра масс ЛА в скоростной (воздушной) системе осей координат и при измерении и в связанной системе осей координат ЛА, угол атаки можно было бы вычислить из уравнения (9) При небольших углах атаки (a = 10°) при cosa = 1 и sina = a угол атаки определяется по формуле (10) В качестве средства непосредственного измерения может быть использован акселерометр, смонтированный на флюгарке. ПРИЛОЖЕНИЕ 4Справочное СИСТЕМА ОСЕЙ КООРДИНАТ, ИХ ВЗАИМНОЕ ПОЛОЖЕНИЕ1. Для оценки исследования движения ЛА и оценки работы акселерометра применяются прямоугольные правые системы координат: связанная система осей координат ЛА, нормальная земная система осей координат, скоростная система осей координат, траекторная система осей координат и др. (ГОСТ 20058-80). 2. Пространственное положение систем осей координат характеризуется соответствующими углами поворота одной системы осей координат относительно другой (ГОСТ 20058-80). 2.1. Положение связанной системы осей координат ЛА относительно нормальной земной системы осей координат определяется тремя углами: углом тангажа , углом крена g, углом рыскания (углом курса) ψ. 2.2. Положение скоростной системы осей координат относительно нормальной земной системы осей координат определяется углами: скоростным углом тангажа , скоростным углом крена ga, скоростным углом рыскания (курса) ψа. 2.3. Положение траекторной системы осей координат относительно нормальной земной системы координат характеризуется углом наклона траектории к горизонту θ и углом пути ψ. 2.4. Угол атаки a и угол скольжения b характеризуют положение скоростной системы осей координат относительно осей связанной системы координат ЛА. Взаимное положение связанной и скоростной системы осей координат относительно нормальной земной системы осей координат и относительно друг друга показано на чертеже.
ПРИЛОЖЕНИЕ 5Обязательное РАСЧЕТ УСКОРЕНИЯ СВОБОДНОГО ПАДЕНИЯ ЛЕТАТЕЛЬНОГО АППАРАТА1. Ускорение свободного падения ЛА в полете определяется из условий высоты полета h, скорости полета Vn, широты места j, направления полета ψ, скорости вращения Земли и гравитационного поля Земли и может быть выражено формулой (1) 2. Ускорение свободного падения тела , неподвижно расположенного на высоте h в гравитационном поле вращающейся Земли, определяется по формуле (2) где j - широта, град; h - высота; ω - угловая скорость вращения Земли, равная 0,728·10-4 рад/с; r3 - радиус Земли равный 6360·103 м; - постоянный коэффициент, выбранный из условия аппроксимации принятого закона изменения ускорения свободного падения с высотой и равный 9,8235 м/с2. Первый член этой формулы описывает изменение гравитационной составляющей ускорения свободного падения в зависимости от высоты полета, а второй - изменение gL за счет вращения Земли. 3. При движении ЛА относительно Земли с путевой скоростью Vn в абсолютном движении будет происходить сложение путевой скорости с круговой скоростью Земли, что в свою очередь вызывает появление дополнительного ускорения (3) 4. Для условий полета, представленных на черт. 1 (определение ускорения свободного падения ЛА), формулу определения суммарного центробежного ускорения можно вывести, исходя из следующих условий: - при полете с ψ = const меридиональная составляющая скорости , а широтная составляющая скорости - суммарная угловая скорость вращения в широтном и меридиональном направлении составят соответственно (4) (5) - составляющие центробежного ускорения за счет составляющих суммарной скорости вращения в том и в другом направлениях будут равны (6) (7) - составляющая ускорения Dgg направлена по радиусу Земли, а Dgш - по перпендикуляру к оси ее вращения. Спроектировав Vn на радиус, сложив обе составляющие и имея в виду, после некоторых преобразований получим выражение для суммарного центробежного ускорения (8) а далее - формулу для вычисления ускорения свободного падения ЛА, движущегося на данной высоте h, широте j, с постоянным ψ в гравитационном поле Земли, т.е. (9) Сравнив полученное выражение с выражением для ускорения свободного падения неподвижного объекта, получим, что изменение ускорения свободного падения только за счет путевой скорости Vn и высоты полета h равно (10) Черт. 1 На черт, 2, 3, 4 приведены графики, позволяющие оценить ускорение свободного падения и его составляющих для заданных условий полета от стандартного значения ускорения свободного падения, равного gc = 9,80665 м/с2. Черт. 2 Черт. 3 Черт. 4 Пользуясь приведенным материалом (см. черт. 1 - 4) ускорение свободного падения для заданных условий полета можно определить по формуле (11) Пример. Значение ускорения свободного падения ЛА при полете на высоте h = 5 км, широте j = 45°, с курсом ψ = 90° и скоростью полета М = 0,5 в безветрие будет равно: gла = 9,80665(1 + 0,025 - 0,035) = 9,80665 · 0,994 = 9,74781 м/с2, т.е. уменьшение ускорения свободного падения ЛА от gc составит 0,6 %. ПРИЛОЖЕНИЕ 6Рекомендуемое ТИПЫ И ОСНОВНЫЕ ХАРАКТЕРИСТИКИ ПОВЕРОЧНЫХ СРЕДСТВ
ПРИЛОЖЕНИЕ 7Обязательное
СОДЕРЖАНИЕ
|