有限元分析报告

有限元分析报告
有限元分析报告

有限元法在工程领域的发展现状和应用

有限元法(Finite Element Method,FEM),是计算力学中的一种重要的方法,它是20世纪50年代末60年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。有限元法最初应用在工程科学技术中,用于模拟并且解决工程力学、热学、电磁学等物理问题。对于过去用解析方法无法求解的问题和边界条件及结构形状都不规则的复杂问题,有限元法则是一种有效的分析方法。

近年来随着计算机技术的普及和计算速度的不断提高,有限元分析在工程设计和分析中得到了越来越广泛的重视,已经成为解决复杂的工程分析计算问题的有效途径,现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算,其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器,国防军工,船舶,铁道,石化,能源,科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃,主要表现在以下几个方面:

(1)增加产品和工程的可靠性

(2)在产品的设计阶段发现潜在的问题

(3)经过分析计算,采用优化设计方案,降低原材料成本

(4)模拟试验方案,减少试验次数,从而减少试验经费

一、有限元法的基本思想

有限元法的基本思想是先将研究对象的连续求解区域离散为一组有限个且

按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合,且单元本身又可以有不同形状,因此可以模拟成不同几何形状的求解小区域;

然后对单元(小区域)进行力学分析,最后再整体分析。这种化整为零,集零为整的方法就是有限元的基本思路。

有限元法分析计算的思路和做法可归纳如下:

1物体离散化

将某个工程结构离散为由各种单元组成的计算模型,这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来;单元节点的设置、性质、数目等应视问题的性质,描述变形形态的需要和计算进度而定(一般情况单元划分越细则描述变形情况越精确,即越接近实际变形,但计算量越大)。所以有限元中分析的结构已不是原有的物体或结构物,而是同新材料的由众多单元以一定方式连接成的离散物体。这样,用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理,则所获得的结果就与实际情况相符合。

2单元特性分析

(1)选择位移模式

在有限单元法中,选择节点位移作为基本未知量时称为位移法;选择节点力作为基本未知量时称为力法;取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化,所以,在有限单元法中位移法应用范围最广。

当采用位移法时,物体或结构物离散化之后,就可把单元中的一些物理量如

位移,应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常,有限元法我们就将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数。

(2)分析单元的力学性质

根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等,找出单元节点力和节点位移的关系式,这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式,从而导出单元刚度矩阵,这是有限元法的基本步骤之一。

(3)计算等效节点力

物体离散化后,假定力是通过节点从一个单元传递到另一个单元。但是,对于实际的连续体,力是从单元的公共边传递到另一个单元中去的。因而,这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去,也就是用等效的节点力来代替所有作用在单元上的力。

3单元组集

利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来,形成整体的有限元方程1-1:

kq

f

式中k—整体结构的刚度矩阵

q—节点位移列阵

f—载荷列阵

4求解未知节点位移

求解有限元方程式(1-1)得出位移。这里,可以根据方程组的具体特点来选择合适的计算方法。

通过上述分析,可以看出,有限单元法的基本思想是“一分一合”,分是为了进行单元分析,合则是为了对整体结构进行综合分析。

二、有限元法的发展现状

有限元法是R.Courant于1943年首先提出的。自从提出有限元概念以来,有限元理论及其应用得到了迅速发展。过去不能解决或能解决但求解精度不高的问题,都得到了新的解决方案。传统的FEM假设:分析域是无限的;材料是同质的,甚至在大部分的分析中认为材料是各向同性的;对边界条件简化处理。但实际问题往往是分析域有限、材料各向异性或边界条件难以确定等。为解决这类问题,美国学者提出用GFEM (Gener-alized Finite Element Method)解决分析域内含有大量孔洞特征的问题;比利时学者提出用HSM(the Hybrid metis Singular element of Membraneplate)解决实际开裂问题。

FEM理论研究成果为其应用奠定了基础,计算机技术的发展为其提供了条

件。20世纪70年代以来,相继出现了一些通用的有限元分析(FEA: Finite Element Analysis)系统,如SAP、ASKA、NASTRAN等,这些FEA系统可进行航空航天领域的结构强度、刚度分析,从而推动了FEM在工程中的实际应用。

20世纪80年代以来,随着工程工作站的出现和广泛应用,原来运行于大中型机上的FEA系统得以在其上运行,同时也出现了一批通用的FEA系统,如ANSYS-PC、NISA,SUPERSAP等。20世纪90年代以来,随着微机性能的显著提高,大批FEA系统纷纷向微机移植,出现了基于Windows的微机版FEA系统。

在FEM应用领域不断扩展、求解精度不断提高的同时,FEM也从分析比较向优化设计方向发展。印度Mahanty博士用ANSYS对拖拉机前桥进行优化设计,结果不但降低了约40%的前桥自重,还避免了在制造过程中的大量焊接工艺,降低了生产成本。

FEM在国内的应用也十分广泛。自从我国成功开发了国内第一个通用有限元程序系统JIGFEX后,有限元法渗透到工程分析的各个领域中,从大型的三峡工程到微米级器件都采用FEM进行分析,在我国经济发展中拥有广阔的发展前景。

目前在进行大型复杂工程结构中的物理场分析时,为了估计并控制误差,常用基于后验误差估计的自适应有限元法。基于后处理法计算误差,与传统算法不同,将网格自适应过程分成均匀化和变密度化2个迭代过程。在均匀化迭代过程中,采用均匀网格尺寸对整体区域进行网格划分,以便得到一个合适的起始均匀网格;在变密度化迭代过程中只进行网格的细化操作,并充分利用上一次迭代的结果,在单元所在的曲边三角形区域内部进行局部网格细化,保证了全局网格尺寸分布的合理性,使得不同尺寸的网格能光滑衔接,从而提高网格质量。整个方案简单易行,稳定可靠,数次迭代即可快速收敛,生成的网格布局合理,质量高。

经过半个多世纪的发展,FEM已从弹性力学平面问题扩展到空间问题、板壳问题;从静力问题扩展到动力问题、稳定问题和波动问题;从线性问题扩展到非线性问题;从固体力学领域扩展到流体力学、传热学、电磁学等其他连续介质领域;从单一物理场计算扩展到多物理场的耦合计算。它经历了从低级到高级、从简单到复杂的发展过程,目前已成为工程计算最有效的方法之一。

三、有限元应用的有关问题

1 FEM的应用过程

FEM应用于实际问题须经历以下过程,如图1所示。

(1)问题的数学描述。

对问题客观规律的数学描述(通常是微分方程及边界条件)是建立有限元方程的前提。单元特性矩阵和整体有限元方程都是基于数学模型建立的。常见的弹性力学基本方程、运动方程、热传导方程等都是对客观现象的数学描述。

(2)有限元方程的建立。

利用变分原理,通过离散、单元分析、整体分析等过程,建立数学模型的有限元方程,它通常是一组易于用数值方法求解的代数方程。

(3)算法研究。

有限元方程的计算量庞大,须有有效的算法来保证计算效率和精度,同时考虑

对计算条件的要求。如求解大型线性方程组的带宽法、波前法,求解大型特征值问题的分块Lanczos 法等。

(4)程序开发。

数值计算依赖于计算机,因此求解算法需用相应的计算程序来实现。

(5)有限元建模。

对应于FEA 系统的前处理(Pre-pro-cessing)。它为数值计算提供所有原始输入数据(节点数据、单元数据和边界条件数据)。因为模型形式直接决定计算精度和规模,且建模所需时间约占整个FEA 的70%左右,所以建模质量和效率是FEA 的关键。图2列出了有限元建模中的关键技术。

图1 FEM 的应用过程

图2 有限元建模的关键技术

(6)数值计算。

理 论

研 究 系 统

开 发

对应于FEA系统的计算(Solving)。它由一系列计算程序组成,计算程序又称求解器(solver)。每个求解器完成特定类型的计算。因此求解器越多,系统功能越强。

(7)结果处理。

对应于FEA系统的后处理(Post-pro-cessing)。它对计算结果进行处理、显示、运算和列表等。

若按照(1)~(7)过程,问题得以解决,则FEM应用结束;反之,则需根据求解结果提出改进方案,循环执行(5)~(7)过程,直至问题解决或得到最佳设计。

对于一个全新的问题,必须从第一步开始。而对已知的问题,可从第(5)步开始,即直接利用已有的FEA系统,建立有限元模型。在实际应用中,绝大多数问题都属于第二类问题。

2 FEM的应用领域

有限元法最初应用在求解结构的平面问题上,发展至今,已由二维问题扩展到三维问题、板壳问题,由静力学问题扩展到动力学问题、稳定性问题,由结构力学扩展到流体力学、电磁学、传热学等学科,由线性问题扩展到非线性问题,由弹性材料扩展到弹塑性、塑性、粘弹性、粘塑性和复合材料,从航空技术领域扩展到航天、土木建筑、机械制造、水利工程、造船、电子技术及原子能等,由单一物理场的求解扩展到多物理场的耦合,其应用的深度和广度都得到了极大的拓展。

2.1 静力分析

包括线性非线性静力分析。线性静力分析研究线弹性结构的变形和应力,它是工程结构分析和设计中最基本的方法。非线性结构静力分析主要研究外载作用下引起的非线性响应,其中非线性来源主要是材料非线性、几何非线性和边界条件非线性3大类。

2.2 动力分析

主要包括以下分析类型:

(1)模态分析。用于求解多自由度系统的模态参数。图3为计算得到的计算机主板的前三阶振型。

(2)瞬态响应分析。求解在时域内结构承受随时间变化的载荷和速度作用时的动力响应。

(3)简谐响应分析。对简谐激励结构在其平衡位置的振动进行分析。

(4)频谱响应分析和随机振动分析。用于分析结构受已知频率激励时的最大响应。

(5)屈曲和失稳分析。分析考察结构的极限承载能力,研究结构总体或局部的稳定性,获得结构失稳形态和失稳路径。

(6)自动接触分析。用于接触边界定义和摩擦分析。

2.3 失效和破坏分析

包括断裂分析(线弹性断裂分析和弹塑性断裂分析)、裂纹萌生与扩展分析、跌落分析和疲劳失效分析。图4是对电视机进行的跌落分析。

2.4 热传导分析

包括稳态热传导分析、瞬态热传导分析、热辐射、强迫对流及温度的耦合分析。图5是一个铸造过程中的热传导分析,目的是追踪固化过程中铸件和模具的温度分布。

2.5 电磁场分析

它用于对电磁场中电感、电容、磁通量密度、涡流、电场分布、磁力线分布、能量损失等物理量进行分析。图6为E型电机的磁场分布和吸力特性图。

2.6 声场分析

它用来研究在含有流体介质中声波的传播问题,或分析浸在流体中的固体结构的动态特性。

2.7 流体分析

研究流体速度、压强、密度变化规律和粘滞流体的运动规律及粘滞流体中运动物体所受阻力及其它热力学性质。图7是离心泵叶轮叶片表面相对速度和压力变化曲线。

2.8 耦合场分析

考虑两种或两种以上物理场的交叉作用和相互影响(耦合)。图8是双压电晶片梁在结构场和电场共同作用下的变形。

四、有限元的发展趋势和研究热点

4.1 有限元的发展趋势

(1)与CAD软件的无缝集合

当今有限元分析软件的一个发展趋势是与通用CAD软件的集成使用,即在用CAD软件完成部件和零件的造型设计后,能直接将模型传送到CAE软件中进行有限元网格划分并进行分析计算,如果分析的结果不满足设计要求则重新进行设计和分析,直到满意为止,从而极大地提高了设计水平和效率。为了满足工程师快捷地解决复杂工程问题的要求,许多商业化有限元分析软件都开发了和著名的CAD软件(例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley和AutoCAD等)的接口。有些CAE软件为了实现和CAD软件的无缝集成而采用了CAD的建模技术,如ADINA软件由于采用了基于Parasolid 内核的实体建模技术,能和以Parasolid为核心的CAD软件(如Unigraphics、SolidEdge、SolidWorks)实现真正无缝的双向数据交换。

(2)更强大的网格处理能力

有限元法求解问题的基本过程主要包括:分析对象的离散化、有限元求解、计算结果的后处理三部分。由于结构离散后的网格质量直接影响到求解时间及求解结果的正确性与否,近年来各软件开发商都加大了其在网格处理方面的投入,使网格生成的质量和效率都有了很大的提高,但在有些方面却一直没有得到改进,如对三维实体模型进行自动六面体网格划分和根据求解结果对模型进行自适应网格划分,除了个别商业软件做得较好外,大多数分析软件仍然没有此功能。自动六面体网格划分是指对三维实体模型程序能自动的划分出六面体网格单元,现在大多数软件都能采用映射、拖拉、扫略等功能生成六面体单元,但这些功能都只能对简单规则模型适用,对于复杂的三维模型则只能采用自动四面体网格划分技术生成四面体单元。对于四面体单元,如果不使用中间节点,在很多问题中将会产生不正确的结果,如果使用中间节点将会引起求解时间、收敛速度等方面的一系列问题,因此人们迫切的希望自动六面体网格功能的出现。自适应性网格划分是指在现有网格基础上,根据有限元计算结果估计计算误差、重新划分网格和再计算的一个循环过程。对于许多工程实际问题,在整个求解过程中,模型的某些区域将会产生很大的应变,引起单元畸变,从而导致求解不能进行下去或求解结果不正确,因此必须进行网格自动重划分。自适应网格往往是许多工程问题如裂纹扩展、薄板成形等大应变分析的必要条件。

(3)由求解线性问题发展到求解非线性问题

随着科学技术的发展,线性理论已经远远不能满足设计的要求,许多工程问题如材料的破坏与失效、裂纹扩展等仅靠线性理论根本不能解决,必须进行非线性分析求解,例如薄板成形就要求同时考虑结构的大位移、大应变(几何非线性)和塑性(材料非线性);而对塑料、橡胶、陶瓷、混凝土及岩土等材料进行分析或需考虑材料的塑性、蠕变效应时则必须考虑材料非线性。众所周知,非线性问题的求解是很复杂的,它不仅涉及到很多专门的数学问题,还必须掌握一定的理论知识和求解技巧,学习起来也较为困难。为此国外一些公司花费了大量的人力和物

力开发非线性求解分析软件,如ADINA、ABAQUS等。它们的共同特点是具有高效的非线性求解器、丰富而实用的非线性材料库,ADINA还同时具有隐式和显式两种时间积分方法。

(4)由单一结构场求解发展到耦合场问题的求解

有限元分析方法最早应用于航空航天领域,主要用来求解线性结构问题,实践证明这是一种非常有效的数值分析方法。而且从理论上也已经证明,只要用于离散求解对象的单元足够小,所得的解就可足够逼近于精确值。现在用于求解结构线性问题的有限元方法和软件已经比较成熟,发展方向是结构非线性、流体动力学和耦合场问题的求解。例如由于摩擦接触而产生的热问题,金属成形时由于塑性功而产生的热问题,需要结构场和温度场的有限元分析结果交叉迭代求解,即"热力耦合"的问题。当流体在弯管中流动时,流体压力会使弯管产生变形,而管的变形又反过来影响到流体的流动……这就需要对结构场和流场的有限元分析结果交叉迭代求解,即所谓"流固耦合"的问题。由于有限元的应用越来越深入,人们关注的问题越来越复杂,耦合场的求解必定成为CAE软件的发展方向。

(5)程序面向用户的开放性

随着商业化的提高,各软件开发商为了扩大自己的市场份额,满足用户的需求,在软件的功能、易用性等方面花费了大量的投资,但由于用户的要求千差万别,不管他们怎样努力也不可能满足所有用户的要求,因此必须给用户一个开放的环境,允许用户根据自己的实际情况对软件进行扩充,包括用户自定义单元特性、用户自定义材料本构(结构本构、热本构、流体本构)、用户自定义流场边界条件、用户自定义结构断裂判据和裂纹扩展规律等等。

4.2 有限元的研究热点

(1)超收敛应力计算

在超收敛应力计算方面,陈传淼教授是中国学派的代表,他的“有限元超收敛构造理论”成果属于国际领先水平。超收敛应力(及位移)计算是以FEM为代表的各种数值方法争相追求的目标,所以它一直是FEM研究的难点与热点问题。

(2)有限元模型修正技术

有限元模型修正技术(或试验/分析模型相关)是要充分利用结构试验和FEA

两者的优点,用少量的结构试验数据对有限元模型进行修正,获得比较准确的有限元模型。

参考文献

[1]张洪信.有限元基础理论与ANSYS应用.北京:机械工业出版社.2006.2.

[2]高秀华.张小江.王欢.有限单元法原理及应用简明教程.北京:化学工业出版社.2008.6

[3]阎照文. ANSYS10.0工程电磁分析技术与实例详解. 北京:中国水利水电出版社.2010.8

[4]张朝晖. ANSYS12.0热分析工程应用实战手册. 北京:中国铁道出版

社.2010.8

[5]张美超.钟世镇. 国内生物力学中有限元的应用研究进展. 解剖科学进展杂志2003

[6]陈锡栋.杨婕.赵晓栋.范细秋. 有限元法的发展现状及应用. 中国制造业信息化杂志2010.6

[7]于亚婷.杜平安.王振伟. 有限元法的应用现状研究. 机械设计学报2005.3

[8]Daniel S Pipkinsay,Satya N Atlurib. Applications of the three di-mensional method finite element alternating method[J].Finite El-ements in Analysis and Design,1996 (23):133-153.

[9]Theofanis Strouboulis, Lin Zhang,Ivo Babu ska. Generalized fi-nite element method using mesh-based handbooks: applicationto problems in domains with many https://www.360docs.net/doc/ed10477963.html,put[J]. MethodsAppl. Mech. Engrg, 2003 (192) :3 109-3 161.

[10]Birgersson F,Finnveden S. Application of the Spectral Finite Ele-ment Method to Turbulent Boundary Layer Induced Vibretion ofPlates[J].Journal of Sound and Vibration ,2003,259(4):873-891.

[11]韩西,钟厉,等.有限元分析在结构分析和计算机仿真中的应用.[J].重庆交通学院学报,2001,(20):124-126.

有限元分析实验报告

武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析 开课学院机电工程学院 指导老师姓名 学生姓名 学生专业班级机电研 1502班 2015—2016 学年第2学期

实验一方形截面悬臂梁的弯曲的应力与变形分析 钢制方形悬臂梁左端固联在墙壁,另一端悬空。工作时对梁右端施加垂直向下的30KN的载荷与60kN的载荷,分析两种集中力作用下该悬臂梁的应力与应变,其中梁的尺寸为10mmX10mmX100mm的方形梁。 1.1方形截面悬臂梁模型建立 建模环境:DesignModeler 15.0。 定义计算类型:选择为结构分析。 定义材料属性:弹性模量为2.1Gpa,泊松比为0.3。 建立悬臂式连接环模型。 (1)绘制方形截面草图:在DesignModeler中定义XY平面为视图平面,并正视改平面,点击sketching下的矩形图标,在视图中绘制10mmX10mm的矩形。(2)拉伸:沿着Z方向将上一步得到的矩阵拉伸100mm,即可得到梁的三维模型,建模完毕,模型如下图1.1所示。 图1.1 方形截面梁模型 1.2 定义单元类型: 选用6面体20节点186号结构单元。 网格划分:通过选定边界和整体结构,在边界单元划分数量不变的情况下,通过分别改变节点数和载荷大小,对同一结构进行分析,划分网格如下图1.2所示:

图1.2 网格划分 1.21 定义边界条件并求解 本次实验中,讲梁的左端固定,将载荷施加在右端,施以垂直向下的集中力,集中力的大小为30kN观察变形情况,再将力改为50kN,观察变形情况,给出应力应变云图,并分析。 (1)给左端施加固定约束; (2)给悬臂梁右端施加垂直向下的集中力; 1.22定义边界条件如图1.3所示: 图1.3 定义边界条件 1.23 应力分布如下图1.4所示: 定义完边界条件之后进行求解。

有限元分析大作业报告

有限元分析大作业报告 试题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) 三节点常应变单元 三节点常应变单元的位移分布图 三节点常应变单元的应力分布图

有限元分析实验报告

学生学号1049721501301实验课成绩 武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析机电工程学院开课学院 指导老师姓名

学生姓名 学生专业班级机电研1502班 学年第学期2016—20152 实验一方形截面悬臂梁的弯曲的应力与变形分析 钢制方形悬臂梁左端固联在墙壁,另一端悬空。工作时对梁右端施加垂直 向下的30KN的载荷与60kN的载荷,分析两种集中力作用下该悬臂梁的应力与应变,其中梁的尺寸为10mmX10mmX100mm的方形梁。 方形截面悬臂梁模型建立1.1 建模环境:DesignModeler15.0。 定义计算类型:选择为结构分析。 定义材料属性:弹性模量为 2.1Gpa,泊松比为0.3。 建立悬臂式连接环模型。 (1)绘制方形截面草图:在DesignModeler中定义XY平面为视图平面,并正 视改平面,点击sketching下的矩形图标,在视图中绘制10mmX10mm的矩形。 (2)拉伸:沿着Z方向将上一步得到的矩阵拉伸100mm,即可得到梁的三维模型,建模完毕,模型如下图 1.1所示。

图1.1方形截面梁模型 :定义单元类型1.2 选用6面体20节点186号结构单元。 网格划分:通过选定边界和整体结构,在边界单元划分数量不变的情况下,通过分别改变节点数和载荷大小,对同一结构进行分析,划分网格如下图 1.2

所示: 图1.2网格划分 1.21定义边界条件并求解 本次实验中,讲梁的左端固定,将载荷施加在右端,施以垂直向下的集中 力,集中力的大小为30kN观察变形情况,再将力改为50kN,观察变形情况,给出应力应变云图,并分析。 (1)给左端施加固定约束; (2)给悬臂梁右端施加垂直向下的集中力; 1.22定义边界条件如图1.3所示:

华科大有限元分析题及大作业题答案——船海专业(DOC)

姓名:学号:班级:

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

二、几何建模与分析 图1-2力学模型 由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况及方向如图1-2所示,建立几何模型,进行求解。 假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3 三、第1问的有限元建模 本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。 1)设置计算类型:两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural 2)选择单元类型:三节点常应变单元选择的类型是PLANE42(Quad 4node42),该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类型是PLANE183(Quad 8node183),该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题,故对Element behavior(K3)设置为plane strain。 3)定义材料参数 4)生成几何模 a. 生成特征点 b.生成坝体截面 5)网格化分:划分网格时,拾取所有线段设定input NDIV 为10,选择网格划分方式为Tri+Mapped,最后得到200个单元。 6)模型施加约束: 约束采用的是对底面BC全约束。 大坝所受载荷形式为Pressure,作用在AB面上,分析时施加在L AB上,方向水平向右,载荷大小沿L AB由小到大均匀分布(见图1-2)。以B为坐标原点,BA方向为纵轴y,则沿着y方向的受力大小可表示为: ρ(1) = gh P- =ρ g = - 10 {* } 98000 98000 (Y ) y

有限元分析报告

《有限元基础理论》报告 学院: 班级: 姓名: 学号: 任课老师: 二〇一一年十二月

题目一:三维托架实体受力分析 题目:1、三维托架实体受力分析:托架顶面承受50psi的均匀分布载荷。托架通过有孔的表面固定在墙上,托架是钢制的,弹性模量E=29×106psi,泊松比v=0.3.试通过ANSYS输出其变形图及其托架的von Mises应力分布。 题目1的分析:先进行建模,此建模的难点在对V3的构建(既图中的红色部分)。要想构建V3,首先应将A15做出来,然后执行Main Menu>Preprocessor>Modeling>Operate>Booleans>Add>V olumes命令,将所有的实体合并为一个整体。建模后,就对模型进行网格的划分,实行Main Menu>Preprocessor>Meshing>MeshTool,先对网格尺寸进行编辑,选0.1,然后点Meshing,Pick all进行网格划分,所得结果如图1.1。划分网格后,就可以对模型施加约束并进行加载求解了。施加约束时要注意,由于三维托架只是通过两个孔进行固定,故施加约束应该只是针对两孔的内表面,执行Main Menu>Solution>Define Loads>Apply>Structrual>Displacement>Symmetry B.C>On Areas命令,然后拾取两孔的内表面,单击OK就行了。施加约束后,就可以对实体进行加载求解了,载荷是施加在三维托架的最顶上的表面的,加载后求解运算,托架的变形图如图1.2。

图1.1、托架网格图 图1.2输出的是原型托架和施加载荷后托架变形图的对比,虚线部分即为托架的原型,从图1.2可看出,由于载荷的作用,托架上面板明显变形了,变形最严重的就是红色部分,这是因为其离托板就远,没有任何物体与其分担载荷,故其较容易变形甚至折断。这是我们 在应用托架的时候应当注意的。

有限元上机实验报告

有限元上机实验报告结构数值分析与程序设计 上机实验 院系: 土木工程与力学学院专业: 土木工程 班级: 姓名: 学号: 指导教师: 1、调试教材P26-30程序FEM1。 1.1、输入数据文件为: 6,4,12,6,1.0E0,0.0,1.0,0.0,1 3,1,2 5,2,4 3,2,5 6,3,5 0.0,2.0 0.0,1.0 1.0,1.0 0.0,0.0 1.0,0.0 2.0,0.0

1,3,7,8,10,12 1.2、输出数据文件为: NN NE ND NFIX E ANU T GM NTYPE 6 4 12 60.1000E+01 0.000 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 -0.100000E+01 2 0.000000E+00 0.000000E+00 3 0.000000E+00 0.000000E+00 4 0.000000E+00 0.000000E+00 5 0.000000E+00 0.000000E+00 6 0.000000E+00 0.000000E+00 NODE X-DISP Y-DISP 1 -0.879121E-15 -0.325275E+01 2 0.879121E-16 -0.125275E+01 3 -0.879121E-01 -0.373626E+00 4 0.117216E-1 5 -0.835165E-15 5 0.175824E+00 -0.293040E-15 6 0.175824E+00 0.263736E-15 ELEMENT X-STR Y-STR XY-STR 1 -0.879121E-01 -0.200000E+01 0.439560E+00 2 0.175824E+00 -0.125275E+01 0.256410E-15 3 -0.879121E-01 -0.373626E+00 0.307692E+00 4 0.000000E+00 -0.373626E+00 -0.131868E+00 2、修改FEM1,计算P31例2-2。

有限元分析大作业试题

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

习题及上机指南:(试题见上机指南) 例题1 坝体的有限元建模与受力分析 例题2 平板的有限元建模与变形分析 例题1:平板的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: plane 0.5 m 0.5 m 0.5 m 0.5 m 板承受均布载荷:1.0e 5 P a 图1-1 受均布载荷作用的平板计算分析模型 1.1 进入ANSYS 程序 →ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: plane →Run 1.2设置计算类型 ANSYS Main Menu : Preferences →select Structural → OK 1.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 → OK 1.5定义实常数 ANSYS Main Menu: Preprocessor →Real Constant s… →Add … →select T ype 1→ OK →input THK:1 →OK →Close (the Real Constants Window)

有限元分析报告样本

《有限元分析》报告基本要求: 1. 以个人为单位完成有限元分析计算,并将计算结果上交;(不允许出现相同的分析模型,如相 同两人均为不及格) 2. 以个人为单位撰写计算分析报告; 3. 按下列模板格式完成分析报告; 4. 计算结果要求提交电子版,报告要求提交电子版和纸质版。(以上文字在报告中可删除) 《有限元分析》报告 一、问题描述 (要求:应结合图对问题进行详细描述,同时应清楚阐述所研究问题的受力状况和约束情况。图应清楚、明晰,且有必要的尺寸数据。) 一个平面刚架右端固定,在左端施加一个y 方向的-3000N 的力P1,中间施加一个Y 方向的-1000N 的力P2,试以静力来分析,求解各接点的位移。已知组成刚架的各梁除梁长外,其余的几何特性相同。 横截面积:A=0.0072 m2 横截高度:H=0.42m 惯性矩:I=0.0021028m4x 弹性模量: E=2.06x10n/ m2/ 泊松比: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”如图

机械零件有限元分析——实验报告

中南林业科技大学机械零件有限元分析 实验报告 专业:机械设计制造及其自动化 年级: 2013级 班级:机械一班 姓名:杨政 学号:20131461 I

一、实验目的 通过实验了解和掌握机械零件有限元分析的基本步骤;掌握在ANSYS 系统环境下,有限元模型的几何建模、单元属性的设置、有限元网格的划分、约束与载荷的施加、问题的求解、后处理及各种察看分析结果的方法。体会有限元分析方法的强大功能及其在机械设计领域中的作用。 二、实验内容 实验内容分为两个部分:一个是受内压作用的球体的有限元建模与分析,可从中学习如何处理轴对称问题的有限元求解;第二个是轴承座的实体建模、网格划分、加载、求解及后处理的综合练习,可以较全面地锻炼利用有限元分析软件对机械零件进行分析的能力。

实验一、受内压作用的球体的有限元建模与分析 对一承受均匀内压的空心球体进行线性静力学分析,球体承受的内压为 1.0×108Pa ,空 心球体的内径为 0.3m ,外径为 0.5m ,空心球体材料的属性:弹性模量 2.1×1011,泊松比 0.3。 承受内压:1.0×108 Pa 受均匀内压的球体计算分析模型(截面图) 1、进入 ANSYS →change the working directory into yours →input jobname: Sphere 2、选择单元类型 ANSYS Main Menu : Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element Types window)→ Options… →select K3: Axisymmetric →OK →Close (the Element Type window) 3、定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3→ OK 4、生成几何模型生成特征点 ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标:input :1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3)→OK 生成球体截面 ANSYS 命令菜单栏: Work Plane>Change Active CS to>Global Spherical ANSYS Main Menu: Preprocessor →Modeling →Create →Lines →In ActiveCoord → 依次连接 1,2,3,4 点生成 4 条线→OK Preprocessor →Modeling →Create →Areas →Arbitrary →By Lines →依次拾取四条线→OK ANSYS 命令菜单栏: Work Plane>Change Active CS to>Global Cartesian 5、网格划分 ANSYS Main Menu : Preprocessor →Meshing →Mesh Tool →(Size Controls) lines: Set

ansys有限元分析大作业

ansys有限元分析大作业

有限元大作业 设计题目: 单车的设计及ansys有限元分析 专业班级: 姓名: 学号: 指导老师: 完成日期: 2016.11.23

单车的设计及ansys模拟分析 一、单车实体设计与建模 1、总体设计 单车的总体设计三维图如下,采用pro-e进行实体建模。 在建模时修改proe默认单位为国际主单位(米千克秒 mks) Proe》文件》属性》修改

2、车架 车架是构成单车的基体,联接着单车的其余各个部件并承受骑者的体重及单车在行驶时经受各种震动和冲击力量,因此除了强度以外还应有足够的刚度,这是为了在各种行驶条件下,使固定在车架上的各机构的相对位置应保持不变,充分发挥各部位的功能。车架分为前部和后部,前部为转向部分,后部为驱动部分,由于受力较大,所有要对后半部分进行加固。

二、单车有限元模型 1、材料的选择 单车的车身选用铝合金(6061-T6)T6标志表示经过热处理、时效。 其属性如下: 弹性模量:) .6+ 90E (2 N/m 10 泊松比:0.33 质量密度:) 3 2.70E+ N/m (2 抗剪模量:) 60E .2+ N/m (2 10 屈服强度:) .2+ (2 75E 8 N/m 2、单车模型的简化 为了方便单车的模拟分析,提高电脑的运算

效率,可对单车进行初步的简化;单车受到的力的主要由车架承受,因此必须保证车架能够有足够的强度、刚度,抗振的能力,故分析的时候主要对车架进行分析。简化后的车架如下图所示。 3、单元体的选择 单车车架为实体故定义车架的单元类型为实体单元(solid)。查资料可以知道3D实体常用结构实体单元有下表。 单元名称说明 Solid45 三维结构实体单元,单元由8个节点定义,具有塑性、蠕变、应力刚化、 大变形、大应变功能,其高阶单元是 solid95

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

船海1004 黄山 U201012278 有限元分析大作业报告 试题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个单元。

1 船海1004 黄山 U201012278 (6)模型施加约束:约束采用的是对底面BC全约束。大坝所受载荷形式为Pressure,作用在AB面上,分析时施加在L上,方向水平向右,载荷大小沿L 由小到大均匀分布。以ABAB B为坐标原点,BA方向为纵轴y,则沿着y方向的受力大小可表示为: P?gh?gyY}*{?)??98000?9800(10? 2、计算结果及结果分析 (1)三节点常应变单元 三节点常应变单元的位移分布图

有限元分析报告

有限元仿真分析实验 一、实验目的 通过刚性球与薄板的碰撞仿真实验,学习有限元方法的基本思想与建模仿真的实现过程,并以此实践相关有限元软件的使用方法。本实验使用HyperMesh 软件进行建模、网格划分和建立约束及载荷条件,然后使用LS-DYNA软件进行求解计算和结果后处理,计算出钢球与金属板相撞时的运动和受力情况,并对结果进行可视化。 二、实验软件 HyperMesh、LS-DYNA 三、实验基本原理 本实验模拟刚性球撞击薄板的运动和受力情况。仿真分析主要可分为数据前处理、求解计算和结果后处理三个过程。前处理阶段任务包括:建立分析结构的几何模型,划分网格、建立计算模型,确定并施加边界条件。 四、实验步骤 1、按照点-线-面的顺序创建球和板的几何模型 (1)建立球的模型:在坐标(0,0,0)建立临时节点,以临时节点为圆心,画半径为5mm的球体。 (2)建立板的模型:在tool-translate面板下node选择临时节点,选择Y-axis,magnitude输入5.5,然后点击translate+,return;再在2D-planes-square 面板上选择Y-axis,B选择上一步移下来的那个节点,surface only ,size=30。 2、画网格

(1)画球的网格:以球模型为当前part,在2D-atuomesh面板下,surfs 选择前面建好的球面,element size设为0.5mm,mesh type选择quads,选择elems to current comp,first order,interactive。 (2)画板的网格:做法和设置同上。 3、对球和板赋材料和截面属性 (1)给球赋材料属性:在materials面板内选择20号刚体,设置Rho为2.000e-08,E为200000,NU为0.30。 (2)给球赋截面属性:属性选择SectShll,thickness设置为0.1,QR设为0。 (3)给板赋材料属性:材料选择MATL1,其他参数:Rho为2.000e-08,E 为100000,Nu为0.30,选择Do Not Export。 (4)给板赋截面属性:截面选择SectShll,thickness设为0.2。其他参数:SHRE为8.333-01,QR为0,T1为0.2。 (5)给板设置沙漏控制:在Properties-Create面板下Card image选择HourGlass,IHQ为4,QM为0.100。更新平板。 4、加载边界条件 (1)将板上最外面的四行节点分别建成4个set。 (2)建立一个load collector。 (3)Analysis-constraints面板中,设置SIZE为1,nodes通过by sets 选择set_1、set_2、set_3、set_4,然后点击creat即可,边界条件加载完毕。 5、建立载荷条件(给球一个3mm的位移) (1)建立一个plot: post-xy plots-plots-creat plot,然后点击return;

有限元实验报告模板

有限元实验报告 T1013-5 20100130508 蔡孟迪

ANSYS有限元上机报告(一) 班级:T1013-5 学号:20100130508 姓名:蔡孟迪 上机题目: 图示折板上端固定,右侧受力F=1000N,该力均匀分布在边缘各节点上;板厚t=2mm 材料选用低碳钢,弹性模量E=210Gpa,μ=0.33. 一、有限元分析的目的: 1.利用ANSYS构造实体模型; 2.根据结构的特点及所受载荷的情况,确定所用单元类型;正确剖分网格并施加外界条件;3.绘制结构的应力和变形图,给出最大应力和变形的位置及大小;并确定折板角点A处的应力和位移; 4.研究网格密度对A处角点应力的影响; 5.若在A处可用过渡圆角,研究A处圆角半径对A处角点应力的影响。 二、有限元模型的特点: 1.结构类型 本结构属于平面应力类型 2.单位制选择 本作业选择N(牛),mm(毫米),MPa(兆帕)。 3.建模方法 采用自左向右的实体建模方法。 4.定义单元属性及类型 1)材料属性:弹性模量:EX=2.10E5MPa, 泊松比:PRXY=0.33 2)单元类型:在Preferences选Structural,Preprocessor>ElemmentType>Add/Edit/Delete中定义单元类型为:Quad4 node 182,K3设置为:平面薄板问题(Plane strs w/thk) 3)实常数:薄板的厚度THK=2mm 5.划分网格 在MeshTool下选set,然后设置SIZE Element edge length的值,再用Mesh进行网格划分。6.加载和约束过程:在薄板的最上端施加X、Y方向的固定铰链,在薄板的最右端施加1000N 的均匀布置的载荷。

重庆大学研究生有限元大作业教学内容

重庆大学研究生有限 元大作业

课程研究报告 科目:有限元分析技术教师:阎春平姓名:色学号: 2 专业:机械工程类别:学术 上课时间: 2015 年 11 月至 2016 年 1 月 考生成绩: 阅卷评语: 阅卷教师 (签名)

有限元分析技术作业 姓名: 色序号: 是学号: 2 一、题目描述及要求 钢结构的主梁为高160宽100厚14的方钢管,次梁为直径60厚10的圆钢管(单位为毫米),材料均为碳素结构钢Q235;该结构固定支撑点位于左右两端主梁和最中间。主梁和次梁之间是固接。试对在垂直于玻璃平面方向的2kPa 的面载荷(包括玻璃自重、钢结构自重、活载荷(人员与演出器械载荷)、风载荷等)作用下的舞台进行有限元分析。 二、题目分析 根据序号为069,换算得钢结构框架为11列13行。由于每个格子的大小为1×1(单位米),因此框架的外边框应为11000×13000(单位毫米)。 三、具体操作及分析求解 1、准备工作 执行Utility Menu:File → Clear&start new 清除当前数据库并开始新的分析,更改文件名和文件标题,如图1.1。选择GUI filter,执行 Main Menu: Preferences → Structural → OK,如图1.2所示

图1.1清除当前数据库并开始新的分析 图1.2 设置GUI filter 2、选择单元类型。 执行Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add→ select→ BEAM188,如图2.1。之后点击OK(回到Element Types window) →Close

有限元分析报告大作业

有限元分析》大作业基本要求: 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.3 3.4几何建模由于扳手结构比较简单,所以可以直接在ANSYS软件上直接建模,在ANSYS建 立正六 边形,再创立直线,面沿线挤出体,得到扳手几何模型 图4 几何建模

ANSYS实体建模有限元分析-课程设计报告

南京理工大学 课程设计说明书(论文) 作者:学号: 学院(系):理学院 专业:工程力学 题目:ANSYS实体建模有限元分析 指导者: (姓名) (专业技术职务) 评阅者: (姓名) (专业技术职务) 20 年月日

练习题一 要求: 照图利用ANSYS软件建立实体模型和有限元离散模型,说明所用单元种类、单元总数和节点数。 操作步骤: 拟采用自底向上建模方式建模。 1.定义工作文件名和工作标题 1)选择Utility Menu>File>Change Jobname命令,出现Change Jobname对话框,在[/FILNAM ] Enter new jobname文本框中输入工作文件名learning1,单击OK按钮关闭该对话框。 2)选择Utility Menu>File>Change Title命令,出现Change Title对话框,在[/TITLE] Enter new title文本框中输入08dp,单击OK按钮关闭该对话框。 2.定义单元类型 1)选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令,出现Element Types对话框,单击Add按钮,出现 Library of Element Types 对话框。在Library of Element Types 列表框中选择 Structural Solid, Tet 10node 92,在Element type reference number文本框中输入1,单击OK按钮关闭该对话框。 2)单击Element Types对话框上的Close按钮,关闭该对话框。 3.创建几何模型 1)选择Utility Menu>P1otCtrls>Style>Colors>Reverse Video命令,设置显示颜色。 2)选择Utility Menu>P1otCtrls>View Settings>Viewing Direction命令,出现Viewing Direction对话框,在XV,YV,ZV Coords of view point文本框中分别输入1, 1, 1,其余选项采用默认设置,单击OK按钮关闭该对话框。 3)建立支座底块 选择Main Menu>Preprocessor> Modeling>Create>volumes>Block>By Demensios 命令,出现Create Block by Demensios对话框,在X1,X2 X-coor dinates文本框

有限元分析及应用大作业

有限元分析及应用大作业 作业要求: 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.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.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 →OK 2.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) →OK 2.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参数名→OK 2.8 分析计算: ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load

厚壁圆筒有限元分析报告

有限元与CAE分析报告 专业: 班级: 学号: 姓名: 指导教师: 实习时间: 年月日

平面问题的厚壁圆筒问题 一、问题提出 如图所示为一厚壁圆筒,其内半径为r1=50mm,外半径为r2=100mm,作用在内孔上的压力p=10 Mpa,无轴向压力,轴向长度很长可视为无穷,要求对其进行结构静力分析,并计算厚壁圆筒径向应力和切向应力沿半径r方向的分布。弹性模量E=200 Gpa,泊松比μ=0.3。 图1 厚壁圆筒 二、建模步骤 1 定义工作文件名 依次单击Utility Menu>File>Change Jobname,在文本框中输入:1245523229,在“New Log and error files”处选中“yes”,单击“OK”。 2 定义工作标题 依次单击Utility Menu>File>Change Title ,在文本框中输入:1245523229,单击“OK”。依次单击Plot>Replot, 3 定义单元类型 1)依次单击Main Menu>Prefrences,选中“Structural”,单击“OK”。

2)依次单击Main Menu>Preprocessor>Element type>Add/Edit/Delete,出现对话框,单击“Add”,出现一个“Library of Element Type”对话框,。在“Library of Element Type”左面的列表栏中选择“Structural Solid”,在右面的列表栏中选择“Quard 4node 182”,单击“OK”。 2) 单击对话框中的“Options”,在弹出的单元属性对话框中,选择K3关键字element behavior为“Plane strain”,再单击“Close”,完成单元的设置。

相关文档
最新文档