import numpy as np import matplotlib.pyplot as plt import matplotlib.tri as tri from collections import OrderedDict # 一、定义函数 # 函数1：计算形函数矩阵N def Q4_N(kesi, yita): N1 = (1-kesi)*(1-yita)/4 N...

Q4单元介绍平面应力问题及平面应变问题的单元列式基本一直，差别仅在弹性矩阵 $[D]$ 。基本列式基本方程几何方程对于平面应力问题，正应力 $\sigma_z=0$，且剪应力 $\tau_{zx}=\tau_{xz}=0$，$\tau_{zy}=\tau_{yz}=0$，仅考虑单元平面内的三个平面应力分量 $\sigma_x$、$\sigma_y$ 和 $\tau_{xy}=\tau_{yx...

#桁架单元:桁架结构作用集中力荷载 import numpy as np import matplotlib.pyplot as plt #一、输入信息 #模尺材料参数 E=200000 #弹性模量 A=4532 #截面面积 EA=E*A #抗拉刚度 #节点坐标及单元节点序号 NodeCoord=np.array([[-4500,0],[-1500,0],...

问题描述如果所示为 $xy$ 平面内的桁架结构，结构几何信息如图所示。节点 1 处为固定铰支座，节点 4 处为滑动铰支座，节点 5、6、7 处分别受到 $-y$ 方向 $P=100000N$ 的集中力作用；结构中各杆件采用相同的规格，其中弹性模量 $E=200000MPa$，截面面积 $A=4532mm^2$。编程思路单元坐标（7个节点）：序号1234567节点坐标(-4500,0)(-15...

桁架单元介绍桁架单元又称为杆单元，该单元具有以下特性：（1）单元仅承受轴向载荷，不承受弯矩；（2）单元节点仅有平动自由度，没有转动自由度，单元本身只发生轴向变形。桁架单元又可分为平面桁架和空间桁架，平面桁架单元每个节点又 2 个自由度，空间桁架单元每个节点有 3 个自由度。基本方程桁架单元仅在单元端点收到轴向力作用，对于单元中任意长度为 $dx$ 的隔离体，满足以下方程：几何方程：$$ \v...