有限元分析大作业

合集下载

有限元分析大作业 展示

有限元分析大作业 展示

二 不同单元对计算结果的影响
实例:一简支梁,各项数据如表所示:
q / kN l / m b / m h / m E / GPa

0.25
Hale Waihona Puke 100161
3
200
分析其变形和 受力情况。 (分析均取梁 结构的一半进 行计算)
对本例分别采用3节点三角形、4节点矩形、8节点六 面体、20节点六面体进行计算分析,且用平面单元分析时, 还考虑网格的划分不同对结果的影响,分别对梁高划分15 和30等份进行计算, 比较结果。 用平面plane 42 单元建立模型,建模图如下:
谱分析整体结果列表
之后读取结果,在File>Read Input from 里面选择 “.MCOM” 文件,查看各个结果。
谱分析6个节点的反作用力表
读取第一组结果的变形和应力图,如下图 所示:
第一组结果的变形图
第一组结果的应力图
3. 时程分析 实例:上题的三层框架结构,施加了两个x方向大小随时 间变化的力(如下图所示),进行时程分析,最后选择节 点,绘制时程图。
梁高划分30等份三角形3节点建模图
梁高划分15等份三角形3节点建模图
梁高划分30等份4节点矩形建模 图
梁高划分15等份4节点矩形建模图
用六面体8节点单元solid 45和六面体20节点单元solid 95进行分析,建模图如下:
六面体8节点建模图
六面体20节点建模图
各个分析数据均是结构的x方向的应力图,最后进过 汇总,所得的结果如下:
计算分以下几歩完成: 1. 设置为模态分析,进行模态分析求解; 2. 设置为谱分析,输入单点响应谱的频率和 谱值之后(阻尼为0),获得谱解; 3. 再设置为模态分析,扩展模态,扩展5阶模 态,进行扩展模态求解; 4. 最后设置成谱分析,进行模态叠加,选择 SRSS方法,输出类型为Displacement , 最后求解,得出结果,整体结果如下图:

有限元分析大作业

有限元分析大作业

超静定梁的有限元分析本文分别通过材料力学解法和有限元解法,求出了超静定梁的支反力、最大位移及最大位移出现位置,并对两者进行了比较和误差分析。

一、超静定梁的材料力学解法梁的约束反力数目超过了有效平衡方程数,单纯使用静力平衡不能确定全部未知力的梁称为超静定梁。

超静定梁比静定梁有许多优点,如可用较少材料获得较大的刚度和强度,个别约束破坏后仍可工作等。

因而超静定梁在工程中得到较多的应用。

超静定梁的解法有很多种,本文采用力法的一种——变形比较法求解未知量。

图1图2选取C 点的支座为多余约束,Rc 为多余支座反力,则相应的基本静定梁为一外伸梁,如图2所示,其上受集中载荷P 、均布载荷q 和多余支座反力Rc 的作用。

相应的变形条件为:c cP cq cRc f f f f =++=其中316cP B Pl f l EI θ=⨯= 4724cq ql f EI =-323c cRc R l f EI =则316Pl EI 4724ql EI -+323c R l EI=0 将已知数据带入可求得 6.25c R =- 负号表示c R 的方向与假设的方向相反。

再列出平衡方程:0X =∑AX R =0A M =∑ 232022B C ql Pl R l R l ---=0C M =∑ 232022AY B ql PllR R l +--=带入已知条件求得:AX R = 393.75AY R = 812.5B R =用叠加法求最大位移:最大的向下位移在A 与B 两点中间:334410.7910481632C R l Pl ql f EI EI EI -=-++=-⨯最大的向上位移在B 与C 两点中间:3344213490.22525103248512C R l Pl ql f EI EI -=--=⨯二、超静定梁的有限元解法在ANSYS 平台上,求解超静定梁。

建模、单元划分、加载后结果如图3所示。

图3求解后可以通过图形和列表两种方式查看结果。

有限元分析大作业报告

有限元分析大作业报告

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

有限元分析大作业模板

有限元分析大作业模板

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

《有限元分析》大作业小组成员:Job name:完成日期:一、问题描述(要求:应结合图对问题进行详细描述,同时应清楚阐述所研究问题的受力状况和约束情况。

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

)二、数学模型(要求:针对问题描述给出相应的数学模型,应包含示意图,示意图中应有必要的尺寸数据;如进行了简化等处理,此处还应给出文字说明。

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

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

)3.3 材料模型(要求:指出选择的材料模型,包括必要的参数数据。

)3.4 网格划分方案(要求:指出网格划分方法,网格控制参数,最终生成的单元总数和节点总数,此外还应附上最终划分好的网格截图。

)3.5 载荷及边界条件处理(要求:指出约束条件和载荷条件。

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

)五、多方案计算比较(要求:节点规模增减对计算精度的影响分析、单元改变对计算精度的影响分析、不同网格划分方案对计算结果的影响分析等,至少应选择其一进行分析,此外还应附上相应截图及数据。

)附件1:小组成员工作说明(要求:明确说明小组各个成员在本次大作业中所做的工作,工作内容将作为口试提问的依据之一,同时也作为成绩评定的依据之一。

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

)附件2:详细的计算过程说明(按照上机指导的格式撰写)。

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

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

平面刚架问题如图示,一个平面刚架右端固定,在左端施加一个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]按钮,完成第五点的生成,系统显示生成的五个点的位置,如图校验所输入的节点坐标的正确与否,可以选择菜单路径:弹出如图所示一个文本窗口,如图,列出了所有的节点及其坐标。

有限元分析大作业

有限元分析大作业

一、有限元方法的手工计算结果与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. 定义分析类型
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列的关键点。

华科大有限元分析题及大作业题答案――船海专业_图文.

华科大有限元分析题及大作业题答案――船海专业_图文.

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

二、几何建模与分析图1-2力学模型由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。

因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况及方向如图1-2所示,建立几何模型,进行求解。

假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3三、第1问的有限元建模本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。

1设置计算类型:两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural2选择单元类型:三节点常应变单元选择的类型是PLANE42(Quad 4node42,该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类型是PLANE183(Quad 8node183,该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。

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

3定义材料参数4生成几何模a. 生成特征点b.生成坝体截面5网格化分:划分网格时,拾取所有线段设定input NDIV 为10,选择网格划分方式为Tri+Mapped,最后得到200个单元。

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

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

有限元分析及应用作业报告试题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种请况可知随着网格的细分,所得的应变以及应力的结果是收敛的。

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

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

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

二、几何建模与分析图1-2 力学模型由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。

因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况及方向如图1-2所示,建立几何模型,进行求解。

假设大坝的材料为钢,则其材料参数:弹性模量E=,泊松比σ=¥三、第1问的有限元建模本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。

1)设置计算类型:两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural2)选择单元类型:三节点常应变单元选择的类型是PLANE42(Quad 4node42),该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类型是PLANE183(Quad 8node183),该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。

因研究的问题为平面应变问题,故对Element behavior(K3)设置为plane strain。

3)定义材料参数4)生成几何模a. 生成特征点b.生成坝体截面*5)网格化分:划分网格时,拾取lineAB和lineBC进行Size Conrotls,设定input NDIV 为15;拾取lineAC,设定input NDIV 为20,选择网格划分方式为Tri+Mapped,最后得到600个单元。

有限元分析报告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)个人按上机指南步骤至少选择习题中3个习题独立完成,并将计算结果上交;也可根据自己科研工作给出计算实例。

2)以小组为单位完成有限元分析计算;3)以小组为单位编写计算分析报告;4)计算分析报告应包括以下部分:A、问题描述及数学建模;B、有限元建模(单元选择、结点布置及规模、网格划分方案、载荷及边界条件处理、求解控制)C、计算结果及结果分析(位移分析、应力分析、正确性分析评判)D、多方案计算比较(结点规模增减对精度的影响分析、单元改变对精度的影响分析、不同网格划分方案对结果的影响分析等)E、建议与体会4)1月8日前必须完成,并递交计算分析报告(报告要求打印)。

习题及上机指南:(试题见上机指南)例题1 坝体的有限元建模与受力分析例题2 平板的有限元建模与变形分析例题1:平板的有限元建模与变形分析计算分析模型如图1-1 所示, 习题文件名: plane图1-1 受均布载荷作用的平板计算分析模型1.1进入ANSYS程序→ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: plane→Run1.2设置计算类型ANSYS Main Menu: Preferences →select Structural →OK1.3选择单元类型ANSYS Main Menu: Preprocessor →Element T ype→Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element T ypes window)→Options… →select K3: Plane stress w/thk →OK→Close (the Element T ype window)1.4定义材料参数ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OK1.5定义实常数ANSYS Main Menu: Preprocessor →Real Constants… →Add… →select T ype 1→OK→input THK:1 →OK →Close (the Real Constants Window)1.6生成几何模型✓生成特征点ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS→依次输入五个点的坐标:input:1(0,0),2(1,0), 3(1,1),4(0,1),5(0.5,0.5) →OK✓生成平板ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS→连接特征点1,2,5 →Apply →连接特征点2,3,5 →Apply →连接特征点3,4,5 →Apply →连接特征点4,1,5 →OK1.7网格划分ANSYS Main Menu: Preprocessor →Meshing→Mesh T ool →(Size Controls) lines: Set →Pick All(in Picking Menu) →input NDIV:1→OK→(back to the mesh tool window)Mesh: Areas, Shape: Tri, Free →Mesh →Pick All (in Picking Menu) →Close( the Mesh T ool window)1.8模型施加约束✓给模型施加x方向约束ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement→On Lines→拾取模型左部的竖直边:Lab2: UX →OK✓施加y方向约束ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement→On Keypoints→拾取4# 特征点:Lab2: UY →OK1.9 分析计算ANSYS Main Menu: Solution →Solve →Current LS→OK(to close the solve Current Load Step window) →OK1.10 结果显示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, UX,UY, Def + Undeformed →OK1.11 退出系统ANSYS Utility Menu: File→Exit →Save Everything→OK例题2:坝体的有限元建模与应力应变分析计算分析模型如图2-1 所示, 习题文件名: dam。

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

《有限元分析及应用》大作业——齿根弯曲应力计算报告
班级:无可奉告
姓名:无可奉告
学号:无可奉告
指导老师:无可奉告
目录
目录 (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问题分析
大小齿轮啮合,小齿轮受载荷情况较为严峻,故分析对象应当为小齿轮。

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

小齿轮齿宽b1=48mm,远大于单齿齿高,所以问题又可以简化为平面应变问题。

最终,研究单个齿,划分平面网格,成了我们解决问题的思路。

2.建模过程
本次大作业使用的三维建模软件为Solidworks 2016,前处理软件为HyperMesh (版本为10.0),求解器为系统自带求解器,后处理软件为HyperView。

2.1几何建模
由于齿轮曲线为渐开线,直接在CAE软件中建模比较复杂,而CAD软件中已有成熟的齿轮
工具,故采用CAD建模导入CAE分析的方式。

首先进行几何建模,由于solidworks软件自带齿轮工具箱,故直接设置参数生成需要的齿轮。

只取一个齿。

齿根加入圆角后导入Hypermesh。

2.2CAE网格划分与计算取中性面。

划分网格。

受载荷和应力集中处网格划分细密一些。

质量检查。

注:这里的零件尺寸过小,按最小网格尺寸检查全部不合格,同时其他选项全部合格,最小网格尺寸又没办法改到我们需要的值,故不考虑网格尺寸的影响。

建立材料,认为材料各向同性。

查百度百科,20CrMnTi的弹性模量E=207GPa,泊松比0.25,密度7.8e3kg/m3。

建立厚度属性。

加载荷和约束。

这里将下面三条边线约束住,靠近齿顶处加力。

根据题目条件,工作转矩为T1=457876N·mm
齿顶圆直径为
Da=(Z+2ha*)m=(29+2)*2.5mm=77.5mm
Da:齿顶圆直径,m:模数,z:齿数,ha*:齿顶高系数
故受力为
F=2*T1/Da=457876*2/77.5N=11816N
均匀加在齿顶8个节点上,每个节点力为1477N
进行计算。

2.3后处理
max=5.185e2=518.5MPa,位置在齿根,与经验一致。

hypermesh默认的单位为吨、秒、毫米,算出来的应力单位为兆帕。

3.多方案比较与结果分析
3.1多方案比较
直接生成的齿轮没有齿根圆角,直接对其进行网格划分与分析,结果如下:
极限应力只有478MPa,齿根应力集中处的网格划分密度虽然比带圆角的模型更加细密,但是精度仍然不够,可以看到,应力集中的部位被完全局限在了某一个网格中,如果将那个网格再度细分,算出来的最大应力还会继续提高。

因为边缘存在突变,理论上突变处的内部应力将会非常大,所以这种建模与划分方式不精确。

实际情况下,也要在齿根处增加过渡圆角。

对原模型齿根处的网格进一步细分,结果如下:
极限应力仍在齿根处,其值增大到了647MPa,说明最开始的网格划分在应力集中处仍旧不够细密。

进一步细分,结果如下:
极限应力为675MPa。

再进一步细分,齿根最大应力不会继续增大。

说明齿根的极限应力在675MPa左右。

说明这一结果是比较可靠的。

3.2结果分析
该算例中,通过有限元方法计算得到的齿根弯曲极限应力为675MPa左右,与书本上按照经验公式计算得出的极限应力613.3MPa较为接近。

最后我们依然会得出齿根抗弯强度略显不足的结论。

但是这略微不足的强度是否能满足工程上的要求,仍需实验验证。

只是我们目前还没有这样的条件。

相关文档
最新文档