网站设计学的科目,无锡锡山区建设局网站,科技进步法,品牌型网站制目录 1、feon框架结构
2. 支持的单元类型
3、实例 1、feon框架结构
包含三个包#xff1a;
sa#xff1a;结构分析包
ffa#xff1a;流体分析包
derivation#xff1a;刚度矩阵包
2. 支持的单元类型
Spring1D11 - 一维弹簧单元
Spring2D11 - 二维弹簧单元
Spring…目录 1、feon框架结构
2. 支持的单元类型
3、实例 1、feon框架结构
包含三个包
sa结构分析包
ffa流体分析包
derivation刚度矩阵包
2. 支持的单元类型
Spring1D11 - 一维弹簧单元
Spring2D11 - 二维弹簧单元
Spring3D11 - 三维弹簧单元
Link1D11 - 一维杆单元
Link2D11 - 二维弹簧单元
Link3D11 - 三维弹簧单元
Beam1D11- 一维梁单元
Beam2D11- 二维弹簧单元
Beam3D11- 三维弹簧单元
Tri2d11S---- 平面应力三角形单元
Tri2D11 ---- 平面应变三角形单元
Tetra3D11 ---- 四面体单元
Quad2D11S ---- 平面应力四边形单元
Quad2D11 ---- 平面应变四边形单元
Plate3D11 ---Midline 板单元
Brick3D11 --- 六面体单元 单元命名方式为 Name dimension order type, type 1 means elastic
import feon.sa as fsa
#创建二维节点12为坐标
nfsa.Node(1,2)
#访问x坐标
print(n.x)
#访问y坐标
print(n.y)
#访问坐标
print(n.coord)
#访问节点位移keys
print(n.nAk)
#访问节点位移
print(n.disp)
#访问节点力
print(n.force)
#访问节力移keys
print(n.nBk)
#设置节点力
n.set_force(Fx100,Fy100)
#设置节点位移
n.set_disp(Ux0.1)
1.0
2.0
(1.0, 2.0)
(Ux, Uy, Phz)
{Ux: 0.0, Uy: 0.0, Phz: 0.0}
{Fx: 0.0, Fy: 0.0, Mz: 0.0}
(Fx, Fy, Mz)
3、实例 # 从feon.sa结构分析子包中导入所有的对象
from feon.sa import *
if __name__ __main__:# 定义材料参数E 210e6 # 弹性模量A 0.005 # 截面面积K 79e3 # 弹簧刚度
# 通过Node对象创建节点需要输入节点坐标n0 Node(0, 0)n1 Node(5, 0)n2 Node(10, 0)n3 Node(15, 0)n4 Node(5, 7)n5 Node(10, 7)n6 Node(15, -1)
# 通过节点创建单元Link2D11为二维杆单元输入参数为单元节点弹性模量和截面面积e0 Link2D11((n0, n1), E, A)e1 Link2D11((n1, n2), E, A)e2 Link2D11((n2, n3), E, A)e3 Link2D11((n4, n0), E, A)e4 Link2D11((n4, n1), E, A)e5 Link2D11((n4, n2), E, A)e6 Link2D11((n4, n5), E, A)e7 Link2D11((n5, n2), E, A)e8 Link2D11((n5, n3), E, A)
# 创建弹簧单元e9 Spring2D11((n3, n6), K)
# 创建一个有限元系统s System()
# 将节点和单元添加到系统中s.add_nodes(n0, n1, n2, n3, n4, n5, n6)s.add_elements(e0, e1, e2, e3, e4, e5, e6, e7, e8, e9)
# 添加力边界条件s.add_node_force(4, Fx30)
# 添加固定支座输入节点的编号0和6s.add_fixed_sup(0, 6)
# 求解系统s.solve()
#查看节点的水平和竖向位移。
print([nd.disp[Ux] for nd in s.get_nodes()])
print([nd.disp[Uy] for nd in s.get_nodes()])
#查看1号节点的位移
print(n1.disp)
#查看单元轴力
print([el.force[N] for el in s.get_elements()])
#查看1号单元的轴力
print(e1.force)
[0.0, 9.523809523809543e-05, 0.00019047619047619086, 0.00023809523809523858, 0.0005274897818499404, 0.0004798707342308927, 0.0]
[0.0, -0.00020356152587617366, -0.0002710686299836396, -0.00017721518987341803, -0.00020356152587617366, -0.00017773529665030604, 0.0]
{Ux: 9.523809523809543e-05, Uy: -0.00020356152587617366, Phz: 0.0}
[array([[-20.], [ 20.]]), array([[-20.], [ 20.]]), array([[-10.], [ 10.]]), array([[-17.20465053], [ 17.20465053]]), array([[0.], [0.]]), array([[ 17.20465053], [-17.20465053]]), array([[ 10.], [-10.]]), array([[-14.], [ 14.]]), array([[ 17.20465053], [-17.20465053]]), array([[ 14.], [-14.]])]
{N: array([[-20.], [ 20.]])} # 从feon.sa结构分析子包中导入所有的对象
from feon.sa import *
if __name__ __main__:# 定义材料参数E 210e6 # 弹性模量A 0.005 # 截面面积K 79e3 # 弹簧刚度
# 通过Node对象创建节点需要输入节点坐标n0 Node(0, 0)n1 Node(0.4, 0)n2 Node(0.4, 0.3)n3 Node(0, 0.3)
# 通过节点创建单元Link2D11为二维杆单元输入参数为单元节点弹性模量和截面面积e0 Link2D11((n0, n1), E,A)e1 Link2D11((n1, n2), E,A)e2 Link2D11((n2, n3), E,A)e3 Link2D11((n1, n3), E,A)e4 Link2D11((n0, n2), E,A)
# 创建一个有限元系统s System()
# 将节点和单元添加到系统中s.add_nodes(n0, n1, n2, n3)s.add_elements(e0, e1, e2, e3, e4)
# 添加力边界条件s.add_node_force(1, Fx20)s.add_node_force(2, Fy-25)
# 添加固定支座输入节点的编号0和6s.add_fixed_sup(0, 3)
# 求解系统s.solve()
#查看节点的水平和竖向位移。
print([nd.disp[Ux] for nd in s.get_nodes()])
print([nd.disp[Uy] for nd in s.get_nodes()])
#查看1号节点的位移
print(n1.disp)
#查看单元轴力
print([el.force[N] for el in s.get_elements()])
#查看1号单元的轴力
print(e1.force)
[0.0, 4.891240446796015e-07, 5.5684891240446785e-06, 0.0]
[0.0, -1.791534391534391e-05, -2.1925925925925923e-05, 0.0]
{Ux: 4.891240446796015e-07, Uy: -1.791534391534391e-05, Phz: 0.0}
[array([[-1.28395062], [ 1.28395062]]), array([[ 14.03703704], [-14.03703704]]), array([[-14.61728395], [ 14.61728395]]), array([[-23.39506173], [ 23.39506173]]), array([[ 18.27160494], [-18.27160494]])]
{N: array([[ 14.03703704], [-14.03703704]])}
