悬臂梁MATLAB有限元算例注释

用有限元法对悬臂梁分析的算例

算例：

如下图所示的悬臂梁，受均布载荷q＝1N／mm2作用。E＝2．1×105N／mm2, μ＝0.3厚度h＝10mm。现用有限元法分析其位移及应力。

梁可视为平面应力状态，先按图示尺寸划分为均匀的三角形网格，共有8×10＝80个单元，5×ll＝55个节点，坐标轴以及单元与节点的编号如图。将均布载荷分配到各相应节点上，把有约束的节点5l、52、53、54、55视作固定铰链，建立如图所示的离散化计算模型。

程序计算框图：

（续左）

程序中的函数功能介绍及源代码

1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym)――该函数用于计算平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)、第二个节点坐标为(xj,yj)、第三个节点坐标为(xm，ym)时的线性三角形元的单元刚度矩阵.该函数返回6×6的单位刚度矩阵k.

2. LinearTriangleAssemble(K,k,i,j,m)――该函数将连接节点i,j,m的线性三角形元的单元刚度矩阵k集成到整体刚度矩阵K。每集成一个单元，该函数都将返回2N×2N的整体刚度矩阵K.

3. LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u)-- 该函数计算在平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)第二个节点坐标为(xj,yj)、第三个节点坐标为(xm，ym)以及单元位移矢量为u时的单元应力。该函数返回单元应力矢量。

函数源代码：

function y = LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym)

A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2;%三角形单元面积，单元节点应该按逆时针排序，保证每个三角形单元的面积都为正值（也可作为一个小函数：LinearTriangleElementArea）

betai = yj-ym;betaj = ym-yi;betam = yi-yj;gammai = xm-xj;

gammaj = xi-xm;gammam = xj-xi;

B = [betai 0 betaj 0 betam 0 ;

0 gammai 0 gammaj 0 gammam ;

gammai betai gammaj betaj gammam betam]/(2*A);

%B为应变矩阵，其中betai=yi-ym,betaj=ym-yi,betam=yi-yj.

gammai=xm-xj, gammaj=xi-xm, gammam=xj-xi.

D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2];

%D为弹性矩阵，分为平面应力问题和平面应变问题

对于平面应力问题D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2];

对于平面应变问题E1=E/（1-NU*NU），NU1=NU/（1-NU）

y = t*A*B'*D*B;%单元刚度矩阵

function y = LinearTriangleAssemble(K,k,i,j,m)

K(2*i-1,2*i-1) = K(2*i-1,2*i-1) + k(1,1); K(2*i-1,2*i) = K(2*i-1,2*i) + k(1,2);

K(2*i-1,2*j-1) = K(2*i-1,2*j-1) + k(1,3); K(2*i-1,2*j) = K(2*i-1,2*j) + k(1,4);

K(2*i-1,2*m-1) = K(2*i-1,2*m-1) + k(1,5); K(2*i-1,2*m) = K(2*i-1,2*m) + k(1,6);

K(2*i,2*i-1) = K(2*i,2*i-1) + k(2,1); K(2*i,2*i) = K(2*i,2*i) + k(2,2);

K(2*i,2*j-1) = K(2*i,2*j-1) + k(2,3); K(2*i,2*j) = K(2*i,2*j) + k(2,4);

K(2*i,2*m-1) = K(2*i,2*m-1) + k(2,5); K(2*i,2*m) = K(2*i,2*m) + k(2,6);

K(2*j-1,2*i-1) = K(2*j-1,2*i-1) + k(3,1); K(2*j-1,2*i) = K(2*j-1,2*i) + k(3,2);

K(2*j-1,2*j-1) = K(2*j-1,2*j-1) + k(3,3); K(2*j-1,2*j) = K(2*j-1,2*j) + k(3,4);

K(2*j-1,2*m-1) = K(2*j-1,2*m-1) + k(3,5); K(2*j-1,2*m) = K(2*j-1,2*m) + k(3,6);

K(2*j,2*i-1) = K(2*j,2*i-1) + k(4,1); K(2*j,2*i) = K(2*j,2*i) + k(4,2);

K(2*j,2*j-1) = K(2*j,2*j-1) + k(4,3); K(2*j,2*j) = K(2*j,2*j) + k(4,4);

K(2*j,2*m-1) = K(2*j,2*m-1) + k(4,5); K(2*j,2*m) = K(2*j,2*m) + k(4,6);

K(2*m-1,2*i-1) = K(2*m-1,2*i-1) + k(5,1); K(2*m-1,2*i) = K(2*m-1,2*i) + k(5,2);

K(2*m-1,2*j-1) = K(2*m-1,2*j-1) + k(5,3); K(2*m-1,2*j) = K(2*m-1,2*j) + k(5,4);

K(2*m-1,2*m-1) = K(2*m-1,2*m-1) + k(5,5); K(2*m-1,2*m) = K(2*m-1,2*m) + k(5,6);

K(2*m,2*i-1) = K(2*m,2*i-1) + k(6,1); K(2*m,2*i) = K(2*m,2*i) + k(6,2);

K(2*m,2*j-1) = K(2*m,2*j-1) + k(6,3); K(2*m,2*j) = K(2*m,2*j) + k(6,4);

K(2*m,2*m-1) = K(2*m,2*m-1) + k(6,5); K(2*m,2*m) = K(2*m,2*m) + k(6,6);

y = K;%对号入座，如前所述，每集成一次都将返回2N×2N的整体刚度矩阵K.此题为110×110

function y = LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u)

A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2;

betai = yj-ym;betaj = ym-yi;betam = yi-yj;gammai = xm-xj;gammaj = xi-xm;gammam = xj-xi;

B = [betai 0 betaj 0 betam 0 ;

0 gammai 0 gammaj 0 gammam ;

gammai betai gammaj betaj gammam betam]/(2*A);

D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2];%平面应力和平面应变问题两种情况

y = D*B*u;%单元应力计算

主程序源代码

E=21e7;NU=0.3;t=0.01;

stifflike5=LinearTriangleElementStiffness(E,NU,t,0.4,0.08,0.36,0.08,0.36,0.06,1)%选取2个基本单元，调用M文件

stifflike1=LinearTriangleElementStiffness(E,NU,t,0.4,0.08,0.36,0.06,0.4,0.06,1)

K=sparse(110,110); %creat a xishu matrix for total stiff创建一个稀疏矩阵

for i=1:49

if rem(i,5)%模取余， bool型变量，非零即为真

j=i;

K=LinearTriangleAssemble(K,stifflike5,j,j+5,j+6);%节点编号

K=LinearTriangleAssemble(K,stifflike1,j,j+6,j+1);

end

end%将每个单元刚度矩阵集成到总刚中

K=full(K);%转化稀疏矩阵 k=K(1:100,1:100);

k=[K,zeros(100,10);zeros(10,100),eye(10)];

k=sparse(k);%利用边界条件简化基本方程

Q=sparse(2:10:92,1,[-200,-400,-400,-400,-400,-400,-400,-400,-400,-400,],110,1);%外部荷载，此处不包括约束条件，通过形函数确定，是不是可以理解为梁的两端为中间的一半呢？

d=k\Q;%高斯消元法，比克莱姆法则在计算速度上有绝对的优势！

x=0:0.04:0.4;

plot(x,d(106:-10:6))%基本绘图命令

grid%带网格

y=zeros(80,3);

q=0;

for i=1:49

switch rem(i,5)

case 1

j=2*i;

u=[d(j-1) d(j) d(j+11) d(j+12) d(j+1) d(j+2)];

u=u';

xl=0.4;yl=0.08;xm=0.36;ym=0.06;xn=0.4;yn=0.06;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

case 2

j=2*i;

u=[d(j-1) d(j) d(j+11) d(j+12) d(j+1) d(j+2)];

u=u';

xl=0.4;yl=0.06;xm=0.36;ym=0.04;xn=0.4;yn=0.04;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

case 3

j=2*i;

u=[d(j-1) d(j) d(j+11) d(j+12) d(j+1) d(j+2)];

u=u';

xl=0.4;yl=0.04;xm=0.36;ym=0.02;xn=0.4;yn=0.02;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

case 4

j=2*i;

u=[d(j-1) d(j) d(j+11) d(j+12) d(j+1) d(j+2)];

u=u';

xl=0.4;yl=0.02;xm=0.36;ym=0;xn=0.4;yn=0;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

otherwise

q=q+3;

end

q=4;

for i=1:49

switch rem(i,5)

case 1

j=2*i;

u=[d(j-1) d(j) d(j+9) d(j+10) d(j+11) d(j+12)];

u=u';

xl=0.4;yl=0.08;xm=0.36;ym=0.08;xn=0.36;yn=0.06;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

case 2

j=2*i;

u=[d(j-1) d(j) d(j+9) d(j+10) d(j+11) d(j+12)];

u=u';

xl=0.4;yl=0.06;xm=0.36;ym=0.06;xn=0.36;yn=0.04;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

case 3

j=2*i;

u=[d(j-1) d(j) d(j+9) d(j+10) d(j+11) d(j+12)];

u=u';

xl=0.4;yl=0.04;xm=0.36;ym=0.04;xn=0.36;yn=0.02;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

case 4

j=2*i;

u=[d(j-1) d(j) d(j+9) d(j+10) d(j+11) d(j+12)];

u=u';

xl=0.4;yl=0.02;xm=0.36;ym=0.02;xn=0.36;yn=0;

y(i+q,:)=LinearTriangleElementStresses(E,NU,t,xl,yl,xm,ym,xn,yn,u)';

xl=xl-0.04;xm=xm-0.04;xn=xn-0.04;

otherwise

q=q+3;

end

end% y(i+q,:)这是实现什么的？没见过这种用法，算法上应该就是通过节点位移实现指定单元的内力，这部分本人看的也晕晕的，望高人指点

N=y(73:80,1)

结果图及数据输出

悬臂梁轴线挠度图：

一单元的单元刚阵

1.0e+006 *

0.8077 0 0 -0.4038 -0.8077 0.4038 0 2.3077 -0.3462 0 0.3462 -2.3077 0 -0.3462 0.5769 0 -0.5769 0.3462 -0.4038 0 0 0.2019 0.4038 -0.2019 -0.8077 0.3462 -0.5769 0.4038 1.3846 -0.7500

0.4038 -2.3077 0.3462 -0.2019 -0.7500 2.5096 五单元的单元刚阵

1.0e+006 *

0.5769 0 -0.5769 0.3462 0 -0.3462 0 0.2019 0.4038 -0.2019 -0.4038 0 -0.5769 0.4038 1.3846 -0.7500 -0.8077 0.3462 0.3462 -0.2019 -0.7500 2.5096 0.4038 -2.3077 0 -0.4038 -0.8077 0.4038 0.8077 0 00.050.10.150.20.250.30.350.4-0.35-0.3-0.25-0.2-0.15

-0.1-0.05

x/m

w /m

-0.3462 0 0.3462 -2.3077 0 2.3077

根部73-80各单元应力计算结果如下（n/m2）：

1.0e+007 *

2.1119 -0.0621 -2.2816 -4.8824 5.0479 2.4065 0.0352 -2.3753

matlab有限元分析实例

MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件，用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人，控制系统等领域。 MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室），软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。 MATLAB有限元分析与应用:

《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书，作者是卡坦，译者是韩来彬。内容简介：《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用，书中的每个示例都以交互的方式求解，使读者很容易就能把MATLAB用于有限分析和应用。另外，《MATLAB有限元分析与应用》还提供了大量免费资源。《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂，循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题，具有很高的工程应用价值。

悬臂梁分析报告

悬臂梁受力分析报告高一博 2016.11.13 西安理工大学机械与精密仪器工程学院

摘要利用ANSYS对悬臂梁进行有限元静力学分析,得到悬臂梁的最大应力和挠度位移。从而校验结构强度和尺寸定义，从而对结构进行最优化设计修正。关键词：悬臂梁，变形分析，应力分析

目录一.问题描述： (4) 二.分析的目的和内容： (4) 三.分析方案和有限元建模方法： (4) 四.几何模型 (4) 五．有限元模型 (4) 六．计算结果： (5) 七.结果合理性的讨论、分析 (8) 八．结论 (8) 参考文献 (8)

一.问题描述：现有一悬臂梁，长500MM，一端固定，另外一端施加一个竖直向下的集中力200N。其截面20MMX20MM的矩形，现在要分析该梁的在集中力作用下产生的位移，应力和局部应力。二.分析的目的和内容： 1.观察悬臂梁的变形情况； 2.观察分析悬臂梁的应力变化； 3.找出其最大变形和最大应力点，分析形成原因；三.分析方案和有限元建模方法： 1.使用ANSYS-modeling-create-volumes-block建模， 2.对梁进行材料定义，网格划分。 3.一端固定，另外一端施加一个向下的200N的力。 4.后处理中查看梁的应力和变形情况。四.几何模型 500X20X20的梁在在ANSYS中进行绘制.由于结构简单规则，无需简化。五．有限元模型单元类型：solid brick8node45 材料参数：弹性模量2e+11pa,泊松比0.3 边界条件：一端固定，一端施加载荷载荷：F=200N 划分网格后的悬臂梁模型

Matlab有限元分析操作基础

Matlab 有限元分析20140226 为了用Matlab 进行有限元分析，首先要学会Matlab 基本操作，还要学会使用Matlab 进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵 11221 21200k k k k k k k k -????-????--+??

2. Matlab有限元分析的基本操作（1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…）（3）组装系统刚度矩阵（集成整体刚度矩阵）（4）引入边界条件（消除冗余方程）（5）解方程（6）后处理（扩展计算）

3. Matlab有限元分析实战【实例1】

分析：步骤一：单元划分

步骤二：构造单元刚度矩阵 >>k1=SpringElementStiffness(100) >>…?

步骤三：构造系统刚度矩阵 a) 分析SpringAssemble库函数function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？今天的系统刚度矩阵是什么？因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ?? ?? -- ?? ？

ANSYS悬臂梁的自由端受力的有限元计算[1]

悬臂梁自由端受力的有限元计算任柳杰10110290005 一、计算目的 1、掌握ANSYS软件的基本几何形体构造、网格划分、边界条件施加等方法。 2、熟悉有限元建模、求解及结果分析步骤和方法。 3、利用ANSYS软件对梁结构进行有限元计算。 4、梁的变形、挠曲线等情况的分析。 5、一维梁单元，二维壳单元，三维实体单元对计算结果的影响。 6、载荷施加在不同的节点上对结果的影响。二、计算设备 PC，ANSYS软件（版本为11.0）三、计算内容悬臂梁受力模型如上图所示，一段长100[mm]的梁，一端固定，另一段受到平行于梁截面的集中力F的作用，F=100[N]。梁的截面为正方形，边长为10[mm]。梁所用的材料：弹性模量E=2.0 105[MPa]，泊松比0.3。四、计算步骤（以梁单元为例） 1、分析问题。分析该物理模型可知，截面边长/梁长度=0.1是一个较小的值，我们可以用梁单元来分析这样的模型。当然，建立合适的壳单元模型和实体单元模型也是可以的。故拟采用这三种不同的方式建立模型。以下主要阐述采用梁单元的模型的计算步骤。 2、建立有限元模型。 a)创建工作文件夹并添加标题；在个人的工作目录下创建一个文件夹，命名为beam，用于保存分析过程中生成的各种文件。启动ANSYS后，使用菜单“File”——“Change Directory…”将工作目录指向beam 文件夹；使用/FILNAME,BEAM命令将文件名改为BEAM，这样分析过程中生成的文件均以BEAM为前缀。偏好设定为结构分析，操作如下： GUI: Main Menu > Preferences > Structural b)选择单元；进入单元类型库，操作如下： GUI: Main Menu > Preprocessor > Element Type > Add/Edit/Delete > Add… 对话框左侧选择Beam选项，在右侧列表中选择2D elastic 3选项，然后单击OK按钮。

有限元的MATLAB解法

有限元的MATLAB解法 1.打开MATLAB。 2.输入“pdetool”再回车，会跳出PDE Toolbox的窗口（PDE意为偏微分方程，是partial differential equations的缩写），需要的话可点击Options菜单下Grid命令，打开栅格。 3.完成平面几何模型：在PDE Toolbox的窗口中，点击工具栏下的矩形几何模型进行制作模型，可画矩形R，椭圆E，圆C，然后在Set formula栏进行编辑并（如双脊波导R1+R2+R3改为RI-R2-R3,设定a、b、s/a、d/b的值从而方便下步设定坐标）用算术运算符将图形对象名称连接起来，若还需要，可进行储存，形成M文件。 4.用左键双击矩形进行坐标设置：将大的矩形left和bottom都设为0，width是矩形波导的X轴的长度，height是矩形波导的y轴的长度，以大的矩形左下角点为原点坐标为参考设置其他矩形坐标。 5.进行边界设置：点击“Boundary”中的“Boundary Mode”,再点击

“Boundary”中的“Specify Boundary Conditions”,选择符合的边界条件，Neumann为诺曼条件,Dirichlet为狄利克雷条件，边界颜色显示为红色。 6.进入PDE模式：点击"PDE"菜单下“PDE Mode”命令，进入PDE 模式，单击“PDE Specification”,设置方程类型，“Elliptic”为椭圆型，“Parabolic”为抛物型，“Hyperbolic”为双曲型，“Eigenmodes”为特征值问题。 7.对模型进行剖分：点击“Mesh”中“Initialize Mesh”进行初次剖分，若要剖的更细，再点击“Refine Mesh”进行网格加密。 8.进行计算：点击“Solve”中“Solve PDE”,解偏微分方程并显示图形解，u值即为Hz或者Ez。 9.单击“Plot”菜单下“Parameters”选项，打开“Plot Selection”对话框。选中Color,Height(3-D plot)和Show mesh三项，然后单击“Plot”按钮，显示三维图形解。 10.如果要画等值线图和矢量场图，单击“Plot”菜单下“Parameters”选项,打开“Plot Selection”对话框。选中Contour和Arrows两项，然后单击Plot按钮，可显示解的等值线图和矢量场图。 11.将计算结果条件和边界导入MATLAB中：点击“Export Solution”,再点击“Mesh”中“Export Mesh”。

有限元悬臂梁仿真

有限元方法大作业课程设计题目：若干个质量不等的仪器要安装在均匀悬臂梁（或板）不同位置上，仪器间要有预留安全距离，试确定一种安装方法，使梁（或板）的变形最小或第一阶固有频率最高。题目分析： 1 题目中没有给定梁的材料和形状、仪器的数量和质量，以及仪器的安全距离。在这里不妨假定，梁的材料为结构钢，其密度为、杨氏模量为Pa、泊松比为0.3，梁的形状为。仪器的数量为3个，均匀的安装在梁上，其质量及其组合如表1所示。表1 仪器的质量、及其组合 2 本次采用solidworks建立梁的实体模型，并导入ansys workbench软件中进行计算。梁模型左端固定，仪器安装顺序依次从左到右。 3 在题目中，需要找到一种安装组合使得梁的变形最小或第一阶固有频率最高，这分别是静力学分析问题和模态分析问题。在静力学分析中，如图2-1所示，在梁上安装仪器的位置上，加上一个加力面（半径为20mm的圆）。在加力面上可以施加均布载荷，这里将仪器的质量换算成相应的均布载荷，施加到相应的加力面上，如图2-2所示。图2-1 ansys workbench实体梁的静力学分析

图2-2 加力面和加力面上的均布载荷采用solid187单元对模型进行网格划分，solid187单元是一个高阶3维10节点固体结构单元，如图2-3所示，单元通过10个节点来定义，每个节点有3个沿着xyz方向平移的自由度。并对加力面附近进行加密，如图2-4所示。进而进行静力学分析，得到梁的总体变形量（total-Deformation）。图2-3 solid187单元图2-4 梁模型网格划分和加力面加密 4 梁的固有频率可由无阻尼自由振动方程求解：令：得到：当：从而求的梁的自振频率。在ansys workbench中，将仪器的质量用质量点代替，并安置在相应的位置上，如图2-5所示。采用solid186单元对模型进行网格划分，其结果如图2-6所示，solid186是一个高阶3维20节点固体结构单元，如图2-7所示，单元通过20个节点来定义，每个节点有3个沿着xyz方向平移的自由度。然后，求解梁模型的前6阶的固有频率。

《有限元基础教程》_【MATLAB算例】3.3.7(2)__三梁平面框架结构的有限元分析(Beam2D2Node)

【MA TLAB 算例】3.3.7(2) 三梁平面框架结构的有限元分析 (Beam2D2Node) 如图3-19所示的框架结构，其顶端受均布力作用，结构中各个截面的参数都为：113.010Pa E =?，746.510I m -=?，426.810A m -=?。试基于MA TLAB 平台求解该结构的节点位移以及支反力。图3-19 框架结构受一均布力作用解答：对该问题进行有限元分析的过程如下。（1）结构的离散化与编号将该结构离散为3个单元，节点位移及单元编号如图3-20所示，有关节点和单元的信息见表3-5。（a ）节点位移及单元编号

（b）等效在节点上的外力图3-20 单元划分、节点位移及节点上的外载（2）各个单元的描述首先在MA TLAB环境下，输入弹性模量E、横截面积A、惯性矩I和长度L，然后针对单元1，单元2和单元3，分别二次调用函数Beam2D2Node_ElementStiffness，就可以得到单元的刚度矩阵k1(6×6)和k2(6×6)，且单元2和单元3的刚度矩阵相同。 >> E=3E11; >> I=6.5E-7; >> A=6.8E-4; >> L1=1.44; >> L2=0.96; >> k1=Beam2D2Node_Stiffness(E,I,A,L1); >> k2=Beam2D2Node_Stiffness(E,I,A,L2); （3）建立整体刚度方程将单元2和单元3的刚度矩阵转换成整体坐标下的形式。由于该结构共有4个节点，则总共的自由度数为12，因此，结构总的刚度矩阵为KK(12×12)，对KK清零，然后两次调用函数Beam2D2Node_Assemble进行刚度矩阵的组装。 >> T=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1] ; >> k3=T'*k2*T; >> KK=zeros(12,12); >> KK=Beam2D2Node_Assemble(KK,k1,1,2);

基于Matlab语言的按平面三角形单元划分的结构有限元程序设计模板

基于Matlab语言的按平面三角形单元划分的结构有限元程序设计专业：建筑与土木工程班级：建工研12-2 姓名：韩志强学号： 471220580

基于Matlab语言的按平面三角形单元划分结构有限元程序设计一、有限单元发及Matlab语言概述 1. 有限单元法随着现代工业、生产技术的发展，不断要求设计高质量、高水平的大型、复杂和精密的机械及工程结构。为此目的，人们必须预先通过有效的计算手段，确切的预测即将诞生的机械和工程结构，在未来工作时所发生的应力、应变和位移因此，需要寻求一种简单而又精确的数值分析方法。有限单元法正是适应这种要求而产生和发展起来的一种十分有效的数值计算方法。有限元法把一个复杂的结构分解成相对简单的“单元”，各单元之间通过结点相互连接。单元内的物理量由单元结点上的物理量按一定的假设内插得到，这样就把一个复杂结构从无限多个自由度简化为有限个单元组成的结构。我们只要分析每个单元的力学特性，然后按照有限元法的规则把这些单元“拼装”成整体，就能够得到整体结构的力学特性。有限单元法基本步骤如下： (1)结构离散：结构离散就是建立结构的有限元模型，又称为网格划分或单元划分，即将结构离散为由有限个单元组成的有限元模型。在该步骤中，需要根据结构的几何特性、载荷情况等确定单元体内任意一点的位移插值函数。 (2)单元分析：根据弹性力学的几何方程以及物理方程确定单元的刚度矩阵。 (3)整体分析:把各个单元按原来的结构重新连接起来，并在单元刚度矩阵的基础上确定结构的总刚度矩阵，形成如下式所示的整体有限元线性方程： {}[]{}δ F=① K 式中，{}F是载荷矩阵,[]K是整体结构的刚度矩阵,{}δ是节点位移矩阵。 (4)载荷移置:根据静力等效原理，将载荷移置到相应的节点上，形成节点载荷矩阵。 (5)边界条件处理:对式①所示的有限元线性方程进行边界条件处理。 (6)求解线性方程:求解式①所示的有限元线性方程，得到节点的位移。在该步骤中，若有限元模型的节点越多，则线性方程的数量就越多，随之有限元分析的计算量也将越大。 (7)求解单元应力及应变根据求出的节点位移求解单元的应力和应变。

有限元分析及应用报告-利用ANSYS软件分析带孔悬臂梁

有限元分析及应用报告题目：利用ANSY软件分析带孔悬臂梁姓名：xxx 学号：xxx 班级：机械xxx 学院: 机械学院指导老师：xxx 二零一五年一月

问题概述图示为一隧道断面，其内受均布水压力q,外受土壤均布压力p；试采用不同单元计算断面内的位移及应力，并分别分析q=0或p=0时的位移和应力分布情况。（材料为钢，隧道几何尺寸和压力大小自行确定）本例假定内圆半径为1m，外圆半径为2m，外受均布压力p=10000pa ，内受均布压力为q=20000pa 。问题分析由题目可知，隧道的的长度尺寸远远大于截面尺寸，并且压力在长度方向上均匀分布，因此本问题可以看作为平面应变问题。由于在一个截面内，压力沿截面四周均匀分布，且截面是对称的圆环，所以可以只取截面1/4进行有限元建模分析，这样不仅简化了建模分析过程，也能保证得到精确的结果。由以上分析，可以选取单元类型plane42进行有限元分析，在option中选择K3 为plane strain。

三.有限元建模 1.设置计算类型由问题分析可知本问题属于平面静应力问题，所以选择preferences 为structure 。 2.单元类型选定选取平面四节点常应变单元plane42,来计算分析隧道截面的位移和应力。由于此问题为平面应变问题，在设置element type的K3时将其设置为plane strain。 3.材料参数隧道的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比(T =0.3 4.几何建模按照题目所给尺寸利用ansys的modeling依次建立keypoint : 1(0,0),2(1,0),3(2,0),4(0,2),5(0,1) , create LINES 依次连接keypoint 2、3和4、5即可创建两条直线，使用create article 的By cent & radius 创建两条圆弧。create AREAS依次选择四条线即建立了所需的1/4截面。 5.网格划分

(完整版)有限元大作业matlab---课程设计例子

有限元大作业程序设计学校：天津大学院系：建筑工程与力学学院专业：01级工程力学姓名：刘秀学号：\\\\\\\\\\\ 指导老师：

连续体平面问题的有限元程序分析 [题目]：如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]：由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。

[程序原理及实现]：用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。该程序的特点如下：问题类型：可用于计算弹性力学平面问题和平面应变问题单元类型：采用常应变三角形单元位移模式：用用线性位移模式载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷材料性质：弹性体由单一的均匀材料组成约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束方程求解：针对半带宽刚度方程的Gauss消元法

输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：

悬臂梁ansys有限元分析求最大挠度

（一）悬臂梁ansys 有限元分析求最大挠度问题：悬臂梁长1000mm ，宽50mm ，高10mm ，左端固定，求其在自重作用下的最大挠度？解：弯矩方程： 221) ()(x l q x M --= 微分方程： 22 1'')(x l q y EI z -= 积分求解：D Cx qx qlx x ql y EI C qx qlx x ql y EI z z +++-=++-=4322322'24 1 6125.06 1 5.05.0 由边界条件：0; 0, 0' ' ====A A A y y x θ 得：C=0， D=0 I=1/12*h^3*b,h 为梁截面的高，b 为梁截面的宽。 q=ρ*g*a*h*l 材料力学公式求：Y=EI 85 gahl^ρ=5.733mm L

ANSYS 模拟求：Y=5.5392mm，详细见下步骤 ANSYS 软件设置及其具体过程如下：步骤1：建立一个模型，在model下creat一个长1，宽0.05，高0.01的长方体实体。（单位默认为m）步骤2：材料属性设置。密度：7800，杨氏模量：2E11，泊松比0.3。

步骤3：划分网格。设置网格单元为structure solid brick 8node 185，mesh tool中设置网格大小为0.002，HEX下点击mesh。

步骤4：施加载荷；在preprocessor中inertia中设置重力加速度Y方向为9.8。在左面施加固定约束（三个方向固定）

步骤5:：求解。在solve下solve current LS。步骤6：后处理查看。在result中plot result，查看nodes displacement。List查看文本，观察nodes的最大位移点。

Matlab有限元分析操作基础共11页

Matlab有限元分析20140226 为了用Matlab进行有限元分析，首先要学会Matlab基本操作，还要学会使用Matlab进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵

3. Matlab有限元分析实战【实例1】

分析：步骤一：单元划分

>>k1=SpringElementStiffness(100)

a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？今天的系统刚度矩阵是什么？因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??？

ansys-二维悬臂梁有限元分析

2 软件知识学习 2.1 软件的使用与介绍软件介绍： ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer, NASTRAN, Alogor, I－DEAS, AutoCAD等，是现代产品设计中的高级CAE工具之一。 ANSYS有限元软件包是一个多用途的有限元法计算机设计程序，可以用来求解结构、流体、电力、电磁场及碰撞等问题。因此它可应用于以下工业领域：航空航天、汽车工业、生物医学、桥梁、建筑、电子产品、重型机械、微机电系统、运动器械等。软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型；分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。软件提供了100种以上的单元类型，用来模拟工程中的各种结构和材料。该软件有多种不同版本，可以运行在从个人机到大型机的多种计算机设备上，如PC，SGI，HP，SUN，DEC，IBM，CRAY等。

基于matlab的有限元法分析平面应力应变问题刘刚

姓名：刘刚学号：15 平面应力应变分析有限元法 Abstruct:本文通过对平面应力/应变问题的简要理论阐述，使读者对要分析的问题有大致的印象，然后结合两个实例，通过MATLAB软件的计算，将有限元分析平面应力/应变问题的过程形象的展示给读者，让人一目了然，快速了解有限元解决这类问题的方法和步骤！一.基本理论有限元法的基本思路和基本原则以结构力学中的位移法为基础，把复杂的结构或连续体看成有限个单元的组合，各单元彼此在节点出连接而组成整体。把连续体分成有限个单元和节点，称为离散化。先对单元进行特性分析，然后根据节点处的平衡和协调条件建立方程，综合后做整体分析。这样一分一合，先离散再综合的过程，就是把复杂结构或连续体的计算问题转化简单单元分析与综合问题。因此，一般的有限揭发包括三个主要步骤：离散化单元分析整体分析。二.用到的函数 1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym,p) (K k I f) (k u) (k u A) (E NU t) 三.实例例1.考虑如图所示的受均布载荷作用的薄平板结构。将平板离散化成两个线性三角元，假定E=200GPa，v=，t=0.025m,w=3000kN/m. 1.离散化 2.写出单元刚度矩阵

通过matlab 的LinearTriangleElementStiffness 函数，得到两个单元刚度矩阵1k 和2k ，每个矩阵都是6 6的。 >> E=210e6 E = >> k1=LinearTriangleElementStiffness(E,NU,t,0,0,,,0,,1) k1 = +006 * Columns 1 through 5 0 0 0 0 0 0 0 0 Column 6 >> NU= NU = >> t= t = >> k2=LinearTriangleElementStiffness(E,NU,t,0,0,,0,,,1)

悬臂梁有限元模拟分析步骤

Introduction to Simulation I-DEAS Tutorials: Simulation Projects Simulation involves three major steps: Pre-processing (modeling, applying boundary conditions, meshing); solving the model; and post-processing (displaying the results). Learn how to: ?create a finite element model ?apply boundary conditions ?mesh the FE model ?solve the FE model ?display the results

Before you begin... Prerequisite tutorials:?Introducing the I-DEAS Interface Quick Tips to Using I-DEAS –and– Creating Parts ?Extruding and Revolving Features

If you didn’t start I-DEAS with a new (empty) model file, open a new one now and give it a unique name. File Open Open Model File form Model File name: any unique name OK Simulation Master Modeler Set your units to mm. Options Units mm (milli newton)

专业课设,悬臂梁有限元分析

1 研究目的与问题阐述 1.1 基本研究目的 (1) 掌握ANSYS软件的基本几何形体构造、网格划分、边界条件施加等方法。 (2) 熟悉有限元建模、求解及结果分析步骤和方法。 (3) 利用ANSYS软件对梁结构进行有限元计算。 (4) 研究不同泊松比对同一位置应力的影响。 1.2 基本问题提出图1.1 模型示意图如图1.1所示，当EX=3.01e6，F=5000N，悬臂梁杆一端固定，另一端为自由端。当悬臂梁的泊松比u为：0.2、0.25、0.3、0.35、0.4时，确定同一位置的应力分布，得出分布云图。二维模型，3*0.09m。 2 软件的介绍与使用 2.1 ANSYS 简介 ANSYS程序是一个功能强大的灵活的设计分析及优化、融结构、流体、电场、磁场、声场分析于一体的大型通用有限元商用分析软件，可广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等一般工业及科学研究。该软件提供了一个不断改进的功能清单，集体包括：结构高度非线性分析、电磁分析、计算流体动力分析、设计优化、接触分析、自适应网格划分、大应变/有限转动工功能一接利用ANSYS参数设计的扩展宏命令功能。 ANSYS由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数

系统下生成的集合数据传入ANSYS,如Pro/Engineer, NASTRAN, Alogor, I－DEAS, AutoCAD等，并通过必要的修补可准确地在该模型上划分网格并求解。 2.2 ANSYS软件的功能介绍 ANSYS软件含有多种有限元分析的能力，包括从简单线性静态分析到复杂非线性动态分析。一个典型的ANSYS分析过程可分为以下三个步骤：创建有限元模型；施加载荷进行求解；查看分析结果；在有限元的分析过程中，程序通常使用以下三个部分：前处理模块，分析求解模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，通过这个模块用户可以建立自己想要的工程有限模型。分析求解模块即是对已建立好的模型在一定的载荷和边界条件下进行有限元计算，求解平衡微分方程。包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析（热-应力耦合、流-固耦合以及电-磁-热-应力耦合）等，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。下面对ANSYS软件的三种模块的功能进行简要介绍： 1.前处理模块 ANSYS软件的前处理模块主要实现三种功能：参数定义、实体建模和网格划分。（1）参数定义 ANSYS程序在进行结构建模的过程中，首先要对所有被建模型的材料进行参数定义。包括定义使用单位制，定义所有使用单元的类型，定义单元的实常数，定义材料的特性以及使用材料库文件等。（2）实体建模

matlab有限元分析实例

1.物理现象：这个对工程师来说是直观的物理现象和物理量，温度多少度，载荷是多大等等。通常来说，用户界面中呈现的、用户对工程问题进行设置时输入的都是此类信息。 2.数学方程：将物理现象翻译成相应的数学方程，例如流体对应的是NS方程，传热对应的是传热方程等等；大部分描述这些现象的方程在空间上都是偏微分方程，偶尔也有ODE（如粒子轨迹、化学反应等）。在这个层面，软件把物理现象“翻译” 为以解析式表示的数学模型。 3.数值模型：在定义了数学模型，并执行了网格剖分后，商业软件会将数学模型离散化，利用有限元方法、边界元法、有限差分法、不连续伽辽金法等方法生成数值模型。软件会组装并计算方程组雅可比矩阵，并利用求解器求解方程组。这个层面的计算通常是隐藏在后台的，用户只能通过一些求解器的参数来干预求解。有限元是一种数值求解偏微分方程的方法。基本过程大致是设置形函数，离散，形成求解矩阵，数值解矩阵，后处理之类的。 MATLAB要把这些过程均自己实现，不过在数值求解矩阵时可以调用已有函数。可以理解为MATLAB是一个通用的计算器，当然它的功能远不止如此。

而ANSYS之类的叫做通用有限元软件，针对不同行业已经将上述过程封装，前后处理也比较漂亮，甚至不太了解有限元理论的人也能算些简单的东西，当然结果可靠性又另说了。比较两者，ANSYS之类的用起来容易得多，但灵活性不如MATLAB。MATLAB用起来很困难，也有人做了一些模块，但大多数只能解决一些相对简单的问题。对于大多数工程问题，以及某些领域的物理问题，一般都用通用有限元软件，这些软件还能添加一些函数块，用以解决一些需要额外设置的东西。但是对于非常特殊的问题，以及一般性方程的有限元解，那只能用MATLAB或C，Fortran之类的了。

Matlab-PDE工具箱有限元法求解偏微分方程

在科学技术各领域中，有很多问题都可以归结为偏微分方程问题。在物理专业的力学、热学、电学、光学、近代物理课程中都可遇见偏微分方程。偏微分方程，再加上边界条件、初始条件构成的数学模型，只有在很特殊情况下才可求得解析解。随着计算机技术的发展，采用数值计算方法，可以得到其数值解。偏微分方程基本形式而以上的偏微分方程都能利用PDE工具箱求解。 PDE工具箱 PDE工具箱的使用步骤体现了有限元法求解问题的基本思路，包括如下基本步骤： 1) 建立几何模型 2) 定义边界条件 3) 定义PDE类型和PDE系数 4) 三角形网格划分

5) 有限元求解 6) 解的图形表达以上步骤充分体现在PDE工具箱的菜单栏和工具栏顺序上，如下具体实现如下。打开工具箱输入pdetool可以打开偏微分方程求解工具箱，如下首先需要选择应用模式，工具箱根据实际问题的不同提供了很多应用模式，用户可以基于适

当的模式进行建模和分析。在Options菜单的Application菜单项下可以做选择，如下或者直接在工具栏上选择，如下列表框中各应用模式的意义为： ① Generic Scalar：一般标量模式（为默认选项）。 ② Generic System：一般系统模式。 ③ Structural Mech.，Plane Stress：结构力学平面应力。 ④ Structural Mech.，Plane Strain：结构力学平面应变。

⑤ Electrostatics：静电学。 ⑥ Magnetostatics：电磁学。 ⑦ Ac Power Electromagnetics：交流电电磁学。 ⑧ Conductive Media DC：直流导电介质。 ⑨ Heat Tranfer：热传导。 ⑩ Diffusion：扩散。可以根据自己的具体问题做相应的选择，这里要求解偏微分方程，故使用默认值。此外，对于其他具体的工程应用模式，此工具箱已经发展到了Comsol Multiphysics软件，它提供了更强大的建模、求解功能。另外，可以在菜单Options下做一些全局的设置，如下 l Grid：显示网格 l Grid Spacing…：控制网格的显示位置 l Snap：建模时捕捉网格节点，建模时可以打开 l Axes Limits…：设置坐标系围 l Axes Equal：同Matlab的命令axes equal命令建立几何模型使用菜单Draw的命令或使用工具箱命令可以实现简单几何模型的建立，如下各项代表的意义分别为

悬臂梁—有限元ABAQUS线性静力学分析实例

线性静力学分析实例——以悬臂梁为例线性静力学问题是简单且常见的有限元分析类型，不涉及任何非线性（材料非线性、几何非线性、接触等），也不考虑惯性及时间相关的材料属性。在ABAQUS中，该类问题通常采用静态通用（Static，General）分析步或静态线性摄动（Static，Linear perturbation）分析步进行分析。线性静力学问题很容易求解，往往用户更关系的是计算效率和求解效率，希望在获得较高精度的前提下尽量缩短计算时间，特别是大型模型。这主要取决于网格的划分，包括种子的设置、网格控制和单元类型的选取。在一般的分析中，应尽量选用精度和效率都较高的二次四边形/六面体单元，在主要的分析部位设置较密的种子；若主要分析部位的网格没有大的扭曲，使用非协调单元（如CPS4I、C3D8I）的性价比很高。对于复杂模型，可以采用分割模型的方法划分二次四边形/六面体单元；有时分割过程过于繁琐，用户可以采用精度较高的二次三角形/四面体单元进行网格划分。悬臂梁的线性静力学分析 1.1 问题的描述一悬臂梁左端受固定约束，右端自由，结构尺寸如图1-1所示，求梁受载后的Mises应力、位移分布。 ν 材料性质：弹性模量3 = E=，泊松比3.0 2e 均布载荷：F=103N 图1-1 悬臂梁受均布载荷图 1.2 启动ABAQUS 启动ABAQUS有两种方法，用户可以任选一种。（1）在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 --

ABAQUS/CAE。（2）在操作系统的DOS窗口中输入命令：abaqus cae。启动ABAQUS/CAE后，在出现的Start Section（开始任务）对话框中选择Create Model Database。 1.3 创建部件在ABAQUS/CAE顶部的环境栏中，可以看到模块列表：Module：Part，这表示当前处在Part（部件）模块，在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建悬臂梁的几何模型。（1）创建部件。对于如图1-1所示的悬臂梁模型，可以先画出梁结构的二维截面（矩形），再通过拉伸得到。单击左侧工具区中的（Create Part）按钮，或者在主菜单里面选择Part--Create，弹出如图1-2所示的Create Part对话框。图1-2 Create Part对话框在Name（部件名称）后面输入Beam，Modeling Space（模型所在空间）设

有限元钢架结构分析手算matlabansys模拟

有限元大作业——钢架结构分析选题人：日期：2016年6月2日

目录：第一章：问题重述 (2) 一、题目内容： (3) 二、题目要求： (3) 第二章：有限元法手工求解 (3) 一、平面两单元离散化 (4) 二、单元分析 (4) 三、单元组装 (6) 四、边界条件引入及组装总体方程 (7) 五、求解整体刚度方程，计算节点2的位移和转角 (7) 六、求节点1、3支撑反力 (8) 七、设定数据，求解结果 (8) 八、绘制轴力图、弯矩图、剪力图 (9) 第三章、matlab编程求解： (11) 一、总体流程图绘制： (11) 二、输入数据： (12) 三、计算单元刚度矩阵： (12) 四、建立总体刚度矩阵： (13) 五、计算未约束点位移： (13) 六、计算支反力： (13) 七、输出数据： (13) 八、编程： (13) 第四章有限元求解 (13) 一、预处理 (13) 二、模型建立： (15) 二、分析计算 (17) 三、求解结果 (18) 四、绘制图像 (19) 第五章结果比较 (22) 第六章心得体会 (22) 第七章附录 (23) 一、matlab程序 (24) 第一章：问题重述

一、题目内容：图示平面钢架结构图题目内容二、题目要求：（1）采用平面梁单元进行有限元法手工求解，要求写出完整的求解步骤，包括： a）离散化：单元编号、节点编号； b）单元分析：单元刚度矩阵，单元节点等效载荷向量； c）单元组长：总体刚度矩阵，总体位移向量，总体节点等效载荷； d）边界条件的引入及总体刚度方程的求解； e）B点的位移，A、C处支撑反力，并绘制该结构的弯矩图、剪力图和轴力图。（2）编制通用平面钢架分析有限元Matlab程序，并计算盖提，与手工结果进行比较；（3）利用Ansys求解，表格列出B点的位移，A、C处支反力，绘制弯矩图、剪力图和轴力图，并与手算和Matlab程序计算结果比较。（4）攥写报告，利用A4纸打印；（5）心得体会，并简要说明各成员主要负责完成的工作。第二章：有限元法手工求解