nodes.add(2, [-b, 0.]) #Left support nodes.add(3, [b, 0.]) #Right support nodes.add(4, [0., h]) #Loading point #ElementSet from pyfem.fem.ElementSet import ElementSet from pyfem.elements.Truss import Truss from pyfem.elements.Spring import Spring elements = ElementSet(nodes) elements.add(1, Truss([2, 4])) elements.add(2, Truss([3, 4])) elements.add(3, Spring([1, 4])) #Add groups elements.addGroup('Trusses', [1, 2]) elements.addGroup('Spring', [3]) #DofSpace from pyfem.fem.DofSpace import DofSpace dofs = DofSpace(elements) dofs.constrain(1, ['u', 'v']) dofs.constrain(2, ['u', 'v']) dofs.constrain(3, ['u', 'v']) ################################### # Store in global data dictionary # ###################################