File CRSM.F4

Directory of image this file is from
This file as a plain text file

	REAL IL,ICM,IR
	DIMENSION D(2),RI(2),SKL(2),X(3,3),ST(3,3),FRC(3),PM(3,3)
	2,SIG(2,4),PER(2),GAM(7),PHI(3)
	PI=3.1415926
	DTR=PI/180.
	RTD=1.0/DTR
	A=3.25
	B=1.0
	ICM=.0039
	D(1)=.375
	D(2)=.375
	RI(1)=(PI*D(1)**4)/64.
	RI(2)=(PI*D(2)**4)/64.
	E=.285E8
	RL=14.
	AM1=.286*PI*RL*D(1)**2/(16.*386)
	AM2=.286*PI*RL*D(2)**2/(16.*386)
	IL=AM1+AM2+.615/386.
C------COMPUTE THE COMPONENTS OF THE SOLENOID FORCES
	FR=.25
	FX=0.0
	FY=0.0
	PHI(1)=14.24
	PHI(2)=24.26
	PHI(3)=0.0
	GAM(1)=31.5781*DTR
	GAM(2)=330.9301*DTR
	GAM(6)=91.0496*DTR
	GAM(7)=270.0*DTR
	GAM(3)=63.44*DTR
	GAM(5)=299.1013*DTR
	GAM(4)=0.0
	II=1
	DO 50 I=1,7
	IF(I.GT.4) II=2
	IF(I.GT.6) II=3
	FX=FX+FR*COS(PHI(II)*DTR)
50	FY=FY+FR*SIN(PHI(II)*DTR)*COS(GAM(I))
	WRITE(3,104)
104	FORMAT(1H1,'  CLAMPED-PINNED',//)
C------LOOP ON POSITION
	DO 1 IX=30,70,10
	XX=.1*FLOAT(IX)
	DO 51 JJ=1,2
	PK=(6.*RL*E*RI(JJ))/((RL-XX)*(2.*RL*XX**2 - 2.*XX**3))
	CK=(3.*E*RI(JJ)*RL**3)/((XX**3)*((RL-XX)**3))
51	SKL(JJ)=.5*(PK+CK)
	WNL=SQRT(SKL(2)/IL)
	PER(1)=1.0/WNL
	X1=(A-B)*SKL(2)/(SKL(1)+SKL(2))
	X2=(A-B)*SKL(1)/(SKL(1)+SKL(2))
	R0=B+X2
	WRITE(3,100) XX,(D(I),RI(I),SKL(I),I=1,2)
100	FORMAT(//,'  CARR. POSITION : ',F10.4,/,T6,'DIA(1)',T18,
	2'MOM(1)',T30,'K(1)',T42,'DIA(2)',T54,'MOM(2)',T66,'K(2)',
	3/,6F12.5)
	WRITE(3,101)
101	FORMAT(//,T10,'SPEED',T25,'PERIOD',T40,'MAX X',T55,'MAX Y',
	2T70,'MAX ANG',/,T10,'(IPS)',T25,'(SEC)',T40,'(IN)',T55,
	3'(IN)',T70,'(DEG)',/)
C-------LOOP ON CARRIAGE SPEED
	DO 2 ISP=4,8,2
	SP=FLOAT(ISP)
	PR=.01/SP
C------INITIALIZE VARIABLES
	N=0
	DO 31 I=1,4
	SIG(2,I)=0.0
	DO 31 J=1,3
	PM(I,J)=0.0
	X(I,J)=0.0
31	ST(I,J)=0.0
	PM(1,1)=1.0
	PM(2,2)=1.0
	PM(2,3)=X1+.5
	PM(3,3)=RTD
	R=R0
	DT=.05*PER(1)
C------LOOP OVER TIME
	DO 3 IT=1,2001
	T=FLOAT(IT-1)*DT
C------DETERMINE THE MAGNITUDE OF THE FORCE
	FCTR=0.0
	TP=1000.*(T-N*PR)
	IF(TP) 6,6,7
7	IF(TP-1.2) 4,4,5
4	FCTR=(1.2-TP)/.2
	IF(TP.LT.1.) FCTR=(1.-(1.-TP)/.2)
	IF(TP.LT..8) FCTR=0.
	IF(TP.LT..6) FCTR=(.6-TP)/.1
	IF(TP.LT..5) FCTR=1.0
	IF(TP.LT..3) FCTR=(1.-(.3-TP)/.3)
	GO TO 6
5	N=N+1
C------DETERMINE THE INSTANTANEOUS ROT. INERTIA, STIFFNESS, FORCE
6	IR=ICM+IL*R**2+AM1*X1**2+AM2*X2**2
	SKA=SKL(1)*X1**2 + SKL(2)*X2**2
	ST(1,1)=-SKL(2)/IL
	ST(2,2)=-(SKL(1)+SKL(2))/IL
	ST(3,3)=-SKA/IR
	FRC(1)=-FCTR*FX/IL
	FRC(2)=-FCTR*FY/IL
	FRC(3)=-FCTR*FX*.8/IR
	PM(2,3)=X1+.5
C------PERFORM OPERATION D2XDT2 = (-K * X + F)/I
	DO 40 I=1,3
	X(I,1)=0.0
	DO 41 J=1,3
41	X(I,1)=X(I,1)+ST(I,J)*X(J,3)
40	X(I,1)=X(I,1)+FRC(I)
	DO 42 I=1,3
	DO 42 J=1,2
42	X(I,J+1)=X(I,J+1)+X(I,J)*DT
	DO 45 I=1,3
	SIG(1,I)=0.0
	DO 46 J=1,3
46	SIG(1,I)=SIG(1,I)+PM(I,J)*X(J,3)
45	IF(ABS(SIG(1,I)).GT.SIG(2,I)) SIG(2,I)=ABS(SIG(1,I))
	R=0.0
	IF(X(3,3)) 43,44,43
43	R=-X(2,3)/X(3,3)
44	IF(ABS(R).GT.R0) R=SIGN(R0,R)
	X1=A-R
	X2=B-R
3	CONTINUE
2	WRITE(3,47) SP,PR,(SIG(2,I),I=1,3)
47	FORMAT(5F15.5,/)
1	CONTINUE
	END



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search