Code:
!------------------------------------------------------------------------------
! SOLVER FOR TRANSIENT-DYNAMIC CALCULATION
!------------------------------------------------------------------------------FINISH
/SOLU ! start solver
!-----Belastung-------
!------Voreinstellungen Dynamik (Zeit)------
!=======================================================================
dTime=(EndTime-IniTime)/nrTimeSteps
EndFX = -Axialkraft ! Final axial load
EndFY = Radialkraft ! Final radial load
EndRotat = Drehwinkel ! Final rotation angle
!------FØlle Spaltenvektoren mit Zeitwerten und Lastschritten aus abgefahrener Kurve------
!==============================================================================================
*DO, ii, 1, nrTimeSteps+1
valsTime(ii) = IniTime + (ii-1)*dTime
valsFX(ii) = -FVorspAx !Axialkraft
valsFY(ii) = FVorspRad+deltaFRad*SIN(2*Pi*freL*valsTime(ii)) !Radialkraft
valsRotX(ii) = (EndRotat)*SIN(2*Pi*freD*valsTime(ii)) !Rotation
*ENDDO
!-----Berechnung Vorspannung--------
TIMINT, OFF !statische Berechnung
write_log, 'INFO','Vorspannung::Start','--','--','--','--',0
TIME, 0.0001
F, PilotNode, FX, -FVorspAX/100
F, PilotNode, FY, FVorspRad/100
NSUBST, 20, 20, 10
SOLVE
TIME, 0.0002
F, PilotNode, FX, -FVorspAx
F, PilotNode, FY, FVorspRad
NSUBST, 40, 50, 20
SOLVE
save
!/EOF
!-----Berechnung Belastung--------
TIMINT, ON !Turns on transient effects
!*GET, UX_reac, NODE, PilotNode, U, X
!*GET, UY_reac, NODE, PilotNode, U, Y
!Time, 0.0003
!D, PilotNode, UX, UX_reac
!D, PilotNode, UY, UY_reac
!NSUBST, 1, 2, 1
!SOLVE
*DO, ik, 1, nrTimeSteps+1
TIME, dTime*ik
F, PilotNode, FX, forceFX(dTime*ik)
F, PilotNode, FY, forceFY(dTime*ik)
D, PilotNode, ROTX, rotatX(dTime*ik)
NSUBST, Substeps, nint(Substeps*SubstepFaktor), nint(Substeps/SubstepFaktor), ON
SOLVE
*ENDDO