C****************************************************************
C***
C*** CL MAGISTRALE IN ING. MECCANICA
C***
C*** CORSO DI COSTRUZIONE DI MACCHINE
C***
C*** ANALISI DI TRANSITORIO DI TRAVE INCASTRATA ALLA BASE
C*** SOGGETTA ALLA ROTTURA DI UN CAVO ELETTRICO APPLICATO
C*** ALL'ESTREMITA' SUPERIORE
C***
C***
C*** METODO SOVRAPPOSIZIONE MODALE
C***
C***
C***-----------------------------------------------------------
C***
FINISH
/CLEAR
/CONFIG,NRES,10000	! AUMENTA IL NUMERO DI RISULTATI SCRIVIBILI SUL FILE DI OUTPUT PER LS
/PREP7
/TITLE,ANALISI DI TRANSITORIO DINAMICO DI COLONNA SOGGETTA A COLPO DI VENTO
C***
C*** PARAMETRI
C***
NEL=50
NELB=5
H=20                           ! Altezza colonna
D=1                             ! Diametro tubo
S=0.01                        ! Spessore tubo
B=2				! BRACCIO
FC=15000		! Tiro del cavo [N]
MYOUNG=2.1E11
RO=7.8E3
NRLS=1000                ! numero totale di risultati per load-step da registrare

*ASK,DMP,SMORZAMENTO :,0.005

*ASK,TIMINC,PASSO TEMPORALE DI INTEGRAZIONE : ,0.001

*ASK,NMOD,NUMERO DI MODI DA CONSIDERARE:,8
C***
C*** NODI
C***
N,1
N,NEL+1,0,H
FILL,1,NEL+1
N,NEL+NELB+1,0,H,B
FILL,NEL+1,NEL+NELB+1
N,NEL+NELB+2,0,H,-B/NELB
N,NEL+2*NELB+1,0,H,-B
FILL,NEL+NELB+2,NEL+2*NELB+1

C***
C*** ELEMENTI
C***
ET,1,4
A=3.14*D**2-3.14*(D-2*S)**2
J=3.14/64*D**4-3.14/64*(D-2*S)**4
R,1,A,J,J,2*J
E,1,2
EGEN,NEL+NELB,1,1,1
E,NEL+1,NEL+NELB+2
E,NEL+NELB+2,NEL+NELB+3
EGEN,NELB-1,1,-1

C***
C*** PROPRIETA' MATERIALE
C***
MP,EX,1,MYOUNG
MP,DENS,1,RO       	! densità
MP,DAMP,1,DMP	! smorzamento

C***
C*** VINCOLI
C***
D,1,ALL,0
FINISH
C*** 
C*** SOLUZIONE
C***
/SOLU
ANTYPE,MODAL            ! ANALISI MODALE PRELIMINARE
MODOPT,QRDAMP,NMOD 	! OPZIONE NECESSARIA PER INCLUDERE LO SMORZAMENTO
SOLVE
SAVE
FINISH
/eof

/SOLU

= = = DA COMPLETARE = = =

C***
C*** POST-PROCESSAMENTO
C***
/POST26
NSOL,2,NBRAC,U,X,SPOST
NSOL,3,NBRAC,ROT,Y,ROT
STORE,NEW
PLVAR,2,3
*ASK,IFL,RETURN PER CONTINUARE: ,0

/POST1
SET
/VIEW,,1,1,1
/dscale,,1
PLDISP,1
*ASK,IFL,RETURN PER CONTINUARE: ,0

ANTIME,100,0.01,,2,2

ANIM,10,1,0.01