Результаты расчета влияют на выбор двигателей разгрузки гиростабилизатора
Последовательность расчета возмущающих моментов
С помощью программы Torques 4kurs.mw (автор П.Г. Русанов),
дополняющей рекомендуемое Уч. Пособие [1]
Арсеньев В.Д. Учебное пособие: Расчет и синтез параметров гиростабилизаторов для маневренных объектов. Часть 1. Издательство МГТУ им. Н.Э. Баумана, 2013, 44 с.
Примечания:
Результаты расчета влияют на выбор двигателей разгрузки гиростабилизатора
2. Красным цветом помечены места возможные исправлений в тексте программы
Шаг 1. Назначение исходных данных (ввести числовые значения в операторы присвоения)
| Параметр
| Значение
| Единицы
измерения
| Физ.-мех. характеристика
|
| g
| 10.
| м/с2
| ускорение поля силы тяжести;
|
| Moc
| 2.
| кг
| масса оптики
|
| Mchel
| 0.5
| кг
| масса чувствит. элемента
|
| Mos
| 0.2
| кг
| масса полуосей
|
| Mdop
| 0.1
| кг
| масса доп. элементов
|
Шаг 2. Расчет масс и весов (это результаты автономного выполнения программы)
| Параметр
| Значение
| Единицы
измерения
| Физ.-мех. характеристика
|
| Mpl:= Moc +Mchel+ Mos+Mdop
| 2.8
| кг
| масса платформы
|
| Mpl:=ceil(Mpl)
|
| кг
| округление величины
|
| Ppl:=Mpl*g;
|
| Н
| вес платформы
|
| Mnr:=сeil(0.3*Mpl)
|
| кг
| масса наружной рамы
|
| Pnr:=Mnr*g
|
| Н
| вес наружной рамы
|
| Psum:=ceil(Ppl+Pnr)
|
| Н
| суммарный вес платформы и нар.рамы
|
Шаг 3. Расчет опорных сил реакции (ввести числовые значения для параметров, выделенных цветом)
| Параметр
| Значение
| Ед.
измер.
| Физ.-мех. характеристика
|
| Pi180:=evalf(Pi/180)
| 0.0175
| рад/град
| коэффициент пересчета величин из [град] в [рад]
|
| nx:=ny:=nz:
| 3.
| -
| Перегрузки ЛА по лин. ускорениям
|
| nv
| 1.
| -
| Линейные виброперегрузки
|
| FRnr:= Psum/2*(evalf(sqrt(nx^2+nz^2))+nv)
|
| Н
| Радиальная сила реакции опоры нар.рамы
|
| FAnr:=ceil(Psum*(ny+nv))
|
| Н
| Аксиальная сила реакции опоры нар.рамы
|
| bkGr
|
| град
| угол контакта в ШП в [град]
|
| bk:=evalf(bkGr *Pi180)
| 0.35
| рад
| угол контакта в ШП в [рад]
|
| APNnr:=
ceil(0.1*(1.58*FRnr*tan(bk)+0.5*FAnr))*10
|
| Н
| Аксиальный преднатяг
ШП-ков нар.рамы
|
| FRpl:= ceil(Ppl/2* (evalf(sqrt(nx^2+ny^2+nz^2))+nv))
|
| Н
| Радиальная сила реакции опоры платформы
|
| FApl:=ceil(Ppl*(sqrt(nx^2+nz^2)+nv))
|
| Н
| Аксиальная сила реакции опоры платформы
|
| APNpl:=
ceil(0.1*(1.58*FRpl*tan(bk)+0.5*FApl))*10
|
|
Н
| Аксиальный преднатяг ШП-ков платформы
|
Шаг 4. Расчет моментов трения в опорах (ввести числовые значения для параметров, выделенных цветом)
| Параметр
| Значение
| Единицы
измерения
| Физ.-мех. характеристика
|
|
M0
|
7.2
| сН см =
= Г см-
| Момент страгивания ШП данной серии
при R > F0 = 5Н=500сН
|
| F0
| 500.
| сН
| Пороговое значение радиальной силы реакции в ШП
|
| k1
| 0.0069
| см
| Эмпирический коэффициент
|
| k2
| 0.0035
| см
| Эмпирический коэффициент
|
| k3
| 0.0008
| см
| Эмпирический коэффициент
|
| Mtry:=2*ceil(0.1*evalf(M0+F0*(k1-k2)+
+(k2*FRnr+k3*APNnr)*100))*10
|
| сН см
| Момент трения вокруг оси н. рамы
|
| Mtrx:=2*ceil(0.1* evalf (M0+F0*(k1-k2)+
+(k2*FRpl+k3*APNpl)*100))*10
|
| сН см
| Момент трения вокруг оси платформы
|
Шаг 5. Расчет моментов от остаточного небаланса масс платформы и нар.рамы при балансировке
| Параметр
| Знач.
| Ед.
изм.
| Физ.-мех. характеристика
|
| FRnr0:= Psum/2
|
| Н
| Радиальная сила реакции опоры н.рамы при балансировке
|
| FAnr0:=
|
| Н
| Аксиальная сила реакции опоры н.рамы при балансировке
|
| M00:=M0+F0*(k1-k2)
| 8.9
| сН см
| Вспомогательная константа расчета
|
| Mnby0:=2*ceil(M00+(k2*FRnr0)*100)
|
|
сН см
| Максимальный момент небаланса вокруг оси нар. рамы при балансировке
|
| FRpl0:= Ppl/2
|
| Н
| Радиальная сила реакции опоры платформы при балансировке
|
| FApl0:= 0
|
| Н
| Аксиальная сила реакции опоры платформы при балансировке
|
| Mnbx0:=2*ceil(M00+(k2*FRpl0)*100)
|
| сН см
| Максимальный момент небаланса вокруг оси платформы при балансировке
|
Шаг 6. Расчет моментов от остаточного небаланса масс платформы и н.рамы в рабочем режиме эксплуатации (ввести числовые значения для параметров, выделенных цветом, в соответствии с ТЗ)
| Параметр
| Значение
| Единицы
измерения
| Физ.-мех. характеристика
|
| aM
|
| град
| Амплитуда угла прокачки вокруг оси н.р.
|
| bM
|
| град
| Амплитуда угла прокачки вокруг оси платформы
|
| MnbYmin
MnbХmin
MnbYmax
MnbXmax
|
-10000
-10000
| сН см
| Для последующего цикла вычислений предварительные значения наименьшей (с учетом знака) и наибольшей из величин моментов от небалансов вокруг осей нар. рамы и оси платформы в рабочем режиме
|
| Ниже представлены 2 вложенных цикла для расчета экстремальных значений моментов от небалансов (с индексом nb) методом перебора значений каждого из углов прокачки в цикле через 1 град от мин. до макс значения с фиксацией экстремальных значений моментов от небалансов по двум осям двухосного гиростабилизатора
(при желании сократить время расчетов шаги перебора значений углов, равные 1, помеченные в программе цветом, можно изменить в сторону увеличения)
|
Знач.
|
Ед.
изм.
|
| for na from -aM by 1 to +aM do
alf:=na*Pi/180;ca:=cos(alf);sa:=sin(alf);
for nb from -bM by 1 to +bM do
bet:=nb*Pi/180; cb:=cos(bet);sb:=sin(bet);
MnbY:=evalf(Mnby0*(nx*(sa-ca)+nz*(sa+ca))+
+Mnbx0*cb*(nx*(sa-ca*cb)+
+ny*sb+nz*(sa*cb+ca)));
MnbX:=evalf(Mnbx0*((nx*ca-nz*sa)*(cb+sb)+ny*(sb-cb)));
if MnbY < MnbYmin then MnbYmin:=MnbY:end if;
if MnbX < MnbXmin then MnbXmin:=MnbX:end if;
if MnbY > MnbYmax then MnbYmax:=MnbY:end if;
if MnbX > MnbXmax then MnbXmax:=MnbX:end if;
od:
od:
print (" diapazon MY nb =",ceil(MnbYmin),ceil(MnbYmax));
print (" " diapazon MX nb =",ceil(MnbXmin),ceil(MnbXmax));
|
|
сН см
сН см
|
|
diapazon MY nb
| Расчетные границы диапазонов изменения величин моментов от небаланса вокруг оси н. рамы
| -362
|
сН см
|
|
diapazon MY nb
| Расчетные границы диапазонов изменения величин моментов от небаланса вокруг оси платформы
| -169
|
сН см
|
| MnbY:=
ceil(max(abs(MnbYmin), abs(MnbYmax)))
| Наибольшая по модулю величина момента от небаланса вокруг оси Y н.рамы.
|
|
сН см
|
| MnbX:=
ceil(max(abs(MnbXmin),abs(MnbXmax)))
| Наибольшая по модулю величина момента от небаланса вокруг оси X платформы
|
|
сН см
|
Шаг 7. Расчет наибольшей величины инерционного момента вокруг оси у наружной рамы в двухосном гиростабилизаторе при трехосном вращении подвижного объекта по трем осям с угловой скоростью w и угловым ускорением e
(ввести числовые значения для параметров, выделенных цветом)
| Параметр
| Значение
| Ед.
изм.
| Физ.-мех. характеристика
|
|
apl ´ bpl ´ cpl
| 0.22
0.22
0.02
|
м
| Габаритные размеры платформы вдоль ее осей xyz
|
| JplZ:= Mpl*(apl^2+bpl^2)/12;
| 0.0242
| кг м2
| момент инерции массы платформы вокруг оси z
|
|
anr ´ bnr ´ cnr
| 0.24
0.24
0.08
|
м
| Габаритные размеры нар.рамы вдоль ее осей xyz
|
| t
| 0.01
| м
| Толщина стенки рамы
|
| V1:= anr*bnr*cnr
| 0.001152
| м3
| Объем нар. рамы вместе с полостью
|
| V2:= (anr- 2 *t) * (bnr- 2 *t) *cnr
| 0.000968
| м3
| Объем полости нар. рамы
|
| V:= V1 - V2;
| 0.000184
| м3
| Объем наружной рамы
|
| plotn:=Mnr/V
| 5434.
| кг/м3
| Удельная плотность материала нар. рамы
|
| JnrX:=plotn*(Vnr1*(bnr^2+cnr^2) -Vnr2*((bnr-t)^2+cnr^2))/12
| 0.00887
| кг м2
| момент инерции массы нар. рамы вокруг оси х
|
| JnrY:=plotn*(Vnr1*(anr^2+cnr^2) -Vnr2*((anr-t)^2+cnr^2))/12
| 0.00887
| кг м2
| момент инерции массы нар. рамы вокруг оси y
|
| JnrZ:=plotn*(Vnr1*(anr^2+bnr^2) -Vnr2*((anr-t)^2+(bnr-t)^2))/12
| 0.0177
| кг м2
| момент инерции массы нар. рамы вокруг оси z
|
| B:=JnrY+JplZ;
| 0.0330
| кг м2
| Вспомогательная величина
|
| Pi180:=evalf(Pi/180)
| 0.01745
| -
| Вспомогательная математ.величина
|
|
w100:=100*Pi180
|
1.745
|
рад/с
| Угловые скорости ЛА по трем осям ЛА, равные 100 град/с, в [рад/с]
|
| wxc:=w100; wyc:=w100; wzc:=w100;
|
1.745
|
рад/с
| Проекции угловой скорости ЛА на оси ЛА
|
|
wt140:=140*Pi180
|
2.443
|
рад/с2
| Угловые ускорения ЛА по трем осям ЛА, равные 140 град/с2, в [рад/с2]
|
| wtxc:=wt140; wtyc:=wt140;
wtzc:=wt140;
|
2.443
|
рад/с2
| Проекции углового ускорения ЛА на оси ЛА
|
|
wwMYmin
wwMYmax
|
-10000
|
Н м
| Для последующего цикла вычислений предварительные значения наименьшей (с учетом знака) и наибольшей величин инерционного момента от вращения ЛА вокруг трех осей
|
| Ниже представлены 2 вложенных цикла для расчета экстремальных значений инерционного момента вокруг оси y нар. рамы от угловых скоростей ЛА при изменении значений каждого из углов прокачки в цикле через 1 град от мин. до макс значения с фиксацией экстремальных значений инерционного момента по оси y нар. рамы.
(при желании сократить время расчетов помеченные в программе цветом шаги перебора значений углов, равные 1, можно изменить в сторону увеличения)
|
Знач.
| Ед.
изм.
|
| for nb from -bM by 1 to bM do
bet:=evalf(nb*Pi180); cb:=cos(bet);sb:=sin(bet);tng:=evalf(tan(bet));
A:=(JnrX+JnrY/cb^2+(B+JplZ)*tng^2-JnrZ);
for na from -aM by 1 to aM do
alf:=evalf(na*Pi180); sa:=evalf(sin(alf)); ca:=evalf(cos(alf));
s2a:=evalf(sin(2*alf)); c2a:=evalf(cos(2*alf));
wwinMY:=evalf(B*tng*(wtxc*ca-wtzc*sa+wyc*(wxc*sa+wzc*ca))+
+A*((wxc^2-wzc^2)*s2a +2*wxc*wzc*c2a)/2);
if wwinMY < wwMYmin then wwMYmin:=wwinMY:end if;
if wwinMY > wwMYmax then wwMYmax:=wwinMY:end if;
od:
od:
!! Внимание далее инерционные моменты оцениваются в [сН см]
print ("wwMYmin =",ceil(wwMYmin*10^4));
print ("wwMYmax =",ceil(wwMYmax*10^4));
wwinMY:=ceil(max(abs(wwMYmin),abs(wwMYmax))*10^4);
|
| -
-
кг м2
-
-
Н м
сН см
сН см
сН см
|
|
wwMYmin
| Это наименьшее значение инерционного момента вокруг оси y нар. рамы
|
-6030
|
сН см
|
|
wwMYmax
| Это наибольшее значение инерционного момента вокруг оси y нар. рамы
|
|
сН см
|
| max(abs(wwMYmin), abs(wwMYmax))
| Наибольшее значение из двух
|
|
сН см
|
Шаг 8. Расчет наибольшей величины инерционного момента вокруг оси y наружной рамы в двухосном гиростабилизаторе при угловых колебаниях подвижного объекта по трем осям с амплитудой g0= gamGr и циклической частотой f
(ввести числовые значения для параметров, выделенных цветом, в соответствии с ТЗ)
| Параметр
| Значение
| Единицы
измерения
| Физ.-мех. характеристика
|
| gamGr
|
| град
| Амплитуды углов угловых колебаний в [град]
|
| gam= gamGr*Pi180
| 0.0345
| рад
| Амплитуды углов угловых колебаний в [рад]
|
| f
|
| гц
| Циклическая частота угловых колебаний ЛА вокруг осей xyz
|
| > wf:=2*Pi*f
| 18.85
| с-1
| Круговая частота угловых колебаний ЛА вокруг осей xyz
|
|
gamMYmin
gamMYmах
|
-10000
|
Н м
| Для последующего цикла вычислений предварительные значения наименьшей (с учетом знака) и наибольшей величин инерционного момента от угловых колебаний ЛА вокруг трех осей
|
| Ниже представлены 3 вложенных цикла для расчета диапазона значений инерционного момента вокруг оси y нар. рамы от угловых колебаний ЛА
g(t)= g0 sin(2π f t) c циклической частотой f при изменении значений каждого из углов прокачки платформы в цикле через 1 град от мин. до макс значения, а также изменения времени от 0 до времени 2t - двух периодов угловых колебаний с шагом ‘ ht’, равным t/’ nt’ (‘ nt’ =30) от периода t, для отбора мини- максных значений инерционного момента вокруг оси y нар. рамы.
(при желании сократить время расчетов помеченные в программе цветом шаги перебора значений углов, равные 1, можно изменить в сторону увеличения и увеличить ‘ ht’ - шаг изменения времени путем уменьшения параметра ’ nt’)
|
Знач.
|
Ед.
изм.
|
|
| nt
ht:=2*Pi/(f*nt);
for nb from -bM by 1 to bM do
bet:=evalf(nb*Pi180); cb:=cos(bet);sb:=sin(bet);tng:=evalf(tan(bet));
A:=evalf(JnrX+JnrY/cb^2+(B+JplZ)*tng^2-JnrZ);
for na from -aM by 1 to aM do
alf:=evalf(na*Pi/180); sa:=evalf(sin(alf)); ca:=evalf(cos(alf));
s2a:=evalf(sin(2*alf)); c2a:=evalf(cos(2*alf));
t:=0;
for kt from 0 to 2*nt do
wxc:=-gam*wf*cos(wf*t); wyc:=wxc; wzc:=wxc;
wtxc:=-gam*wf^2*sin(wf*t); wtyc:=wtxc; wtzc:=wtxc;
gamMY:=evalf(B*tng*(wtxc*ca-wtzc*sa+wyc*(wxc*sa+wzc*ca)+ A*((wxc^2- wzc^2)*s2a +2*wxc*wzc*c2a)/2));
if gamMY < wwMYmin then gamMYmin:=gamMY:end if;
if gamMY > wwMYmax then gamMYmax:=gamMY:end if;
t:=t+ht;
od;
od;:od:
!! Внимание:далее инерционные моменты оцениваются в [сН см]
print ("gamMYmin =", ceil(gamMYmin*10^4));
print ("gamMYmax =", ceil(gamMYmax*10^4));
gamMY:=ceil(max(abs(gamMYmin), abs(gamMYmax))*10^4);
|
0.0698
| -
сек
-
кг м2
-
-
сек
1/с
1/с2
Н м
сек
сН см
сН см
сН см
|
|
|
gamMYmin
| Наименьшее отрицательное значение инерционного момента вокруг оси y нар. рамы
|
-6040
|
сН см
|
|
|
gamMYmax
| Наибольшее положительное значение инерционного момента вокруг оси y нар. рамы
|
|
сН см
|
|
| gamMY:= ceil(max(abs(gamMYmin), abs(gamMYmax))*10^4)
| Максимальное значение инерционного момента от угловых колебаний вокруг оси y нар. рамы
|
|
сН см
|
|
| | | | | | | |
Шаг 9. Расчет наибольшей величины момента вокруг оси y наружной рамы в двухосном гиростабилизаторе вследствие непропорциональных упругих деформаций конструкции (неравножесткости) при линейных перегрузках ЛА
(ввести числовые значения для параметров, выделенных цветом)
| Параметр
| Знач.
| Ед. изм
| Физ.-мех. характеристика
|
| СBradHP
|
| Н/мкм
| Коэффициент жесткости ШП-ка нар. рамы в радиальных направлениях
|
| СBaxeHP
|
| Н/мкм
| Коэффициент жесткости ШП-ка нар. рамы в осевом направлении
|
| СBradPL
|
| Н/мкм
| Коэффициент жесткости ШП-ка платформы в радиальных направлениях
|
| СBaxePL
|
| Н/мкм
| Коэффициент жесткости ШП-ка платформы в осевом направлении
|
| Последующие коэффициенты жесткости для наружной рамы из алюминиевого сплава с габаритными размерами 240 мм ´240 мм ´20 мм и прямоугольным поперечным сечением с размерами 10 мм ´20 мм получены с помощью компьютерной программы ЛИРА.
Оси Х, Y – лежат в «плоскости» нар. рамы и являются осями подвеса платформы
|
|
| CXHP
| 1.45
| Н/мкм
| Коэффициент жесткости Нар. Рамы при ее нагружении силой вдоль оси Х
|
|
| CYHP
| 1.58
| Н/мкм
| Коэффициент жесткости Нар. Рамы при ее нагружении силой вдоль оси Y
|
|
| CZHP
| 1.39
| Н/мкм
| Коэффициент жесткости Нар. Рамы при ее нагружении силой вдоль оси Z
|
|
| CqXHP
| 7.8
| Н/мкм
| Коэффициент жесткости Нар. Рамы при ее нагружении распределенными объемными силами вдоль оси Х
|
|
| CqYHP
| 21.6
| Н/мкм
| Коэффициент жесткости Нар. Рамы при ее нагружении распределенными объемными силами вдоль оси Y
|
|
| CqZHP
| 10.4
| Н/мкм
| Коэффициент жесткости Нар. Рамы при ее нагружении распределенными объемными силами вдоль оси Z
|
|
|
|
|
| cHP0X:=1/(1/CXHP +1/CqXHP+ 1/СBradHP);
| 1.217
| Н/мкм
| Эквивалентный коэффициент жесткости HP вдоль оси Х
|
|
| cHP0Z:=1/(1/CZHP +1/CqZHP+ 1/СBradHP);
| 1.220
| Н/мкм
| Эквивалентный коэффициент жесткости HP вдоль оси Z
|
|
| cPL0Z:=СBradPL;
|
| Н/мкм
| Эквивалентный коэффициент жесткости PL вдоль оси Z
|
|
| cPL0Y:=СBradPL;
|
| Н/мкм
| Эквивалентный коэффициент жесткости PL вдоль оси Y
|
|
| McHPy:=ceil(evalf((nx^2+nz^2)*
(Pnr^2*abs(1/cHP0X -1/cHP0Z) +
Pnr*Ppl*(1/cHP0X -1/ cPL0Z)+
Ppl^2*/cPL0Z))/4/100);
|
|
сН см
| Максимальный момент переносных сил инерции вокруг оси Y наружной рамы у неравножесткого подвеса
|
|
|
McPLx:=ceil(evalf((nx^2+nz^2)*
Ppl^2*abs(1/cPL0Y -1/cPL0Z))/4/100);
|
|
сН см
| Максимальный момент переносных сил инерции вокруг оси Х платформы у неравножесткого подвеса
|
|
| | | | | | | | | | |
Шаг 10. Расчет наибольшей величины момента вокруг оси Y наружной рамы в двухосном гиростабилизаторе вследствие упругих деформаций конструкции (неравножесткости) при линейной вибрации основания гиростабилизатора
(ввести числовые значения для параметров, выделенных цветом, в соответствии с условиями ТЗ)
| Параметр
| Знач.
| Ед. изм
| Физ.-мех. характеристика
|
|
lambdaHP
|
| -
| Коэффициент динамичности механической системы гиростабилизатора поперечном направлении к Y - оси вращения нар. рамы
|
|
lambdaPL
|
| -
| Коэффициент динамичности механической системы гиростабилизатора поперечном направлении к X - оси вращения платформы
|
| McvHPy:=ceil(evalf(nv^2* lambdaHP* lambdaHP*
(Pnr^2*abs(1/cHP0X -1/cHP0Z)+
Pnr*Ppl*(1/cHP0X -1/ cPL0Z)+
Ppl^2*/cPL0Z))/2/100);
|
|
сН см
| Максимальный момент у неравножесткого подвеса вокруг оси Y наружной рамы от переносных сил инерции при вибрации основания
|
Шаг 11. Расчет наибольшей величины момента тяжения токоподводов от основания к наружной раме и далее к платформе для наиболее рационального варианта реализации
(ввести предельные числовые значения для параметров, выделенных цветом, в соответствии с ТЗ)
| Параметр
| Знач.
| Ед. изм
| Физ.-мех. характеристика
|
|
| Для коллекторных токоподводов
|
|
| Nknr
|
| -
| Количество токопроводящих колец в коллекторе между корпусом и нар. рамой
|
|
| Nkpl
|
| -
| Количество токопроводящих колец в коллекторе между нар. рамой и платформой
|
|
| Nsh
|
| -
| Количество щеток, контактирующих с одним токопроводящим кольцом
|
|
| Rk
| 0.5
| см
| Радиус кольца
|
|
| Fk
| 0.1
| H
| Сила прижатия щетки к кольцу
|
|
| f
| 0.2
| -
| Коэффициент трения скольжения между щеткой и кольцом
|
|
| Mknr:=f*Fk*Rk*Nknr*Nsh;
|
| сН см
| момент тяжения коллекторного токоподвода к нар. раме
|
|
| Mkpl:=f*Fk*Rk*Nkpl*Nsh;
|
| сН см
| момент тяжения коллекторного токоподвода к платформе
|
|
| Для гибких токоподводов
|
|
| Nnr
|
| -
| Количество электропроводов в жиле токоподвода к нар. раме
|
| Nnr
|
| -
| Количество электропроводов в жиле токоподвода к платформе
|
|
tetaGrnr
|
|
град
| Максимальный угол изгиба и скручивания жилы токоподвода к нар. раме в [град]
|
| tetanr:= tetaGrnr*Pi180;
| 1.22
| рад
| Максимальный угол изгиба или скручивания жилы токоподвода к нар. раме в [рад]
|
|
tetaGrpl
|
|
град
| Максимальный угол изгиба и скручивания жилы токоподвода к платформе в [град]
|
| tetapl:= tetaGrpl*Pi180;
| 1.05
| рад
| Максимальный угол изгиба или скручивания жилы токоподвода к платформе в [рад]
|
| LtpHP
|
| см
| Длина жилы между местами ее закрепления на корпусе и наружной раме
|
| LtpPL
|
| см
| Длина жилы между местами ее закрепления на наружной раме и платформой
|
| MtrwHPy:=ceil(6.5*Nnr*tetanr/
LtpHP);
|
| сН см
| Момент тяжения гибкого токоподвода к нар. раме
|
| MtrwPLx:=ceil(6.5*Npl*tetapl/
LtpPL);
|
| сН см
| Момент тяжения гибкого токоподвода к платформе
|
| Сравнение моментов тяжения для двух вариантов исполнения токоподводов
(результат сравнения указывает на преимущество применения гибких токоподводов)
|
| Mtoky:=min(Mknr,MtrwHPy);
|
| сН см
|
|
| Mtokx:=min(Mkpl,MtrwPLx);
|
| сН см
|
|
| | | | | | | | |
Шаг 12. Итоговый расчет наибольших суммарных возмущающих моментов
вокруг осей нар. рамы и платформы
| №
| Возмущающие моменты вокруг оси нар. рамы
| Знач.
|
|
| момент трения в опорах
|
|
|
| момент от остаточного небаланса
|
|
|
| инерционный момент от вращения ЛА
|
|
|
| инерционный момент от угловой вибрации основания
|
|
|
| момент от неравножесткости карданова подвеса при линейных перегрузках ЛА
|
|
|
| момент от неравножесткости при линейной вибрации основания
|
|
|
| момент тяжения токоподводов
|
|
|
| инерционный момент привода [1]
|
|
|
| демпфирующий момент привода [1]
|
|
|
| Суммарный момент
| 13200 сН см
|
| №
| Возмущающие моменты вокруг оси платформы
| Знач.
|
|
| момент трения в опорах
|
|
|
| момент от остаточного небаланса
|
|
|
| момент тяжения токоподводов
|
|
|
| инерционный момент привода [1]
|
|
|
| демпфирующий момент привода [1]
|
|
|
| Суммарный момент
| 630 сН см
|
Рекомендуемый порядок выполнения расчета.
1. Скопировать исходный файл с текстом программы Torques 4kurs.mw и сохранить в новом файле под новым именем, например: VMom.mw.
2. В текстовом редакторе внести необходимые исправления в текст VMom.mw программы и сохранить его.
3. Переписать файл VMom.mw в каталог пакета Maple.
4. Войти в этот каталог и отправить кликом на счет файл VMom.mw.
5. После появления текста программы на экране в среде Maple кликнуть на иконку !!!, чтобы инициализировать процесс вычислений.
(Следует иметь ввиду!!, что время расчета по программе Torques 4kurs.mw с исходными установками составляет примерно 15 минут. Поэтому для первого пробного просчета можно укрупнить шаги перебора значений углов na, nb с 1 на 10, а параметр nt уменьшить с 30 на 10.)
6. После завершения расчета с помощью команд File à Save As следует сохранить результаты расчета под новым именем, например VMom REZ.mw. Впоследствии из сохраненного файла в среде Maple можно копировать фрагменты результатов в файлы WORD.
Примечание. Внимание!! Параметры форматирования исходного текстового файла Maple VMom.mw и сохраненного VMom REZ.mw различны, что осложняет дополнительную корректировку текста программы VMom REZ.mw текстовым редактором. Поэтому не следует сохранять результаты расчетов в исходный текстовый файл VMom.mw.
Воспользуйтесь поиском по сайту: