python有限元传热求解_Python进行有限元编程-平面应力问题（三节点三角形单元）...

参考书籍是：《有限元方法基础教程》(国际单位制)(第五版)

章节为：第6章 建立平面应力和平面应变刚度方程

重要的事情继续强调(有限元的基本计算流程)：

Step 1: 选择单元类型。

Step 2: 选择位移函数。

Step 3: 定义应变/位移和 应力/应变关系。

Step 4: 推导单元刚度矩阵和方程。

Step 5: 组合单元刚度方程得出总体方程并引进边界条件。

Step 6: 求解未知自由度。

Step 7:求解单元应变和应力。

Step 8: 解释结果。

==> 关于三节点三角形单元的刚度矩阵的推导过程，详见《有限元方法基础教程》第6章 建立平面应力和平面应变刚度方程

==> 本次主要计算的问题见上图。

==> 首先创建一个单元类。

==> 关于python中的构造函数(__init__方法)(因为之前学习java程序，java程序有一个构造函数的说法)传入的参数包括：三节点的坐标，和杨氏模量，泊松比，以及单元的厚度。(单元的面积是根据节点的坐标计算出来的)

import numpy as np

class T3(object):

# 接下来首先创建类的构造函数 __init__

# 对传入的参数进行说明，需要传入单元节点的坐标数据 coordinates

# 传入材料的杨氏模量E； 材料的泊松比 v

# 单元的厚度 t.

def __init__(self, coors, E, v,