有限元作业

合集下载

有限元的大作业报告示例

1.题目概况
矩形板尺寸如下图1，板厚为5mm。

材料弹性模量为
松比μ= 0.27 。

施加约束和载荷并讨论：
图
1 计算简图
1.1基本数据
E = 2⨯105N/mm2，泊
序号载荷约束备注42 向下集中载荷F=800N, 作用于cd 边3/4 处（近d） c d 点简支
1.2分析任务/分析工况
讨论板上开孔、切槽等对于应力分布的影响。

（载荷约束组合不变）。

提示：各种圆孔，椭圆孔随大小、形状、数量，分布位置变化引起的应力分布变化；各种形状，大小的切槽及不同位置引起应力分布的变化等，选择二至三种情况讨论，并思考其与机械零部件的构型的相对应关系。

2.模型建立
2.1单元选择及其分析
由于平板长宽分别为300x100，故可取网格单元大小为1。

如图：
2.2模型建立及网格划分
模型按单元为1 划分后的网格大小如图所示：
2.3载荷处理
向下集中载荷F=800N, 作用于cd 边3/4 处（近d） c d 点简支
3.计算分析
3.1位移分布及其分析
（1）位移分布如图：。

现代设计方法(关于有限元)作业

《现代设计方法》作业关于有限元法的研究学院：机械工程学院专业：机械制造及其自动化0.有限元法有限元法分析起源于50年代初杆系结构矩阵的分析。

随后，Clough于1960年第一次提出了“有限元法”的概念。

其基本思想是利用结构离散化的概念，将连续介质体或复杂结构体划分成许多有限大小的子区域的集合体，每一个子区域称为单元（或元素），单元的集合称为网格，实际的连续介质体（或结构体）可以看成是这些单元在它们的节点上相互连接而组成的等效集合体；通过对每个单元力学特性的分析，再将各个单元的特性矩阵组集成可以建立整体结构的力学方程式，即力学计算模型；按照所选用计算程序的要求，输入所需的数据和信息，运用计算机进行求解。

当前，有限元方法/理论已经发展的相当成熟和完善，而计算机技术的不断革新，又在很大程度上推进了有限元法分析在工程技术领域的应用。

然而，如此快速地推广和应用使得人们很容易忽视一个前提，即有限元分析软件提供的计算结果是否可靠、满足使用精度的前提，是合理地使用软件和专业的工程分析。

有限元法分析一般包括四个步骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算结果的分析。

每一个步骤在操作过程中都或多或少地引入了误差，这些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我们的认识和判断。

1.受内压空心圆筒的轴对称有限元分析例图1.1所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径为120mm的刚性圆孔中，试求圆筒内径处的位移。

结构的材料参数为：200=，0.3E GPaμ=。

图1 结构图对该问题进行有限元分析的过程如下。

（1）结构的离散化与编号由于该圆筒为无限长，取出中间一段（20mm高），采用两个三角形轴对称单元，如图1.2所示。

对该系统进行离散，单元编号及结点编号如图1.3所示，有关结点和单元的信息见表1.1。

图1.2 有限元模型图1.3 节点位移编号及单元编号表1.1 单元编号及结点编号单元编号结点编号 ①②1 2 3 2 3 4结构的结点位移列阵为11223344[]T r r r r u w u w u w u w δ= (1.1) 结构的结点外载列阵12[000000]T r r F F F = （1.2）1r F 和2r F 为由内压作用而等效在结点1和结点2上的载荷，其大小为1122240202//502622r r r h p F N F N ππ-⨯⨯⨯==== （1. 3）约束的支反力矩阵123344[00T z z r z r z R R R R R R R = ] (1.4)其中1z R 和2z R 为结点1和结点2在Z 方向的约束支反力，（3r R ，3z R ）和（4r R ，4z R ）为结点3和结点4在r 方向和Z 方向的约束支反力。

有限元大作业

1．推导有限元计算格式，理解有限元原理：建立图示受拉直杆在自重（设单位长度重度为q ，截面积为A ）和外力P 作用下的拉伸问题的微分方程，并分别利用不同的原理（变分求极值（最小势能或虚功原理）、加权残值法）推导有限元计算格式（取两个单元）。

手工求出端点的位移（自己给定参数值）。

设杆长为L ，截面面积为A(x)，弹性模数为E,单位长重量q ，受拉杆x 处的位移为u(x)。

取微元dx 的力平衡，建立受拉杆位移所满足的微分方程()du x dx ε=，()du x E E dxσε== dx 上下截面内力与微元自重相等得()*()()*()A x dx x dx A x x dx qdx σσ++-+=-(()())dA x x q dxσ∴=- (())d duEA x q dx dx=- 0x L << ()0u x = 0x =()duEA x p dx= x L = 得解析解：2()2q x P u Lx x EA EA=-+将其分为两个单元，节点为1,2,3，得22382qL PL u EA EA=+232qL PL u EA EA=+有限元法：1)位移函数01u α= 2111u u l α-=得1211(1)x x u u u l l =-+ 令11(1)x N l =-21x N l = 11122122u u N u N u N N u⎧⎫⎪⎪⎡⎤=+=⎨⎬⎣⎦⎪⎪⎩⎭{}1u N d ⎡⎤=⎣⎦ 2)应变、应力表达{}{}111211du dN d d dx dx l l ε⎡⎤⎡⎤===-⎢⎥⎣⎦⎣⎦{}1B d ε⎡⎤=⎣⎦ {}1E E B d σε⎡⎤==⎣⎦ {}1S d σ⎡⎤=⎣⎦3)势能表示{}{}(){}{}(){}{}{}{}{}1111''112211''121112210111111111111111121221222T V ll T T T T T U W D dV F u F u qdx u u d B E d Adx F u F u ql EA EA ql l l d d d F d EA EA ql l l εε⎡⎤=-=-+-⎣⎦+⎡⎤=-+-⎣⎦⎡⎤⎡⎤-⎢⎥⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦⎢⎥⎣⎦∏⎰⎰⎰4)单元平衡方程 a)最小势能原理110u ∂=∂∏120u ∂=∂∏111111212112112ql F u AE l u ql F ⎧⎫-⎪⎪⎧⎫⎡⎤-⎪⎪⎪⎪=⎨⎬⎨⎬⎢⎥-⎣⎦⎪⎪⎪⎪⎩⎭+⎪⎪⎩⎭b)虚位移原理{}(){}(){}TeTdd F qdx d δδδεσΩ+=Ω⎰⎰{}{}1B d σεδ⎡⎤=⎣⎦ {}1E E B d σεδ⎡⎤==⎣⎦{}(){}{}(){}111111TTT l d F d B E B d Adxδδ⎡⎤=⎣⎦⎰ 由虚位移任意性得，{}{}1111T lF B E B Adxd ⎡⎤=⎣⎦⎰ 积分得111111212112112ql F u AE l u ql F ⎧⎫-⎪⎪⎧⎫⎡⎤-⎪⎪⎪⎪=⎨⎬⎨⎬⎢⎥-⎣⎦⎪⎪⎪⎪⎩⎭+⎪⎪⎩⎭ 记为{}{}111k d F ⎡⎤=⎣⎦ 同理222212323112112ql F u AE l u ql F ⎧⎫-⎪⎪⎧⎫⎡⎤-⎪⎪⎪⎪=⎨⎬⎨⎬⎢⎥-⎣⎦⎪⎪⎪⎪⎩⎭+⎪⎪⎩⎭{}{}222k d F ⎡⎤=⎣⎦ {}{}ei i eF R =∑ 12220F F += 23F P =11111112211223222022202EAEAql F l l u ql ql EA EA EA EA u l l l l u ql EAEA P l l ⎡⎤⎧⎫-⎢⎥+⎪⎪⎢⎥⎧⎫⎪⎪⎢⎥⎪⎪⎪⎪⎢-+-⎥=+⎨⎬⎨⎬⎢⎥⎪⎪⎪⎪⎢⎥⎩⎭⎪⎪+⎢⎥⎪⎪--⎢⎥⎩⎭⎣⎦可得：22382qL PLu EA EA=+232qL PL u EA EA=+与解析解结果一致。

有限元分析大作业报告

有限元分析大作业报告一、引言有限元分析是工程领域中常用的数值模拟方法，通过将连续的物理问题离散为有限个子区域，然后利用数学方法求解，最终得到数值解。

有限元分析的快速发展和广泛应用，为工程领域提供了一种强大的工具。

本报告将介绍在大作业中所进行的有限元分析工作及结果。

二、有限元模型建立本次大作业的研究对象是工程结构的应力分析。

首先，通过对结构进行几何建模，确定了结构的尺寸和形状。

然后，将结构离散为有限个单元，每个单元又可以看作一个小的子区域。

接下来，为了求解结构的应力分布，需要为每个单元确定适当的单元类型和单元属性。

最后，根据结构的边界条件，建立整个有限元模型。

三、材料属性和加载条件在建立有限元模型的过程中，需要为材料和加载条件确定适当的参数。

本次大作业中，通过实验获得了结构材料的弹性模量、泊松比等参数，并将其输入到有限元模型中。

对于加载条件，我们选取了其中一种常见的加载方式，并将其施加到有限元模型中。

四、数值计算和结果分析为了求解结构的应力分布，需要进行数值计算。

在本次大作业中，我们选用了一种常见的有限元求解器进行计算。

通过输入模型的几何形状、材料属性和加载条件，求解器可以根据有限元方法进行计算，并得到结构的应力分布。

最后，我们通过对计算结果进行分析，得出了结论。

五、结果讨论和改进方法根据计算结果，我们可以对结构的应力分布进行分析和讨论。

根据分析结果，我们可以得出结论是否满足设计要求以及结构的强度情况。

同时，根据分析结果，我们还可以提出改进方法，针对结构的特点和问题进行相应的优化设计。

六、结论通过对工程结构进行有限元分析，我们得到了结构的应力分布，并根据分析结果进行了讨论和改进方法的提出。

有限元分析为工程领域提供了一种有效的数值模拟方法，可以帮助工程师进行结构设计和分析工作，提高设计效率和设计质量。

【1】XXX，XXXX。

【2】XXX，XXXX。

以上是本次大作业的有限元分析报告，总结了在建立有限元模型、确定材料属性和加载条件、数值计算和结果分析等方面的工作，并对计算结果进行讨论和改进方法的提出。

平面有限元法作业

在单元的边界 x = ±a 及 y = ±b 上，位移是按线性变化的，而在公共边界上有两个节点相
连，这两个公共节点有共同的节点位移值，从而保证了两个相邻单元在其公共边界上位移的连续性。故四节点矩形单元满足位移连续性条件。#
{ } 3-7：求以下受力单元的等效节点载荷 R 。已知：lij、lim、lmj 、
⎢⎣ 0 0 0
0 0.5 0
0
0 − 0.5 − 0.5 0 0.5 ⎥⎦12×12
利用矩阵的运算关系
[ ] [ ] [k]T =
B]T [D][B]tA T
= [B]T [D]T
[B]T
T
tA
由于 [D]是对称矩阵， [D]T = [D]
所以 [k]T = [B]T [D] [B]tA = [k]，即 [k]为对称矩阵。#
3-5：图示平面等腰三角形单元，若 μ = 0.3 ，弹性模量为 E，厚度为 t，求形函数矩阵 [N ]、应变矩阵 [B] 及单元刚度矩阵 [K ]。（补充题意：平面应力情况）
q、P，厚度 t，P 点作用在 jm 中点处，沿 x 方向，三角形分布载荷垂直于 ij 边。
4
解：q 的单元 N/m2 ，设厚度为 t，如图示
Xi
=
−
1 3
qlij
t
cos
30°
=
−
3 6
qlij
t
Yi
=
−
1 3
qlij
t
sin
30°
=
−
1 6
qlij
t
等效节点载荷
X
j
=
−
1 6
qlijt cos30° +
第三章作业

有限元作业整理版

一、写出下图所示的三结点三角形单元的插值函数(形函数)Ni,Nj,Nm及插值函数矩阵[N],应变矩阵[B]。

二、如下图所示的三结点三角形单元,厚度为t,弹性模量是E,泊松比μ.试求:插值函数(形函数)矩阵[N],应变矩阵[B],应力矩阵[S],单元刚度矩阵[K]e.
三、下图所示的三结点三角形单元在jm边作用有线性分布的面载荷（x方向）,试求等效结点载荷向量.
四、如下图所示,一正方形平板,厚度为t,边长为a,弹性模量E,泊松比μ.划分为两个三角形单元.求:1,3点的位移.
五、下图所示矩形板，分成四个三结点三角形单元.要求:
(1)写出由单元刚度矩阵组装总体刚度矩阵的表达式;
(2)如1234就是一个矩形单元,试求此单元的单元刚度矩阵.。

(完整word版)有限元分析大作业报告要点

有限元分析大作业报告试题1：一、问题描述及数学建模图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较：（1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；（2）分别采用不同数量的三节点常应变单元计算；（3）当选常应变三角单元时，分别采用不同划分方案计算。

该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图所示。

二、采用相同单元数目的三节点常应变单元和六节点三角形单元计算1、有限元建模（1）设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取Preferences 为Structural（2）选择单元类型：三节点常应变单元选择的类型是Solid Quad 4 node182；六节点三角形单元选择的类型是Solid Quad 8 node183。

因研究的问题为平面应变问题，故对Element behavior（K3）设置为plane strain。

（3）定义材料参数：弹性模量E=2.1e11，泊松比σ=0.3（4）建几何模型：生成特征点；生成坝体截面（5）网格化分：划分网格时，拾取lineAB和lineBC，设定input NDIV 为15；拾取lineAC，设定input NDIV 为20，选择网格划分方式为Tri+Mapped，最后得到600个单元。

（6)模型施加约束：约束采用的是对底面BC 全约束。

大坝所受载荷形式为Pressure ，作用在AB 面上，分析时施加在L AB 上，方向水平向右，载荷大小沿L AB 由小到大均匀分布。

以B 为坐标原点，BA 方向为纵轴y ，则沿着y 方向的受力大小可表示为：}{*980098000)10(Y y g gh P -=-==ρρ2、计算结果及结果分析（1）三节点常应变单元三节点常应变单元的位移分布图三节点常应变单元的应力分布图（2）六节点三角形单元六节点三角形单元的变形分布图六节点三角形单元的应力分布图①最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近，即坝底和水的交界处，且整体应力和位移变化分布趋势相似，符合实际情况；②结果显示三节点和六节点单元分析出来的最大应力值相差较大，原因可能是B点产生了虚假应力，造成了最大应力值的不准确性。

有限元作业

有限元作业一、名称:平面结构的静力学分析二、问题描述如图所示不规则板，板厚0.5，孔直径为1，倒角半径为0.4。

约束条件：在左边的孔约束全部自由度，载荷施加在右边圆的下半部分，两端施加最小值（50）逐渐均匀变大到底部的最大值（500）。

三、过程:1.选择单元格类型：新建材料类型Main Menu> Preprocessor > Material Props> Material Models> Structural > Linear > Elastic >Isotropic.并设置属性四、建立模型1.运行Main Menu> Preprocessor> Modeling> Create> Areas> Rectangle> By Dimensions分别输入：X1 = 0 X2 = 6 Y1 = -1 Y2 = 1X1 = 4 X2 = 6 Y1 = -1 Y2 = -3（1）在两个矩形的基础上依次建立两个圆形1运行Utility Menu> PlotCtrls> Pan, Zoom, Rotate2单击small dot3运行Utility Menu> WorkPlane> Display Working Plane4运行Utility Menu> WorkPlane> WP Settings点Polar .Grid. Triad5运行Main Menu> Preprocessor> Modeling> Create> Areas> Circle> Solid Circle6输入WP X = 0 WP Y = 07建立半径为1的圆.8同样的方法建立第二个圆（2）进行布尔运算运行Main Menu> Preprocessor> Modeling> Operate> Booleans> Add>Areas（3）倒圆角1.运行Utility Menu> PlotCtrls> Numbering点line numbering.2.运行Utility Menu> WorkPlane> Display Working Plane Main Menu>Preprocessor> Modeling> Create> Lines> Line Fillet选择17 和 8.3.输入半径0.44.运行Utility Menu> Plot> Lines5.运行Utility Menu> PlotCtrls> Pan, Zoom, Rotate点Zoom button.放大该区域6.运行Main Menu> Preprocessor> Modeling> Create> Areas> Arbitrary> ByLines选线4.5和17.运行Utility Menu> Plot> Areas分别建立两个孔1.运行Utility Menu> WorkPlane> Display Working Plane2.运行Main Menu> Preprocessor> Modeling> Create> Areas> Circle> SolidCircle输入WP X = 0WP Y = 03.建立半径为0.4的孔4.运行Utility Menu> WorkPlane> Offset WP to> Global Origin5.运行Main Menu> Preprocessor> Modeling> Create> Areas> Circle> SolidCircle选取点WP X = 0WP Y = 0建立第二个孔6.运行Utility Menu> WorkPlane> Display Working Plane7.运行Utility Menu> Plot> Replot和Utility Menu> Plot> Lines进行布尔运算运行Main Menu> Preprocessor> Modeling> Operate> Booleans> Add> Areas【一】设置材料属性运行Main Menu> Preprocessor> Material Props> Material Models如下图设置【二】定义实常数运行Main Menu> Preprocessor> Real Constants> Add/Edit/Delete输入【三】划分网格1.运行Main Menu> Preprocessor> Meshing> Mesh Tool Set Global Sizecontrol. 输入0.5。

有限元编程作业

1,
*Elset, elset=__PickedSurf9_S4, internal, instance=Ball-1
7, 149, 162, 292, 485, 487,……,1098, 1192, 1218, 1260
*Elset, elset=__PickedSurf9_S2, internal, instance=Ball-1
*Element, type=C3D4
1, 163, 164, 165, 166
……
1357, 297, 49, 47, 51
*Nset, nset=_PickedSet2, internal, generate
1, 318, 1
*Elset, elset=_PickedSet2, internal, generate
**定义材料Mat-Ball和Mat-Plate
*Material, name=Mat-Ball
*Density
7800.,
*Elastic
2.068e+11, 0.3
*Material, name=Mat-Plate
*Density
7800.,
*Elastic
2.078e+11, 0.3
*Element, type=C3D8R
1, 243, 244, 17, 16, 1561, 1562, 1335, 1334
……
6135, 6327, 5480, 6359, 7645, 6798, 6797, 7677
**内部节点集
*Nset, nset=_PickedSet2, internal, generate
** STEP: Step-1
**定义一般静态分析步

有限元作业

Bierenzuode,kanbudong作业1：有一个等截面两节点二力杆，杆长为L ，截面积为A ，材料弹性模量为E 。

每个节点只考虑一个水平位移，对于图 (a)、(b) 所示的坐标系统和位移插值函数，分别求相应的[B]矩阵和单元刚度矩阵[K]。

解：（a ）、212()u x x αα=+，由边界条件确定常数1α、2α：当0x =时，1i u α=；当x l =时，212j u l αα=+，可得2222()1i j x x u x u u l l ⎛⎫=-+ ⎪⎝⎭因每个节点只考虑一个水平位移故以矩阵形式表示的单元位移函数为：{}{}{}2211122222()1u u x x f x u N N u u ll ⎧⎫⎧⎫⎧⎫==-=⎨⎬⎨⎬⎨⎬⎩⎭⎩⎭⎩⎭单元的几何矩阵：[]''122222x x B N N l l ⎡⎤⎡⎤==-⎣⎦⎢⎥⎣⎦{}{}{}12x u E E B u σε⎧⎫=⨯=⎨⎬⎩⎭，即[][]D E = 对于矩形截面梁单元，积分：yzd dA =⎰⎰为单元横截面面积。

梁单元刚度矩阵[]0leT EA K B B dx ⎡⎤⎡⎤⎡⎤⎣⎦⎣⎦⎣⎦=⎰222202222l x x x l EA x ll l dx ⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦-=-⎰44334433EAEA l l EAEA ll --⎡⎤⎢⎥⎢⎥⎣⎦=（b ）、212()u x x x αα=+，由边界条件确定常数1α、2α：当2l x =-时，21224i a l a l u =-+；当2l x =时，21224j a l a l u =+可得222222()i j x lx x lxu x u u l l -+=+ 因每个节点只考虑一个水平位移故以矩阵形式表示的单元位移函数为：{}{}{}221112222222()u u x lxx lx f x u N N u u ll ⎧⎫⎧⎫⎧⎫-+===⎨⎬⎨⎬⎨⎬⎩⎭⎩⎭⎩⎭单元的几何矩阵：[]''122244x l x l B N N l l -+⎡⎤⎡⎤==⎣⎦⎢⎥⎣⎦{}{}{}12x u E E B u σε⎧⎫=⨯=⎨⎬⎩⎭，即[][]D E =对于矩形截面梁单元，积分：yzd dA =⎰⎰为单元横截面面积。

有限元分析大作业

一、有限元方法的手工计算结果与ansys分析结果的对比1分析的问题描述如图1所示，桁架的杆截面面积为8，由钢制成（E=200GPa）。

用有限元法计算出每个节点的位移以及反作用力。

(1)(2)(3)图1对于上述问题，本文将用手工计算和ansys软件分别计算出结果，对计算出来的结果进行对比。

2手工计算2.1桁架结构的有限元计算方法对于桁架结构，每个单元的刚度矩阵为，（2-1）YX图2其中，为桁架单元在整体坐标系中与X轴的夹角；，A为桁架的截面积，E 为弹性模量，L为桁架长度。

在固体力学问题中，有限元公式通常由如下的一般形式，Ku=F（2-2）其中，K为刚度矩阵，u为位移矩阵，F为载荷矩阵。

运用公式（2-3），就能求出反作用力，R=Ku-F（2-3）其中，R为反作用力矩阵。

2.2计算过程计算每个桁架单元的刚度，用公式（2-1）计算每个每个桁架单元的刚度矩阵，将每个单元放入总刚度矩阵，他们的位置分别为：10-100000 00000000 -10100000 00000000 00000000 00000000 00000000 0000000000000000 00000000 0010-1000 0000000000-101000 00000000 00000000 000000003.9-4.90000-3.9 4.9-4.9 6.10000 4.9-6.1 00000000 00000000 00000000 00000000-3.9 4.90000 3.9-4.9 4.9-6.10000-4.9 6.100000000 00000000 00000000 000 1.28000-1.28 00000000 00000000 00000000 000-1.28000 1.2800000000 00000000 00000000 00000000 0000 3.9 4.9-3.9-4.9 0000 4.9 6.1-4.9-6.1 0000-3.9-4.9 3.9 4.9 0000-4.9-6.1 4.9 6.1将个刚度矩阵相加得到总刚度矩阵为，19.9-4.90-16000-3.9 4.9 -4.9 6.100000 4.9-6.1 -160320-16000 00012.8000-12.8 00-16019.9 4.9-3.9-4.9 0000 4.9 6.1-4.9-6.1 -3.9 4.900-3.9-4.97.80 4.9-6.10-12.8-4.9-6.1025应用边界条件施加载荷，将总刚度矩阵带入式（2-2）得：19.9-4.90-16000-3.9 4.9Ux1 -4.9 6.100000 4.9-6.1Uy1 -160320-16000Ux2 00012.8000-12.8Uy200-16019.9 4.9-3.9-4.9Ux3 0000 4.9 6.1-4.9-6.1Uy3 -3.9 4.900-3.9-4.97.80Ux4 4.9-6.10-12.8-4.9-6.1025Uy4带入边界条件解得：将结果带入（2-3）得：=Fx1Fy1Fx2Fy2Fx3Fy3Fx4Fy43用ansys软件求解（单位统一N,mm，Mpa）（1）选择单元（图3）图3（2）附材料属性（图4）图4（3）创建模型（图5）图5（4）施加载荷（图6）图6（5）求解每个节点的位移（图7）图7节点的反力（图8）图8（6）模型变形图（7）位移等值线分布图4结果对比及分析手算结果ansys 计算结果位移(mm)Ux100Uy100Ux2-0.0016-0.0016Uy2-0.0468-0.0468Ux300Uy300Ux4-0.0066-0.0066Uy4-0.0317-0.0317表1手算结果ansys计算结果节点反力(N)Fx1-1027.8-1027.8 Fy11608.31608.3 Fx2 5.60 Fy2-100 Fx32066.72063.1 Fy32257.12255.4 Fx48.80 Fy4-2.70表2由表1和表2可以看出，手工计算的结果与ansys计算的结果基本一致。

有限元大作业

一、试题一1、问题描述图示无限长刚性地基上的三角形大坝，受齐顶的水压力，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下集中方案进行比较：1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；2）分别采用不同数量的三节点常应变单元计算；3）当选常应变三角形单元时，分别采用不同划分方案计算。

图 1.12、数学建模及有限元建模2.1数学建模将无限长的三维问题转化为二维问题，取出大坝的一个截面，建立模型如图2.1.1所示，图 2.1.12.2单元选择分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算，分别采用不同数量的三节点常应变单元计算，当选常应变三角形单元时，分别采用不同划分方案计算。

由于是平面应力问题，故Element Behavior K3选择Plane stress。

E=210Gpa,u=0.3。

2.3网格划分按照题意，采用不同数量的三节点常应变单元计算，划分具体的方案见第三部分计算结果。

2.4载荷及边界条件处理边界条件：底端将X，Y方向全约束。

载荷：由于水的压强是随着深度越来越大的P=密度*重力加速度*深度，故对大坝的梯梯形面施加载荷Result = 9800*(10-{Y})3、计算结果及结果分析3.1三节点常应变单元网格划分及结果图 1.3.1.1 图 1.3.1.2图 1.3.1.3 图1.3.1.4 由图可知：DMX=0.143E-4MSMN=49520PASMX=150581PA3.2六节点常应变单元网格划分及结果图1.3.2.1 图 1.3.2.2图 1.3.2.3 图 1.3.2.4由图可知：DMX=0.264E-4MSMN=1893PASMX=306986PA3.3三节点常应变单元网格划分及结果（24个单元）图 1.3.3.1 图 1.3.3.2图 1.3.3.3 图 1.3.3.4由图可知：DMX=0.203E-4MSMN=24617PASMX=199290PA3.4三节点常应变单元网格划分及结果（48个单元）图 1.3.4.1 图 1.3.4.2图 1.3.4.3 图 1.3.4.4由图可知：DMX=0.228E-4MSMN=16409PASMX=243234PA3.5三节点划分方案A图 1.3.5.1 图 1.3.5.2图 1.3.5.3 图 1.3.5.4由图可知：DMX=0.106E-4MSMN=50110PASMX=18194PA3.6三节点划分方案B图 1.3.6.1 图 1.3.6.2图 1.3.6.3 图 1.3.6.44、多方案分析比较4.1相同数目三节点和六节点单元比较DMX(mm) SMN（pa）SMX(pa) 方案数值三节点三角形单元0.143E-4 49520 150581六节点三角形单元0.264E-4 1893 306986表 1.4.1结论：相同的单元划分方案和单元规模，采用不同的阶次的相同形状的单元，分析结果不同。

有限元课程大作业

金属坯料挤压过程有限元分析一、前言：金属挤压是将放在挤压模具内的金属锭坯从一端施加外力，强迫其从特定的模孔中流出，获得所需要的断面形状和尺寸的制品的技术。

冷挤压时由于材料是在冷态下成形，而且变形量一般都很大，挤压过程中作用在模具上的单位压力很大,此时模具有开裂破坏的可能，对压力机也构成威胁，金属坯料在通过模具过程中，坯料与模具之间产生相当大的应力，这就要求模具需要有相当大的强度、硬度、以及耐磨性，因此冷挤压时要进行挤压力的计算。

挤压力的计算是模具设计的重要依据，也是选择挤压设备的依据。

模具角度、接触表面的摩擦系数、坯料变形量都会影响应力变化，在保证加工要求的前提下，应当通过适当方式降低坯料及模具之间的应力。

通过有限元分析，得出应力分布图，分析变形区域、死区，对模具进行优化改进。

二、有限元介绍：ANSYS概述ANSYS软件是融结构、热、流体、电磁、声学于一体的大型通用有限元软件，可广泛地用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、水利、日用家电等一般工业及科学研究。

该软件提供了不断改进的功能清单，具体包括：结构高度非线性分析、电磁分析、计算流体力学分析、设计优化、接触分析、自适应网格划分及利用ANSYS参数设计语言扩展宏命令功能。

ANSYS软件功能强大，主要特点有：实现多场及多场耦合分析；实现前后处理、求解及多场分析统一数据库的一体化；具有多物理场优化功能；强大的非线性分析功能;多种求解器分别适用于不同的问题及不同的硬件设备；支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容；强大的并行计算功能支持分布式并行及共享内存式并行；多种自动网格划分技术；良好的用户开发环境。

ANSYS不仅支持用户直接创建模型，也支持与其他CAD软件进行图形传递，其支持的图形传递有：SAT、Parasolid、STEP。

相应地，可以进行接口的常用CAD 软件有：Unigraphics、Pro/Engineer、I-Deas、Catia、CADDS、SolidEdge、SolidWorks等。

中科大工程中的有限元作业答案（最新）

整个物体也无转动，可推出 b1 = b2 = b3 = 0
、最小势能原理法
2
1 1 12 1 2 22 1 3 32 1 1 2 2
P = 2 k d + 2 k d + 2 k d - Fq - F q
2的位移。
d1 , d 2 , d 3分别代表弹簧 (1)，
( 2 )，( 3)的伸长量，q1和q2表示节点1，
= êêk21(1) k22(1) 0úú + êê0 k11(2) k12(2) úú
êë 0 0 0úû êë0 k21(2) k22(2) úû
k12(1)
0ù
ék11(1)
ê
(1)
(1)
(2)
= êk21 k22 + k11 k12(2) úú
êë 0
k21(2)
k22(2) úû
EA
é E1A1
ë
û
2
2
sinq cos
-cos q
-sinq cos
q
qù
é cos q
2
2
ê
AE sinq cosq
sin q
-sinq2cosq
-sin q úú
2
K=
ê
L ê -cos q
cos q
sinq cos
-sinq cos
q
qú
2
2
ê-sin cos
sinq cosq
cos q úû
ë q q -sin q
3
Þ RB = ql
8
3
R Bl
ql
, y BR =
= 3EI
8EI
4
3
ql

有限元分析题及大作业题答案

有限元分析及应用作业报告试题10一、问题描述确定图示扳手中的应力, E=210Gpa，μ=0.3, 假设厚度为10mm；并讨论采用何种处理可降低最大应力或改善应力分布。

图1为扳手的基本形状和基本尺寸图二、数学建模与分析由图1及问题描述可知，板手的长宽尺寸远远大于厚度，研究结构为一很薄的等厚度薄板，满足平面应力的几何条件；作用于薄板上的载荷平行于板平面且作用在沿厚度方向均匀分布在办手柄的左边缘线，而在两板面上无外力作用，满足平面应力的载荷条件。

故该问题属于平面应力问题，薄板所受的载荷为面载荷，分布情况及方向如图1所示，建立几何模型，并进行求解。

薄板的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3三、有限元建模1、单元选择：选取三节点常应变单元来计算分析薄板扳手的位移和应力。

由于此问题为平面应力问题，：三节点常应变单元选择的类型是PLANE42（Quad 4node42），该单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三节点单元。

2、定义材料参数：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK3、生成几何模型：a.创建关键点点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入16个点的坐标→OKb、将这16个关键点有直线依次连起来，成为线性模型4、生成实体模型：ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →连接特征点→生成两个area→Operate→Subtract→拾取整个扳手区域→OK→生成扳手模型5、结点布置及规模6、网格划分方案ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →Mesh: Areas, Shape: Tri,Free →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)7、载荷及边界条件处理8、求解控制A、模型施加约束给模型施加x方向约束ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Lines →拾取模型左部的竖直边→OKB、给模型施加载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →force→on keypoints→拾取上面左端关键点→700N/mm→okC、分析计算：ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK6)结果显示：ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… →select Def + Undeformed →OK (back to Plot Results window) →Contour Plot →Nodal Solu →select: DOF solution →displacement vector sum，von mises stress→OK四、计算结果及结果分析1、三节点常应变单元1）三节点单元的网格划分图2 常应变三节点单元的网格划分平面图图3 常应变三节点单元的网格划分立体图2）三节点单元的约束受载情况图4 常应变三节点单元的约束受载图3)三节点单元的位移分析图5 常应变三节点单元的位移分布图4)三节点单元的应力分析图6 常应变三节点单元的应力分布图2、六节点三角形单元1)六节点三角形单元网格划分图7 六节点三角形单元网格划分图2)六节点三角形单元约束和受载情况分析图8 六节点三角形单元约束受载图3）六节点三角形单元位移分析图9 六节点三角形单元的变形分布图4) 六节点三角形单元的应力分析图9 六节点三角形单元的应力分布图图10 六节点三角形单元的局部应力分布图根据以上位移和应力图，可以得出常应变三节点单元和六节点三角形单元的最小最大位移应力如表1-1所示。

北科大研究生有限元作业分析报告

作业二：有限元分析报告（1）在机械运动中通常遇到轴与轴承，轴与齿轮，以及盘与轴连接问题，在很多情况下都是采用过盈配合来连接。

本报告主要对转子轴和圆盘的连接处应力进行分析。

简图如下其他具体参数如下表：部分圆环轴圆盘内径/mm 25 34 外径/mm 35 200 厚度/mm 15025轴与盘左边距 10 弹性模量/Mpa 2.01*105 接触摩擦系数 0.2 泊松比0.3载荷：在过盈配合时轴与轴盘间有接触压力，无外力。

而在分析轴从圆盘中抽出时，对轴施加位移载荷，位移载荷大小为40。

（2）题目主要是对接触的非线性进行分析。

设置了两个载荷步分别对过盈配合进行分析，以及对轴从盘中拔出过程的接触应力进行分析。

由于模型为轴对称的，可以用轴对称方法进行分析，即为了便于后处理观察，在建模时采用1/4部分进行建模，最后将其扩展，来观察整个结构变形及应力。

圆盘圆环轴本题主要分析接触处应力，在建模时不考虑倒角与圆角，这样便于计算与划分网格。

（3）几何模型的建立。

考虑到模型为轴对称问题，在建模时不需要把整个模型都建立出来，只需要建立模型的1/4，在分析完后再映射为整个模型。

报告中模型数据已给出。

在ANSYS中通过前处理模块建模功能创建圆环实体。

利用Create-V olumes-Cylinder-Partial-Cylinder命令赋予圆盘与圆环轴的参数，先后创建圆盘与圆环轴。

进行等轴测显示Pan Zoom Rotate 中Iso，并将圆环轴移动，使其有正确位置（和圆盘Z方向有位移10）。

最终得到模型图如下。

（4）本题目选用的8节点单元，具体为Structural Solid-Brick 8 node185单元。

因为分析内容属于实体范畴，模型形状规则，简单，可以比较方便的划分为六面体单元，计算量小。

故优先选用8 node185。

（5）两接触固体均为规则的几何平面，进行网格划分时采用的线控方式。

对圆环轴以及圆盘分别对各边线进行等分，具体在size controls下Lines中Set 中进行设置。

中科大有限元作业答案(第一次到第五次)

有限元作业一1、对图示杆结构，已知节点 3 的位移 10mm，试用有限元法求出各节点的节点力和节点 2 的位移。

已知两段杆的长度为 l1, l2 截面积为 A 1, A 2 ，弹性模量为 E 1, E2 。

2.求解如图所示桁架节点 1 处的水平位移和垂直位移分量以及每一杆单元的应力。

已知所有单元 A  5106m2 , E  200GPa, L 1m。

有限元作业一1、对图示杆结构，已知节点 3 的位移 10mm，试用有限元法求出各节点的节点力和节点 2 的位移。

已知两段杆的长度为 l1, l2 截面积为 A 1, A 2 ，弹性模量为 E 1, E2 。

解：将整个杆件可以划分 3 个节点，2 个单元。

如图所示：设图中力为 F，则 1 节点所受反作用力为－F。

对单元①，单元刚度矩阵为：ke(1) 同理，单元②的单元刚度矩阵为：1 1 AE 1 1  L1 1 1  ke(2) 将单元刚度矩阵集成，有整体刚度矩阵为：A2E2  1 1  L2  1 1 k  ke(1) [k]e(2)(1) (1) k11 k12 0 0 0 0    (1) (1) (2) (2)   k21 k22 0   0 k11 k12  (2) (2)  0 k21 k22    0 0 0   (1) (1) k11 k12 0  (1) (1) (2) (2)   k21 k22  k11 k12  (2) (2)   0 k k 21 22   E1A  EA 1  1 1 0   L L1  1   E1A EA E A EA   1 1 1  2 2  2 2  L1 L2 L2   L1  EA E2 A2   2 2  0  L2 L2  F 1x  F      整体节点载荷矩阵为 R  F 2x   0  ， F  F   3x   1x  0      整体节点位移矩阵为：   2x   2x 。

有限元分析大作业报告

该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图所示。

因研究的问题为平面应变问题，故对Element behavior（K3）设置为plane strain。

（6)模型施加约束：约束采用的是对底面BC 全约束。

大坝所受载荷形式为Pressure ，作用在AB 面上，分析时施加在L AB 上，方向水平向右，载荷大小沿L AB 由小到大均匀分布。

以B 为坐标原点，BA 方向为纵轴y ，则沿着y 方向的受力大小可表示为：}{*980098000)10(Y y g gh P -=-==ρρ2、计算结果及结果分析（1）三节点常应变单元三节点常应变单元的位移分布图三节点常应变单元的应力分布图（2）六节点三角形单元六节点三角形单元的变形分布图六节点三角形单元的应力分布图单元类型最小位移（mm）最大位移（mm）最小应力（Pa）最大应力（Pa）三节点0 0.0284 5460.7 392364六节点0 0.0292 0.001385 607043①最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近，即坝底和水的交界处，且整体应力和位移变化分布趋势相似，符合实际情况；②结果显示三节点和六节点单元分析出来的最大应力值相差较大，原因可能是B点产生了虚假应力，造成了最大应力值的不准确性。

有限元作业

1 有限元的基础理论有限单元法是20世纪50年代以来随着电子计算机的广泛应用而发展起来的有一种数值解法。

有限元分析（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题有限元分析后再求解。

它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。

这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。

有限元法的基本思路和基本原则以结构力学中的位移法为基础，把复杂的结构或连续体看成有限个单元的组合，各单元彼此在节点处连接而组成整体。

把连续体分成有限个单元和节点，称为离散化。

先对单元进行特性分析，然后根据各节点处的平衡和协调条件建立方程，综合后作整体分析。

这样一分一合，先离散再综合的过程，就是把复杂结构或连续体的计算问题转化为简单单元的分析与综合的问题。

因此，一般的有限元解法包括三个主要步骤：离散化、单元分析、整体分析。

(1)离散化一个复杂的弹性体可以看作由无限个质点组成的连续体。

为了进行解算，可以将此弹性体简化为有限个单元组成的集合体，这些单元只在有限个节点上铰接，因此，这集合体只具有有限个自由度，这就为解算提供了可能。

有无限个质点的连续体转化为有限个单元的集合体，就称为离散化。

(2)单元分析单元分析首先要进行单元划分。

在工程结构中，一般采用四种类型的基本单元，即标量单元、线单元（杆、梁单元）、面单元和体单元。

四种基本单元的若干例子及各单元节点自由度（节点位移）表示在图（1-1）中。

而单元划分一般注意下面几点：一、从有限元本身来看，单元划分的越细，节点布置得越多，计算的结果越精确。

但计算时间和计算费用的增加。

所以在划分单元时对应兼顾这两个方面。

二、在边界比较曲折，应力比较集中，应力变化较大的地方，单元应划分的细点，而在应力变化平缓处单元划分的大些。

有限元作业

有限元的大作业报告示例

现代设计方法(关于有限元)作业

有限元大作业

有限元分析大作业报告

平面有限元法作业

有限元作业整理版

(完整word版)有限元分析大作业报告要点

有限元作业

有限元编程作业

有限元作业

有限元分析大作业

有限元大作业

有限元课程大作业

中科大工程中的有限元作业答案（最新）

有限元分析题及大作业题答案

北科大研究生有限元作业 分析报告

中科大有限元作业答案(第一次到第五次)

有限元分析大作业报告

有限元作业

北科大研究生有限元作业分析报告