Installation 安装
支持的操作系统包括：Windows，Linux和MacOS。
Recommend 推荐
Use the package manager pip to install pyfem:
使用pip命令安装：
pip install -U pyfem
If you have no root access on Linux/MacOS, please try:
如果您在Linux/MacOS上没有root访问权限，请尝试以下操作：
python -m pip install -U pyfem
Users in China can install pyfem from mirrors such as:
中国用户可以使用以下镜像：
From Source 基于源代码
git clone https://github.com/sunwhale/pyfem.git
cd pyfem
pip install .
or 或者
git clone https://github.com/sunwhale/pyfem.git
cd pyfem
python install.py
采用第二种方法需要将可执行文件或批处理文件写入环境变量。
Quickstart 快速开始
Run in command line 在命令行运行
pyfem --help
Run the first example 执行第一个算例
当前算例文件存储目录 examples/tutorial ，该算例定义了一个二维平面应变模型，材料为塑性随动强化，载荷为y方向的循环拉伸-压缩。
cd examples/tutorial
pyfem -i Job-1.toml
算例配置文件 Job-1.toml
title = "Job-1"
[mesh] # 前处理网格文件
type = "gmsh"
file = 'mesh.msh'
[dof] # 自由度
names = ["u1", "u2"]
order = 1
family = "LAGRANGE"
[[amplitudes]] # 幅值列表
name = "Amp-1"
type = "TabularAmplitude"
start = 0.0
data = [
[0.0, 0.0],
[0.5, 1.0],
[1.0, 0.0],
[1.5, -1.0],
[2.0, 0.0],
[2.5, 1.0],
[3.0, 0.0],
[3.5, -1.0],
[4.0, 0.0],
[4.5, 1.0],
[5.0, 0.0],
]
[[bcs]] # 边界条件列表
name = "BC-1"
category = "DirichletBC"
type = ""
dof = ["u2"]
node_sets = ['bottom']
element_sets = []
value = 0.0
[[bcs]] # 边界条件列表
name = "BC-2"
category = "DirichletBC"
type = ""
dof = ["u1"]
node_sets = ['left']
element_sets = []
value = 0.0
[[bcs]] # 边界条件列表
name = "BC-3"
category = "DirichletBC"
type = ""
dof = ["u2"]
node_sets = ['top']
element_sets = []
value = 0.01
amplitude_name = "Amp-1"
[solver] # 求解器属性
type = "NonlinearSolver"
option = "NewtonRaphson"
total_time = 5.0
start_time = 0.0
max_increment = 1000000
initial_dtime = 0.05
max_dtime = 0.05
min_dtime = 0.001
[[materials]] # 材料列表
name = "Material-1"
category = "Plastic"
type = "KinematicHardening"
data = [100000.0, 0.25, 400.0, 1000.0]
[[sections]] # 截面列表
name = "Section-1"
category = "Solid"
type = "PlaneStrain"
option = "SmallStrain"
element_sets = ["rectangle"]
material_names = ["Material-1"]
data = []
[[outputs]] # 输出列表
type = "vtk"
field_outputs = ['S11', 'S22', 'S12', 'E11', 'E22', 'E12']
is_save = true
gmsh格式的网格文件 mesh.msh
$MeshFormat
4.1 0 8
$EndMeshFormat
$PhysicalNames
5
1 5 "left"
1 6 "right"
1 7 "top"
1 8 "bottom"
2 9 "rectangle"
$EndPhysicalNames
$Entities
4 4 1 0
1 0 0 0 0
2 1 0 0 0
3 1 1 0 0
4 0 1 0 0
1 0 0 0 1 0 0 1 8 2 1 -2
2 1 0 0 1 1 0 1 6 2 2 -3
3 0 1 0 1 1 0 1 7 2 3 -4
4 0 0 0 0 1 0 1 5 2 4 -1
1 0 0 0 1 1 0 1 9 4 3 4 1 2
$EndEntities
$Nodes
9 9 1 9
0 1 0 1
1
0 0 0
0 2 0 1
2
1 -0 0
0 3 0 1
3
1 1 0
0 4 0 1
4
0 1 0
1 1 0 1
5
0.4999999999986921 0 0
1 2 0 1
6
1 0.4999999999986921 0
1 3 0 1
7
0.5000000000020595 1 0
1 4 0 1
8
0 0.5000000000020595 0
2 1 0 1
9
0.5000000000003758 0.5000000000003758 0
$EndNodes
$Elements
5 12 1 12
1 1 1 2
1 1 5
2 5 2
1 2 1 2
3 2 6
4 6 3
1 3 1 2
5 3 7
6 7 4
1 4 1 2
7 4 8
8 8 1
2 1 3 4
9 3 7 9 6
10 6 9 5 2
11 7 4 8 9
12 9 8 1 5
$EndElements
Postprocess 后处理
算例计算完成后将在配置文件所在目录下生成 .pvd 或 .vtu 文件，可以使用开源可视化软件 paraview 进行查看。
Preprocess 前处理
本项目暂不提供前处理模块，基于 meshio 库，可以识别gmsh、abaqus 和 ansys等有限元软件的网格文件。