Simple beam bending

[oataga]

Hi,
I have a problem with writing a simple simulation, where a loaded .stl beam is supposed to bend in a Y direction.
Did anyone try to write something like this and could help?

I get the following error:

[ERROR]   [BeamFEMForceField(force_field)] Topology is empty.

And this is the code:

import Sofa

def createScene(node):

    node.addObject('RequiredPlugin', name='constraint', pluginName='Sofa.Component.Constraint.Projective')
    node.addObject('RequiredPlugin', name='mesh', pluginName='Sofa.Component.IO.Mesh')
    node.addObject('RequiredPlugin', name='load', pluginName='Sofa.Component.MechanicalLoad')
    node.addObject('RequiredPlugin', name='fem', pluginName='Sofa.Component.SolidMechanics.FEM.Elastic')
    node.addObject('RequiredPlugin', name='state', pluginName='Sofa.Component.StateContainer')
    node.addObject('RequiredPlugin', name='topology_constant', pluginName='Sofa.Component.Topology.Container.Constant')
    node.addObject('RequiredPlugin', name='topology_dynamic', pluginName='Sofa.Component.Topology.Container.Dynamic')
    node.addObject('RequiredPlugin', name='visual', pluginName='Sofa.Component.Visual')
    node.addObject('RequiredPlugin', name='ogl', pluginName='Sofa.GL.Component.Rendering3D')
    node.addObject('RequiredPlugin', name='mass', pluginName='Sofa.Component.Mass')
    node.addObject('RequiredPlugin', name='loop', pluginName='Sofa.Component.AnimationLoop')
    node.addObject('RequiredPlugin', name='sol', pluginName='Sofa.Component.Constraint.Lagrangian.Solver')
    node.addObject('RequiredPlugin', name='euler', pluginName='Sofa.Component.ODESolver.Backward')
    node.addObject('RequiredPlugin', name='linearowy', pluginName='Sofa.Component.LinearSolver.Iterative')
    node.addObject('FreeMotionAnimationLoop')
    node.addObject('LCPConstraintSolver', tolerance='1e-3', maxIt='10000')
    

    beam = node.addChild("beam")
    beam.addObject('EulerImplicitSolver', name='cg_solver', printLog=False, rayleighMass='0.1')
    beam.addObject('CGLinearSolver', iterations='200', name='linear solver', tolerance='1e-09', threshold='1e-09')

    beam.addObject('MeshSTLLoader', name='loader', filename='bela.stl')
    beam.addObject('MeshTopology', name='mshld',src='@loader')
    beam.addObject('PointSetTopologyContainer', src='@mshld')


    beam.addObject('MechanicalObject', name='mechanical', template='Rigid3d', position='@loader.position')
    beam.addObject('UniformMass', name='masa', totalMass='50')


    beam.addObject('BeamFEMForceField', name='force_field', template='Rigid3d', youngModulus=6.9e4, poissonRatio=0.33)

    beam.addObject('ConstantForceField', name='force', indices='0', forces='0 100 0  0 1 0 0') 
    

    beam.addObject('FixedConstraint', name='fixed_dofs', indices='0')


    beam.addObject('OglModel', name='visual', template='Vec3d', src='@loader', color='1 0 0 1')


    beam.addObject('LineAxis', name='axis', size='0.1')


    beam.addObject('DefaultAnimationLoop', name='loop')
    beam.addObject('DefaultVisualManagerLoop', name='vizloop')

    node.dt = 1 
    node.gravity = [0, 0, 0,] 

    return node

if __name__ == '__main__':
    Sofa.GUI.Application().initialize()
    Sofa.GUI.main_loop()
    root = Sofa.SimulationGraph.getRoot(node)
    root.addObject('DefaultAnimationLoop', name='animation')
    createScene(root)
    root.animate(True).run()

[hugtalbot]

Hi @oataga

Does your bela.stl file contain edges?
If so, you should use a EdgeSetTopologyContainer instead of a PointSetTopologyContainer
Cheers

[hugtalbot]

Any update @oataga ?

[github-actions[bot]]

⚠️ ⚠️ ⚠️
@oataga
We are setting a new forum management policy: topics like this one, in which the last reply is more than 4 months old, will automatically be closed. Without further news in the coming weeks, our bot will proceed to automatic closure, thus keeping this forum clean and fresh 🌱 Thank you for your understanding

[github-actions[bot]]

⚠️ ⚠️ ⚠️
@oataga
In accordance with our forum management policy, the last reply is more than 4 months old and is therefore closed. Our objective is to keep the forum up to date and offer the best support experience.

Please feel free to reopen it if the topic is still active by providing us with an update. Please feel free to open a new thread at any time - we'll be happy to help where and when we can.

[github-actions[bot]]

⚠️ ⚠️ ⚠️
@oataga
In accordance with our forum management policy, the last reply is more than 4 months old and is therefore closed. Our objective is to keep the forum up to date and offer the best support experience.

Please feel free to reopen it if the topic is still active by providing us with an update. Please feel free to open a new thread at any time - we'll be happy to help where and when we can.