A finite element package for learning
Project description
pyfem
pyfem是一个完全基于python语言实现的极简有限元求解器。基于numpy、scipy和meshio等第三方库,主要用于有限元方法的学习、有限元算法验证和快速建立材料本构模型的程序原型。
Installation 安装
Use the package manager pip to install pyfem. 使用pip命令安装。
pip install pyfem
Tutorial 指南
Run in command line 在命令行运行:
pyfem --help
Run the first example 执行第一个算例:
pyfem -i Job-1.toml
算例配置文件 Job-1.toml:
title = "Job-1"
[mesh]
type = "abaqus"
file = "hex8.inp"
[dof]
names = ["u1", "u2", "u3"]
order = 1
family = "LAGRANGE"
[[bcs]]
name = "BC-1"
category = "DirichletBC"
type = ""
dof = ["u1"]
node_sets = ['Set-X0']
element_sets = []
value = 0.0
[[bcs]]
name = "BC-2"
category = "DirichletBC"
type = ""
dof = ["u2"]
node_sets = ['Set-Y0']
element_sets = []
value = 0.0
[[bcs]]
name = "BC-3"
category = "DirichletBC"
type = ""
dof = ["u3"]
node_sets = ['Set-Z0']
element_sets = []
value = 0.0
[[bcs]]
name = "BC-4"
category = "NeumannBC"
type = "Distributed"
dof = ["u1"]
node_sets = ['Set-X1']
element_sets = ['Set-X1']
value = 1.0
[solver]
type = "NonlinearSolver"
option = "NewtonRaphson"
total_time = 1.0
start_time = 0.0
max_increment = 1000000
initial_dtime = 0.1
max_dtime = 1.0
min_dtime = 0.001
[[materials]]
name = "Material-1"
category = "Plastic"
type = "KinematicHardening"
data = [100000.0, 0.25, 400.0, 1000.0]
[[amplitudes]]
name = "Amp-1"
type = "TabularAmplitude"
data = [
[0.0, 0.0],
[1.0, 1.0]
]
[[sections]]
name = "Section-1"
category = "Solid"
type = "Volume"
option = "SmallStrain"
element_sets = ["Set-All"]
material_names = ["Material-1"]
data = []
[[outputs]]
type = "vtk"
field_outputs = ['S11', 'S22', 'S33', 'S12', 'S13', 'S23', 'E11', 'E22', 'E33', 'E12', 'E13', 'E23']
on_screen = false
采用abaqus格式的网格文件 hex8.inp:
*Heading
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
1, 1., 1., 1.
2, 1., 0., 1.
3, 1., 1., 0.
4, 1., 0., 0.
5, 0., 1., 1.
6, 0., 0., 1.
7, 0., 1., 0.
8, 0., 0., 0.
*Element, type=C3D8R
1, 5, 6, 8, 7, 1, 2, 4, 3
*Nset, nset=Set-X0, generate
5, 8, 1
*Elset, elset=Set-X0
1,
*Nset, nset=Set-X1, generate
1, 4, 1
*Elset, elset=Set-X1
1,
*Nset, nset=Set-Y0, generate
2, 8, 2
*Elset, elset=Set-Y0
1,
*Nset, nset=Set-Y1, generate
1, 7, 2
*Elset, elset=Set-Y1
1,
*Nset, nset=Set-Z0
3, 4, 7, 8
*Elset, elset=Set-Z0
1,
*Nset, nset=Set-Z1
1, 2, 5, 6
*Elset, elset=Set-Z1
1,
*Nset, nset=Set-All, generate
1, 8, 1
*Elset, elset=Set-All
1,
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*End Instance
**
*End Assembly
Documents 帮助文档
https://pyfem-doc.readthedocs.io/
Development
ToDo list
- 增加Neumann边界条件,支持concentrated force,distributed和pressure定义方式
- 增加hdf5计算结果输出格式
- 完善帮助文档
- 完善输入文件的校检
- 增加测试模块
- 增加日志模块
- 增加后台运行模式
- 处理平面应力状态的面外应力平衡
- 增加粘弹性力学本构模型
- 增加晶体塑性力学本构模型
- 增加温度场求解单元
- 增加温度场-位移场耦合求解单元
- 增加相场-位移场耦合求解单元
- 增加内聚区单元
- 增加动力学求解器
- 建立前处理界面
Bug list
- 采用abaqus网格文件时,如果存在node不属于任何element则在计算时会导致全局刚度矩阵奇异。
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
pyfem-0.1.3.tar.gz
(55.4 kB
view hashes)
Built Distribution
pyfem-0.1.3-py3-none-any.whl
(90.6 kB
view hashes)