Влияние жёсткости и длины боковых стоек

Точность и быстродействие гибридной модели

Параметры эксперимента и влияние жёсткости стоек

Упругий увод платформы по xxx:

IKFK​

, а не длиной: иначе сила натяга растёт пропорционально жёсткости и серия смешивает два эффекта.

Девять конфигураций: kx∈{20,200,2000}k_x \in \{20, 200, 2000\}

Онлайн

В них отражение механизма (y,φ)→(y,−φ)(y, \varphi) \to (y, -\varphi)(y,φ)→(y,−φ) действует как (p,s)→(p,−s)(p, s) \to (p, -s)(p,s)→(p,−s); отсюда чётность: yyy -- чётная функция sss, φ\varphiφ -- нечётная, и коэффициенты при нечётных (соответственно чётных) степенях sss заведомо нулевые. Базис с фиксированной чётностью содержит шесть коэффициентов для yyy и четыре для φ\varphiφ вместо десяти и даёт те же ошибки; свободная подгонка чётность воспроизводит сама.

В них отражение механизма (y,φ)→(y,−φ)(y, \varphi) \to (y, -\varphi)(y,φ)→(y,−φ) действует как (p,s)→(p,−s)(p, s) \to (p, -s)(p,s)→(p,−s); отсюда чётность: yyy -- чётная функция sss, φ\varphiφ -- нечётная, и коэффициенты при нечётных (соответственно чётных) степенях sss заведомо нулевые. Базис с фиксированной чётностью содержит шесть коэффициентов для yyy и четыре для φ\varphiφ вместо десяти и даёт те же ошибки; свободная подгонка чётность воспроизводит сама.

Вклад работы: количественный анализ влияния жёсткости и длины боковых стоек kxk_xkx​, LLL на упругий увод платформы Δx\Delta xΔx и на расхождение прямой и обратной задач (девять конфигураций, kx∈{20,200,2000}k_x \in \{20, 200, 2000\}kx​∈{20,200,2000} Н/м); гибридная модель для обеих задач: IKhybrid=IKideal+P3\mathrm{IK}_{\mathrm{hybrid}} = \mathrm{IK}_{\mathrm{ideal}} + \mathbf{P}_3IKhybrid​=IKideal​+P3​ и FKhybrid=FKideal+P3′\mathrm{FK}_{\mathrm{hybrid}} = \mathrm{FK}_{\mathrm{ideal}} + \mathbf{P}'_3FKhybrid​=FKideal​+P3′​ в повёрнутых координатах Q~=(p,s)\tilde{Q} = (p, s)Q~​=(p,s), коэффициенты которых идентифицируются МНК по кинетостатическим невязкам при активных моментах на приводах; верификация точности и быстродействия во фреймворке Flexia: ошибки отображения (0.330.330.33 мкм и 0.0240.0240.024 мрад в прямой задаче, 0.00950.00950.0095 мрад в обратной) при 0.13…0.140.13\dots0.140.13…0.14 мкс на вызов -- на четыре порядка быстрее расчёта численного равновесия.

Кинетостатика во Flexia и расхождение прямой и обратной задач

Именно это делает возможным сравнение приближённых моделей с эталоном: аналитика IKideal\mathrm{IK}_{\mathrm{ideal}}IKideal​ отличается от эталона на поправку ΔΘ\Delta\ThetaΔΘ, которая и подлежит аппроксимации.

а не как отдельная постановка с иными граничными условиями.

(реакции связей)

Модель

управляется

Время расчёта

перебор всех четырёх сборок для диагностики стоит 0.60.60.6 мкс и 211221122112 Б, то есть O(1)O(1)O(1) без итераций и с детерминированным временем шага.

(считается только запрошенный вариант сборки, без выделения памяти)

полиномиальная поправка

Величина x=12(xC+L2cos⁡α+xF+L2cos⁡γ)x = \tfrac{1}{2}(x_C + L_2\cos\alpha + x_F + L_2\cos\gamma)x=21​(xC​+L2​cosα+xF​+L2​cosγ) тождественно равна нулю -- это встроенный контроль формул; если ∣cos⁡α∣>1|\cos\alpha| > 1∣cosα∣>1, ∣cos⁡γ∣>1|\cos\gamma| > 1∣cosγ∣>1 или ∣sin⁡φ∣>1|\sin\varphi| > 1∣sinφ∣>1, поза недостижима, и функция возвращает признак отсутствия решения.

Отображение IKideal:(y,φ)↦Θ\mathrm{IK}_{\mathrm{ideal}}: (y, \varphi) \mapsto \Theta

правило «минимум ∣φ∣|\varphi|∣φ∣» её не выделяет (поза у всех четырёх общая), поэтому сборка фиксируется один раз по номинальной позе и далее ведётся продолжением.

Рабочая сборка -- (+1,−1)(+1, -1)(+1,−1)

попарное расстояние по (θ1,θ2)(\theta_1, \theta_2)(θ1​,θ2​) -- не менее 0.590.590.59 рад при медиане 0.770.770.77 рад (400 точек рабочего поля, во всех по четыре решения).

0.28+0.10+0.28

Исследовательский вопрос

Топология и параметры PRBM-модели

тяг обязателен: трос не воспринимает сжатие. При таком натяге область значений ОЗК целиком остаётся в области растяжения (+396+396+396…+712+712+712 мкм),

тяг обязателен: трос не воспринимает сжатие. При таком натяге область значений ОЗК целиком остаётся в области растяжения (+396+396+396…+712+712+712 мкм),

жёсткость каждой kx=20k_x = 20kx​=20 Н/м, заделки в x=±0.33x = \pm 0.33x=±0.33 м, предварительный натяг P=535P = 535P=535 мкм (усилие 10.710.710.7 мН). Н

Смещение по xxx удерживает не кинематическая связь, а пара упругих растяжек -- тросов с линейными пружинами

: основание, два входных кривошипа L1L_1L1​, два шатуна L2L_2L2​ и выходное звено (платформа) L3L_3L3​. Оси кривошипов AAA и BBB разнесены на D=0.312D = 0.312D=0.312 м; длины L1=L2=0.15L_1 = L_2 = 0.15L1​=L2​=0.15 м, L3=0.10L_3 = 0.10L3​=0.10 м;

геометрия

Расчётная схема --

Статья

данным: сплайны [21] и регрессия по выборке

Альтернатива -- аппроксимация

Для параллельных роботов ту же роль играют аналитические решения:

-- поэтому для манипуляторов общего вида строят геометрические решения в замкнутой форме

, итерационно,

Кинетостатика плоских податливых цепей с большими перемещениями решается методом нескольких сегментов [13], а для параллельных структур силовой анализ строится на тех же уравнениях равновесия [14]. Для податливых шестизвенников типа Уатта исследованы кинематические и динамические отклики [15], однако задача быстрого отображения «поза → углы приводов» для них не ставилась.

но не отвечают на вопрос «какие углы приводов соответствуют заданной позе»

Границы применимости матричных методов -- малые деформации, отсутствие контакта -- исследованы отдельно [10], а сводка всего семейства (PRBM, матричные, энергетические и МКЭ-методы) дана в обзоре [11] [11] M. Ling, L. L. Howell, J. Cao, и G. Chen, «Kinetostatic and dynamic modeling of flexure-based compliant mechanisms: A survey», ASME Applied Mechanics Reviews, т. 72, вып. 3, с. 030802, 2020, doi: 10.1115/1.4045679. .

заменяется вращательной парой с

Вклад работы

Обратная задача в контуре реального времени

разделяются:

Доминирующий подход -- метод псевдотвёрдого тела (PRBM)

Обзор методов моделирования податливости

Контекст и актуальность

Список литературы

параметры

Определение модели

Численное интегрирование

Вычисление правой части

Приведение к системе первого порядка

Упругие опоры. Матрицы демпфирования C и жёсткости K

Матрица масс M

Введение

параметры

При равных дисбаланса

Суммарный кинетический момент двух роторов из

ля роторов с тройной симметрией лопастей тензоры Jup и Jdown не зависят от угла поворота.

Матрица масс зависит от времени только через фазу вращения роторов α(t), п

Для учёта различного положения центров масс тел относительно точки, к которой приводится тензор, используется теорема Гюйгенса–Штейнера.

блочную структуру

Для роторов rrot,i=pi+ci, где pi — положение системы координат ротора, а ci — положение его центра масс в этой системе.

Динамика системы описывается линейной матричной моделью второго порядка

бщая для обоих роторов угловая скорость вращения по модулю в момент времени t,

величина дисбалансирующей массы

диагональными по осям X, Y, Z:

упруго-вязких

еремещения центра вращения qt и малые угловые отклонения системы qr

φ, θ, ψ

Координаты вектора x, y, z

𝐂

where

The scheme investigated here is single-stage and does not require inner iterations, so the cost of a step is fixed and predictable, which makes it especially interesting for hardware-in-the-loop (HIL) applications where execution time is more important than accuracy.

A- and L-stability. The L-stability of the scheme means that the infinitely stiff decaying modes are suppressed entirely, so that the stiffness of the constraints does not restrict the step size from the standpoint of stability.

The physical meaning of the multipliers was checked statically: at a constant prescribed crank angle the mechanism must be at rest, and the actuator multiplier must equal the derivative of the potential energy with respect to the crank angle, dU/dθ, taken along the constraint manifold (gravity plus the energy of the spring). The configuration for each angle is taken from the exact analytical solution of the loop geometry, and dU/dθ is computed by automatic differentiation of the same potential energy along that solution; the states are indeed static — the maximal acceleration in them does not exceed 2.8⋅10−14 m/s2. The relative discrepancy between the actuator multiplier and dU/dθ is 2.9⋅10−15 over the range θ∈[0.4,2.2] rad.

This reference shares no code with the CROS core: it uses analytically eliminated coordinates rather than numerical projection onto the null space of the constraint Jacobian.

Constraint residual. The closure error of a hinge — the distance between the hinge point computed from one link and from the other — is evaluated from the raw state without using the multipliers, over t∈[0,20] s. The attachment to ground is satisfied to machine precision: its constraint rows contain only identity blocks, so the solution of the linear system reproduces them exactly. On halving the step the closure error of the hinges decreases by factors of 3.44, 3.83, 3.91 and 3.99, that is, the ratios converge to 4 — the signature of second order. The closure error is not equal to zero, so the scheme does not enforce the constraint exactly; it enforces it to the accuracy of the discretization, and the achievable geometric tolerance is determined by the step size.

and the resulting index-2 or index-1 system is discretized directly.

In the figures these variants are denoted as index 1, index 1 + stabilization and index 1 + projection, and the direct form as index 3.

using StaticArrays.jl; the figures were produced with Makie.jl.

index 3: {Mq¨=Q(q,q˙,t)+CT(q)λ,Φ(q,t)=0, where

The scheme is A- and L-stable, so it tolerates the stiff constraint dynamics that force explicit index-reduced integrators to very small steps.