![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
有限元程序
Smith Eric
网络工程师、有限元软件开发、计算机视觉、结构动力学、机械臂
展开
-
等参元:平面四节点四边形等参元的刚度矩阵的计算
补充: 其中:--------------------------------------------------------------------------------------例题:如图为一个平面3节点四边形等参元,采用4点高斯积分计算该单元刚度矩阵。[1有限元基础教程-曾攀...原创 2022-05-09 17:49:32 · 5669 阅读 · 10 评论 -
CST常应变三角形单元及其Matlab代码
function k=Triangle2DD3Node_Stiffness(E,NU,t,xi,yi,xj,yj,xm,ym,ID)%该函数计算单元的刚度矩阵%输入弹性模量E、泊松比NU和厚度t%输入3个节点i,j,m的坐标xi,yi,xj,yj,xm,ym%输入平面问题性质指示参数ID(1为平面应力,2为平面应变)%输出单元刚度矩阵k(6x6)%--------------------------------------------------A=(xi*(yj-ym)+x...原创 2022-05-08 19:02:26 · 1125 阅读 · 2 评论 -
余因子法(伴随法)
[1] 有限元方法基础教程(国际单位制版)(第五版)原创 2022-05-08 16:12:37 · 418 阅读 · 0 评论 -
有限元分析之边界条件
原创 2022-05-04 17:21:49 · 2378 阅读 · 1 评论 -
Matlab-有限单元法-2D梁单元的刚度矩阵组装(曾攀)
funtion z=Bar2D2Node_Assembly(KK,k,i,j)%input:单元刚度矩阵k,单元的节点编号i,j%output:总刚度矩阵 %A点,3个自由度DOF(1)=3*i-2;DOF(2)=3*i-1;DOF(3)=3*i;%B点,3个自由度DOF(4)=3*i-2;DOF(5)=3*i-1;DOF(6)=3*i;%一个梁单元具有6个自由度for n1=1:6 for n2=1:6 KK(DOF(n1),DOF(n2))=KK(DO.原创 2022-04-30 22:52:48 · 1913 阅读 · 0 评论 -
Matlab-四杆桁架结构的有限元分析(有限元分析基础-曾攀)
求解步骤如下:程序见最后Step1:结构离散和编号Step2:计算各单元的刚度矩阵 输入弹性模量、横截面积、坐标位置、夹角等参数 代入程序分别求得四个单元刚度矩阵(k1.k2,k3,k4)Step3:建立整体刚度方程%四个节点,共8个自由度>>KK=zeros(8,8);%分别对号入座,代入总刚矩阵>>KK=Bar2D2Node_Assembly(KK,k1,1,2);>>KK=Bar2D2Nod...原创 2022-04-29 23:08:49 · 2396 阅读 · 0 评论 -
Matlab-梁单元有限元分析(有限元基础-曾攀)
3.3.1梁件分析的基本力学原理M1:采用一般的建模和分析方法,即从对象中取出dxdy微元体进行分析,建立最一般的方程。这样所用的力学变量较多,方程复杂,关键是未考虑的“细长梁”这一问题的具体特征。M2: 针对细长梁用“特征建模”的简化方法来推导3大方程,其基本思想是采用工程宏观量来进行问题的描述。可以看出:该问题的的主要特征为: 1.梁为细长梁(long beam),因此可主要采用x坐标来刻画; 2.主要变形为垂直于x对的挠度,可只用挠度(deflection...原创 2022-04-29 21:00:53 · 4484 阅读 · 0 评论 -
Matlab-杆单元整体刚度矩阵组装(有限元基础-曾攀)
一维杆单元的组装:funtion z=Bar1D2Node_Assembly(KK,k,i,j)%input:单元刚度矩阵k,单元的节点编号i,j%output:总刚度矩阵DOF(1)=i;DOF(2)=j;for n1=1:2 for n2=1:2 KK(DOF(n1),DOF(n2))=KK(DOF(n1),(n2))+k(n1,n2); endendz=KK;二维杆单元组装funtion z=Bar2D2Node_Assembly(K原创 2022-04-29 17:31:17 · 8948 阅读 · 9 评论