Matlab有限元

%%%%% 真解 u=sin(pi*x)*sin(pi*y)*sin(t) %%%%% 方程 diff(u,t)-Laplace(u)=f

%%%%% f=sin(pi*x)*sin(pi*y)*cos(t)+2*pi^2*sin(pi*x)*sin(pi*y)*sin(t) %clear all % clc

%%%%finite element code for parabolic equation with constant coefficient %%%mesh%%

node=[0,0;1,0;1,1;0,1]; elem=[2,3,1;4,1,3]; T=1;

bdEdge=setboundary(node,elem,’Dirichlet’); n=input(‘Please input initial mesh:’); M=input(‘M=’); for i=1:n

[node,elem,bdEdge]=uniformrefine(node,elem,bdEdge); end

N=size(node,1); NT=size(elem,1); S=1/NT; r=1/M;

A=zeros(N,N); u=zeros(N,M+1)

有 限

元 大 作 业 程 序 设

学校：天津大学

院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\\\\\\\\\\\\ 指导老师：

计

连续体平面问题的有限元程序分析

[题目]：

如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界

上受正向分布压力，

p?1kNm，同时在沿对角线y轴上受一对集中压

力，载荷为2KN，若取板厚t?1，泊松比v?0。

2kN 1kN/m 2kN

[分析过程]：

由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。

[程序原理及实现]：

用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下：

问题类型：可用于计算弹性力学平面问题和平面应变问题 单元类型：采用常应变三角形单元 位移模

%%%%% 真解 u=sin(pi*x)*sin(pi*y)*sin(t) %%%%% 方程 diff(u,t)-Laplace(u)=f

%%%%% f=sin(pi*x)*sin(pi*y)*cos(t)+2*pi^2*sin(pi*x)*sin(pi*y)*sin(t) %clear all % clc

%%%%finite element code for parabolic equation with constant coefficient %%%mesh%%

node=[0,0;1,0;1,1;0,1]; elem=[2,3,1;4,1,3]; T=1;

bdEdge=setboundary(node,elem,’Dirichlet’); n=input(‘Please input initial mesh:’); M=input(‘M=’); for i=1:n

[node,elem,bdEdge]=uniformrefine(node,elem,bdEdge); end

N=size(node,1); NT=size(elem,1); S=1/NT; r=1/M;

A=zeros(N,N); u=zeros(N,M+1)

风vrvb

有限元部分实验报告

F0805102班 5080519046 王江

一、问题描述

一个带圆孔平板如图，内孔半径1mm，平板为方形，其边长为20mm。两侧受均布拉 伸载荷q=1000N/mm。平板材料性能参数包括：泊松比0.3，弹性模量E=200GPa。试分析平 板内部应力场。扩展讨论：当小孔直径变化时，孔边上的应力将会如何变化。

二、模型描述

2.1 模型简化

利用对称性原理，我们可以只对平板

的四分之一进行研究。

如右图所示，考虑第一象限中的平板：

对于X轴上的分应力fxx及fxy，由于对称性

可知fxy=0，且X轴上的质点在Y方向应没有

位移。 同理对于Y轴上的分应力fyx及fyy，

可由对称性推出 fyx=0，且Y轴上的质点在

X方向应没有位移。 因此可将该部分平板看

做只有一边受外载荷q，且在X轴上受Y=0，

Y轴上受X=0的边界约束。 而由对称性可知，

二、三、四象限中的平板受载荷及边界条件

情况与第一象限完全一致。因此只研究1/4

平板是合理的，与研究整体平板结果相同。

2.2、实验模型

模型单元如右上图所示，建立以（0 0 0）为圆心，（1 0 0）和（0 1 0）为边界的圆弧，再以（10 0 0）及（10 10 0）、（10 1

齿轮弯曲应力的有限元分析

摘 要：本文对有限元的概念和分析方法做了介绍，利用有限元分析软件ANSYS对UG建模的齿轮进行了分析，得出了齿轮在不同载荷下，弯曲应力的变化情况，对齿轮的设计提供了理论依据。

关键词：ANSYS；有限元；齿轮

有限元分析是使用有限元方法来分析静态或动态的物理物体或物理系统。在这种方法中一个物体或系统被分解为由多个相互联结的、简单、独立的点组成的几何模型。在这种方法中这些独立的点的数量是有限的，因此被称为有限元。由实际的物理模型中推导出来得平衡方程式被使用到每个点上，由此产生了一个方程组。这个方程组可以用线性代数的方法来求解。有限元分析的精确度无法无限提高。元的数目到达一定高度后解的精确度不再提高，只有计算时间不断提高。

有限元分析可被用来分析比较复杂的、用一般地说代数方法无法足够精确地分析的系统，它可以提供使用其它方法无法提供的结果。在实践中一般使用电脑来解决在分析时出现的巨量的数和方程组。

在分析一个物体或系统中的压力和变形时有限元分析是一种常用的手段，此外它还被用来分析许多其它问题如热传导、流体力学和电力学。

通用有限元分析软件有：德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQU

基于UG的有限元分析

1. 模型的建立

利用UG8.0/ Modeling 模块建立模型，如图1所示：

图1 模型

2. 新建有限元模型

1) 单击【开始】→【高级仿真】命令，在【仿真导航器】窗口中右击单击【Rocker.prt】节点，在出的快捷菜单中单击【新建FEM】命令，弹出【新建部件文件】对话框，默认名称、文件夹，单击【确定】按钮。

2) 弹出【新建FEM】对话框，设置求解器为 NX NASTRAN。分析类型为结构分析。单击【确定】按钮，进入了创建有限元模型的环境。 3) 单击工具栏的【材料属性】

图标，弹出【指派材料】

对话框，选择好实体模型，在【材料】列表框中单击【Steel】，

再单击【确定】按钮即完成部件材料属性设置。 4) 单击工具栏中的【物理属性】

图标，弹出图2所示的

【物理属性表管理器】对话框，单击【创建】按钮，弹出【PSOLID】（体单元）对话框，如图2所示，在【材料】列表框中选取【Steel】选项，其他选项默认，单击【确定】按钮。返回到【物理属性表管理器】对话框。单击【关闭】按钮退出。

图2 【PSOLID】对话框

5) 单击工具栏中的【网格捕集器】

图标，弹出图3所示

的【网格捕集器】对话框，在【实体属性】列表框中选取上述设置的

1. Geom即为几何体，是我们分析对象的真实模型，实际物体的三维表现形式；elems即

为网格单元，是我们分析对象的力学模型，是对实际物体的一种近似模拟，是把实际物体转换成可计算的力学和数学模型，它不是简单的线和面，是带有数据的线和面。 在HyperMesh中，我们把geom和elems统称为comps，comps可以理解为图层，这里的图层和CAD的图层的概念不同。这里comps是以后赋予模型材料和几何性质的一个最小单元，或者说对于不同材料性质和不同几何性质的elems要处于不同的comps中。每个comps都会有个名字，所以同一个名字的comps包含两个部分，即XXX（名字）geom和XXX（名字）elems。当然几何体和力学模型是两个完全独立的部分，所以两者完全可以放在不同的comps。

窗口下方是主菜单，共分7类，分别是Geom、1D、2D、3D、BCs、Tool、Post，每一类中有一些重复的比较经常使用的命令。

Geom:主要是对模型的修改和操作。 1d:主要是对线单元的修改和操作。 2d:是对面单元的修改和操作。 3d:是对固体单元的修改和操作。 BCS:边界条件。 TOOL:使用的方法。 POST:后处理的命令。

窗口右下方是对

有限元大作业

学院：机械工程学院专业：机械工程及自动化班级：姓名：学号：日期：

题号：135

题目：如图1所示的托架，其顶面承受100N的均匀分布载荷，托架通过有孔的表面固定在墙上，托架是钢制的，杨氏模量

E=1*1011Pa，泊松比v=0.3，试通过ANSYS输入其变形图及托架的von Misese应力分布。（题中in单位换成mm，如1in= 25.4mm,图中所示每个尺寸分别乘以：11）

该题的实际模型及尺寸如图2所示。

图1 显示原始尺寸的模型 图2 转化尺寸后的模型

1前处理

1.1改变文件名。单击File\\Change Jobname，弹出Change Jobname对话框，在jobname框中输入Entity Analysis，单击OK，退出Change Jobname对话框。

1.2改变工作目录。单击File\\Change directory,浏览至桌面AnsysCourse文件夹，单击确定，退出浏览文件夹对话框。

1.3改变文件标题。单击File\\change Title,弹出Change Title对话框，在Title框中输入Entity Analysis。单击OK，退出Change

1. Geom即为几何体，是我们分析对象的真实模型，实际物体的三维表现形式；elems即

为网格单元，是我们分析对象的力学模型，是对实际物体的一种近似模拟，是把实际物体转换成可计算的力学和数学模型，它不是简单的线和面，是带有数据的线和面。 在HyperMesh中，我们把geom和elems统称为comps，comps可以理解为图层，这里的图层和CAD的图层的概念不同。这里comps是以后赋予模型材料和几何性质的一个最小单元，或者说对于不同材料性质和不同几何性质的elems要处于不同的comps中。每个comps都会有个名字，所以同一个名字的comps包含两个部分，即XXX（名字）geom和XXX（名字）elems。当然几何体和力学模型是两个完全独立的部分，所以两者完全可以放在不同的comps。

窗口下方是主菜单，共分7类，分别是Geom、1D、2D、3D、BCs、Tool、Post，每一类中有一些重复的比较经常使用的命令。

Geom:主要是对模型的修改和操作。 1d:主要是对线单元的修改和操作。 2d:是对面单元的修改和操作。 3d:是对固体单元的修改和操作。 BCS:边界条件。 TOOL:使用的方法。 POST:后处理的命令。

窗口右下方是对

基于UG的有限元分析

1. 模型的建立

利用UG8.0/ Modeling 模块建立模型，如图1所示：

图1 模型

2. 新建有限元模型

1) 单击【开始】→【高级仿真】命令，在【仿真导航器】窗口中右击单击【Rocker.prt】节点，在出的快捷菜单中单击【新建FEM】命令，弹出【新建部件文件】对话框，默认名称、文件夹，单击【确定】按钮。

2) 弹出【新建FEM】对话框，设置求解器为 NX NASTRAN。分析类型为结构分析。单击【确定】按钮，进入了创建有限元模型的环境。 3) 单击工具栏的【材料属性】

图标，弹出【指派材料】

对话框，选择好实体模型，在【材料】列表框中单击【Steel】，

再单击【确定】按钮即完成部件材料属性设置。 4) 单击工具栏中的【物理属性】

图标，弹出图2所示的

【物理属性表管理器】对话框，单击【创建】按钮，弹出【PSOLID】（体单元）对话框，如图2所示，在【材料】列表框中选取【Steel】选项，其他选项默认，单击【确定】按钮。返回到【物理属性表管理器】对话框。单击【关闭】按钮退出。

图2 【PSOLID】对话框

5) 单击工具栏中的【网格捕集器】

图标，弹出图3所示

的【网格捕集器】对话框，在【实体属性】列表框中选取上述设置的