<!-- GENERATED by tools/generators/examplesDocsEmitter.py from python/ - do not edit -->
(examples-ancfcablecantilevered)=
# ANCFcableCantilevered.py

You can view and download this file on Github: [ANCFcableCantilevered.py](https://github.com/jgerstmayr/EXUDYN/tree/master/python/Examples/ANCFcableCantilevered.py)

```{code-block} python
:linenos:

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# This is an EXUDYN example
#
# Details:  ANCFCable test with many cable elements
#
# Author:   Johannes Gerstmayr
# Date:     2023-10-15
#
# Copyright:This file is part of Exudyn. Exudyn is free software. You can redistribute it and/or modify it under the terms of the Exudyn license. See 'LICENSE.txt' for more details.
#
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

import exudyn as exu
from exudyn.utilities import * #includes itemInterface and rigidBodyUtilities
import exudyn.graphics as graphics
from exudyn.beams import *

import numpy as np

#create an environment for mini example
SC = exu.SystemContainer()
mbs = SC.AddSystem()

oGround=mbs.AddObject(ObjectGround(referencePosition= [0,0,0]))
nGround = mbs.AddNode(NodePointGround(referenceCoordinates=[0,0,0]))

np.random.seed(0) #reproducible results (also for SolutionViewer when reloading ...)

rhoA = 78.
EA = 100000.
EI = 200

nCables = 10000
for i in range(nCables):
    p0 = np.random.rand(3)
    if i < nCables/2:
        p0[0]=0.5
        p1 = p0 + [2+np.random.rand()*0.25,0,0]
    else:
        p0[0]=-0.5
        p1 = p0 - [2+np.random.rand()*0.25,0,0]
    

    cable = ObjectANCFCable(massPerLength=rhoA, 
                  bendingStiffness=EI*(0.5+np.random.rand()*0.25), 
                  bendingDamping = EI,
                  axialStiffness=EA, 
                  )

    ancf=GenerateBeamElementsAlongLine(mbs=mbs,
                  positionStart=p0, positionEnd=p1,
                  numberOfElements=12, #converged to 4 digits
                  beamTemplate=cable, #this defines the beam element properties
                  gravity = [0,-9.81,0],
                  groundConstraintsStart = [1,1,1, 0,1,1], #add constraints for pos and rot (r'_y,r'_z)
                  )

#assemble and solve system for default parameters
mbs.Assemble()

endTime=10
stepSize = 0.5e-3

simulationSettings = exu.SimulationSettings()

#simulationSettings.solution.file.write = False
simulationSettings.solution.file.writePeriod = 0.02 #data not used
simulationSettings.solution.file.binary = True
simulationSettings.solution.precision = 6 #float
simulationSettings.solution.sensors.writePeriod = 0.002 #data not used
simulationSettings.timeIntegration.verboseMode = 1 #turn off, because of lots of output
simulationSettings.linearSolver.solverType = exu.LinearSolverType.EigenSparse
simulationSettings.parallel.numberOfThreads = 8
simulationSettings.show.computationTime = True
simulationSettings.show.statistics = True

simulationSettings.timeIntegration.numberOfSteps = int(endTime/stepSize)
simulationSettings.timeIntegration.endTime = endTime
simulationSettings.timeIntegration.newton.useModifiedNewton = True

#simulationSettings.timeIntegration.realtime.active = True
#simulationSettings.timeIntegration.realtime.factor = 0.5

SC.visualizationSettings.general.graphicsUpdateInterval = 0.02
SC.visualizationSettings.view0.window.renderWindowSize=[1200,1024]
SC.visualizationSettings.nodes.show = False
SC.visualizationSettings.loads.show = False
SC.visualizationSettings.connectors.show = False
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++

SC.renderer.Start()
#SC.renderer.DoIdleTasks()

mbs.SolveDynamic(simulationSettings)

SC.renderer.DoIdleTasks()
SC.renderer.Stop() #safely close rendering window!

#mbs.SolutionViewer()
```