摘要：
以下面的问题为例，讨论FEM进行梁的模态分析：单跨梁受约束如图所示，梁的密度ρ=3e3kg/m³，试求梁的自振频率及振型。设梁的弹性模量为E=200GPa,横截面积A=1.0e-2㎡，长度L=2，界面惯性矩为I=1E-4m4。从题意，不考虑轴向拉伸，因此，刚度矩阵和质量矩阵为4*4的矩阵（对于6*6矩阵，见文章有限元法中，求解2D梁问题（http://www.cnblogs.com/heaventian/articles/2766227.html）下面给出其计算的程序：首先，给定单元划分，单元编号，节点的单元及整体编号的子程序coord2DBeam.mView Code % THE 阅读全文
摘要：
以下面的问题为例：对于如图所示的平面传热问题，若上端有给定的热流-2W/m2，即从下往上传输热量，结构下端有确定的温度100°，周围介质温度为20°，在两侧有换热，换热系数为α=100W/㎡/K，热导率200W/m/K,试分析其稳定温度场。使用下面的程序包进行求解：首先：将区域划分为4个单元，各单元包含的节点在element3.dat中显示：1 2 42 5 42 6 52 3 6每个节点的横纵坐标在coordinates.dat文件中显示：-10 00 010 0-5 100 105 10边界包含三个：恒温边界的节点编号，在dirichelet.dat文件中显示：1 22 阅读全文