Composipy

What it is

Composipy is a Python-based library designed to address the challenges of composite plate analysis and optimization in the aerospace industry. The library offers tools for laminate stress-strain, plate buckling analysis and lamination parameter optimization, helping engineers during the design process. Utilizing object-oriented programming and native Python structures, Composipy ensures a intuitive workflow and easy integration into existing engineering practices, such as defining material, defining properties, so on.

The library is built using the most powerful Python numerical libraries, so users can think of using Composipy in the same way as Pandas, NumPy, and SciPy. In fact, Composipy is built using their objects and structures.

<img src="https://github.com/rafaelpsilva07/composipy/blob/main/doc/images/composipy_features.PNG" width="700">

It is especially useful for leveraging the data-driven culture in companies and can be used to build response surfaces, generate samples, and much more. Refer to the Composipy Examples

Main features

• ABD Matrix
• Laminate Stress and Strain
• Lamination Parameters
• Plate Buckling
• Plate Optimization

How to install

Directly From PYPI

pip install composipy

Built from source

python setup.py install

Quick start

Create the Material Properties.

from composipy import OrthotropicMaterial
 
E1 = 60800
E2 = 58250
v12 = 0.07
G12 = 4550
t = 0.21

mat_1 = OrthotropicMaterial(E1, E2, v12, G12, t)

See OrthotropicMaterial for reference.

Define the Laminate.

from composipy import LaminateProperty
          #[Ply1 --------------------------- PlyN]
          #[BOTTOM  ------------------------ TOP]
stacking = [-45, 45, 90, 0, 0, 0, 0, 90, 45, -45]
laminate1 = LaminateProperty(stacking, mat_1)

See LaminateProperty for reference.

Calculate Stiffnnes Matrix and Lamination Parameters

print(laminate1.ABD) # prints the ABD matrix as a np.ndarray
print(laminate1.xiA) # prints lamination parameters of extension as a np.ndarray
print(laminate1.xiD) # prints lamination parameters of bending as a np.ndarray

Create a Laminate Strength Analysis.

from composipy import LaminateStrength

laminate_strength = LaminateStrength(laminate1, Nxx=100, Mxx=10)
laminate_strength.epsilon0() #strains at the midplane
laminate_strength.calculate_strain() #strain ply by ply
laminate_strength.calculate_strain() #stress ply by ply

See LaminateStrength for reference.

Also, check the Stress Strain Calculation of a Laminate to see a complete example.

Create a Plate Structure for Buckling Analysis

from composipy import PlateStructure
 
constraints = {    
     'x0' : 'pinned',
     'xa' : 'pinned',
     'y0' : 'free',
     'yb' : 'pinned'
}
# Edges among 'x0', 'xa', 'y0', 'yb' can be 'clamped', 'pinned' or 'free'.
panel = PlateStructure(laminate1, 360, 360, m=10, n=10, Nxx=-1, constraints=constraints)

See PlateStructure for reference.

Calculate Buckling

print(panel.buckling_analysis()) # solve the eigenvalue problem.

Plot Buckling shape mode

print(panel.plot_eigenvalue())

Composipy is able to perform buckling calculation considering different boundary conditions and in-plane load applications. Check the Critical Buckling Examples With Varying Boundary Conditions to see a complete example.

If you are interesting in plate optimization, you may want to check the Optimization of a Plate Subjected to Buckling Loads

Theoretical References

The implementation of composipy is based on the following reference:

• SILVA, Rafael Pereira da. Composite Plate optimization combining semi-analytical model, Lamination Parameters and a Gradient-Based Optimizer. 2023. 82f. Dissertation of Master of Science – Instituto Tecnológico de Aeronáutica, São José dos Campos.
• rafaelpsilva07. (2024). rafaelpsilva07/rafaelmscdissertation: v1.0.0 (1.0.0). Zenodo. https://doi.org/10.5281/zenodo.10546621
• Application repository: https://github.com/rafaelpsilva07/rafaelmscdissertation