Entities creation from scratch
Here we present the interface in order to create new entities in an empty model from Python.
Creating nodes
We can create a node by doing. If we Try to create a new node with the same Id and different coordinates we would get an error.
main_model_part.CreateNewNode(1, 1.0,0.0,0.0)
#main_model_part.CreateNewNode(1, 0.0,0.0,0.0) # Here an error is thrown
However if we try to create a node with the same coordinates twice nothing is actually done (and no error is thrown)
main_model_part.CreateNewNode(1, 1.00,0.00,0.00)
print(main_model_part.NumberOfNodes()) # This still returns 1!!
We can now access the node as needed, for example:
print(main_model_part.GetNode(1).Id) # Gives 1
print(main_model_part.GetNode(1,0).X) # Gives 1.0
Nodes
can be created in every order
main_model_part.CreateNewNode(2000, 2.00,0.00,0.00)
main_model_part.CreateNewNode(2, 2.00,0.00,0.00)
We could then loop over all the nodes:
for node in main_model_part.Nodes:
print(node.Id, node.X, node.Y, node.Z)
Or eventually remove nodes one by one by doing:
main_model_part.RemoveNode(2000)
Let’s now see what happens if we add a node to a submodelpart. Here the node will be both in root ModelPart
and “BC”, but for example not in derived submodelparts from “BC” or root ModelPart
.
bc_model_part = main_model_part.GetSubModelPart("BC")
bc_model_part.CreateNewNode(3, 3.00,0.00,0.00)
Multiple nodes can be removed at once (and from all levels) by flagging them:
for node in main_model_part.Nodes:
if node.Id < 3:
node.Set(KratosMultiphysics.TO_ERASE,True)
main_model_part.RemoveNodesFromAllLevels(KratosMultiphysics.TO_ERASE)
One could call simply the function RemoveNodes
and remove them from the current level down.
Creating Elements and Conditions
Elements
and Conditions
can be created from the python interface by providing their connectivity as well as the Properties
to be employed in the creation. The string to be provided is the name by which the element is registered in Kratos. An error is thrown if i try to create an element with the same Id
main_model_part.AddProperties(KratosMultiphysics.Properties(1))
main_model_part.CreateNewElement("Element2D3N", 1, [1,2,3], main_model_part.GetProperties()[1])
An identical interface is provided for Conditions, as well as functions equivalent to the nodes for removing from one level or from all the levels.
Creating MasterSlaveConstraints
The case of the constraints, like it requires the existance of DoF
we will create a new model part in order to avoid problems. This example will advance some methods that will be introduced in future sections, as AddNodalSolutionStepVariable
. The current interface for python works directly with nodes, it is not possible to manipulate DoF
manually.
mp = this_model.CreateModelPart("constraint_example")
mp.AddNodalSolutionStepVariable(KratosMultiphysics.DISPLACEMENT)
mp.AddNodalSolutionStepVariable(KratosMultiphysics.REACTION)
mp.CreateNewNode(1, 0.00000, 0.00000, 0.00000)
mp.CreateNewNode(2, 0.00000, 1.00000, 0.00000)
KratosMultiphysics.VariableUtils().AddDof(KratosMultiphysics.DISPLACEMENT_X, KratosMultiphysics.REACTION_X, mp)
mp.CreateNewMasterSlaveConstraint("LinearMasterSlaveConstraint", 1, mp.Nodes[1], KratosMultiphysics.DISPLACEMENT_X, mp.Nodes[2], KratosMultiphysics.DISPLACEMENT_X, 1.0, 0)