![](https://img-blog.csdnimg.cn/20201014180756918.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
有限元分析
力学猿
在路上的小白一枚
展开
-
LU分解解线性方程组python实现
本人是湖南大学的一名力学学生,2020年上数值分析课时写下的作业,希望对大家有所帮助,一起进步。1、下面是自己整理的LU解线性方程组的基本理论。2、接下来是python实现代码import numpy as npdef LU(A): ''' 生成值全位0的U矩阵,和单位矩阵L ''' L = np.eye(len(A)) U = np.zeros(np.shape(A)) for r in range(1, len(A)): # 求U的第一行和原创 2020-06-09 18:28:32 · 6190 阅读 · 0 评论 -
使用Python简单实现jacobi迭代
def Jacobi(mx,mr,n=100,c=0.0001): if len(mx) == len(mr): # mx为系数矩阵,mr x = [] for i in range(len(mr)): x.append([0]) # x矩阵初始化 count = 0 # 迭代次数的计算 while count < n: nx = [] #保存单次迭代后的值的集合原创 2020-07-12 11:48:33 · 1313 阅读 · 0 评论 -
简单有限元分析python实现——二维1单元4节点刚度矩阵求解然后得到4个节点的位移和力
2020数值分析作业,已成功实现功能,可直接复制代码就能实现1、先给出大概的理论推导过程这里取高斯点坐标为Ɛ=-Ƞ=1/,权重为W1=W2=12、边界条件和解题思路如下3、源代码双手奉上from math import *import numpy as np# from my_inv import invmat # 自己在网上看的矩阵求逆方法# from LU import my_LUsolve # LU解线性方程# 1个单元四个节点def shapefunction原创 2020-06-09 21:41:32 · 6574 阅读 · 6 评论 -
用python进行简单有限元分析——二维4单元9节点总刚度组装以及9个节点的力和位移求解
2020年数值分析作业,已成功实现,可直接复制代码运行!!1、大概的理论讲解。这里主要讲总刚度矩阵的组装原理,如何得到单元刚度矩阵请看原创 2020-06-11 00:29:30 · 6062 阅读 · 0 评论