C*************************************************************
!*
!*  MODELLO VOLANTE CON  ELEMENTI TRAVE (razze in acciaio, corona in alluminio)
!* 
!*  FEATURES : SMORZAMENTO PROPORZIONALE
!*             USO DELLA VARIABILI DIMENSIONATE
!*             CONFONTO DATI NUMERICI E SPERIMENTALI
!* 
!* 
FINISH
/CLEAR
/PREP7
C*****************************************************************
C***
C*** PARAMETRI
C***
NFS=1		! PASSO DI FREQUENZA
F1=15		! FREQUENZE CONTENUTE NELL'ECCITAZIONE
F2=29
F3=57
F4=70		! FREQUENZA FINALE DELL'ANALISI
*ASK,BET,MOLTIPLICATORE MATRICE RIGIDEZZA :,0.001
*ASK,ALF,MOLTIPLICATORE MATRICE MASSA     :,1

C**** SPOSTAMENTO IN INGRESSO PUNTO 1

*DIM,SP1,ARRAY,3,3,2
SP1(1,1,1)=-6E-5	! 15 Hz
SP1(1,1,2)= 1.4E-5
SP1(1,2,1)=4E-6
SP1(1,2,2)=1.8E-5
SP1(1,3,1)=-1.3E-5
SP1(1,3,2)=-4.087E-6

SP1(2,1,1)=-6.287E-8	! 29 Hz
SP1(2,1,2)= -9.014E-7
SP1(2,2,1)=-4.128E-6
SP1(2,2,2)=-2.486E-6
SP1(2,3,1)=1.092E-6
SP1(2,3,2)=5.097E-7

SP1(3,1,1)=7.425E-8	! 57 Hz
SP1(3,1,2)= -2.378E-8
SP1(3,2,1)=4.668E-7
SP1(3,2,2)=-2.827E-7
SP1(3,3,1)=-2.132E-7
SP1(3,3,2)=-9.625E-8


C**** SPOSTAMENTO IN USCITA PUNTO 4 E PUNTO 5

*DIM,SP4X,ARRAY,F4/NFS+3
*DIM,SP4Y,ARRAY,F4/NFS+3
*DIM,SP4Z,ARRAY,F4/NFS+3
*DIM,SP5X,ARRAY,F4/NFS+3
*DIM,SP5Y,ARRAY,F4/NFS+3
*DIM,SP5Z,ARRAY,F4/NFS+3
*DO,IJK,1,F4/NFS+3,1
		SP4X(IJK)=0
		SP5X(IJK)=0
		SP4Y(IJK)=0
		SP5Y(IJK)=0
		SP4Z(IJK)=0
		SP5Z(IJK)=0
*ENDDO

SP4X(15)=8.2E-5
SP4X(29)=3.4E-5
SP4Y(15)=1.8E-5
SP4Y(29)=8.4E-6
SP4Z(57)=7.0E-7
SP4Z(57)=2.6E-6

SP4X(14)=8.2E-5
SP4X(28)=3.4E-5
SP4Y(14)=1.8E-5
SP4Y(28)=8.4E-6
SP4Z(56)=7.0E-7
SP4Z(56)=2.6E-6

SP4X(16)=8.2E-5
SP4X(30)=3.4E-5
SP4Y(16)=1.8E-5
SP4Y(30)=8.4E-6
SP4Z(58)=7.0E-7
SP4Z(58)=2.6E-6


SP5X(14)=7.4E-5
SP5Y(14)=2.1E-5
SP5Y(28)=7.5E-6
SP5Y(56)=1.0E-6
SP5Z(14)=1.7E-5
SP5Z(28)=4.8E-6
SP5Z(56)=1.5E-6

SP5X(15)=7.4E-5
SP5Y(15)=2.1E-5
SP5Y(29)=7.5E-6
SP5Y(57)=1.0E-6
SP5Z(15)=1.7E-5
SP5Z(29)=4.8E-6
SP5Z(57)=1.5E-6

SP5X(16)=7.4E-5
SP5Y(16)=2.1E-5
SP5Y(30)=7.5E-6
SP5Y(58)=1.0E-6
SP5Z(16)=1.7E-5
SP5Z(30)=4.8E-6
SP5Z(58)=1.5E-6


C*** ==========================================================================================
C***
C*** INIZIO DEFINIZIONE MODELLO
C***
!*
!* definiz elementi
!*   
ET,1,BEAM4  
!*  
!* definiz real costants 
!* 
C*** CORONA
R,1,92E-6,1.55e3*1E-12,1.04e3*1E-12,10E-3,14E-3, , 
RMORE, ,291E-12, , , , ,
!*  
R,2,140E-6,2287E-12,1167E-12,10E-3,14E-3,0,
RMORE,0,2658E-12,0,0,0,0, 

R,3,140E-6,2287E-12,1167E-12,10E-3,14E-3,0,
RMORE,0,2658E-12,0,0,0,0,   
!*
C*** RAZZE LATERALI  
R,4,342E-6,2.109e5*1E-12,456E-12,4E-3,90E-3,0,   
RMORE,0,1824E-12,0,0,0,0,   
!*  
R,5,272E-6,1.146e5*1E-12,362.67E-12,4E-3,80E-3,0, 
RMORE,0,1451E-12,0,0,0,0,   
!*  
R,6,200E-6,4.347e4*1E-12,266.67E-12,4E-3,54E-3,0, 
RMORE,0,1067E-12,0,0,0,0,   
!*  
R,7,148E-6,1.87E-12,197.33E-12,4E-3,44E-3,0,  
RMORE,0,789.33E-12,0,0,0,0,
!*  
R,8,128E-6,1.105e4*1E-12,170.67E-12,4E-3,34E-3,0, 
RMORE,0,682.67E-12,0,0,0,0,
!*  
R,9,188E-6,4.255e4*1E-12,250.67E-12,4E-3,60E-3,0, 
RMORE,0,1003E-12,0,0,0,0,   
!*  
R,10,1E-6,7.2e4*1E-12,320*1E-12,4E-3,60E-3,0, 
RMORE,0,1280E-12,0,0,0,0, 

C*** RAZZA CENTRALE
!*  
R,11,420E-6,3.885e5*1E-12,560E-12,4E-3,110E-3, ,
RMORE, ,2.24e3*1e-12, , , , , 
!*  
R,12,320E-6,2.027e5*1E-12,426.67E-12,4E-3,100E-3,0,
RMORE,0,1707E-12,0,0,0,0,   
!*  
R,13,184E-6,4.146e4*1E-12,245.33E-12,4E-3,60E-3,0,
RMORE,0,981.33E-12,0,0,0,0,
!*  
R,14,124E-6,9961E-12,165.33E-12,4E-3,32E-3,0,   
RMORE,0,661.33E-12,0,0,0,0,
!*  
R,15,156E-6,2.293e4*1E-12,208E-12,4E-3,48E-3,0,
RMORE,0,832E-12,0,0,0,0,
!*  
R,16,216E-6,5.443e4*1E-12,288*1E-12,4E-3,60E-3,0,
RMORE,0,1.152e3*1E-12,0,0,0,0,
!*  
R,17,1E-6,7.2e4*1E-12,320*1E-12,4E-3,60E-3,0,
RMORE,0,1280e-12,0,0,0,0, 

C*** TELAIO RIGIDO SIMULANTE IL PIANTONE

R,18,1e-1,1e-3,1e-3,1e-3,1e-3,0,
RMORE,0,1e-3,0,0,0,0,   

C*** CORONA IN POLIMERO

A1=3.14*(0.03**2)/4
J1=3.14*(0.03**4)/64
R,19,A1,J1,J1,0.03,0.03,0,
RMORE,0,2*J1,0,0,0,0,   
!*  
!*  definiz proprietà materiali: 
!*  materiale 1 = materiale corona = alluminio
!*  materiale 2 = materiale razze = acciaio
!* 
MP,DENS,1,2.7e3   
MP,EX,1,7E10 
MP,GXY,1,7E10/2.6
MP,EX,2,2.1E11
MP,GXY,2,2.1E11/2.6 
MP,DENS,2,7.8e3

MP,EX,3,5.2E6
MP,GXY,3,5.2E6/2.6 
MP,DENS,3,4.9e2

!*  
!*  definiz geometria
!* 
K,1,55E-3,0,0, 
K,2,63E-3,0,0, 
K,3,87E-3,0,23.1E-3,  
K,4,103E-3,0,38.5E-3, 
K,5,135E-3,0,69.3E-3,  
K,6,143E-3,0,77E-3,   
K,7,156E-3,0,77E-3,   
K,8,168E-3,0,77E-3,   
K,9,-55E-3,0,0,
K,10,-63E-3,0,0,   
K,11,-87E-3,0,23.1E-3,
K,12,-103E-3,0,38.5E-3,   
K,13,-135E-3,0,69.3E-3,   
K,14,-143E-3,0,77E-3, 
K,15,-156E-3,0,77E-3, 
K,16,-168E-3,0,77E-3,
K,17,0,-52E-3,0,   
K,18,0,-58E-3,0,   
K,19,0,-82E-3,18E-3,  
K,20,0,-114E-3,42E-3, 
K,21,0,-134E-3,77E-3, 
K,22,0,-150E-3,77E-3, 
K,23,0,-156E-3,77E-3, 
K,24,0,-168E-3,77E-3, 
K,25,0,168E-3,77E-3,  
K,26,164.5E-3,34E-3,77E-3,   
K,27,160E-3,-52E-3,77E-3,
K,28,-164.5E-3,34E-3,77E-3,  
K,29,-160E-3,-52E-3,77E-3,   
K,30,34E-3,-164.5E-3,77E-3,  
K,31,-34E-3,-164.5E-3,77E-3,
K,32,,,77E-3, 
!*  
   

K,109,-40.5E-3
K,110,-40.5E-3,17.5E-3
K,111,40.5E-3,17.5E-3
K,112,40.5E-3,
K,113,40.5E-3,-17.5E-3
K,114,0,-17.5E-3
K,115,-40.5E-3,-17.5E-3

!*  
LARC,26,25,32,168E-3,  
LARC,25,28,32,168E-3,
LARC,29,31,32,168E-3, 
LARC,30,27,32,168E-3, 
CM,CORONA1,LINES

LSEL,NONE
LARC,26,25,32,168E-3,  
LARC,25,28,32,168E-3,
LARC,29,31,32,168E-3, 
LARC,30,27,32,168E-3, 
CM,CORPOL1,LINES

LSEL,NONE
LARC,28,16,32,168E-3,  
LARC,16,29,32,168E-3,  
LARC,8,26,32,168E-3, 
LARC,27,8,32,168E-3,
CM,CORONA2,LINES

LSEL,NONE
LARC,28,16,32,168E-3,  
LARC,16,29,32,168E-3,  
LARC,8,26,32,168E-3, 
LARC,27,8,32,168E-3,
CM,CORPOL2,LINES

!*  
LSEL,NONE
LARC,31,24,32,168E-3,  
LARC,24,30,32,168E-3, 
CM,CORONA3,LINES 

LSEL,NONE
LARC,31,24,32,168E-3,  
LARC,24,30,32,168E-3, 
CM,CORPOL3,LINES 
!*  

*DO,IJK,1,7
	LSEL,NONE
	L,IJK,IJK+1
	L,IJK+8,IJK+9
	CM,RAZZA%IJK%,LINES
*ENDDO
*DO,IJK,17,23
	LSEL,NONE
	L,IJK,IJK+1
	CM,RAZZA%IJK%,LINES
*ENDDO

LSEL,NONE
*DO,IJK,109,114
	L,IJK,IJK+1
*ENDDO
L,115,109
L,9,109
L,112,1
L,114,17
CM,PIANTONE,LINES

C*** INTRODUZIONE ELEMENTI
*DO,IJK,1,3
	CMSEL,,CORONA%IJK%
	REAL,IJK
	MAT,1
	LESIZE,ALL,,,3
	LMESH,ALL
*ENDDO
*GET,NEMX,ELEM,0,NUM,MAX
EGEN,2,0,ALL
ESEL,,ELEM,,NEMX+1,999
EMODIF,ALL,MAT,3
EMODIF,ALL,REAL,19
ALLS

*DO,IJK,1,7
	CMSEL,,RAZZA%IJK%
	REAL,IJK+3
	MAT,2
	LESIZE,ALL,,,1
	LMESH,ALL
*ENDDO
*DO,IJK,17,23
	CMSEL,,RAZZA%IJK%
	REAL,IJK-6
	MAT,2
	LESIZE,ALL,,,1
	LMESH,ALL
*ENDDO

CMSEL,,PIANTONE
REAL,18
LESIZE,ALL,,,1
LMESH,ALL


ALLS
C***
C*** FINE DEFINIZIONE MODELLO
C***
C*** ===============================================================================

!* 
!* ricerca soluzione
!*  
  
/SOLU  



!*  
ANTYPE,HARMIC
HROPT,FULL
ALPHAD,ALF
BETAD,BET

KBC,1

C***
C*** 0-15 HZ
C***
DK,111,ALL,0
DK,110,ALL,0
DK,113,ALL,0
DK,115,ALL,0

HARFRQ,0,F1
NSUBST,(F1-1)/NFS
SOLVE

C***
C*** 15 HZ
C***
DKDELE,ALL,ALL
KSEL,,KP,,110,111
KSEL,A,KP,,113 
KSEL,A,KP,,115  
DK,ALL,UX,SP1(1,1,1),SP1(1,1,2)
DK,ALL,UY,SP1(1,2,1),SP1(1,2,2)
DK,ALL,UZ,SP1(1,3,1),SP1(1,3,2)
KSEL,ALL

HARFRQ,F1,F1
NSUBST,1
SOLVE


= = = = = DA COMPLETARE = = = = =


/post26


NSEL,,LOC,X,164.499E-3,164.501E-3
NSEL,R,LOC,Y,33.999E-3,34.001E-3
NSEL,R,LOC,Z,76.999E-3,77.001E-3
*GET,N4,NODE,,NUM,MAX

NSEL,,LOC,X,-0.001E-3,0.001E-3
NSEL,R,LOC,Y,-168.001E-3,167.999E-3
NSEL,R,LOC,Z,76.999E-3,77.001E-3
*GET,N5,NODE,,NUM,MAX

alls
numvar,200
nsol,2,N4,u,X,X4
nsol,3,N4,u,Y,Y4
nsol,4,N4,u,Z,Z4

nsol,5,N5,u,X,X5
nsol,6,N5,u,Y,Y5
nsol,7,N5,u,Z,Z5

store

VPUT,SP4X,8,,0,X4EXP	! TRASFERISCE NELLE VARIABILI I CONTENUTI DEI VETTORI
VPUT,SP4Y,9,,0,Y4EXP
VPUT,SP4Z,10,,0,Z4EXP

VPUT,SP5X,11,,0,X5EXP
VPUT,SP5Y,12,,0,Y5EXP
VPUT,SP5Z,13,,0,Z5EXP

PLCPLX,0		! sceglie tra rappresentare il modulo (0)
			! o la fase (1) o altro
PLVAR,2,8
*ASK,DD,RETURN PER PROSEGUIRE,0
PLVAR,3,9
*ASK,DD,RETURN PER PROSEGUIRE,0
PLVAR,4,10
*ASK,DD,RETURN PER PROSEGUIRE,0
PLVAR,5,11
*ASK,DD,RETURN PER PROSEGUIRE,0
PLVAR,6,12
*ASK,DD,RETURN PER PROSEGUIRE,0
PLVAR,7,13