有限元分析报告大作业

合集下载

有限元的大作业报告示例

有限元的大作业报告示例

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)位移分布如图:。

有限元分析报告大作业

有限元分析报告大作业

《有限元分析》大作业基本要求:1. 以小组为单位完成有限元分析计算,并将计算结果上交;2. 以小组为单位撰写计算分析报告;3. 按下列模板格式完成分析报告;4. 计算结果要求提交电子版,一个算例对应一个文件夹,报告要求提交电子版和纸质版。

《有限元分析》大作业小组成员:储成峰李凡张晓东朱臻极高彬月Job name :banshou完成日期:2016-11-22一、问题描述(要求:应结合图对问题进行详细描述,同时应清楚阐述所研究问题的受力状况和约束情况。

图应清楚、明晰,且有必要的尺寸数据。

)如图所示,为一内六角螺栓扳手,其轴线形状和尺寸如图,横截面为一外接圆半径为0.01m的正六边形,拧紧力F为600N,计算扳手拧紧时的应力分布<图1扳手的几何结构(要求:针对问题描述给出相应的数学模型,应包含示意图,示意图中应有必要的尺寸数据;数学模型图2数学模型如图二所示,扳手结构简单,直接按其结构进行有限元分析。

三、有限元建模3.1单元选择(要求:给出单元类型,并结合图对单元类型进行必要阐述, 包括节点、自由度、 实常数等。

)图3单元类型扳手截面为六边形,采用4节点182单元,182单元可用来对固体结构进行如进行了简化等处理,此处还应给出文字说明。

)Figure 1B2.1: PLANE1S2 Geometry二维建模。

182单元可以当作一个平面单元,或者一个轴对称单元。

它由4个结点组成,每个结点有2个自由度,分别在x,y方向。

扳手为规则三维实体,选择8节点185单元,它由8个节点组成,每个节点有3个自由度,分别在x,y,z方向。

3.2实常数(要求:给出实常数的具体数值,如无需定义实常数,需明确指出对于本问题选择的单元类型,无需定义实常数。

)因为该单元类型无实常数,所以无需定义实常数3.3材料模型(要求:指出选择的材料模型,包括必要的参数数据。

)对于三维结构静力学,应力主要满足广义虎克定律,因此对应ANSYS^的线性,弹性,各项同性,弹性模量EX 2e11 Pa,泊松比PRXY=0.33.4几何建模由于扳手结构比较简单,所以可以直接在ANSYS软件上直接建模,在ANSYS建立正六边形,再创立直线,面沿线挤出体,得到扳手几何模型VULUHES 1TYPE NUMEZ图4几何建模3.5网格划分方案(要求:指出网格划分方法,网格控制参数,最终生成的单元总数和节点总数,此外还应附上最终划分好的网格截图。

有限元分析大作业报告

有限元分析大作业报告

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

【1】XXX,XXXX。

【2】XXX,XXXX。

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

有限元分析大作业

有限元分析大作业

有限元大作业一题目要求:图1所示为一悬臂梁,在端部承受载荷,材料弹性模量为E,泊松比为1/3,悬臂梁的厚度(板厚)为t,若该粱被划分为两个单元,单元和节点编号如图所示,试按平面应力问题计算各个节点位移计支反力。

一、单元划分1.计算简图及单元划分如下所示:2.进行节点及单元编号节点i j m单元① 2 3 4② 3 2 13.节点坐标值节点号1 2 3 4坐标值X 2 2 0 0Y 1 0 1 0二、计算单元刚度矩阵1、计算每个单元面积△以及i b ,i c (m j i i ,,=) ①②单元的面积相等,即12121=⨯⨯=∆ 单元①的i b ,i c⎩⎨⎧=--==-=0)(1m j i m j i y x c y y b ⎩⎨⎧=--==-=2)(0i m ji m j x x c y y b ⎩⎨⎧-=--=-=-=2)(1j i mj i m y x c y y b 对平面应力问题,其表达式为[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+-+-+-+∆-=s r s r sr s r s r s r s r s r b b uc c cb u b uc b c u c ub c c u b b u Et Krs 21212121)1(42 然后对单元①求解单元刚度子矩阵2==i r 2==i s []⎥⎦⎤⎢⎣⎡=3/1001329)1(22Et K 2==i r 3==j s []⎥⎦⎤⎢⎣⎡=03/23/20329)1(23Et K2==i r 4==m s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)1(24Et K 3==j r 3==j s []⎥⎦⎤⎢⎣⎡=4003/4329)1(33Et K 3==j r 2==i s []⎥⎦⎤⎢⎣⎡=03/23/20329)1(32Et K 3==j r 4==m s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)1(34Et K 4==m r 4==m s []⎥⎦⎤⎢⎣⎡=3/133/43/43/7329)1(44Et K 4==m r 2==i s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)1(42Et K 4==m r 3==j s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)1(43Et K由子矩阵[]e rs K 合成单元刚度矩阵[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=3/133/443/23/13/23/43/73/23/43/2143/24003/23/23/403/43/203/13/203/23/103/213/2001329)1(Et K将单元①的单元刚度矩阵补零升阶变为单元刚度矩阵,其在总体刚度矩阵中的位置为:节点号→单元②的i b ,i c⎩⎨⎧=--=-=-=0)(1m j im j i y x c y y b ⎩⎨⎧-=--==-=2)(0i m ji m j x x c y y b ⎩⎨⎧=--==-=2)(1j i mj i m y x c y y b 然后对单元 求解单元刚度子矩阵:3==i r 3==i s []⎥⎦⎤⎢⎣⎡=3/1001329)2(33Et K 3==i r 2==j s []⎥⎦⎤⎢⎣⎡=03/23/20329)2(32Et K 3==i r 1==m s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)2(31Et K 1 2 3 412[])1(22K[])1(23K[])1(24K3[])1(32K[])1(33K[])1(34K4[])1(42K[])1(43K[])1(44K2==j r 2==j s []⎥⎦⎤⎢⎣⎡=4003/4329)2(22Et K 2==j r 3==i s []⎥⎦⎤⎢⎣⎡=03/23/20329)2(23Et K 2==j r 1==m s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)2(21Et K 1==m r 1==m s []⎥⎦⎤⎢⎣⎡=3/133/43/43/7329)2(11Et K 1==m r 3==i s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)2(13Et K 1==m r 2==j s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)2(12Et K 由子矩阵[]e rs K 合成单元刚度矩阵[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=3/133/443/23/13/23/43/73/23/43/2143/24003/23/23/403/43/203/13/203/23/103/213/2001329)2(Et K将单元②的单元刚度矩阵补零升阶变为单元贡献矩阵,其在总体刚度矩阵中的位置为:节点号→1 2 3 41 [])2(11K[])2(12K[])2(13K2 [])2(21K[])2(22K[])2(23K3 [])2(31K [])2(32K [])2(33K 4三、计算总体刚度矩阵总体刚度矩阵是由各单元的贡献矩阵迭加而成)2()1(][][][][K K K K e +==∑四、进行节点约束处理根据节点约束情况,在总刚矩阵中可采用划行划列处理约束的方法,由题目易知,节点3和4的已知水平位移和垂直位移都为零,划去其相对应的行和列,则总刚矩阵由8阶变为4阶,矩阵如下:⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------2/02/03/13043/203/73/23/443/23/133/43/23/43/43/73292211p p v u v u Et329][Et K =1 2 3 413/133/43/43/743/23/23/4----3/13/23/21----000243/23/23/4----3/13003/73/43/403/13/23/21----33/13/23/21----3/43/403/13003/743/23/23/4----40003/13/23/21----43/23/23/4----3/133/43/43/7化简⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------Et p Et p v u v u 3/1603/160130122072412213424472211 五、求解线性方程组方法:采用LU 分解法 1.求解矩阵[]U 各元素⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------75/10775/640075/6475/353007/767/27/7502447~7/877/87/7607/87/337/207/767/27/7502447~13012207241221342447⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----353/44900075/6475/353007/767/27/7502447~ 得到的[]U 矩阵如下:[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=353/44900075/6475/353007/767/27/7502447U 2.求解矩阵[]L 各元素⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----13012207241221342447353/44900075/6475/353007/767/27/75024471353/6475/767/20175/27/40017/40001 得到的[]L 矩阵如下:[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=13012207241221342447L3.进行求解⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⇒⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=Et p Et p Et p y Et p Et p Ly 79425/850800225/323/1603/1603/160⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⇒=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡Et p Et p Et p v u v u y v u v u U 79425/850800225/323/160353/44900075/6475/353007/7675/27/750244722112211 解得Et p v /422.82-= Et p u /497.12-= Et p v /028.91-= Et p u /897.11=于是求得各节点的位移为:⎩⎨⎧-==Etp v Etp u /028.9/897.111 ⎩⎨⎧-=-=Etp v Etp u /422.8/497.122 ⎩⎨⎧==033v u ⎩⎨⎧==044v u 六、求解相应的支反力(运用静力学的平衡方程进行求解)3号节点和4号节点的支反力如下图所示:。

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

(完整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点产生了虚假应力,造成了最大应力值的不准确性。

有限元分析大作业

有限元分析大作业

《有限元分析及应用》大作业——齿根弯曲应力计算报告班级:无可奉告姓名:无可奉告学号:无可奉告指导老师:无可奉告目录目录 (2)1.概述 (3)1.1工程问题描述 (3)1.2问题分析 (3)2.建模过程 (4)2.1几何建模 (4)2.2CAE网格划分与计算 (5)2.3后处理 (8)3.多方案比较与结果分析 (9)3.1多方案比较 (9)3.2结果分析 (11)1.概述1.1工程问题描述我在本次作业中的选题为齿根弯曲应力的计算与校核。

通过对机械设计的学习,我们可以知道,齿轮的失效形式主要是齿面接触疲劳和齿根弯曲断裂,而闭式传动硬齿面齿轮的失效形式以齿根弯曲断裂,这个时候进行齿根弯曲应力的校核才比较有意义,在设计问题的时候应当选取这种类型的算例。

设计计算的另一个主要思路是将有限元计算的结果与传统机械设计的结算结果进行对比,以从多方面验证计算结果的准确性。

综上,我们最终选取了《机械原理》(第三版)P50例3-1中的问题进行校核计算。

已知起重机械用的一对闭式直齿圆柱齿轮,传动,输入转速n1=730r/min,输入功率P1=35kW,每天工作16小时,使用寿命5年,齿轮为非对称布置,轴的刚性较大,原动机为电动机,工作机载荷为中等冲击。

z1=29,z2=129,m=2.5mm,b1=48mm,b2=42mm,大、小齿轮均为20CrMnTi,渗碳淬火,齿面硬度为58~62HRC,齿轮精度为7级,试验算齿轮强度。

齿面为硬齿面,传动方式为闭式传动。

根据设计手册查出的许用接触应力为1363.6Mpa,计算结果为1260Mpa,强度合格。

根据设计手册查出的许用弯曲应力为613.3MPa,计算结果为619Mpa,强度略显不够。

1.2问题分析大小齿轮啮合,小齿轮受载荷情况较为严峻,故分析对象应当为小齿轮。

可以看出,由于齿轮单侧受载荷,传动过程中每个齿上载荷的变化过程是相同的,故问题可被简化为反对称问题,仅需研究单个齿。

有限元分析报告

有限元分析报告

有限元分析大作业计算分析报告A、问题描述及数学建模;B、有限元建模(单元选择、结点布置及规模、网格划分方案、载荷及边界条件处理、求解控制)C、计算结果及结果分析(位移分析、应力分析、正确性分析评判)D、多方案计算比较(结点规模增减对精度的影响分析、单元改变对精度的影响分析、不同网格划分方案对结果的影响分析等)E、建议与体会试题1图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析,并对以下几种计算方案进行比较:1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算;(X)2)分别采用不同数量的三节点常应变单元计算;(V)3)当选常应变三角单元时,分别采用不同划分方案计算。

(V)1 •有限元建模单元选择、结点布置及规模、网格划分方案、载荷及边界条件处理、求解控制1)单元选择:由于ANSYS没有提供三角形常应变单元,故采用六节点三角形单元进行计算。

2)结点布置:(0,0)(6,0)(10,0)(3,0)(0,5)(3,5)3)单元数目:44)网格划分方案万案1 万案25)边界条件底边加上UX , UY的约束6)载荷受齐顶的水压力作用,呈阶梯状分布,载荷函数为F=1000 (10-Y)7)求解控制:默认,单一载荷步2 •计算结果及结果分析(位移分析、应力分析、正确性分析评判)1)位移结果与应力分析方1 :最大位移为0.109X10-05最大应力为15936最小应力为5181万案2:最大位移为0.130X10-05最大应力为15058最小应力为78342)正确性分析从应力分布图中,我们比较这两种网格划分方案的优劣,方案1的应力最大位置在(0,0)处,方案2的应力最大位置在(6,0)处,显然,方案1更贴近实际情况,因为其左下角的单元与整体单元的受力情况相似,而方案2则有较大的差别。

但是,由于这两种网格划分都非常粗燥,根据常识,在坝顶处,其受力为0,应力应接近于0,而这两种情况的最小应力分别为5181和7834。

ansys有限元分析报告大作业

ansys有限元分析报告大作业

有限元大作业设计题目: 单车的设计及ansys有限元分析专业班级:姓名:学号:指导老师:完成日期: 2016.11.23单车的设计及ansys模拟分析一、单车实体设计与建模1、总体设计单车的总体设计三维图如下,采用pro-e进行实体建模。

在建模时修改proe默认单位为国际主单位(米千克秒 mks)Proe》文件》属性》修改2、车架车架是构成单车的基体,联接着单车的其余各个部件并承受骑者的体重及单车在行驶时经受各种震动和冲击力量,因此除了强度以外还应有足够的刚度,这是为了在各种行驶条件下,使固定在车架上的各机构的相对位置应保持不变,充分发挥各部位的功能。

车架分为前部和后部,前部为转向部分,后部为驱动部分,由于受力较大,所有要对后半部分进行加固。

二、单车有限元模型 1、材料的选择单车的车身选用铝合金(6061-T6)T6标志表示经过热处理、时效。

其属性如下:弹性模量:)(2N/m 1090E .6泊松比:0.33质量密度:)(2N/m 32.70E + 抗剪模量:)(2N/m 1060E .2+ 屈服强度:)(2N/m 875E .2+ 2、单车模型的简化为了方便单车的模拟分析,提高电脑的运算效率,可对单车进行初步的简化;单车受到的力的主要由车架承受,因此必须保证车架能够有足够的强度、刚度,抗振的能力,故分析的时候主要对车架进行分析。

简化后的车架如下图所示。

3、单元体的选择单车车架为实体故定义车架的单元类型为实体单元(solid )。

查资料可以知道3D 实体常用结构实体单元有下表。

4、网格的划分1)建模使用proe,导出为IGES,再将IGES格式的车架导入ansys中file→import→IGES→找到格式为IGES的车架文件,文件导入后如下图。

2)选择单元体Main Menu→Preprocessor→Element Type→Add/Edit/Delete→Add(弹出对话框选择下图所示的单元体)3)定义车架材料的属性Main Menu →Preprocessor →Material Props →Material Models (按照下图给车架定义材料属性)定义材料的弹性模量:)(2N/m 1090E .6+、泊松比:0.33定义材料的质量密度:)(2N/m 32.70E +4)画分网格Main Menu→Preprocessor→Meshing →Mesh Tool→弹出的对话框中设置网格的边界长度设为10→点Mesh→pick all(自动生成网格如下图所示)三、单车静强度分析1)导入已画分好的模型网格2)进行菜单过滤Main Menu→Preferences→弹出的对话框中勾选Structural→点OK3)定义分析类型Main Menu→Solution→Analysis Type→New Analysis→勾选Static→点OK4)定义求解控制菜单选择Main Menu→Solution→Analysis Type→Sol’n Controls命令→打开对话框勾选Calculate prestress effects项→OK5)施加约束选择Main Menu→Solution→Define Loads→Apply→Structural→Displacement→OnAreas→选择下边的横杆→点OK→弹出的对话框中选择All DOF→OK完成约束6)施加载荷假设驾驶人的重量为60kg,那么作用在单车上的力的大小就为600N。

有限元分析报告作业(8)

有限元分析报告作业(8)

中原工学院2011-2012-2《有限元基础》大作业姓名:学号:成绩:六角板手在向下载荷作用下应力分析一个六角板手(截面高度10mm),几何和材料特性:截面高度=10mm;形状=六角形;长度=7.5cm;把手长度=20cm;弯曲半径=1cm;弹性模量=2.07x10e11Pa;向下荷载=20 N。

分析板手在这种荷载作用下的应力分布。

1 设置分析标题1、选择“Utility Menu>File>Change Title”。

2、输入“Static Analysis of an Allen Wrench”并按“OK”。

2 设置单位1、点击“ANSYS Input window”右下角。

2、输入“/UNITS,SI”命令然后按回车。

注意在“ANSYS Input window”的输入行上方出现了这个命令。

3、选择“Utility Menu>Parameters>Angular Units”,出现“Angular Units for ParametricFunctions”对话框。

4、在其中选择“Degrees DEG”。

5、按“OK”。

3 定义参数1、选择“Utility Menu>Parameters>Scalar Parameters”。

出现“Scalar Parameters”对话框。

2、输入下表的参数。

在定义每个参数后按“Accept”按钮。

参数值描述EXX 2.07E11 Young‘smodulusis2.07E11PaW_HEX .01 Widthofhexacrossflats=.01mW_FLAT W_HEX*TAN(30) Widthofflat=.0058mL_SHANK .075 Lengthofshank(shortend).075mL_HANDL.2 Lengthofhandle(longend).2mEBENDRAD .01 Bendradius.01mL_ELEM .0075 Elementlength.0075mNO_D_HE2 Numberofdivisionsalonghexflat=2XTOL 25E-6 Toleranceforselectingnode=25E-6m 注意--可以用大写或小写字母,但ANSYS总是用大写显示。

有限元分析大作业精华-实验报告

有限元分析大作业精华-实验报告

平面刚架问题如图示,一个平面刚架右端固定,在左端施加一个y方向的-3000N的力P1,中间施加一个Y方向的-1000N的力P2,试以静力来分析,求解各接点的位移。

已知组成刚架的各梁除梁长外,其余的几何特性相同。

横截面积:A=0.0072 m²横截高度:H=0.42m惯性矩:I=0.0021028m4x弹性模量:E=2.06x10n/ m²/泊松比:u=0.3用ANSYS 分析平面刚架1.设定分析模块选择菜单路径:MainMenu—preference 弹出“PRreferences for GUI Filtering”对话框,如图示,在对话框中选取:Structural”,单击[OK]按钮,完成选择。

2.选择单元类型并定义单元的实常数(1)新建单元类型并定(2)定义单元的实常数在”Real Constants for BEAM3”对话框的AREA中输入“0。

0072”在IZZ中输入“0。

0002108”,在HEIGHT中输入“0.42”。

其他的3个常数不定义。

单击[OK]按钮,完成选择3.定义材料属性在”Define Material Model Behavier”对话框的”Material Models Available”中,依次双击“Structural→Linear→Elastic→Isotropic”如图在如下图的对话框EX中输入“2.06e11”,在PRXY框中输入“0.3”,完成材料模型的定义。

4建立平面刚架节点和单元(1)生成节点选择菜单路,生成节点于目前坐标系统命令,单击以后弹出如图对话框,在对话框的Node number 中输入“1“接着依序输入第一点XYZ的坐标值”2,0,0,然后单击[APPL Y]按钮继续生成第二点,如图接着在对话框的Node number 中输入“2“接着依序输入第一点XYZ的坐标值”6,0,0,然后单击[APPL Y]按钮继续生成第三点,如图接着在对话框的Node number 中输入“3“接着依序输入第一点XYZ的坐标值”0,2,0,然后单击[APPL Y]按钮继续生成第四点,如图接着在对话框的Node number 中输入“4“接着依序输入第一点XYZ的坐标值”4,2,0,然后单击[APPL Y]按钮继续生成第五点,如图接着在对话框的Node number 中输入“5“接着依序输入第一点XYZ的坐标值”8,2,0,然后单击[OK]按钮,完成第五点的生成,系统显示生成的五个点的位置,如图校验所输入的节点坐标的正确与否,可以选择菜单路径:弹出如图所示一个文本窗口,如图,列出了所有的节点及其坐标。

有限元受力分析大作业

有限元受力分析大作业

1. 定义分析类型
Main Menu >Solution > Analysis Type > New Analysis, 选择Static 2. 定义位移约束
Main Menu >Solution> Define Loads >Apply >Structural>Displacement > On Areas, 在 弹出对话框中选Circle,点击OK后,在弹出的对话框中填入:0
后处理
1. 查看总体变形
Main Menu >General Postproc >Plot Results > Deformed shape
(2)
(2)
四. 加载和求解
1. 定义分析类型
Main Menu >Solution > Analysis Type > New Analysis, 选择Static 2. 定义位移约束
Main Menu >Solution> Define Loads >Apply >Structural>Displacement > On Areas, 在 弹出对话框中选Circle,点击OK后,在弹出的对话框中填入:0
七. 后处理
1. 查看总体变形
Main Menu >General Postproc >Plot Results > Deformed shape
S
作业三:轴承座盈利分析
1、 题目
2、建模
3、ansys前处理
导入
ansys前处理
1、定义单元类型 Main Menu>Proprocessor>Element Type >Add/Edit/Delete 弹出对话框 中后,点“Add” 。双弹出对话框,选“Solid”和“10node 92”,点 “OK”,退回到前一个对话框。

有限元分析技术课程大作业

有限元分析技术课程大作业

有限元分析技术课程大作业1 工程介绍现需要对某露天大型玻璃平面舞台的钢结构进行分析,该钢结构布置在xy 平面内。

学生序号为079,分格的列数(x向分格)=0×10+7+5=12,分格的行数(y向分格)=9+4=13,共有156个分格。

每个分格x方向尺寸为1m,y方向尺寸为1m。

钢结构的主梁为高160宽100厚14的方钢管;次梁为直径60厚10的圆钢管(单位为毫米),材料均为碳素结构钢Q235;该结构固定支撑点位于左右两端主梁和最中间的次梁的两端。

玻璃采用四点支撑与钢结构连接(采用四点支撑表明垂直作用于玻璃平面的面载荷将传递作用于玻璃所在钢结构分格四周的节点处,表现为点载荷;试对在垂直于玻璃平面方向的22/KN m的面载荷(包括玻璃自重、钢结构自重、活载荷(人员与演出器械载荷)、风载荷等)作用下的舞台进行有限元分析(每分格面载荷对于每一支撑点的载荷可等效于0.5KN的点载荷)。

作业提交的内容:(1)屏幕截图显示该结构的平面布置结构;(2)该结构每个支座的支座反力;(3)该结构节点的最大位移及其所在位置;(4)对该结构中最危险单元(杆件)进行强度校核。

2有限元模型的建立该钢结构中每一分格x方向尺寸为1m,y方向尺寸为1m,x方向分格数量为12,y方向分格数量为13。

该钢结构由主梁和次梁构成,其中主梁为高160mm、宽100mm、厚14mm的方钢管,次梁为直径60mm、厚10mm的圆钢管。

由于在该结构中所有构件均为梁单元,而Ansys程序中提供了多种梁单元,以模拟不同场合的应用,且对于每种梁单元类型都有特定的算法。

在本次建模过程中,考虑到需要对该结构中的危险单元进行强度校核,因此,选择了BEAM188单元类型来建立本钢架结构,进而对其进行有限元分析。

BEAM188为三维线性有限应变梁单元,该单元基于铁木辛哥的梁结构理论,考虑了剪切变形的影响,能够满足本次分析的需求。

以下为基于ANSYS图形界面(Graphic User Interface , GUI)的菜单操作流程(1) 进入ANSYS(设定工作目录和工作文件)程序→ ANSYS → ANSYS Interactive → Working directory(设置工作目录)→Initial Jobname(设置工作文件名):Analysis → Run → OK(2) 设置计算类型ANSYS Main Menu:Preferences → Structural → OK(3) 定义单元类型ANSYS Main Menu:Preprocessor → Element Type → Add/Edit/Delete... → Add → Beam: 3D 2node 188 → OK(返回到Element Types窗口)→ Close(4) 定义材料参数ANSYS Main Menu: Preprocessor → Material Props → Material Models → Structural → Linear → Elastic → Isotropic → input EX: 2.0E5, PRXY: 0.3(定义泊松比及弹性模量) → OK → Close(关闭材料定义窗口)(5)定义梁单元截面ANSYS Main Menu:Preprocessor →Sections→Beam→Common Sections→Beam Tool(6) 构造梁模型生成舞台几何模型ANSYS Main Menu:Preprocessor → Modeling → Create → Keypoints → In Active CS → NPT Keypoint number:1,X,Y,Z Location in active CS:0,0,0 → Apply 通过复制关键点操作,形成14行13列的关键点。

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

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

有限元分析及应用作业报告试题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所示。

有限元分析报告Ansys大作业

有限元分析报告Ansys大作业

有限元分析作业作业名称扳手静态受力分析姓名学号宁波理工学院班级题目:扳手静态受力分析:扳手的材料参数为:弹性模量E=210GPa,泊松比u=0.3:此模型在左侧内六角施加固定位移约束,在右侧表面竖直方向上施加648 N的集中力。

10模型如下图:1-11.定义工作文件名和文件标题(1)定义工作文件名:执行File-Chang Jobname-3090601048(2)定义工作标题:执行File-Change Tile-3090601048(3)更改工作文件储存路径:执行File-Chang Directory-E:\ANSYS2.定义分析类型、单元类型及材料属性(1)定义分析类型,执行Main Menu-Preferences,如下图所示:2-1(2)定义单元类型,执行Main Menu-Preprocessor-Element Type-Add 弹出Element Type 对话框.如下图所示:2-2(3)定义材料属性执行Main menu-Preprocessor-Material Props-Material models,在Define material model behavior对话框中,双击Structual-Linear-Elastic-Isotropic.如下图所示:2-33.导入几何模型将模型导入到ANSYS,执行File-Import—PRAR…—浏览上述模型,如下图所示:3-13-24. 网格划分执行Main Menu-Preprocessor-meshing-Mesh Tool命令,考虑到零件的复杂性,采用智能网格划分,精度为1,其他选项为默认,如下图所示:4-14-25. 加载以及求解(1)添加位置约束执行Solution-apply-structural-displacement-on areas(对六角内表面进行约束),如下图所示:5-15-2(2)添加载荷,执行Solution-apply-structural-force-on keypoints,如下图所示:5-35-4(3)求解执行Main menu-Solution-Solve-Current LS,求解。

有限元分析大作业报告要点

有限元分析大作业报告要点

有限元分析大作业报告试题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)六节点三角形单元六节点三角形单元的变形分布图六节点三角形单元的应力分布图(3)计算数据表单元类型最小位移(mm)最大位移(mm)最小应力(Pa)最大应力(Pa)三节点0 0.0284 5460.7 392364六节点0 0.0292 0.001385 607043 (4)结果分析①最大位移都发生在A点,即大坝顶端,最大应力发生在B点附近,即坝底和水的交界处,且整体应力和位移变化分布趋势相似,符合实际情况;②结果显示三节点和六节点单元分析出来的最大应力值相差较大,原因可能是B点产生了虚假应力,造成了最大应力值的不准确性。

有限元大作业

有限元大作业

有限元大作业第一篇:有限元大作业有限元应力分析报告大作业机械与运载工程学院车辆四班龙恒 20110402415 2014年8月30日一、问题描述桦木板凳材料参数如图形状参数:长40mm,宽30mm,高45mm(其他详细参数见零件图)通过施加垂直于板凳上表面的均匀载荷600N,分析板凳的应变和应力?二、使用inventor进行建模及应力分析1、通过inventor建立板凳3D模型利用草图拉伸等方法建立与零件图中尺寸一致的三维立体板凳模型2、点选环境下的应力分析开始对板凳进行应力分析3、根据所给条件设置材料等参数、将安全系数设为屈服强度,因为板凳主要受压变形点开“木材(桦木)”根据前面所给参数对其进行参数设置4、固定约束如图板凳的4个脚底面设置为固定约束,使得板凳受载后,脚底面不会沿垂直方向位移,模拟真实情况5、施加载荷在板凳上表面施加大小为600N的垂直均布载荷(这里是模拟一个成人坐上去的重力)6、划分网格通过设置网格的尺寸参数来划分出5种不同网格数量,从而得出5种不同网格数划分得出的应力应变分布图,最后分析划分不同网格数对结果的影响。

(1)网格最大(2)网格较大(3)网格一般大小(4)网格较小(5)网格最小7、求解得出结果得出5组不同网格数所得数据(应力云图,应变云图,所有结果数据)(1)网格数1437根据应力云图可知,红色地方所受的应力最大,最大应力为:15.48Mpa 根据应变云图可知,红色地方的应变最大,最大应变为:0.001434μl(2)网格数8651根据应力云图可知,红色地方所受的应力最大,最大应力为:18.88Mpa 根据应变云图可知,红色地方的应变最大,最大应变为:0.001755μl(3)网格数20484根据应力云图可知,红色地方所受的应力最大,最大应力为:22.62Mpa 根据应变云图可知,红色地方的应变最大,最大应变为:0.002103μl(4)网格数41578根据应力云图可知,红色地方所受的应力最大,最大应力为:23.76Mpa 根据应变云图可知,红色地方的应变最大,最大应变为:0.002206μl(5)网格数68788根据应力云图可知,红色地方所受的应力最大,最大应力为:25.97Mpa 根据应变云图可知,红色地方的应变最大,最大应变为:0.002454μl综合上述5种请况可知随着网格的细分,所得的应变以及应力的结果是收敛的。

有限元分析报告Ansys大作业

有限元分析报告Ansys大作业

有限元分析作业作业名称扳手静态受力分析姓名学号宁波理工学院班级题目:扳手静态受力分析:扳手的材料参数为:弹性模量E=210GPa,泊松比u=0.3:此模型在左侧内六角施加固定位移约束,在右侧表面竖直方向上施加648 N的集中力。

10模型如下图:1-11.定义工作文件名和文件标题(1)定义工作文件名:执行File-Chang Jobname-3090601048(2)定义工作标题:执行File-Change Tile-3090601048(3)更改工作文件储存路径:执行File-Chang Directory-E:\ANSYS2.定义分析类型、单元类型及材料属性(1)定义分析类型,执行Main Menu-Preferences,如下图所示:2-1(2)定义单元类型,执行Main Menu-Preprocessor-Element Type-Add 弹出Element Type 对话框.如下图所示:2-2(3)定义材料属性执行Main menu-Preprocessor-Material Props-Material models,在Define material model behavior对话框中,双击Structual-Linear-Elastic-Isotropic.如下图所示:2-33.导入几何模型将模型导入到ANSYS,执行File-Import—PRAR…—浏览上述模型,如下图所示:3-13-24. 网格划分执行Main Menu-Preprocessor-meshing-Mesh Tool命令,考虑到零件的复杂性,采用智能网格划分,精度为1,其他选项为默认,如下图所示:4-14-25. 加载以及求解(1)添加位置约束执行Solution-apply-structural-displacement-on areas(对六角内表面进行约束),如下图所示:5-15-2(2)添加载荷,执行Solution-apply-structural-force-on keypoints,如下图所示:5-35-4(3)求解执行Main menu-Solution-Solve-Current LS,求解。

有限元分析报告与应用大作业

有限元分析报告与应用大作业

有限元分析及应用大作业课程名称: 有限元分析及应用班级:姓名:试题2:图示薄板左边固定,右边受均布压力P=100Kn/m作用,板厚度为0.3cm;试采用如下方案,对其进行有限元分析,并对结果进行比较。

1)三节点常应变单元;(2个和200个单元)2)四节点矩形单元;(1个和50个单元)3)八节点等参单元。

(1个和20个单元)图2-1 薄板结构及受力图一、建模由图2-1可知,此薄板长和宽分别为2m和1.5m,厚度仅为0.3cm,本题所研究问题为平面应力问题。

经计算,平板右边受均匀载荷P=33.33MPa,而左边被固定,所以要完全约束个方向的自由度,如图2-2所示。

取弹性模量E=2.1×11Pa,泊松比μ=0.3。

P=33.33MPa图2-2 数学模型二、第一问三节点常应变单元(2个和200个单元)三节点单元类型为PLANE42,设置好单元类型后,实常数设置板厚为0.3M。

采用2个单元的网格划分后的结果如图2-3,200个单元的网格划分图如图2-6所示。

约束的施加方式和载荷分布如图2-2中所示。

约束右边线上节点全部自由度。

计算得到的位移云图分别如图2-4、7所示,应力云图如图2-5、8所示。

图2-3 2个三角形单元的网格划分图图2-4 2个三角形单元的位移云图图2-5 2个三角形单元的应力云图图2-6 200个三角形单元的网格划分图图2-7 200个三角形单元的位移云图图2-8 200个三角形单元的应力云图三、第二问四节点矩形单元的计算四节点单元类型为PLANE42,设置好单元类型后,实常数设置板厚为0.3M。

采用1个单元的网格划分后的结果如图2-9,50个单元的网格划分图如图2-12所示。

约束的施加方式和载荷分布如图2-2中所示。

约束右边线上节点全部自由度。

计算得到的位移云图分别如图2-10、11所示,应力云图如图2-13、14所示。

图2-9 1个四边形单元的网格划分图图2-10 1个四边形单元的位移云图图2-11 1个四边形单元的应力云图图2-12 50个四边形单元的网格划分图图2-13 50个四边形单元的位移云图图2-14 50个四边形单元的应力云图四、第三问八节点等参单元的计算四节点单元类型为PLANE82,设置好单元类型后,实常数设置板厚为0.3M。

有限元分析大作业

有限元分析大作业

基于ANSYS软件的有限元分析报告机制1205班杜星宇U201210671一、概述本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变.然后与实际情况进行比较,证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据,并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想,对有限元分析方法的实际应用有一个大致的认识。

二、问题分析已知:桌子几何尺寸如图所示,单位为mm。

假设桌子的四只脚同地面完全固定,桌子上存放物品,物品产生的均匀分布压力作用在桌面,压力大小等于300Pa,其中弹性模量E=9。

3GPa,泊松比μ=0。

35,密度ρ=560kg/m3,分析桌子的变形和应力.将桌脚固定在地面,然后在桌面施加均匀分布的压力,可以看作对进行平面应力分析,桌脚类似于梁单元。

由于所分析的结构比较规整且为实体,所以可以将单元类型设为八节点六面体单元。

操作步骤如下:1、定义工作文件名和工作标题(1)定义工作文件名:执行UtilityMenu/ File/ChangeJobname,在弹出Change Jobname 对话框修改文件名为Table。

选择New log anderrorfiles复选框。

(2)定义工作标题:Utility Menu/File/Change Title,将弹出ChangeTit le对话框修改工作标题名为The analysis of table。

(3)点击:Plot/Replot。

2、设置计算类型(1)点击:Main Menu/Preferences,选择Structural,点击OK。

3、定义单元类型和材料属性(1)点击:Main Menu/Preprocessor/Element Type/Add/Edit/Delete,点击Add,选择Solid〉Brick 8node 185,点击OK,点击Close。

(2)点击Main menu/preprocessor/Material Props/Material Models / Structural/ Linear/ Elastic/Isotropic,设置EX为9.3e9,PRXY为0。

有限元分析大作业报告

有限元分析大作业报告

有限元分析大作业报告试题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)六节点三角形单元六节点三角形单元的变形分布图六节点三角形单元的应力分布图单元类型最小位移(mm)最大位移(mm)最小应力(Pa)最大应力(Pa)三节点0 0.0284 5460.7 392364六节点0 0.0292 0.001385 607043①最大位移都发生在A点,即大坝顶端,最大应力发生在B点附近,即坝底和水的交界处,且整体应力和位移变化分布趋势相似,符合实际情况;②结果显示三节点和六节点单元分析出来的最大应力值相差较大,原因可能是B点产生了虚假应力,造成了最大应力值的不准确性。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

有限元分析》大作业基本要求:1.以小组为单位完成有限元分析计算,并将计算结果上交;2.以小组为单位撰写计算分析报告;3.按下列模板格式完成分析报告;4.计算结果要求提交电子版,一个算例对应一个文件夹,报告要求提交电子版和纸质版。

有限元分析》大作业小组成员:储成峰李凡张晓东朱臻极高彬月Job name :banshou完成日期:2016-11-22一、问题描述(要求:应结合图对问题进行详细描述,同时应清楚阐述所研究问题的受力状况和约束情况。

图应清楚、明晰,且有必要的尺寸数据。

)如图所示,为一内六角螺栓扳手,其轴线形状和尺寸如图,横截面为一外接圆半径为0.01m的正六边形,拧紧力F为600N,计算扳手拧紧时的应力分布图1 扳手的几何结构数学模型要求:针对问题描述给出相应的数学模型,应包含示意图,示意图中应有必要的尺寸数据;图 2 数学模型如图二所示,扳手结构简单,直接按其结构进行有限元分析。

三、有限元建模3.1 单元选择要求:给出单元类型, 并结合图对单元类型进行必要阐述, 包括节点、自由度、 实常数等。

)图 3 单元类型如进行了简化等处理,此处还应给出文字说扳手截面为六边形,采用4 节点182单元,182 单元可用来对固体结构进行二维建模。

182单元可以当作一个平面单元,或者一个轴对称单元。

它由4 个结点组成,每个结点有2 个自由度,分别在x,y 方向。

扳手为规则三维实体,选择8 节点185单元,它由8 个节点组成,每个节点有3 个自由度,分别在x,y,z 方向。

3.2 实常数(要求:给出实常数的具体数值,如无需定义实常数,需明确指出对于本问题选择的单元类型,无需定义实常数。

)因为该单元类型无实常数,所以无需定义实常数3.3材料模型(要求:指出选择的材料模型,包括必要的参数数据。

)对于三维结构静力学,应力主要满足广义虎克定律,因此对应ANSYS中的线性,弹性,各项同性,弹性模量EX:2e11 Pa, 泊松比PRXY=0.33.4几何建模由于扳手结构比较简单,所以可以直接在ANSYS软件上直接建模,在ANSYS建立正六边形,再创立直线,面沿线挤出体,得到扳手几何模型图4 几何建模3.5网格划分方案(要求:指出网格划分方法,网格控制参数,最终生成的单元总数和节点总数,此外还应附上最终划分好的网格截图。

)图5 网格划分设定截面每条边线分割份数为3,单元边长为0.01 ,采用映射的划分单元方法最大节点数为2145,单元总数为1000,网格划分如图4 所示3.5载荷及边界条件处理(要求:指出约束条件和载荷条件。

)四、计算结果及结果分析(要求:此处包括位移分析、应力分析、支反力分析等,应附上相应截图及数据,此外还应对正确性进行分析评判。

)。

变形图:如图7 所示,根据材料力学相关分析初步判断,扳手手柄朝力的方向一定的弯曲变形。

图7 变形图位移分析:如图8 所示,扳手手柄向力的方向位移,位移最大在手柄处节点,最大位移:0.514E-03m图8 位移图应力分析:如图 9 所示,根据材料力学分析,应力主要集中在扳手的弯角处,最大应力:0.146E+09Pa附件 1:小组成员工作说明(要求: 明确说明小组各个成员在本次大作业中所做的工作, 工作内容将作为口试提问的依 据之一, 同时也作为成绩评定的依据之一。

需注意,附件 1 的撰写应由小组成员共同完成。

) 储成峰:启动 ANSYS 程序,定义分析类型,单元类型。

张晓东:材料属性,创建几何模型,划分网格。

朱臻极:添加约束,添加载荷。

李凡:处理变形图,应力图,位移图。

高彬月: 计算结果及结果分析, 文档的整理。

附件 2:详细的计算过程说明(按照上机指导的格式撰写)1. 启动 ANSYS 程序开始 — —所 有程序 — — ANSYS ED — — ANSYS Product Launcher — — SimulationEnvironment 下拉框中选择 ANSYSE D , Product 下拉框中选 ANSYSE D ; File Management 选项图 9 应力图卡:设置Working Directory 和Job Name 命名为banshou ——Run。

2.设定分析类型Main Menu ——Preferences ——Structural ——OK。

3.定义单元类型Main Menu——Preprocessor ——Element Type——Add——Structural Solid ,Quad4 Node182 ——Apply —Brick 8node 185--OK 。

4.定义材料属性Main Menu——Preprocessor ——Material Props ——Material Models ——Structural ——Linear ——Elastic ——Isotropic ——EX=2e11,PRXY=0.3.5.创建正六边形拾取菜单main menu —preprocessor —modeling —creat —areas —polygon —hexagon,在wpx,wpy,radius 文本框中分别输入0,0,0.01 。

6.改变视点拾取菜单uility menu —plotctrls —pan zoom rotate ,单击iso ,fit 。

7.显示关键点,线号拾取菜单uility menu —plotctrls —numbering ,把关键点和线号打开,单击ok。

8.创建关键点拾取菜单main menu —preprocessor —modeling —creat —keypoints —in active cs ,在npt 输入7,x,y,z ,输入0,0,0 ,单击apply ,在npt 中输入8,x ,y ,z,输入0,0,0.05, 在npt 中输入9,x,y,z,输入0,0.1,0.05 ,单击ok。

9.创建直线拾取菜单main menu —preprocessor —modeling —creat —lines —lines —straight line ,选取关键点7,8 和8,9 创建两条直线。

10.创建圆角拾取菜单main menu—preprocessor —modeling —creat —lines —linesfillet ,拾取直线7,8 单击ok ,在rad 输入0.01511.创建直线拾取菜单main menu —preprocessor —modeling —creat —lines —lines —straight line ,选取关键点1,4 ,创建直线。

12.将六边形划分为两部分拾取菜单main menu —preprocessor —modeling —operate —Booleans —divide —area by line ,分别拾取六边形面和1,4 间的直线。

13,划分单元拾取菜单main menu—preprocessor —meshing —meshtool ,单击size controls 的lines 后的set ,拾取直线2,3,4 ,单击ok,在NDIV 文本框中输入3,单击apply ,再次拾取直线实用文档7,9,8 单击ok,删除NDIV文本框中的3,在SIZE 文本框中输入0.01 ,单击ok,在mesh 区域中选择单元形状为quad,选择划分单元的方法为mapped,单击mesh,弹出拾取窗口,拾取六边形面的两部分,单击ok14,由面沿直线挤出体拾取菜单main menu —preprocessor —modeling —operate —extrude —areas —alonglines ,分别拾取六边形面的两部分和直线7,9,815,清除面单元拾取菜单main menu—preprocessor —meshing —clear —areas ,拾取z=0 的两个面,点击ok16 . 施加约束拾取菜单main menu —solution —defineloads —apply —structural —displacement —on areas ,弹出拾取窗口,拾取z=0 的两个平面,单击‘ ok',在列表中选择“ ALLDOF”单击“ ok ”拾取菜单main menu —solution —define loads —apply —structural —force/moment —on keypoints, 弹出拾取窗口,拾取扳手长臂端面的六个顶点,单击ok ,选择lab 为fx ,在value 文本框中输入100,单击ok 。

17,求解拾取菜单main menu —solution —solve —current ls —点击solve current load step 中的ok ,出现solution is ok ,可以查看结果18,后处理(1)变形图Main menu ——General Postproc ——Plot Result ——deformed shape(2)应力图Main menu ——General Postproc ——Plot Result ——contour plot ——Nodal solution ——von mises stress(3)位移图Main Menu → General Postproc → Contour Results → Nodal Solution →(Contour Nodal Solution Data 对话框)DOF Solution → Displacement vector sum→ OK。

标准文案。

相关文档
最新文档