SolidsPy_dyn: 2D Dynamic Finite Element Analysis with Python

This repo contains the dynamic version of the finite element analysis code SOLIDSPy. The dynamic code uses an implicit time integrator scheme.

Features

The code allows the user to find the displacement time histories for all the points in the mesh. The response is written into ready to process Paraview VTKs files.

The repo contains 3 main folders:

1. solidspydyn/ which stores the source code in the following modules:

   • solids_GUI.py: The main program;
   • preprocesor.py: Pre-processing subroutines including Gmsh convertion functions using meshio
   • assemutil.py: Assembly of elemental stiffness matrices ;
   • femutil.py: Shape functions, its derivatives and general finite element method subroutines;
   • uelutil.py: Elemental or local matrix subroutines for different elements; and
   • postprocesor.py: Several results handling subroutines.

2. meshes/ Complete models including its Gmsh representation and a Python script to produce the required (nodes, elements, materials and load) text files ready for input.

3. docs/ Documentation files including a User Manual.

Installation

The code is written in Python and it depends on numpy, scipy and sympy. To use the code just download it and store it in the desired folder.

To run the examples with specification of the folder storing the input files through a GUI you will need to install easygui.

To easily generate the required SolidsPy_dyn text files out of a Gmsh model you will need meshio.

These two can be installed with:

pip install easygui
pip install meshio

How to run a simple model

After downloading, you can run an analysis in 3 easy steps:

• Create the model (i.e., geometry and mesh) using Gmsh.
• Generate the text files (eles.txt, nodes.txt, mater.txt and loads.txt) required by SolidsPydyn using a python script based on meshio.
• Execute the script named ejecucion_basico.py which is available in the main folder.

Youtube videos with complete step-by-step explanations are available in the Applied Mechanics channel.

License

This project is licensed under the MIT license. The documents are licensed under

Creative Commons Attribution License.

Authors

• Juan Gomez, Professor at Universidad EAFIT.
• Nicolás Guarín-Zapata, Researcher at the Applied Mechanics Group at Universidad EAFIT.