Результаты расчета влияют на выбор двигателей разгрузки гиростабилизатора
Последовательность расчета возмущающих моментов
С помощью программы 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.
Воспользуйтесь поиском по сайту: