Проектирование мотоустановки среднемагистрального пассажирского самолета
C if(xt.gt.0.1103.and.(fif.gt.1.2915.and.fif.lt.1.85))
C * goto 3
6 continue
r1ef=dcos(fif+dfi/2.)
r2ef=dcos(fif+dfi/2.)
if(xt.gt.0.1103.and.(fif.gt.0..and.fif.le.pi/2.))then
r1ef=1.
r2ef=1.
c write(*,*) '1 *********',fif
c pause
endif
if(xt.gt.0.1103.and.(fif.gt.pi/2.and.fif.le.pi))then
r1ef=-1.
r2ef=-1.
c write(*,*) '2 *********',fif
c pause
endif
c write(*,*) fif*180./pi,xt,r1ef,r2ef
c PAUSE' '
P(1)=(Po(i)+Py(i)*alf*dsin(fif+dfi/2.)+Pza(i)*alf*r1ef
* -Pzb(i)*bett*r2ef)*q*L*dfi*R(i)*sh
c write(*,*) P/(sh*L*R(i)*dfi),xt,fif
c pause
P(2)=Po(i)*q*L*dfi*R(i)*sh
P(3)=Py(i)*alf*dsin(fif+dfi/2.)*q*L*dfi*R(i)*sh
P(4)=Pza(i)*alf*r1ef*q*L*dfi*R(i)*sh
P(5)=-Pzb(i)*bett*r2ef*q*L*dfi*R(i)*sh
do 1000 ik=1,5
Pyc(ik)=Pyc(ik)+P(ik)*dsin(fif+dfi/2.)
Pzc(ik)=Pzc(ik)+P(ik)*dcos(fif+dfi/2.)
Myx(ik)=Myx(ik)+P(ik)*dsin(fif+dfi/2.)*xt
Mzx(ik)=Mzx(ik)+P(ik)*dcos(fif+dfi/2.)*xt
Pxc(ik)=Pxc(ik)+P(ik)*dtan(gamm(i))
1000 continue
3 continue
2 continue
xt=xt+sh
R(i)=R(i)+sh*L*dtan(gamm(i))
1 continue
do 100 ik=1,5
if(dabs(Pyc(ik)).gt.1.d-10) then
Xc(ik)=Myx(ik)/Pyc(ik)
else
write(*,*)ik,'Myx=',Myx(ik)
endif
if(dabs(Pzc(ik)).gt.1.d-10) then
Xcc(ik)=Mzx(ik)/Pzc(ik)
else
write(*,*)ik,'Mzx=',Mzx(ik)
endif
write(*,*) 'Значения аэр. сил в Ц.Д.'
100 write(*,12) Pxc(ik),Pyc(ik),Pzc(ik),Xc(ik)*L,Xcc(ik)*L
12 format(1x,'Pxc=',f10.2,/,1x,'Pyc=',f10.2/,1x,'Pzc=',f10.2/,
* 1x,'Xyc=',f15.7,/,1x,'Xzc=',f15.7)
open(1,file='aer.res')
write(1,*) 'Начальное значение угла Fi'
write(1,*) fin*180./pi
write(1,*) 'Конечное значение угла Fi'
write(1,*) fi*180./pi
write(1,*) 'Начальная Xнач'
write(1,*) Ln*5.6
write(1,*) 'Конечная Xкон'
write(1,*) Lk*5.6
write(1,*) 'Значение угла BETTA'
write(1,*) bett*180./pi
write(1,*) 'Значение угла ALFA'
write(1,*) alf*180./pi
write(1,*) 'Значение q'
write(1,*) q
write(1,*) 'Значения аэр. сил в Ц.Д.'
do 102 ik=1,5
102 write(1,12) Pxc(ik),Pyc(ik),Pzc(ik),Xc(ik)*L,Xcc(ik)*L
close(1)
stop' '
end