Skip to main content

Simulation and FEA environment for Python built upon Calculix and GMSH

Project description

https://github.com/drlukeparry/pyccx/workflows/Python%20application/badge.svg Documentation Status https://badge.fury.io/py/PyCCX.svg Chat on Gitter

Provides a library for creating and running 3D FEA simulations using the opensource Calculix FEA Package.

The aims of this project was to provide a simple framework for implemented 3D FEA Analysis using the opensource Calculix solver. The analysis is complimented by use of the recent introduction of the GMSH-SDK , an extension to GMSH to provide API bindings for different programming languages by the project authors to provide sophisticated 3D FEA mesh generation outside of the GUI implementation. This project aims to provide an integrated approach for generating full 3D FEA analysis for use in research, development and prototyping in a Python environment. Along with setting up and processing the analysis, convenience functions are included.

The inception of this project was a result of finding no native Python/Matlab package available to perfom full non-linear FEA analysis of 3D CAD models in order to prototype a concept related to 3D printing. The project aims to compliment the work of the PyCalculix project, which currently is limited to providing capabilities to generate 2D Meshes and FEA analysis for 2D planar structures. The potential in the future is to provide a more generic extensible framework compatible with different opensource and commercial FEA solvers (e.g. Abaqus, Marc, Z88, Elmer).

An interface that built upon GMSH was required to avoid the use of the GUI, and the domain specific .geo scripts. Learn more.

Structure

PyCCX framework consists of classes for specifying common components on the pre-processing phase, including the following common operations:

  • Mesh generation

  • Creating and applying boundary conditions

  • Creating load cases

  • Creating and assigning material models

  • Performing the simulation

In addition, a meshing class provides an interface with GMSH for performing the meshing routines and for associating boundary conditions with the elements/faces generated from geometrical CAD entities. The Simulation class assembles the analysis and performs the execution to the Calculix Solver. Results obtained upon completion of the analysis can be processed. Currently the analysis is unit-less, therefore the user should ensure that all constant, material paramters, and geometric lengths are consistent - by default GMSH assumes ‘mm’ units.

Current Features

Meshing:

  • Integration with GMSH for generation 3D FEA Meshes (Tet4, Tet10 currently supported)

  • Merging CAD assemblies using GMSH

  • Attaching boundary conditions to Geometrical CAD entities

FEA Capabilities:

  • Boundary Conditions (Acceleration, Convection, Fixed Displacements, Forces, Fluxes, Pressure, Radiation)

  • Loadcase Types (Structural Static, Thermal, Coupled Thermo-Mechanical)

  • Materials (Non-linear Elasto-Plastic Material)

Results Processing:

  • Element and Nodal Results can be obtained across timesteps

Installation

Installation is currently supported on Windows, all this further support will be added for Linux environments. PyCCX can be installed along with dependencies for GMSH automatically using.

pip install pyccx

Depending on your environment, you will need to install the latest version of Calculix. This can be done through the conda-forge calculix package in the Anaconda distribution,

conda install -c conda-forge calculix

or alternatively downloading the package directly. On Windows platforms the path of the executable needs to be initialised before use.

from pyccx.core import Simulation

# Set the path for Calculix in Windows
Simulation.setCalculixPath('Path')

Usage

The following code excerpt shows an example for creating and running a steady state thermal analysis of model using PyCCX of an existing mesh generated using the pyccx.mesh.mesher class.

from pyccx.core import DOF, ElementSet, NodeSet, SurfaceSet, Simulation
from pyccx.results import ElementResult, NodalResult, ResultProcessor
from pyccx.loadcase import  LoadCase, LoadCaseType
from pyccx.material import ElastoPlasticMaterial

# Set the path for Calculix in Windows
Simulation.setCalculixPath('Path')

# Create a thermal load case and set the timesettings
thermalLoadCase = LoadCase('Thermal Load Case')

# Set the loadcase type to thermal - eventually this will be individual analysis classes with defaults
thermalLoadCase.setLoadCaseType(LoadCaseType.THERMAL)

# Set the thermal analysis to be a steady state simulation
thermalLoadCase.isSteadyState = True

# Attach the nodal and element result options to each loadcase
# Set the nodal and element variables to record in the results (.frd) file
nodeThermalPostResult = NodalResult('VolumeNodeSet')
nodeThermalPostResult.useNodalTemperatures = True

elThermalPostResult = ElementResult('Volume1')
elThermalPostResult.useHeatFlux = True

# Add the result configurations to the loadcase
thermalLoadCase.resultSet = [nodeThermalPostResult, elThermalPostResult]

# Set thermal boundary conditions for the loadcase using specific NodeSets
thermalLoadCase.boundaryConditions.append(
    {'type': 'fixed', 'nodes': 'surface6Nodes', 'dof': [DOF.T], 'value': [60]})

thermalLoadCase.boundaryConditions.append(
    {'type': 'fixed', 'nodes': 'surface1Nodes', 'dof': [DOF.T], 'value': [20]})

# Material
# Add a elastic material and assign it to the volume.
# Note ensure that the units correctly correspond with the geometry length scales
steelMat = ElastoPlasticMaterial('Steel')
steelMat.density = 1.0    # Density
steelMat.cp =  1.0        # Specific Heat
steelMat.k = 1.0          # Thermal Conductivity

analysis.materials.append(steelMat)

# Assign the material the volume (use the part name set for geometry)
analysis.materialAssignments = [('PartA', 'Steel')]

# Set the loadcases used in sequential order
analysis.loadCases = [thermalLoadCase]

# Analysis Run #
# Run the analysis
analysis.run()

# Open the results  file ('input') is currently the file that is generated by PyCCX
results = analysis.results()
results.load()

The basic usage is split between the meshing facilities provided by GMSH and analysing a problem using the Calculix Solver. Documented examples are provided in examples .

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

PyCCX-0.1.2.tar.gz (30.7 kB view details)

Uploaded Source

Built Distribution

PyCCX-0.1.2-py3-none-any.whl (34.5 kB view details)

Uploaded Python 3

File details

Details for the file PyCCX-0.1.2.tar.gz.

File metadata

  • Download URL: PyCCX-0.1.2.tar.gz
  • Upload date:
  • Size: 30.7 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.11.4

File hashes

Hashes for PyCCX-0.1.2.tar.gz
Algorithm Hash digest
SHA256 41e816f3ff8d47da8e55684abb01da41e8d397b57685beaa49ad9eb5999e4156
MD5 0b2d485960907112d25a5fe5fd27e79b
BLAKE2b-256 a691187c659379da8921c361dfee1e6f77ec86cbfd2ddb0f36d997c7eb98758f

See more details on using hashes here.

File details

Details for the file PyCCX-0.1.2-py3-none-any.whl.

File metadata

  • Download URL: PyCCX-0.1.2-py3-none-any.whl
  • Upload date:
  • Size: 34.5 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.11.4

File hashes

Hashes for PyCCX-0.1.2-py3-none-any.whl
Algorithm Hash digest
SHA256 00942b58ee0de19e2352b66a9fa925ceb4dbe91db86e18f2c52412dfbe2ec362
MD5 e0bb6479c6893f3ec23aa40f73d08fe3
BLAKE2b-256 04491a2bf6647fe4dfd86e993ae6bed3b9243639c98a3bceff26a14e451bd85f

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page