有限元分析及应用大作业

合集下载

ansys有限元分析工程实例大作业

ansys有限元分析工程实例大作业

ansys有限元分析工程实例大作业————————————————————————————————作者:————————————————————————————————日期:辽宁工程技术大学有限元软件工程实例分析题目基于ANSYS钢桁架桥的静力分析专业班级建工研16-1班(结构工程)学号 471620445姓名日期 2017年4月15日基于ANSYS钢桁架桥的静力分析摘要:本文采用ANSYS分析程序,对下承式钢桁架桥进行了有限元建模;对桁架桥进行了静力分析,作出了桁架桥在静载下的结构变形图、位移云图、以及各个节点处的结构内力图(轴力图、弯矩图、剪切力图),找出了结构的危险截面。

关键词:ANSYS;钢桁架桥;静力分析;结构分析。

引言:随着现代交通运输的快速发展,桥梁兴建的规模在不断的扩大,尤其是现代铁路行业的快速发展更加促进了铁路桥梁的建设,一些新建的高速铁路桥梁可以达到四线甚至是六线,由于桥面和桥身的材料不同导致其受力情况变得复杂,这就需要桥梁需要有足够的承载力,足够的竖向侧向和扭转刚度,同时还应具有良好的稳定性以及较高的减震降噪性,因此对其应用计算机和求解软件快速进行力学分析了解其受力特性具有重要的意义。

1、工程简介某一下承式简支钢桁架桥由型钢组成,顶梁及侧梁,桥身弦杆,底梁分别采用3种不同型号的型钢,结构参数见表1,材料属性见表2。

桥长32米,桥高5.5米,桥身由8段桁架组成,每个节段4米。

该桥梁可以通行卡车,若只考虑卡车位于桥梁中间位置,假设卡车的质量为4000kg,若取一半的模型,可以将卡车对桥梁的作用力简化为P1,P2,和P3,其中P1=P3=5000N,P2=10000N,见图2,钢桥的形式见图1,其结构简图见图3。

图1钢桥的形式图2桥梁的简化平面模型(取桥梁的一半)图3刚桁架桥简图所用的桁架杆件有三种规格,见表1表1 钢桁架杆件规格杆件截面号形状规格顶梁及侧梁 1 工字形400X400X12X12桥身弦杆 2 工字形400X300X12X12底梁 3 工字形400X400X16X16所用的材料属性,见表2表2 材料属性参数钢材弹性模量EX泊松比PRXY 0.3密度DENS 78002 模型构建将下承式钢桁梁桥的各部分杆件,包括顶梁及侧梁,桥身弦杆,底梁均采用BEAM188单元,此空间梁单元可以考虑所模拟杆件的轴向变形; 定义了一套材料属性,各类杆件为钢材,其对应的参数如表2所示;根据表1中的杆件规格定义了三种梁单元截面,根据表1分别定义在相应的梁上;建模时直接建立节点和单元,在后续按照先建节点再建杆的次序一次建模。

有限元分析大作业 展示

有限元分析大作业 展示

二 不同单元对计算结果的影响
实例:一简支梁,各项数据如表所示:
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】XXX,XXXX。

【2】XXX,XXXX。

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

(完整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问题分析大小齿轮啮合,小齿轮受载荷情况较为严峻,故分析对象应当为小齿轮。

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

有限元法理论及应用参考答案(推荐文档)

有限元法理论及应用参考答案(推荐文档)

有限元法理论及应用大作业1、试简要阐述有限元理论分析的基本步骤主要有哪些?答:有限元分析的主要步骤主要有:(1)结构的离散化,即单元的划分;(2)单元分析,包括选择位移模式、根据几何方程建立应变与位移的关系、根据虚功原理建立节点力与节点位移的关系,最后得到单元刚度方程;(3)等效节点载荷计算;(4)整体分析,建立整体刚度方程;(5)引入约束,求解整体平衡方程。

2、有限元网格划分的基本原则是什么?指出图示网格划分中不合理的地方。

题2图答:一般选用三角形或四边形单元,在满足一定精度情况,尽可能少一些单元。

有限元划分网格的基本原则:1.拓扑正确性原则。

即单元间是靠单元顶点、或单元边、或单元面连接2.几何保持原则。

即网络划分后,单元的集合为原结构近似3.特性一致原则。

即材料相同,厚度相同4.单元形状优良原则。

单元边、角相差尽可能小5.密度可控原则。

即在保证一定精度的前提下,网格尽可能的稀疏一些。

(a)(b)中节点没有有效的连接,且(b)中单元边差相差很大。

(c)中没有考虑对称性,单元边差很大。

3、分别指出图示平面结构划分为什么单元?有多少个节点?多少个自由度?题3图答:(a )划分为杆单元, 8个节点,12个自由度。

(b )划分为平面梁单元,8个节点,15个自由度。

(c )平面四节点四边形单元,8个节点,13个自由度。

(d )平面三角形单元,29个节点,38个自由度。

4、什么是等参数单元?。

答:如果坐标变换和位移插值采用相同的节点,并且单元的形状变换函数与位移插值的形函数一样,则称这种变换为等参变换,这样的单元称为等参单元。

5、在平面三节点三角形单元中,能否选取如下的位移模式,为什么?(1).⎪⎩⎪⎨⎧++=++=26543221),(),(y x y x v yx y x u αααααα (2). ⎪⎩⎪⎨⎧++=++=2652423221),(),(yxy x y x v yxy x y x u αααααα 答:(1)不能,因为位移函数要满足几何各向同性,即单元的位移分布不应与人为选取的 坐标方位有关,即位移函数中的坐标x,y 应该是能够互换的。

有限元分析及应用大作业

有限元分析及应用大作业

有限元分析及应用大作业作业要求:1)个人按上机指南步骤至少选择习题中3个习题独立完成,并将计算结果上交;也可根据自己科研工作给出计算实例。

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

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

因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况P=98000-9800*Y;建立几何模型,进行求解;假设大坝的材料为钢,则其材料参数:弹性模量E=,泊松比σ=;2:有限元建模过程:进入ANSYS :程序→ANSYS APDL设置计算类型 :ANSYS Main Menu: Preferences →select Structural →OK选择单元类型 :ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182,六节点三角形单元选择Solid Quad 8node 183) →OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window)定义材料参数:ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:, PRXY: →OK生成几何模型:生成特征点:ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标:input:1(0,0),2(10,0),3(1,5),4,5) →OK生成坝体截面:ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点,1(0,0),2(6,0),3(0,10) →OK网格划分:ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)模型施加约束:给底边施加x和y方向的约束:ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK给竖直边施加y方向的分布载荷:ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量;2) 在Result窗口中出现{X},写入所施加的载荷函数:98000-9800*{Y}; 3) File>Save(文件扩展名:func) →返回:Parameters →Functions →Read from file:将需要的.func文件打开,参数名取meng,它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边;OK →在下拉列表框中,选择:Existing table →OK →选择需要的载荷为meng参数名→OK分析计算:ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current LoadStep window) →OK结果显示:ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… → select Def + Undeformed →OK (back to Plot Results window)→Contour Plot →Nodal Solution→select:DOF solution, UX,UY,Def + Undeformed , Stress ,SX,SY,SZ, Def + Undeformed →OK:结果分析:(第一小题)图为采用三节点常应变单元,200个单元,左下图为位移变形图;右下图为应力变形图图为采用六节点三角形单元,200个单元,左下图为位移变形图;右下图为应力变形图根据位移和应力图可得:单元类型Min(位移)Max(位移)Min(应力)Max(应力)常应变三节点05461392364六节点三角形0607043对比分析可得:最大位移都发生在A点,即大坝顶端,最大应力发生在B点附近,即坝底和水的交界处,且整体应力和位移变化分布趋势相似,符合实际情况;结果显示三节点和六节点单元分析出来的最大应力值相差较大根据结果显示,最小三节点和六节点单元分析出来的最小应力值相差极为悬殊,结合理论分析,实际上A点不承受载荷,最小应力接近于零,显然六节点三角形单元分析在这一点上更准确。

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

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

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

有限元分析及应用大作业作业要求:1)个人按上机指南步骤至少选择习题中3个习题独立完成,并将计算结果上交;也可根据自己科研工作给出计算实例。

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

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

因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况P=98000-9800*Y;建立几何模型,进行求解;假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3;2:有限元建模过程:2.1 进入ANSYS :程序→ANSYS APDL 15.02.2设置计算类型:ANSYS Main Menu: Preferences →select Structural →OK2.3选择单元类型:ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182,六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window)2.4定义材料参数:ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK2.5生成几何模型:生成特征点:ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标:input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK生成坝体截面:ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点,1(0,0),2(6,0),3(0,10) →OK2.6 网格划分:ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)2.7 模型施加约束:给底边施加x和y方向的约束:ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK给竖直边施加y方向的分布载荷:ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量;2) 在Result窗口中出现{X},写入所施加的载荷函数:98000-9800*{Y};3) File>Save(文件扩展名:func) →返回:Parameters →Functions →Read from file:将需要的.func文件打开,参数名取meng,它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边;OK →在下拉列表框中,选择:Existing table →OK →选择需要的载荷为meng参数名→OK2.8 分析计算:ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current LoadStep window) →OK2.9 结果显示:ANSYS Main Menu: General Postproc→Plot Results →Deformed Shape… →select Def + Undeformed→OK (back to Plot Results window)→Contour Plot →Nodal Solution→select:DOF solution, UX,UY,Def + Undeformed, Stress ,SX,SY,SZ, Def + Undeformed→OK3.1:结果分析:(第一小题)图为采用三节点常应变单元,200个单元,左下图为位移变形图;右下图为应力变形图图为采用六节点三角形单元,200个单元,左下图为位移变形图;右下图为应力变形图根据位移和应力图可得:单元类型Min(位移) Max(位移) Min(应力) Max(应力)常应变三节点0 0.0284 5461 392364六节点三角形0 0.0292 0.001385 607043对比分析可得:最大位移都发生在A点,即大坝顶端,最大应力发生在B点附近,即坝底和水的交界处,且整体应力和位移变化分布趋势相似,符合实际情况;结果显示三节点和六节点单元分析出来的最大应力值相差较大根据结果显示,最小三节点和六节点单元分析出来的最小应力值相差极为悬殊,结合理论分析,实际上A点不承受载荷,最小应力接近于零,显然六节点三角形单元分析在这一点上更准确。

六节点的应力范围较大,所以可判断在单元数目相同的前提下,节点数目越多,分析精度就越大;但是节点数目的增多会带来计算工作量增加和计算效率降低的问题。

3.2结果分析(第二小题):图为采用三节点常应变单元,第一个为96个单元;第二个为3750个单元序号单元数Max(位移) Min(应力) Max(应力)1 96 0.0262mm 13655 2754802 200 0.0284mm 5461 3923643 1750 0.029mm 2183 538858单元数影响分析:随着单元数目的增加,最大位移变化不大,应力变化范围逐步增大;随着单元数目的增加,即网格划分越密,分析的结果准确度将会提高;但是单元数目的增加和节点数目的增加都会造成计算量的增加和计算速度的下降的问题;对于本次计算结果,仍可能存在虚假应力,应力的准确值无法准确得出,只是网格划分越密,计算结果越精确。

所以减少虚假应力影响的措施之一就是增加单元的数目,提高网格划分的密度。

3.3 结果分析(第三小题):方案一方案二方案一的位移图和应力图方案二的位移图和应力图分析:编号Max(位移)Min(应力) Max(应力)方案一0.0107mm 50773 156173方案二0.0128mm 76772 157567由结果可知,方案一,二都只有四个单元,所以在计算应力和位移的时结果的准确度较低。

由应力图可知,方案二得出的最大应力处不符合实际情况,而方案一的最大应力所在位置符合实际情况,所以总体来说,方案一的分析结果优于方案二。

4 总结:由以上结果可知,划分网格时,单元数目也不能太少,太少会导致计算精度较低,单元数目的增加也可以提高计算的精度;但是采用较多节点的单元反而会增加计算的工作量,影响工作效率和经济性,所以选用合适的单元数量至关重要;其次,单元类型的选择也相当重要,直接影响计算的精度。

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

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

(1个和20个单元)解:1 建模:由于薄板的长度为2m,而厚度只有0.3cm,所以长度相对于厚度为无限大,满足平面应变问题模型,而且长度上受均布载荷,上下两边不不受作用力,所以采用平面应变问题进行求解。

假设薄板的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.32 有限元ANSYS建模分析:2.1 进入ANSYS :程序→ANSYS APDL 15.02.2设置计算类型:ANSYS Main Menu: Preferences →select Structural →OK2.3选择单元类型:ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182,四节点矩形单元选择Solid Quad 4node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window)2.4定义材料参数:ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK2.5生成几何模型:生成特征点:ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标:input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK生成坝体截面:ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点,1(0,0),2(6,0),3(0,10) →OK2.6 网格划分:ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)2.7 模型施加约束:给底边施加x和y方向的约束:ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK给竖直边施加y方向的分布载荷:ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量;2) 在Result窗口中出现{X},写入所施加的载荷函数:98000-9800*{Y};3) File>Save(文件扩展名:func) →返回:Parameters →Functions →Read from file:将需要的.func文件打开,参数名取meng,它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边;OK →在下拉列表框中,选择:Existing table →OK →选择需要的载荷为meng参数名→OK2.8 分析计算:ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK2.9 结果显示:ANSYS Main Menu: General Postproc→Plot Results →Deformed Shape… →select Def + Undeformed→OK (back to Plot Results window)→Contour Plot →Nodal Solution→select:DOF solution, UX,UY,Def + Undeformed, Stress ,SX,SY,SZ, Def + Undeformed→OK3.结果分析:3.1 第一小题结果分析:三节点常应变单元(2单元,200单元)2个单元的三节点常应变单元的位移图,应力图200个单元的三节点常应变单元的位移图,应力图编号Max(位移) Min(应力)Max(应力)1个单元0.000783mm 86527 105673200个单元0.000775mm 83039 152695分析:由应力图和位移图可知,随着单元数的增加,ANSYS模拟精度变高,两种不同数量的单元,位移变化相差不大;而应力变化相差较大,200单元的最大应力位置在与固定端接触位置,与实际相符,所以200单元数比2单元数更精确,适当的增加单元的数量可以提高模拟精度。

相关文档
最新文档