EA039267B1 - Переходной участок закругления колеи рельсовой дороги - Google Patents

Переходной участок закругления колеи рельсовой дороги Download PDF

Info

Publication number
EA039267B1
EA039267B1 EA202100109A EA202100109A EA039267B1 EA 039267 B1 EA039267 B1 EA 039267B1 EA 202100109 A EA202100109 A EA 202100109A EA 202100109 A EA202100109 A EA 202100109A EA 039267 B1 EA039267 B1 EA 039267B1
Authority
EA
Eurasian Patent Office
Prior art keywords
track
axis
values
transition section
length
Prior art date
Application number
EA202100109A
Other languages
English (en)
Other versions
EA202100109A1 (ru
Inventor
Геннадий Викторович Величко
Original Assignee
Геннадий Викторович Величко
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Геннадий Викторович Величко filed Critical Геннадий Викторович Величко
Priority to EA202100109A priority Critical patent/EA039267B1/ru
Publication of EA202100109A1 publication Critical patent/EA202100109A1/ru
Publication of EA039267B1 publication Critical patent/EA039267B1/ru
Priority to PCT/BY2022/000003 priority patent/WO2022178612A1/ru
Priority to EP22758653.4A priority patent/EP4324982A1/en

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E01CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01BPERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
    • E01B5/00Rails; Guard rails; Distance-keeping means for them
    • E01B5/02Rails
    • E01B5/14Rails for special parts of the track, e.g. for curves
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61FRAIL VEHICLE SUSPENSIONS, e.g. UNDERFRAMES, BOGIES OR ARRANGEMENTS OF WHEEL AXLES; RAIL VEHICLES FOR USE ON TRACKS OF DIFFERENT WIDTH; PREVENTING DERAILING OF RAIL VEHICLES; WHEEL GUARDS, OBSTRUCTION REMOVERS OR THE LIKE FOR RAIL VEHICLES
    • B61F5/00Constructional details of bogies; Connections between bogies and vehicle underframes; Arrangements or devices for adjusting or allowing self-adjustment of wheel axles or bogies when rounding curves
    • B61F5/38Arrangements or devices for adjusting or allowing self- adjustment of wheel axles or bogies when rounding curves, e.g. sliding axles, swinging axles
    • EFIXED CONSTRUCTIONS
    • E01CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01BPERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
    • E01B2/00General structure of permanent way

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Train Traffic Observation, Control, And Security (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Machines For Laying And Maintaining Railways (AREA)

Abstract

Изобретение относится к устройству и эксплуатации переходного участка (далее по тексту ПУ) закругления рельсовой дороги, высокое качество поворотного движения экипажа по которому с заданной скоростью V при переменных по длине l кривизне оси его колеи k(l) и ее поперечного уклона i(l) обеспечивают новые закономерности его гармонизируемой формы. Форма оси колеи переходного участка определена местоположением N равноудаленных друг от друга точек с координатами x[n] и y[n], которые определены по формулам (1) и (2), переменные в которых, в свою очередь, определены по формулам (3)-(5). Поперечный уклон i[n] колеи в этих точках определен по формуле (6), переменные в которой, в свою очередь, определены по формулам (7)-(9). Наряду с традиционно управляемым параметром L (длина переходного участка) в данное изобретение включена дополнительная возможность варьирования значениями управляемых параметров Z и U, которые оказывают влияние на итоговые геометрические свойства нетождественных функций угла, кривизны и поперечного уклона колеи ПУ.

Description

Данное изобретение относится к устройству и эксплуатации переходного участка (далее по тексту ПУ) закругления рельсовой дороги, высокое качество поворотного движения экипажа по которому с заданной скоростью V при переменных по длине l кривизне оси его колеи k(l) и ее поперечного уклона i(l) обеспечивают новые закономерности его гармонизируемой формы.
План трассы колеи рельсовой дороги (рельсового пути) состоит из чередующейся последовательности прямолинейных участков и закруглений, вписанных в углы их поворота. Каждое закругление состоит из кругового и двух переходных участков длиной L. Ось колеи кругового участка описана дугой окружности с заданной постоянной кривизной K=1/R, а ее наружный рельс возвышен над внутренним также на постоянную высоту D. Это позволяет частично уменьшить поперечное ускорение, действующее на круговой части закругления при расчетной скорости V. Это уменьшение достигается за счет направленного к центру его кривизны поперечного уклона I=D/S. В пределах круговой части закругления этот уклон также постоянен, так как на ее протяжении он зависит только от постоянной величины возвышения D и ширины колеи S. Постоянство значений параметров V, R, D и S в пределах круговой части закругления колеи обуславливает постоянство теоретически максимального на всем протяжении значения непогашенного поперечного ускорения (далее по тексту НПУ) amax. Согласованность значений параметров V, R, D и S предусматривает, что между обусловленным ими значением и его допустимым значением адоп соблюден баланс, который задан условием amaxgоп. При этом НПУ amax вычисляется как результат векторного сложения действующих в горизонтальной плоскости гравитационных G и центробежных С ускорений amax=С-G·tan(α), где α - это угол поперечного наклона колеи, вычисляемый как a=arcsin(I).
В отличие от кругового участка закругления колеи ее переходной участок имеет более сложную форму. Она определяется видом функций кривизны k(χ) и поперечного уклона ί(χ), изменяющихся в диапазоне 0<k(χ)<K и 0<i(χ)<I в зависимости от относительной доли текущей длины колеи x=l/L, непрерывно изменяющейся в диапазоне 0 <χ<1. Функции большинства известных форм ПУ предшествующего уровня техники описаны в публикациях [1-4]. В каждой из них закономерность ее кривизны k(χ) и поперечного уклона ί(χ) задана единообразно по принципу k(χ)=K·f(χ) и ί(χ)=!0χ) c соблюдением требования 0<f(χ)<1. Этот принцип обусловил тождественность свойств графиков кривизны и поперечного уклона у подавляющего большинства известных форм ПУ предшествующего уровня техники. При этом основное различие между ними состоит в Gn порядке геометрической гладкости строго монотонных единичных функций f(χ). Он определяет максимальный порядок отличной от нуля n-й производной, которая в точках χ=0 и χ=1 имеет нулевые значения. У широко распространенной формы ПУ G0-го порядка гладкости функция f(χ)=χ. Поэтому графики закономерности ее кривизны k(χ) и поперечного уклона ί(χ) представлены прямыми линиями k(χ)=χ·K и i(χ)=χ·I. Закономерности относительных прямоугольных координат χ(χ) и y(χ) горизонтальной проекции оси колеи такой формы ПУ соответствуют только одной плоской кривой, определяемой в уровне техники как клотоида, clothoid или spiral.
У всех остальных форм ПУ Gn-го порядка гладкости с n>0 плавные очертания графиков закономерностей их кривизны k(χ) и поперечного уклона ί(χ) подобны первой половине графика синусоиды. Это определило общее для их форм название half-sine, т.е. полусинусоидальные. При одном и том же порядке гладкости полусинусоидальные формы ПУ имеют различные закономерности функции f(χ). Они различаются по плавности, которая оценивается по максимуму их первой производной df(χ)/dχ, вычисляемому в точке χ=0.5.
Устройство half-sine форм ПУ снижало остроту проблемы так называемого бокового рывка (lateral jerk), наблюдаемого в начале и в конце ПУ клотоидной формы. Однако в процессе их эксплуатации выявлялись проблемы качества движения, обусловленные низкочастотными колебаниями экипажей. На это указывают результаты исследований, проводимых в конце прошлого века в Японии на прямолинейных и криволинейных участках колеи скоростных рельсовых магистралей ПУ [5]. По мере повышения скоростей движения и развития высокоскоростных магистралей (далее ВСМ) актуальность проблем поворотного движения экипажа с переменной кривизной оси ПУ и с переменным поперечным уклоном его колеи возрастала. На это указывают современные тенденции к усложнению закономерностей half-sine форм ПУ таким образом, чтобы уменьшить значение максимумов первых производных их функций df(χ)/dχ в их центральной части [3]. В дополнение к этому применяются также методы т.н. приподнятого трассирования [2, 3], при котором закономерность кривизны k(χ)=K·f(χ) определяет не ось колеи ПУ, а траекторию движения расчетной точки экипажа, которая приподнята над уровнем верха головок рельс (далее ВГР) на высоту H.
Как правило, все множество этих мер направлено на решение основной проблемы функционирования системы Путь+Экипаж (далее по тексту ПЭ), которая на ПУ закруглений рельсовых дорог проявляется в неблагоприятной динамике силового взаимодействия ее элементов и в снижении уровня удобства движения пассажиров. Из-за этого повышается темп расстройства геометрии колеи ПУ и смежных с ним участков, а также преждевременно изнашиваются рельсы и колеса экипажей. Степень деструктивности этих процессов в значительной мере зависит от геометрических свойств форм ПУ. Результаты анализа кинематических показателей системы ПЭ на ПУ с традиционными и наиболее прогрессивными фор
- 1 039267 мами [6] указывают на то, что каждой из них присущи те или иные недостатки, которые не обеспечивают требуемого качества криволинейного движения экипажей. Однако даже в стандарте стран ЕС нет теоретически обоснованных приоритетов выбора и применения какой-либо одной из 5-ти рекомендуемых в нем форм ПУ, которая обеспечивала бы наиболее высокий уровень качества функционирования системы ПЭ.
Существующие проблемы обоснования и практического применения подавляющего числа форм ПУ предшествующего уровня техники обусловлены недостаточно адекватными и крайне упрощенными моделями системы ПЭ и методами оценки качества ее функционирования. Как правило, они основаны на показателях кинематики абстрактной точки, расположенной на уровне ВГР. В отличие от этого положения данного изобретения основаны на показателях кинематики точек, расположенных на функционально значимых уровнях экипажа. В зависимости от характера перевозок и приоритетных аспектов качества расчетные уровни, на которых они могут быть расположены, учитывают положение центра масс экипажа (далее ЦМ) и/или положение вестибулярного аппарата его пассажира.
Такой подход согласуется с целью обеспечения того качества криволинейного движения экипажа, под которыми традиционно подразумевают приемлемый уровень комфорта пассажиров и плавное изменение сил взаимодействия колес экипажа с рельсами колеи. Однако для его обеспечения с применением любой из форм ПУ предшествующего уровня техники предоставлялась единственная возможность варьирования значениями только их длин L. К тому же упрощенная модель системы ПЭ не позволяла объективно оценивать влияние этого фактора на качество ее функционирования на этих участках.
Проведенные на железных дорогах Японии исследования показали, что дискомфорт и неблагоприятные последствия для здоровья пассажиров обусловлены низкочастотными колебаниями поперечных ускорений с частотой ниже 1 Гц [5]. Как показали результаты расчетов, полученные с применением описанной в статье [6] многофакторной детерминированной кинематической модели (далее МДКМ), такой характер колебаний поперечных ускорений, действующих на расчетном уровне экипажа Н>0, в различной мере присущ всем half-sine формам ПУ предшествующего уровня техники. На это указывают значительные амплитуды осцилляций значений свойственных им функций НПУ а(l) и скоростей его изменения ψ=da/dt, переменных по длине l и времени t движения по ПУ. Эти же колебания приводят к неравномерному изменению сил реакции рельс, в результате чего возрастает динамическая составляющая действующей на них нагрузки.
Значения и закономерности этих показателей обусловлены сложным взаимодействием влияющих на них геометрических и физических факторов, наиболее важный результат которого проявляется и оценивается на расчетном уровне экипажа Н. Этим обусловлен интегративный характер этих показателей, который существенно отличается от локальных показателей кинематики абстрактной точки, традиционно рассчитываемых и оцениваемых ранее на уровне ВГР=0. Сложное взаимодействие свойств системы ПЭ учтено в ее многофакторной детерминированной кинематической модели (далее МДКМ), описанной в статье [6]. Рассчитываемые в МДКМ интегративные показатели качества ее функционирования зависят от значений следующих параметров:
скорости движения экипажа V;
радиуса круговой части закругления R;
расчетного возвышения наружного рельса над внутренним D;
расстояния между осями рельс колеи S;
возвышения H расчетной точки над уровнем ВГР;
длины ПУ L;
закономерности угла e(l) касательной к оси колеи и ее кривизны k(l);
закономерности отвода поперечного уклона колеи i(l);
закономерности кривизны проекции траектории движения расчетной точки kH (l) на горизонтальную плоскость;
значения булевой переменной устанавливаемого в зависимости от симметричного (true) или асимметричного (false) отвода возвышения наружного рельса над внутренним;
закономерности кривизны проекции оси колеи на вертикальную плоскость закономерности кривизны проекции продольной оси наружного (левого) рельса на вертикальную плоскость kL(l);
закономерности кривизны проекции продольной оси внутреннего (правого) рельса на вертикальную плоскость kR(l);
Как показал анализ результатов, полученных с применением МДКМ [6], варьирование в приемлемом диапазоне значениями только длин L ПУ предшествующего уровня техники не позволяет согласовывать взаимодействие всех перечисленных выше свойств их half-sine форм таким образом, чтобы устранить или уменьшить амплитуды осцилляций свойственных им функций НПУ а(l) и скоростей их изменения ψ=da/dt. Вследствие этого не представляется возможным также обеспечить равномерный и плавный рост разности сил реакции наружного и внутреннего рельсов колеи AF=FL-FR.
Характер этих трудно устранимых недостатков, которые присущи half-sine формам ПУ предшествующего уровня техники, иллюстрирует пример оценки интегральных показателей качества одной из
- 2 039267 наилучших форм ПУ [3], описанной функцией Order (3,7). Из перечня всех предложенных в [3] альтернатив она обладает наиболее высоким G4-m порядком гладкости. На графике 1-й производной функции отвода поперечного уклона колеи и кривизны ее оси присутствует характерное плато (см. фиг. 9). Это явно указывало на стремление автора этого решения устранить амплитуды осцилляций НПУ а(l) и уменьшить максимум постоянной скорости его изменения ψmax на значительном протяжении центральной части ПУ. Однако, как показал кинематический анализ функционирования системы ПЭ на ПУ такой формы, достижение этой цели приводит к существенной дисперсии этих показателей и других интегративных показателей качества движения в его начале и в конце (см. фиг. 10-13).
Таким образом, основными для оценки ПУ, полученных тем или иным методом, являются графики эпюр НПУ а(l) и скорости его изменения ψ(l)=da/dt, действующих на расчетном уровне Н. Квазилинейный характер эпюры НПУ а(l) является общепризнанным признаком высокого уровня качества конструкции ПУ, обеспечиваемого при условии гладкой и плавной эпюры скорости его изменения ψ(l)=da/dt на подходах к ее центральному платообразному участку. Из этого следует, что комфорт движения пассажиров и положительная динамика силового взаимодействия элементов СПЭ должна быть обеспечена гладким и плавным изменением эпюр а(l) и ψ(1), действующих на тех расчетных уровнях экипажа, которые критичны для показателей всех аспектов качества его движения. Приоритет комфорта пассажиров обуславливает необходимость линеаризации эпюры ψ(1) на большем протяжении центрального участка ПУ с минимальным отклонением постоянных максимальных значений ψmax ее эпюры от абсолютного, но ~Vp2_ β) практически недостижимого минимума J
По совокупности общих технических признаков наиболее близкими к заявляемой форме ПУ закругления колеи рельсовой дороги является форма ПУ, описанная в изобретении Railroad curve transition spiral design method based on control of vehicle banking motion закономерностью Order (3, 7) [3]. В рамках предложенного там способа проектирования переходной кривой железнодорожного пути осуществляют выбор одного из 18-ти математических выражений, определяющих значение ускорения (т.е. вторую производную) угла поперечного уклона (крена) пути в зависимости от переменного расстояния, изменяющегося от нуля до заданной длины переходной кривой. В зависимости от выбранного математического выражения определяются геометрические свойства формы оси колеи ПУ с длиной L и ее поперечного уклона. Из-за отсутствия методики и функционально обоснованных критериев выбора необходимого математического выражения, назначения длины L ПУ, а также необходимых в ряде случаев значений дополнительных параметров, получаемая в рамках указанного изобретения форма переходного участка не обеспечивает исключение известных из уровня техники недостатков.
Для выявления этих недостатков и поиска путей их устранения за счет новой формы ПУ с соответствующими этой цели геометрическими свойствами была применена более совершенная кинематическая модель криволинейного движения экипажа с переменной кривизной оси колеи и ее поперечного уклона [6].
Результаты анализа форм ПУ, известных из предшествующего уровня техники [6], позволили сформулировать задачу настоящего изобретения в установлении таких закономерностей формы колеи и длины L ПУ закругления рельсовой дороги, соблюдение которых при ее проектировании и строительстве обеспечит достижение наивысшего уровня качества функционирования системы ПЭ при заданных значениях следующих ее параметров: скорости движения V, радиуса закругления R, возвышении наружного рельса над внутренним D, высоты Н расположения расчетной точки экипажа над уровнем ВГ и расстояния между осями рельс S. Таким образом, технический результат заявляемого изобретения должен состоять в обеспечении наивысшего уровня качества функционирования системы ПЭ на ПУ закругления колеи рельсовой дороги при соблюдении расчетных значений ее параметров.
Поставленная задача решается, и указанные технические результаты достигаются заявляемым переходным участком закругления колеи рельсовой дороги, сопрягающим на длине L смежные с ним прямолинейный участок и круговую часть закругления с постоянным радиусом R, содержащим основание и формирующую колею рельсошпальную решетку, уложенную на нем в соответствии с координатами проекции точек оси переходного участка на горизонтальную плоскость с переменной по длине L величиной поперечного уклона i нормальной к оси линии, касательной к поверхности верха головок рельс, при этом форма оси колеи переходного участка определена местоположением N равноудаленных друг от друга точек с координатами x[n] и y[n], а поперечный уклон i[n] колеи в этих точках изменяется вдоль оси переходного участка в диапазоне 0<i[n]<D/S, где D - расчетное возвышение наружного рельса колеи над внутренним в пределах круговой части оси закругления рельсовой дороги, a S - расстояние между вертикальными осями поперечного сечения рельс колеи, по мере равномерного наращивания на постоянную величину Al=L/(N-1) расстояния l[n] до каждой ее n-й точки от l[0]=0.0 м в начале оси переходного участка и до l[N-1]=L в его конце. Поставленная задача решается, и указанные технические результаты достигаются за счет того, что относительные прямоугольные координаты точек проекции оси колеи x[n] и y[n] определены по формулам
- 3 039267 χ[0] = 0.0; χ[«] = χ[λ-1] + Δζ-cos^ (#) + Дд(£)) (1) τ[θ] = θ·θ: y[w] = ^E-l] + Az-sin(^(J) + A^(J)) (2) где n - порядковый номер точки оси переходного участка, изменяющийся в диапазоне l<n<(N-1), x[n] - координата, характеризующая расстояние от начала переходного участка до места проекции n-й точки оси на линию, касательную к ней в точках x[0] и y[0], y[n] - координата, характеризующая смещение в сторону центра кривизны оси ее n-й точки, измеряемое по нормали от линии, касательной к ней в точках x[0] и y[0],
Δ1 - постоянная величина приращения текущей длины l оси колеи переходного участка, измеряемой по ее дуге от ее начала до каждой n-й точки, δ - относительная доля текущей длины 1+Δ1/2 оси колеи переходного участка, измеряемой от ее начала до середины отрезка оси, заключенного между ее смежными точками n-1 и n, вычисляемая по фор муле βb(δ) - значение функции основного угла касательной к оси колеи переходного участка, вычисляемого в зависимости от относительной доли δ текущей длины 1+Δ1/2 до точки касания по формуле
Δβ(δ) - значение функции дополнения к основному углу касательной к оси переходного участка, вычисляемой в зависимости от относительной доли δ текущей длины 1+Δ1/2 до точки касания по формуле
где
U - параметр, определяющий величину дополнений к основным значениям угла касательной к оси колеи переходного участка в каждой ее точке с относительной долей δ текущей длины 1+Δ1/2.
Поперечный уклон i[n] колеи переходного участка в каждой n-й точке ее оси, изменяющейся в интервале 0<n<(N-1), определен по формуле /Η = 4(ζ)+Δ,(ζ) (6) где χ- относительная доля текущей длины l оси колеи до точки n, вычисляемая как <7’ ib(x) - значение функции основного поперечного уклона колеи переходного участка в точке n, вычисляемое в зависимости от относительной доли χ по формуле /,(/)=^/^37-/(102-/(140-/(172-/(198-/(154-/(66-12/))))))) (8)
Δi(χ) - значение функции дополнения к основному поперечному уклону колеи в точке n, вычисляемое в зависимости от относительной доли χ по формуле
Д|(/) = -2“/4 (^/(9-/(35-/(76-/(99-/(77-/(33-6/))))))) (9) где
Z - параметр, определяющий величину дополнений к основным значениям поперечного уклона колеи в каждой ее точке с относительной долей χ.
Востребованных практикой вариантов сочетания предопределенных значений параметров V, R, D, S и Н системы ПЭ может быть очень много. Это существенно усложняет задачу обеспечения стабильно высокого уровня качества функционирования системы ПЭ при существенном различии вариантов сочетания их значений. Как показали результаты исследований [7], достижению наилучшей согласованности взаимодействия элементов системы ПЭ способствует корректировка итоговых значений геометрических ΪΑ л ГА л ГА свойств формы ПУ w за счет учитываемых в них дополнений ХА Величина этих дополнений зависит от разницы между текущим вариантом сочетания проектных значений параметров V, R, D, S и Н и тем вариантом, который был учтен при обосновании функций основных (базовых) значений в ГА и i ГА свойств формы ПУ л Υи Ч
Поэтому наряду с традиционно управляемым параметром L в данное изобретение включена дополнительная возможность варьирования значениями управляемых параметров Z и U, которые оказывают влияние на итоговые геометрические свойства нетождественных функций угла, кривизны и поперечного уклона колеи ПУ. В каждой точке его оси, удаленной от ее начала на расстояние l, это влияние реализуется за счет зависимости ее координат x(l) и y(l) от функции угла касательной к ней состоящего из
- 4 039267 д - , вычисляемого по формуле (4) основного его значения и зависящего от параметра U дополнения к дД нему которое вычисляется по формуле (5), а также за счет зависимости соответствующего этой точке текущего значения превышения d(l) наружного рельса колеи над внутренним от функции поперечΡΊ · ГЙ ного уклона колеи VL) состоящего из вычисляемого по формуле (8) основного его значения \L) и зависящего от параметра Z дополнения к нему 33’ которое вычисляется по формуле (9).
Минимальные и максимальные степени членов полиномиальных функций (4), (5), (8) и (9), их коэффициенты, а также диапазоны варьирования значений управляемых параметров Z и U обоснованы с учетом обеспечения строгой монотонности G4-гладкой функции кривизны оси колеи, описываемой первой производной функции угла касательной к ней по длине ПУ \LJ < и нетождественной _ 3 . _ .
ей G -гладкой функции поперечного уклона колеи При x=l/L функция значений кривизны оси ко леи k(x) в каждой текущей координате длины l ПУ определена как
Ц/)=М/)+М*)( где ^/) = ^7^15-/(2043^67^-/^4085-/(1^5-/(185^-/^178-/(3900-600/^^
AU) = U^/(-z(74-^(340-z(885-z(1425-z(1452-j(914-^(325-50Z)))))^,(12)
U - параметр, определяющий максимальную величину дополнения к основному значению кривизны оси колеи переходного участка.
Такое математическое описание функций (4), (5), (8), (9), (11), (12) существенно расширяет возможности согласования взаимодействия всех свойств системы ПЭ в целях обеспечения стабильно высокого уровня качества ее функционирования при любом целесообразном варианте сочетания значений ее предопределенных V, R, D, H, S и управляемых L, Z, U параметров. Далее по тексту процесс достижения такой согласованности будет называться гармонизацией, получаемая в итоге форма ПУ будет называться гармонизированной, а соответствующие ей значения управляемых параметров L, Z и U будут называться оптимальными.
Квазилинейной G1 гладкой и плавной эпюре НПУ а(l), ассоциируемой с желаемым уровнем качества функционирования системы ПЭ на ПУ гармонизированной формы, должна соответствовать G0 гладкая и плавная эпюра скорости его изменения ψ(l)=da/dt. Соблюдение этого требования исключает риск так называемого бокового рывка (англ.: lateral jerk), наблюдаемого в начале и в конце ПУ. При этом максимальное значение этой скорости на существенном протяжении центрального участка ПУ должно быть постоянным. При идеальном соблюдении этих требований очертания эпюры функции ψ(1) будут подобны очертаниям равнобедренной трапеции со сглаженными углами и с достаточно протяженным плато в ее центральной части. Такие очертания в наибольшей мере способствуют наибольшему приближению линии плато с постоянными значениями ψmax к абсолютному, но практически недостижимому минимуму ψmin. Это соответствует общепринятому мнению о достижении наивысшего уровня качества движения.
Известным из предшествующего уровня техники half-sine формам ПУ присущи низкочастотные осцилляции интегративных показателей кинематики на расчетном уровне экипажей H. На это указывают результаты их оценки с применением МДКМ [6]. Величина и частота амплитуд этих осцилляций обусловлены несогласованным взаимодействием тождественных свойств нелинейных закономерностей отвода поперечного уклона колеи и кривизны ее оси. Устройство ПУ с заявляемой в данном изобретении формой позволяет существенно уменьшить амплитуды этих осцилляций до практически незначимого уровня. При этом постоянство и степень приближения ψmax к ψmin в пределах плато эпюры ψ(1) необходимо оценивать по дисперсии W скоростей изменения значений dψ/dt от соответствующей постоянству ψmax нулевой скорости. Из этого следует, что синергический эффект, достигаемый при разных вариантах сочетания значений предопределенных V, R, D, Н, S и управляемых L, Z, U параметров системы ПЭ, обратно пропорционален дисперсии W. С учетом этого цель поиска оптимальных значений параметров L, Z и U формализована как
где δb - относительная доля начала плато эпюры функционала ψ(1);
δe - относительная доля конца плато эпюры функционала ψ(1);
M - количество равномерно распределенных на интервале от δb до δe точек с относительными долями длины δ, в которых численно оценивают амплитуду осцилляций значений функционала dψ/dt.
По объективным причинам решение этой задачи не может компенсировать то негативное влияние,
- 5 039267 которое может оказать закономерность вертикальной кривизны оси ПУ и/или асимметричный (см. фиг. 16, 17) отвод возвышения D. Поэтому функцию цели (13) следует минимизировать без учета этих факторов, т.е. при гипотетически нулевой вертикальной кривизне оси ПУ и симметричном методе отвода (см. фиг. 14, 15) возвышения D. Как показывают аналитические зависимости МДКМ [6], при симметричном отводе возвышения D в два раза уменьшается вертикальная кривизна эпюр высот по наружному и внутреннему рельсам, а ее противоположные знаки исключают влияние этой кривизны на ускорения, действующие по вертикальной оси экипажа. С учетом этих условий значения функционала dψ/dt в каждой относительной координате δb<δ<δe длины ПУ можно вычислять по упрощенной формуле где kH - функция кривизны траектории расчетной точки экипажа на расчетном уровне H.
Переменные значения параметров L, Z, U и предопределенные значения параметров V, R, D, Н и S входят в уравнения функций, являющихся членами функционала (14). Это обеспечивает целенаправленность поиска оптимального варианта сочетания значений параметров L, Z и U, при котором дисперсии W амплитуд осцилляций, вычисленных по формуле (13), будут стремиться к минимуму. В соответствии с изложенными ранее обоснованиями это обеспечит наивысшее качество криволинейного движения экипажа, оцениваемое на уровне ЦМ или на каком-либо другом функционально значимом уровне экипажа.
В ходе проверки надежности и качества такой гармонизации форм ПУ было установлено, что при расчетных значениях НПУ amax>0.4 м/с2 минимизируемую дисперсию амплитуд осцилляций значений функционала dψ/dt следует оценивать в пределах плато эпюры ψ(1), ограниченного точками с относительными долями длины ПУ δb=1/4 и δe=3/4. При значениях НПУ в диапазоне 0.1<amax<0.4 м/с2 эти относительные доли границы следует вычислять по эмпирическим формулам л(15) ь 4 юЦ 0.35 ) V 0.35J J £ °-4~ЛЙ Ϊ(16) е 4 10Ц 0.35 ) ( 0.35) J
В ходе этой проверки также было установлено, что для закруглений, обеспечивающих так называемый сбалансированный режим движения экипажей при |amax|<0.1 м/с2, требуется существенное увеличение длин ПУ [7]. Ввиду этого целесообразность практической реализации таких конструкций закруглений является крайне сомнительной и далеко не всегда осуществимой. В основном - из-за ограничений по углу поворота или по другим причинам. Поэтому минимальное расчетное значение НПУ amax, при котором целесообразно гармонизировать форму ПУ, не должно быть меньше 0.1 м/с2.
Гармонизация формы ПУ может быть осуществлена путем ручного подбора комбинации близких к оптимуму значений L, Z и U. При этом следует иметь в виду, что при программной гармонизации форм ПУ с применением математически обоснованных методов минимум функционала (13) варьировал в диапазоне 1.OE-7<W<1.OE-4 м/с4. Эти значения W актуальны для расчетной скорости движения экипажа, варьирующей в диапазоне 100<V<400 км/ч, радиусов закругления R, обеспечивающих расчетные значения НПУ amax в диапазоне 0.1<amax<1.0 м/с2 при расчетном возвышении наружного рельса над внутренним, варьирующем в диапазоне 25<D<150 мм, а также при расчетном уровне экипажа Н=2200 мм и ширине колеи S=1520 мм. При других условиях диапазон варьирования минимума функционала (13) может быть иным.
Однако даже для ручной гармонизации форм ПУ крайне желательно автоматизировать рутинные вычисления и построение эпюр, позволяющих контролировать эффективность этого процесса и качество его результатов. В идеале каждый шаг уточнения значений параметров L, Z или U должен приводить к минимизации численных значений функции цели (13), а также делать все более близкими к прямой эпюру НПУ a(l) и участок плато эпюры ψ(1). В ходе работы над данным изобретением достаточно удобной и эффективной оказалась реализация этого процесса в программе Microsoft Excel. Наряду с численной оценкой степени гармонизации формы ПУ ее качество подтверждалось эпюрами всех геометрических и функциональных свойств системы ПЭ. Также в этой программе рассчитывались координаты и другие параметры, необходимые для разбивки закруглений с гармонизированными формами ПУ, соответствующими положениям данного изобретения.
Таким образом, в предпочтительных формах реализации заявляемого переходного участка при заданной величине параметра U, варьируемой в диапазоне от -45 до +55, значение суммы функций основ„ „ о о о о β β Ω ΔΑςΊ о ω Т _ ного угла касательной к оси колеи W и дополнений к нему на всем протяжении длины L переходного участка строго монотонно и непрерывно изменяется от 0 до L/2R, значение ее первой производк\ - | = 4β\- | + δ/- 11 dl, ной по длине l, определяющей закономерность кривизны оси переходного участка Л/ I Л/ \L))I , также строго монотонно и непрерывно изменяется от 0 до 1/R, а для суммы всех последующих 2-х, 3-х,
- 6 039267
4-х и 5-х производных этих функций по этой же длине l обеспечивается непрерывность изменения их значений, а также равенство их нулю на обеих концах интервала 0<l<L. При этом функция кривизны оси колеи k(x) определяется по приведенным выше формулам (10)-(12).
Также предпочтительными являются формы реализации заявляемого переходного участка, в которых при заданной величине параметра Z, варьируемой в диапазоне от -88 до +33, значение суммы функ“ т ций основного поперечного уклона колеи и дополнений к нему на всем протяжении длины L переходного участка строго монотонно и непрерывно изменяются от 0 до D/S, а значения ее 1-х, 2-х и 3-х производных по длине l непрерывно изменяются при соблюдении равенства их нулю на обеих концах интервала 0<l<L.
Следует отметить, что для применения заявляемого изобретения в современных технологиях автоматизированного проектирования рельсовых дорог более эффективна программная гармонизация форм ПУ с применением численного метода Ньютона. Основные положения необходимого для этого алгоритма и результаты его работы, иллюстрирующие преимущества заявляемого переходного участка закругления колеи рельсовой дороги будут рассмотрены далее более подробно на примерах некоторых предпочтительных, но не ограничивающих форм реализации со ссылками на позиции фигур чертежей, на которых схематично представлены фиг. 1 - план закругления колеи рельсовой дороги;
фиг. 2 - схемы расчета действия гравитационных ускорений (G) на уровне ВГР криволинейного участка закругления колеи рельсовой дороги;
фиг. 3 - схемы расчета действия центробежных ускорений (С) на уровне ВГР криволинейного участка закругления колеи рельсовой дороги;
фиг. 4 - гравитационные составляющие ускорений и сил взаимодействия элементов системы ПЭ;
фиг. 5 - центробежные составляющие ускорений и сил взаимодействия элементов системы ПЭ;
фиг. 6 и 7 - суммарные векторы сил реакции рельс FL и FR, определяемые с учетом направления вектора НПУ amax и силы FC;
фиг. 8 - координатный метод описания кинематики расчетной точки экипажа М с учетом геометрических свойств формы ПУ;
фиг. 9 - нормализованные графики функции Order (3, 7) и ее производных;
фиг. 10 - графики кривизны k(l) оси колеи и кривизны kH(l) траекторий движения расчетных точек;
фиг. 11 - графики непогашенных поперечных ускорений а(l), действующих на разных уровнях H расчетного экипажа;
фиг. 12 - графики скоростей изменения непогашенных поперечных ускорений а(l) ψ=da/dt;
фиг. 13 - графики сил реакции рельс и их разности;
фиг. 14 - поперечное сечение колеи ПУ при симметричном отводе возвышения D;
фиг. 15 - диаграммы превышений рельс над продольным профилем оси ПУ при симметричном отводе возвышения D;
фиг. 16 - поперечное сечение колеи ПУ при асимметричном отводе возвышения D;
фиг. 17 - диаграммы превышений наружного рельса и обусловленных ими превышений оси ПУ над проектным профилем внутреннего рельса при асимметричном отводе возвышения D;
фиг. 18 - расчетная схема унифицированного метода вычисления прямоугольных координат n-й точки оси колеи ПУ в локальной системе тангенса закругления;
фиг. 19 - графическое представление процесса минимизации функции цели гармонизируемой формы ПУ в соответствии с изобретением;
фиг. 20-24 - пример геометрических и функциональных свойств системы ПЭ, обеспечиваемых гармонизированной формой ПУ по изобретению;
фиг. 25 - пример зависимости оптимальных длин L гармонизированных форм ПУ от НПУ amax, расчетной скорости V и возвышений D, а также зависимости минимальных длин L клотоидных и рекомендованных [4] half-sine форм ПУ от аналогичной V скорости при D=150 мм и допустимой скорости подъема колеса по возвышению наружного рельса λ=1/10 км/ч;
фиг. 26 - пример инвариантной к скорости V зависимости максимальных значений ψmax гармонизированных форм ПУ от НПУ amax и возвышений D;
фиг. 27 - пример инвариантной к скорости V зависимости сдвижек р круговых участков осей закруглений гармонизированных форм ПУ от НПУ amax и возвышений D.
План трассы рельсовой колеи состоит из чередующейся последовательности прямолинейных участков 1 и закруглений, вписанных в углы их поворота. Каждое закругление (схематично представлено на фиг. 1) состоит из двух переходных участков (ПУ) 2 и расположенного между ними кругового участка 3 рельсовой колеи заданной кривизны K=1/R. При известной ширине колеи S радиус его оси R согласован с расчетными значениями: скорости движения экипажей V, возвышения наружного рельса над внутренним D и величины НПУ amax, действующего на уровне ВГР. Ось ПУ 2 обозначена цифрой 4.
На фиг. 2 и 3 представлены схемы расчета действия гравитационных ускорений (G) и центробеж
- 7 039267 ных ускорений (C), действующих на экипаж на уровне ВГР криволинейного участка закругления колеи рельсовой дороги. Эти схемы соответствуют традиционно применяемым формулам, оценивающим действие векторов гравитационного G=9.81 м/с2 (фиг. 2) и центробежного С=kv2 м/с2 (фиг. 3) ускорений на участках закруглений колеи с постоянной или с переменной кривизной без учета возвышения ЦМ экипажа или вестибулярного аппарата пассажира над уровнем ВГР. На фиг. 3 обозначены позицией 5 называемый далее также как левый наружный и позицией 6 называемый далее также как правый внутренний рельсы.
На фиг. 4 схематично изображены векторы ускорений GL и GR, действующие по нормали к оси колесной пары экипажа 7, значения которых зависят от ускорения силы земного притяжения g=9.81 м/с2, угла α поперечного уклона колеи i и отношения значения параметра Н к S, а также от вертикальных ускорений, обусловленных скоростью движения V, вертикальной кривизной оси kV, а также кривизной kL и kR эпюр отвода возвышения D по левому и правому рельсам колеи (гравитационные составляющие).
На фиг. 5 схематично изображено графическое представление векторов, действующих по нормали к оси колесной пары ускорений CL и CR, значения которых зависят от кривизны траектории расчетной точки экипажа 7 и скорости движения V, а также от угла α поперечного уклона колеи i и отношения значения параметра Н к S (центробежные составляющие).
Фиг. 6 и 7 схематично изображены суммарные векторы силы реакции FL левого 5 (наружного) и силы реакции FR правого 6 (внутреннего) рельсов, определяемые с учетом направления вектора НПУ amax, а также зависящего от этого центростремительного или центробежного направления вектора силы поперечной реакции FC левого (наружного) или правого (внутреннего) рельсов.
На фиг. 8 представлена расчетная схема учета геометрических свойств формы колеи ПУ, а также возвышения Н точки М над уровнем ВГР в координатном методе описания ее кинематики.
На фиг. 9-13 схематично представлены графики нормализованных геометрических свойств наиболее близкой к заявляемой форме ПУ функции Order (3, 7) [3], а также присущие ей амплитуды осцилляций интегративных показателей качества функционирования системы ПЭ, которые представлены графиками: кривизны оси ее колеи и траекторий движения расчетных точек на разных уровнях экипажа (фиг. 10), НПУ а(1) (фиг. 11) и скорости их изменения da/dt (фиг. 12), действующих на этих же уровнях экипажа, а также силы реакции левого FL и правого FR рельсов и их разности (фиг. 13). Все показатели, положенные в основу этих графиков, рассчитаны по зависимостям МДКМ [6] при V=400 км/ч, R=8000 м, D=150 мм, S=1520 мм и L=420 м на расчетном уровне H=2200 мм.
На фиг. 14 схематично представлено поперечное сечение колеи ПУ и соответствующие ему halfsine закономерности эпюр превышений головок рельс над линией продольного профилям оси ПУ при их симметричном отводе на величину ±D/2 (см. фиг. 15).
На фиг. 16 схематично представлено поперечное сечение колеи ПУ и соответствующая ему half-sine закономерность эпюры превышения головки наружного (левого) рельса над линией продольного верха головки внутреннего (правого) рельса при асимметричном его отводе на величину D, а также закономерность эпюры дополнительных превышений оси колеи над ее проектным положением (см. позицию 8), обусловленных данным методом отвода (см. фиг. 17).
На фиг. 18 представлена расчетная схема унифицированного метода вычисления локальных прямоугольных координат n-й точки оси колеи ПУ любой формы при известной функции угла касательной к ней β(δ). На чертеже позициями обозначены ось 9 прямолинейного участка колеи, ось 4 переходного участка закругления колеи и ось 10 кругового участка закругления колеи радиуса R.
На фиг. 19 приведено графическое и аналитическое представление функции цели гармонизации формы ПУ, геометрические свойства которой описаны закономерностями (формулами) (1)-(12). Пунктирными линиями показан характер закономерностей эпюр значений dψ/dt функционала (14) и их осцилляции, дисперсия амплитуд которых последовательно уменьшаются в ходе гармонизации.
На фиг. 20-24 в графическом виде представлен пример геометрических и функциональных свойств системы ПЭ, обеспечиваемых гармонизированной формой ПУ по изобретению (в одной из возможных, но не ограничивающих форм реализации) при оптимальных значениях L=420 м, Z=0.64, U=-0.38 и заданных значениях параметров V=400 км/ч, R=8000 м, D=150 мм, S=1520 мм и H=2200 мм.
На фиг. 25 представлены зависимости оптимальных длин L ПУ закругления колеи S=1520 мм от НПУ amax, полученные в результате их гармонизации при фиксированных значениях Н=2200 мм для 240 вариантов сочетания выборки значений предопределенных параметров V, R и D. Эта выборка включала 4 дискретных значения скорости, варьируемой с шагом 100 км/ч в диапазоне 100<V<400 км/ч, для каждого из которых устанавливалось 6 дискретных значений возвышения наружного рельса над внутренним, варьируемых с шагом 25 мм в диапазоне 25<D<150 мм. Для каждого из образуемых таким образом 24 вариантов сочетания разных значений пар параметров V и D вычислялось 10 значений радиусов R, обеспечивающих соответствующее число дискретных значений НПУ amax, варьируемых с шагом 0.1 м/с2 в диапазоне 0.1<amax<1.0 м/с2. При отображении точек с координатами amax L масштаб оси абсцисс и диапазон отображаемых на ней оптимальных длин L для каждого из 3 значений скорости V=300 км/ч, V=200 км/ч и V=100 км/ч согласовывался с абсциссами оптимальных длин L, вычисленными при расчетной
- 8 039267 скорости V=400 км/ч. В результате этого положение всех 180 точек с одинаковыми ординатами amax и разными абсциссами L, вычисленными при расчетных скоростях V=300 км/ч, V=200 км/ч и V=100 км/ч, совпадало с положением 60 точек с такими же ординатами amax и абсциссами L, вычисленными при расчетной скорости V=400 км/ч. Это позволило представить результаты гармонизации форм ПУ, полученные по всему объему выборки, лишь 6-ю линиями одного графика гипотетически степенных зависимостей оптимальных длин L ПУ от подмножества вариантов сочетания значений параметров V, R, D, Н и S, в каждом из которых радиус R согласован с параметрами V, D и S по расчетной величине НПУ закругления amax, варьируемой в диапазоне 0.1<amax<1.0 м/с2.
Аналогичным образом на фиг. 25 представлена вертикальная линия графика, традиционно ограничивающего минимальные длины L клотоидных форм ПУ по максимально допустимой скорости подъема колеса λ=1/10 км/ч. Также на фиг. 25 представлены линии графиков, которые ограничивают минимальные длины некоторых типов half-sine форм ПУ с такими же значениями параметров V, R и D=150 мм. Эти длины определены с учетом рекомендованных в [4] коэффициентов. При этих параметрах оптимальные длины гармонизированных форм ПУ совпадают с минимальными длинами клотоидных и half-sine форм ПУ лишь при одном из значений НПУ amax, варьируемом в данном случае в диапазоне 0.1<amax<0.3 м/с2. При других расчетных значениях НПУ amax эти длины существенно разнятся. Вполне очевидно, что допускаемое в разных странах варьирование приемлемой скорости подъема колеса в диапазоне ίο - Л - ю [4, 8, 9, 10] не устранит принципиальные отличия в представленных на фиг. 25 зависимостях. Поэтому положения практической реализации данного изобретения предусматривают устройство гармонизированных форм ПУ исключительно только с соответствующими их параметрам оптимальными длинами L.
На фиг. 26 и 27 представлены примеры зависимостей максимальных значений ψmax (фиг. 26) и сдвижек p круговых участков осей закруглений (фиг. 27) от расчетных значений НПУ amax и возвышений D, соответствующих оптимальным значениям параметров L, Z и U тех же 240 гармонизированных форм ПУ по изобретению (в одной из возможных, но не ограничивающих форм реализации), учтенных при построении графиков фиг. 25. В отличие от представленных на фиг. 25 графиков зависимостей оптимальных длин L зависимости максимальных значений ψmax (фиг. 26) и сдвижек р круговых участков осей закруглений (фиг. 27) в данной форме представления инвариантны к расчетной скорости V.
Геометрия заявляемого ПУ закругления колеи рельсовой дороги определяется, и заявляемый ПУ функционирует следующим образом.
Каждое закругление колеи рельсовой дороги состоит из кругового 3 и двух переходных участков (ПУ) 2, связывающих его с соответствующими прямолинейными участками 1. Ось 4 кругового участка колеи описана дугой заданной кривизны K=1/R. При известной ширине колеи S ее радиус R назначают в зависимости от расчетных значений: скорости движения экипажей V, возвышения D наружного рельса 5 над внутренним 6 и величины НПУ amax, действующего на функционально значимом уровне расчетного экипажа H. В соответствии с представленной на фиг. 8 расчетной схемой кинематической модели [6] радиус кривизны горизонтальной проекции траектории точки М, возвышенной над уровнем ВГР на величину Н, на круговом участке траектории с постоянным радиусом кривизны оси колеи R отличается от № него на весьма малую, в сравнении с ним, величину s ’. Поэтому максимальное значение НПУ amax на круговом участке колеи вычисляют с достаточной для практических целей точностью по формуле — , (17)
О
При фиксированном значении расчетной скорости V требуемый баланс постоянных значений amaxдоп на этом участке закругления обеспечивается за счет подбора необходимого сочетания значений параметров R и D, допустимых по условиям проектирования.
В отличие от этого форма заявляемого ПУ закругления колеи рельсовой дороги, предшествующего круговому участку, определяется с учетом двух геометрических и одного функционального требований: соблюдения строго монотонного и плавного изменения поперечного уклона колеи i(l) и кривизны k(l) ее оси на всем протяжении участка изменения значений этих параметров от нулевых значений на прямолинейном участке колеи до максимальных значений на круговом участке закругления;
соответствия геометрических свойств формы ПУ теории изгиба балки на упругом основании;
согласования длины и закономерностей свойств изгиба рельс колеи ПУ в горизонтальной и вертикальной плоскостях проекции в целях достижения наивысшего уровня качества функционирования системы ПЭ.
Соблюдение геометрических требований обеспечивается за счет заявляемых свойств функций изменения угла дуги оси колеи, ее кривизны и ее поперечного уклона. Описываемый ими изгиб рельс обеспечивает G3-й порядок геометрической гладкости в вертикальной плоскости проекции и G6-й - в горизонтальной плоскости проекции.
Соблюдение функционального требования обеспечивается за счет вариативности заявляемых свойств нетождественных функций изменения кривизны оси колеи и ее поперечного уклона. Варьирование длины ПУ и свойств его формы влияет на показатели кинематики расчетной точки на функционально значимом уровне экипажа, возвышенном над уровнем ВГР на величину Н>0. Интегративные показа
- 9 039267 тели качества этого процесса включены в формализованную цель процедуры согласования значений управляемых и предопределенных параметров системы ПЭ, обеспечивающих наибольшую эффективность их взаимодействия. Эта цель достигается при оптимальных значениях управляемых параметров L, Z и U, вычисляемых при заданных значениях параметров V, R, D, Н и S в следующей последовательности:
1. С учетом вычисленного по формуле (17) НПУ amax по формулам (15) и (16) вычисляют соответствующие ему относительные доли δb и 5е длины участка эпюры функционала (14), в пределах которого необходимо оценивать и в соответствии с функцией цели (13) минимизировать дисперсию амплитуд осцилляций W его значений.
2. Устанавливают начальные значения искомых параметров L=20 м, Z=0 и U=0.
3. К оптимальным значения этих параметров приближают в ходе 5-ти циклов, в пределах каждого из которых, в течении 20-ти итераций методом Ньютона последовательно приближают сначала значение L при фиксированных текущих значениях Z и U, затем значение Z при фиксированных текущих значениях L и U и, наконец, значение U при фиксированных текущих значениях L и Z. При этом новые значение каждого из параметров, устанавливаемое в конце каждого 20-го итерационного цикла его приближения, которое далее будет учитываться в процессе приближения значения последующих параметров, выбирают из перечня их корректных значений, обеспечивающих по условию (13) наименьшую дисперсию W;
необходимые для вычисления W значения функционала (14) dψ/dt и его частные производные следует рассчитывать с шагом Δδ=(δeb)/(M-1) при М>1000;
некорректными являются результаты с отрицательной или с превышающей допустимую для данного закругления длиной L;
со значениями параметра Z, выходящими в процессе расчета за пределы допустимых значений 88<Z<33:
со значениями параметра U, выходящими в процессе расчета за пределы его допустимых значений 45<U<55;
кривизну kH траектории расчетной точки экипажа на уровне H, а также ее производные, необходимые для расчета текущих значений функционала (14), следует вычислять с учетом описанных в [6] детерминированных зависимостей МДКМ;
каждое из претендующих на новое приближенное значение параметра Ln+1, Zn+1 и Un+1 вычисляют по формулам dt J Д 4 dt 8t5Ln J уе (/ί6. ye . дУ(<Д
5=¾ I dt J / ( s=gb dt dtdZn J (19) (20) соблюдение этих рекомендаций позволяет минимизировать функцию цели (13) при W<1.OE-O3 или с гораздо лучшим результатом.
4. Вычисленную таким образом оптимальную длину L ПУ округляют до ближайшего большего целого значения, кратного 1 м. Более грубое округление, например, кратно 5, 10 или до 20 и более метров крайне нежелательно. Особенно в меньшую сторону, т.к. приводит к существенному увеличению дисперсии амплитуд осцилляций функции dψ/dt. В конечном итоге это может привести к вибрации, толчкам или раскачиванию кузова экипажа в процессе эксплуатации ПУ при проектных значениях параметров V, R, D, H и S.
5. Значения параметров Z и U могут быть округлены с погрешностью ±0.01.
Соответствующие этим требованиям значения параметров L, Z и U являются гармонизированными со значениями других предопределенных параметров закругления колеи и расчетного экипажа V, R, D, H и S. Они будут актуальны и могут многократно применяться при устройстве ПУ заявляемой формы в других закруглениях рельсовой колеи с аналогичным или с другим углом поворота ±θ при соблюдении условия |θ|>L/R и постоянстве расчетных значений предопределенных параметров V, R, D, H и S. B случае изменения хотя бы одного из них вычисления по пп.1-5 описанной выше процедуры гармонизации значений параметров L, Z и U следует повторить.
При актуальных значениях параметров L, Z и U относительные прямоугольные координаты оси заявляемого ПУ х[п] и у[п] следует вычислять по формулам (1) и (2) при N=L40+1. При соблюдении этого условия разница между длиной дуги Δ1=Ο.1 м и длиной стягивающей ее хорды от точки n-1 и до точки n будет весьма незначительной. Это обеспечит достаточную для практических целей детальность и точность расчета как абсолютных координат оси закругления в установленной на объекте системе, так и других его геометрических характеристик. При этом значения тангенсов, биссектрисы, центра круговой кривой и ее сдвижки, необходимые для устройства закругления колеи с ПУ заявляемой формы, вычис
- 10 039267 ляются по стандартной методике и по общеизвестным формулам, применимым к закруглению с любой другой формой ПУ при известных локальных координатах x[N-1] и y[N-1] точек в конце оси колеи каждого из ПУ и угле β всей дуги его оси.
При традиционной технологии поточечной разбивки оси ПУ заявляемой формы на местности с шагом, большим чем Al, можно использовать координаты только тех точек, кратность номеров которых соответствует этим требованиям. При устройстве колеи ПУ заявляемой формы с применением автоматизированных систем управления рихтово-подбивочными машинами, а также при формировании в современных информационных технологиях их 3D моделей, целесообразно определять местоположение оси колеи ПУ с достаточно подробным дециметровым шагом. Для этого достаточно учесть координаты всего числа N, рассчитанных по формулам (1) и (2) точек.
Конструктивные и функциональные преимущества закруглений колеи с ПУ заявляемой формы подтверждены результатами гармонизации 240 вариантов сочетания значений предопределенных параметров закруглений колеи шириной S=1520 мм, сформированных при 4-х дискретных значениях скорости V, варьируемой с шагом 100 км/ч в диапазоне 100<V<400 км/ч, для каждого из которых устанавливались 6 дискретных значений возвышения наружного рельса над внутренним, варьируемых с шагом 25 мм в диапазоне 25<D<150 мм. Для каждого из 24 вариантов сочетания разных значений пар параметров V и D вычислялось 10 значений радиусов R, обеспечивающих соответствующее число дискретных значений НПУ amax, варьируемых с шагом 0.1 м/с2 в диапазоне 0.1<amax<1.0 м/с2.
Оптимальные значения L, Z и U всех 240 вариантов сочетания значений V, R и D с фиксированным значением возвышения функционально значимого уровня расчетного экипажа Н=2200 мм были вычислены при 1.OE-O7<W<4.OE-O4. Такая степень минимизации функции цели (13) обеспечила высокий уровень качества функционирования системы ПЭ, подтверждаемый G1-гладкими квазилинейными эпюрами НПУ а(Х) (см. фиг. 22) и G0-гладкими трапециевидными эпюрами скорости его изменения ψ=da/dt (см. фиг. 23). При этом отклонения максимальных значений этой скорости ψmax от абсолютного, но практически недостижимого минимума ψmln были существенно меньше, чем у известных форм ПУ предшествующего уровня техники с аналогичными значениями параметров V, R, D и L. Наряду со значениями параметра ψmax результаты гармонизации форм ПУ дополнялись значениями так называемой сдвижки круговой кривой р, по значению которой в комплексе со значениями R и L можно судить о конструктивных преимуществах заявляемой формы ПУ, имеющих важное для трассирования рельсовых дорог значение.
Часть результатов, полученных при гармонизации форм ПУ с расчетным возвышением D=150 мм для расчетных скоростей V=400 км/ч и V=100 км/ч, приведена в таблице. Они указывают на наличие устойчивой их зависимости от расчетных значений НПУ amax. Более наглядно характер и особенности этой зависимости иллюстрируют графики, представленные на фиг. 25-27.
Пример оптимальных значений параметров L, Z и U заявляемых форм ПУ с D=150 мм, гармонизированных для расчетных значений НПУ, варьируемых в диапазоне 0.1<amax<1.0 м/с2 при V=400 км/ч (см. числитель) и V=100 км/ч (см. знаменатель) для H=2200 мм и S=1520 мм
«max, М/С2 R, м L, м Z и ^,м/с3 р, м
0.1 12110/760 930/233 -3.52/-4.52 6.52/8.52 0.014/0.014 0.932 / 0.960
0.2 11030/690 686/ 171 -0.04/-0.31 -0.il/-0.ll 0.035/0.035 0.561 /0.558
0.3 10120/630 580/ 145 0.43 / 0.25 -0.62 / -0.62 0,063 / 0.063 0.437 / 0.436
0.4 9360 / 585 513/128 0.26 / -0.04 -0.15/-0.32 0,095 / 0.095 0.370/0.368
0.5 8700 / 545 464/ 116 0.40/0.16 -0.45 / -0.33 0.132/0.131 0.325/0.325
0.6 8125/510 426/ 106 0.01/-0.26 -0.42/-0.18 0.172 /0.170 0.294/0.295
0.7 7620 / 480 397 / 99 -0.03 / -0.30 -0.18/-0.22 0.216/0.217 0.272/0.271
0.8 7180/450 373 / 93 0.18/-0.53 -0.40 / -0.38 0.263 / 0.26 0.255 / 0.253
0.9 6785 / 425 353/88 0.31/-0.20 -0.59/-0.71 0.313/0.314 0.241 /0.240
1.0 6430 / 400 336 / 84 0.34 / 0.28 -0.43 / -0.39 0.366 / 0.366 0.231 /0.230
Из анализа графиков, представленных на фиг. 25, следует, что степенные зависимости оптимальных длин Lopt гармонизированных ПУ заявляемой формы существенно отличаются от вертикальной линии графика зависимости минимальных длин Lmin клотоидных форм ПУ, традиционно нормируемых по так называемой допустимой скорости подъема наружного колеса по возвышению λ=V·D/L. В данном примере она ограничивает минимальные длины Lmin клотоидных форм ПУ закруглений с расчетным возвышением D=1500 мм при рекомендованной в [8, 9, 10] допустимой скорости подъема наружного колеса λ=1/10 км/ч (—28 мм/с). В случае применения рекомендованных в [4] half-sine форм ПУ их минимальные длины будут еще больше отличаться от оптимальных длин гармонизованных ПУ заявляемой формы. Исходя из явной тенденции к увеличению коэффициентов пропорциональности, коррелирующих с порядком геометрической гладкости соответствующих им типов half-sine форм ПУ, минимальные длины более совершенных и более гладких half-sine форм ПУ могут быть еще больше. Это указывает на явные
- 11 039267 проблемы эмпирических методов нормирования длин ПУ по любому из допустимых значений параметра λ или отношения D/L, не опосредованная причинно-следственная связь которых с интегративными показателями качества движения не обоснована.
Из анализа графиков, представленных на фиг. 26, следует, что гармонизированным ПУ заявляемой формы присущи полиномиальные зависимости максимальной скорости изменения ψmax от расчетных значений НПУ amax, которые в отличии от оптимальных значений длин инвариантны по отношению к расчетной скорости движения V. При этом соответствующий показателю ψmax проектный уровень комфорта движения на гармонизированном ПУ заявляемой формы будет тем выше, чем больше будет расчетное возвышение D. Как следует из графика функции ψmax при D=150 мм, ее значения во всем диапазоне варьирования НПУ 0.1<amax<1.0 м/с2 не превышают даже нижней границы диапазона традиционно рекомендуемых норм 0.4<ymax<0.6 м/с3.
Из анализа графиков, представленных на фиг. 27, следует, что сдвижкам р круговых кривых закруглений с гармонизированными ПУ заявляемой формы также присущи инвариантные по отношению к расчетной скорости движения V полиномиальные зависимости от расчетных значений НПУ amax. При этом значения сдвижек р круговых кривых закруглений с гармонизированными ПУ заявляемой формы в разы меньше сдвижек р круговых кривых весьма распространенных и проблемных закруглений с клотоидными формами ПУ. Как показывают расчеты, кратная разница этих сдвижек способствует переустройству закруглений с клотоидными формами ПУ на закругления с гармонизированными ПУ заявляемой формы в целях повышения скорости и удобства движения с наименьшими параметрами рихтовки существующего пути.
Учет приведенных закономерностей может способствовать принятию оптимальных решений при устройстве новых или реконструкции существующих закруглений рельсовых дорог. Из-за большого числа возможных вариантов сочетания предопределенных значений параметров V, R, D, Н и S, варьирующих в достаточно широком диапазоне, традиционное представление предварительно рассчитанных оптимальных значений L, Z и U гармонизированных форм ПУ в виде таблиц или графиков не представляется целесообразным. Для современных автоматизированных технологий проектирования, строительства и эксплуатации рельсовых дорог более приемлема программная реализация алгоритмов гармонизации форм ПУ и расчета необходимых для их устройства координат в соответствии с положениями данного изобретения.
В отличие от решений предшествующего уровня техники, объективные законы физики и математически и обоснованные методы гармонизации, положенные в основу данного изобретения, обеспечивают поиск оптимальных значений всех необходимых параметров конструкции переходных участков закруглений, геометрические свойства формы которых описаны закономерностями (1)-(12). Они прямо и непосредственно, а не косвенно и опосредовано, обеспечивают интегративное качество функционирования системы Путь+Экипаж. При этом уровень этого качества зависит от согласованности и соизмеримости таких ее параметров, как возвышение расчетной точки экипажа Н, скорость его движения V, радиус закругления R, возвышение наружного рельса D и метода его отвода. Поэтому положения и отличительные особенности данного изобретения следует учитывать также еще и на этапе принятия решений об основных значениях параметров закруглений проектируемых трасс рельсовых магистралей.
Источники информации.
1. Bjorn Kufver, VTI rapport 420A, “Mathematical description of railway alignments and some preliminary comparative studies”, Swedish National Road and Transport Research Institute, 1997. Digitala Vetenskapliga Arkivet (цифровой научный архив). [Электронный ресурс] -17 июля 2020. - Режим доступа:
http://www.diva-portal.org/smash/record. jsf?pid=::diva2%3A675179&dswid=4876.
2. Патент ЕР № 1523597В1, опубл. 16.07.2008 г.
3. Klauder, Louis, Т., Jr. Railroad curve transition spiral design method based on control of vehicle banking motion. Available from:
https://patentscope.wipo.int/search/ru/detail.jsf?docId=WQ2001098938&recNum=l&maxRec= &office=&prevFilter=&sortOption=&queryString=&tab=PCTDescription
4. EN 13803-1:2010: Railway applications -Track-Track alignment design parameters - Track gauges
1435 mm and wider - Part 1 :Plain line [Required by Directive 2008/57/EC]
5. M. Ueno et al.: Motion Sickness Caused by High Curve Speed Railway Vehicles. Jpn J Ind Health
1986; 28:266-274.
6. Velichko, G. 2020. Quality analysis and evaluation technique of railway track + vehicle system performance at railway transition sections with various shape curves, In Transport Means 2020: Proceedings of the 24th International Scientific Conference, Part II: 573-578. Available from: https.7/transportmeans.ktu.edu/wp-content/uploads/sites/307/2018/02/Transport-meansA4-lI-dalis.pdf
7. Velichko, G. 2020. Shape Harmonization of the Railway Track Transition Section & the Kinematics of Vehicle Body Design Point, In Transport Means 2020: Proceedings of the 24th International Scientific Conference, Part II: 910-915. Available from: https://transportmeans.ktu.edu/wp-content/uploads/sites/307/2018/02/Transport-means-A4-lldalis.pdf
8. Шахунянц Г. M. Железнодорожный путь / Г. М. Шахунянц. - Москва : Транспорт, 1987. -
479 с.
9. Теоретичш основи впровадження високошвидк1сного руху повд1в в Украшк монограф1я / М.Б.Курган, Д.М. Курган; Дншропетр. нац. ун-т зал1зн. трансп. ϊμ. акад. В. Лазаряна. - Днтро, 2016. - 283 с.
10. С.В. Шкурников и др. Общине требования к проектированию ВСМ Москва-Казань «Транспорт Российской Федерации» № 2(57) 2015 26-29 с.

Claims (3)

1. Переходный участок закругления колеи рельсовой дороги, сопрягающий на длине L смежные с ним прямолинейный участок и круговую часть закругления с постоянным радиусом R, содержащий основание и формирующую колею рельсошпальную решетку, уложенную на нем в соответствии с координатами проекции точек оси переходного участка на горизонтальную плоскость с переменной по длине L величиной поперечного уклона i нормальной к оси линии, касательной к поверхности верха головок рельс, при этом форма оси колеи переходного участка определена местоположением N равноудаленных друг от друга точек с координатами x[n] и у[п], а поперечный уклон i[n] колеи в этих точках изменяется вдоль оси переходного участка в диапазоне 0<i[n]<D/S, где D - расчетное возвышение наружного рельса колеи над внутренним в пределах круговой части оси закругления рельсовой дороги, a S - расстояние между вертикальными осями поперечного сечения рельс колеи, по мере равномерного наращивания на постоянную величину Al=L/(N-1) расстояния l[n] до каждой ее n-й точки от l[0]=0.0 м в начале оси переходного участка и до l[N-1]=L в его конце, отличается тем, что относительные прямоугольные координаты точек проекции оси колеи х[п] и у[п] определены по формулам χ[θ] = θ.θ; χ[«] = χ[λΐ-1] + Δ/ соб^Д(<5) + Δ/?(^)) (1) у[0] = 0.0; у[«] = у[«-1] + АГ8т(д (£) + ^β (δ)) (2) где n - порядковый номер точки оси переходного участка, изменяющийся в диапазоне l<n<(N-1), x[n] - координата, характеризующая расстояние от начала переходного участка до места проекции n-й точки оси на линию, касательную к ней в точках х[0] и у[0],
У[п] - координата, характеризующая смещение в сторону центра кривизны оси ее n-й точки, измеряемое по нормали от линии, касательной к ней в точках х[0] и у[0],
Al - постоянная величина приращения текущей длины l оси колеи переходного участка, измеряемой по ее дуге от ее начала до каждой n-й точки, δ - относительная доля текущей длины 1+Δ1/2 оси колеи переходного участка, измеряемой от ее начала до середины отрезка оси, заключенного между ее смежными точками n-1 и n, вычисляемая по формуле
2и-1
2(У-1) (3) βb(δ) - значение функции основного угла касательной к оси колеи переходного участка, вычисляемого в зависимости от относительной доли δ текущей длины l+Al/2 до точки касания по формуле
Δβ(δ) - значение функции дополнения к основному углу касательной к оси переходного участка, вычисляемой в зависимости от относительной доли δ текущей длины l+Al/2 до точки касания по формуле
где
U - параметр, определяющий величину дополнений к основным значениям угла касательной к оси колеи переходного участка в каждой ее точке с относительной долей δ текущей длины l+Al/2, а поперечный уклон i[n] колеи переходного участка в каждой n-й точке ее оси, изменяющейся в интервале 0<n<(N-1), определен по формуле + (6) где χ - относительная доля текущей длины l оси колеи до точки n, вычисляемая как (7) ib(x) - значение функции основного поперечного уклона колеи переходного участка в точке n, вычисляемое в зависимости от относительной доли χ по формуле
4(/) = ^/0-/(102-/(140-/(172-/(198-/(154-/(66-12/)^ (8)
Δi(χ) - значение функции дополнения к основному поперечному уклону колеи в точке n, вычисляемое в зависимости от относительной доли χ по формуле \ (z) = -Zy Z* (1 - ζ(β - X (35 - χ (76 - / (99 - /(77 - /(33 - 6/)))))() (9) где
Z - параметр, определяющий величину дополнений к основным значениям поперечного уклона колеи в каждой точке n с ее относительной долей χ.
2. Переходной участок по п.1, отличающийся тем, что при заданной величине параметра U, варьи
- 13 039267 руемой в диапазоне от -45 до +55, значение суммы функций основного угла касательной к оси колеи т
W и дополнении к нему на всем протяжении длины L переходного участка строго монотонно и непрерывно изменяется от 0 до L/2R, значение ее первой производной по длине 1, определяющей зако4ΑΨΑΑΑΪΙΑ номерность кривизны оси переходного участка v y V , также строго монотонно и непрерывно изменяется от 0 до 1/R, а для суммы всех последующих 2-х, 3-х, 4-х и 5-х производных этих функций по этой же длине 1 обеспечивается непрерывность изменения их значений, а также равенство их нулю на обеих концах интервала 0<l<L, при этом функция кривизны оси колеи к(/), где /=1/L, в каждой текущей координате длины 1 ПУ определена как *(z) = ^(/)+A U) 0°) где
Ш) = ^5(з15-ф043-/(б720-/(14085-/(19795-Д18579-Д11178-/(3900-600/))^
A(z) = c|z^7-z[74-Z(340-Z(885-Z(1425-Z(1452-Z(914-z(325-50Z)))))))^ (12)
U - параметр, определяющий максимальную величину дополнения к основному значению кривизны оси колеи переходного участка.
3. Переходной участок по п.1, отличающийся тем, что при заданной величине параметра Z, варьи. г η руемой в диапазоне от -88 до +33, значение суммы функций основного поперечного уклона колеи и
А Л дополнений к нему \ / на всем протяжении длины L переходного участка строго монотонно и непрерывно изменяются от 0 до D/S, а значения ее 1-х, 2-х и 3-х производных по длине 1 непрерывно изменяются при соблюдении равенства их нулю на обеих концах интервала 0<l<L.
EA202100109A 2021-02-23 2021-02-23 Переходной участок закругления колеи рельсовой дороги EA039267B1 (ru)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EA202100109A EA039267B1 (ru) 2021-02-23 2021-02-23 Переходной участок закругления колеи рельсовой дороги
PCT/BY2022/000003 WO2022178612A1 (ru) 2021-02-23 2022-02-23 Переходной участок закругления колеи рельсовой дороги
EP22758653.4A EP4324982A1 (en) 2021-02-23 2022-02-23 Transition section of a railway track curve

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
EA202100109A EA039267B1 (ru) 2021-02-23 2021-02-23 Переходной участок закругления колеи рельсовой дороги

Publications (2)

Publication Number Publication Date
EA202100109A1 EA202100109A1 (ru) 2021-12-22
EA039267B1 true EA039267B1 (ru) 2021-12-24

Family

ID=80631145

Family Applications (1)

Application Number Title Priority Date Filing Date
EA202100109A EA039267B1 (ru) 2021-02-23 2021-02-23 Переходной участок закругления колеи рельсовой дороги

Country Status (3)

Country Link
EP (1) EP4324982A1 (ru)
EA (1) EA039267B1 (ru)
WO (1) WO2022178612A1 (ru)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU40396A1 (ru) * 1932-03-17 1934-12-31 В.Н. Образцов Колейный путь в кривых малого радиуса
US3469784A (en) * 1967-06-06 1969-09-30 Jackson E Campbell Apparatus for reducing shells in outer-curve rails
US4860666A (en) * 1988-02-19 1989-08-29 Utdc Inc. Multi-axle, steered articulated railway vehicle with compensation for transitional spirals
JPH03107002A (ja) * 1989-08-28 1991-05-07 Speno Internatl Sa 鉄道線路のレールの削り直し方法及び装置
AT402211B (de) * 1989-10-10 1997-03-25 Oesterr Bundesbahnen Gleis mit zwei auf einem realen flächenstreifen mit überhöhungswinkel angeordneten schienen

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2001273611A1 (en) * 2000-06-20 2002-01-02 Louis T. Klauder Jr. Railroad curve transition spiral design method based on control of vehicle banking motion
AT412975B (de) 2002-07-23 2005-09-26 Wiener Linien Gmbh & Co Kg Gleis mit übergangsbogen und kräfteminimaler überhöhungsrampe
RU2373318C2 (ru) * 2007-12-21 2009-11-20 Государственное образовательное учреждение высшего профессионального образования "Московский государственный университет путей сообщения" (МИИТ) Способ плавного отвода ширины колеи железнодорожного пути и подрельсовое железобетонное основание для его осуществления (варианты)

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU40396A1 (ru) * 1932-03-17 1934-12-31 В.Н. Образцов Колейный путь в кривых малого радиуса
US3469784A (en) * 1967-06-06 1969-09-30 Jackson E Campbell Apparatus for reducing shells in outer-curve rails
US4860666A (en) * 1988-02-19 1989-08-29 Utdc Inc. Multi-axle, steered articulated railway vehicle with compensation for transitional spirals
JPH03107002A (ja) * 1989-08-28 1991-05-07 Speno Internatl Sa 鉄道線路のレールの削り直し方法及び装置
AT402211B (de) * 1989-10-10 1997-03-25 Oesterr Bundesbahnen Gleis mit zwei auf einem realen flächenstreifen mit überhöhungswinkel angeordneten schienen

Also Published As

Publication number Publication date
EP4324982A1 (en) 2024-02-21
WO2022178612A1 (ru) 2022-09-01
EA202100109A1 (ru) 2021-12-22

Similar Documents

Publication Publication Date Title
CN103132411A (zh) 一种高速铁路轨道静态平顺性分析调整方法
CN109657339A (zh) 一种铁道车辆坡道运行综合性能的评估方法
Long et al. Dynamic analysis of railway transition curves
CN116244841B (zh) 一种轨道交通既有线线路纵断面拟合优化方法
CN1114016C (zh) 校正轨道位置的方法
Iwnicki The results of the Manchester benchmarks
Aceituno et al. An alternative procedure to measure railroad track irregularities. Application to a scaled track
Droździel et al. The effect of railway vehicle-track system geometric deviations on its dynamics in the turnout zone
EA039267B1 (ru) Переходной участок закругления колеи рельсовой дороги
CN114912159A (zh) 一种轨道交通线路平面几何线形的拟合方法
CN113609565B (zh) 铁路桥梁轨道静态长波不平顺的验收方法及装置
CN111832618B (zh) 轨道动、静态检查数据的匹配方法
KR100797055B1 (ko) 철도궤도의 효율적 유지보수를 위한 3차원 데이터 형식의철도 궤도틀림도 측정방법
Loktev et al. Curved turnouts for curves of various radii
Klauder Jr Railroad spiral design and performance
Liu et al. Railway wheel profile optimization design based on nurbs curve
CN112257159B (zh) 一种面向bim设计的桥梁工程线路坐标计算方法
Wan et al. Improvement of train-track interaction in turnouts by optimising the shape of crossing nose
CN113378255A (zh) 一种曲线导轨梁的设计方法和曲线导轨梁
JP2011174351A (ja) 列車走行の安全性を確認する計測管理方法
Li et al. Research on Dynamic Detection and Predictive System of Grooved Rail of Tram Twist Irregularity
RU2276216C2 (ru) Способ измерений горизонтальных неровностей (рихтовки) и кривизны в плане рельсовых нитей
Esveld et al. Shape optimization of a railway wheel profile
Drozdziel et al. Pre-processing of wheel and rail geometry in simulation software
JP3660655B2 (ja) 基準器の平面位置調整量算出装置