当前位置:文档之家› 有限元大作业-黄书文-M201370395

有限元大作业-黄书文-M201370395

有限元大作业-黄书文-M201370395
有限元大作业-黄书文-M201370395

有限元实例分析大作业一

受均匀内压作用的厚壁圆筒: 问题描述: 受均匀内压p=12.5N/mm 2作用的厚壁圆筒。其几何参数为:内径R i =100mm , 外径R e =200mm ,桶壁后h=100mm ,材料参数为:E=8666.67Mpa ,v=0.3, s σ=17.32Mpa ,材料符合Mise 屈服条件。 (a)求理想塑性材料的解,给出应力r σ和θσ沿径向r 的分布曲线,并求完全卸载 后圆筒内的残余应力分布。 (b)求线性强化材料(E 1=0.6E 或E 1=0.6E)的解,即应力r σ和θσ沿径向r 的分布 曲线。 (c)求幂硬化材料的解并绘出当弹塑性比例系数为m=0,1/4,1/2,2/3和m=1.0时, 即应力r σ和θσ沿径向r 的分布曲线。 求解分析: 由于该厚壁筒模型是轴对称模型,所以在求解过程中,我们选取了1/4模型进行了进行建模分析,具体如下图: 建模时取了柱坐标系下厚壁筒从0。~90。范围内的部分,高度取为100mm ,模型完成后进行网格的划分,这里利用了Patran 的Mesh Seed 功能,通过在径向、周向,高度方向撒种生成Mesh 网格,网格划分如上图。 考虑到实体的变形情况,关于模型的边界条件,定义如下: (1)模型的上、下表面为两个平面,在该两平面上限制z 方向的位移为0; (2)对于模型的内外两圆弧面,为了方便定义边界条件,建立了柱坐标,该两平

是延径向变形的,所以ρ坐标是放开的,为了限制模型的刚体移动,这里限制角坐标θ为0。 (3)对于模型两个侧平面,是属于模型的对称面,所以该两平面的单元在垂直于平面的方向上位移为零,这里利用柱坐标,即沿周向的位移为零,所以同样要限制角坐标θ为0。 由于厚壁筒受到均匀内压,所以在施加载荷时选择均布载荷Pressure,大小为p=12.5N/mm2,作用在内圆弧表面上。 对于材料塑性的定义,首先定义样式模量和泊松比,然后在弹塑性对话框里定义屈服载荷和硬化系数或通过在Stress/Strain Curve栏中添加事先定义的材料属性场来表征弹塑性比例系数m。 对于求解分析,求解器选择Nastran进行计算分析,单元属性选择3D Solid 属性,分析类型定义为非线性并设置大变形和跟随力及载荷增量步等,以此来进行弹塑性的非线性求解。 结果分析: (a)对于理性塑性材料,即硬化系数为0,求解结果如下: 该图为100%载荷作用下模型的应力云图及变形情况。观察可知,筒内壁应力较高且首先达到屈服应力发生塑性变形,沿径向方向向外,各层应力逐渐递减,且外层部分属于弹性变形的范畴,模型某一层为弹塑性变形的分界面。

有限元分析大作业报告

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

华科大有限元分析题及大作业题答案——船海专业(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

(完整版)有限元大作业matlab---课程设计例子

有限元大作业程序设计 学校:天津大学 院系:建筑工程与力学学院 专业:01级工程力学 姓名:刘秀 学号:\\\\\\\\\\\ 指导老师:

连续体平面问题的有限元程序分析 [题目]: 如图所示的正方形薄板四周受均匀载荷的作用,该结构在边界 上受正向分布压力, m kN p 1=,同时在沿对角线y 轴上受一对集中压 力,载荷为2KN ,若取板厚1=t ,泊松比0=v 。 [分析过程]: 由于连续平板的对称性,只需要取其在第一象限的四分之一部分参加分析,然后人为作出一些辅助线将平板“分割”成若干部分,再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性,四分之一部分的边界约束,载荷可等效如图所示。

[程序原理及实现]: 用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN,经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下: 问题类型:可用于计算弹性力学平面问题和平面应变问题 单元类型:采用常应变三角形单元 位移模式:用用线性位移模式 载荷类型:节点载荷,非节点载荷应先换算为等效节点载荷 材料性质:弹性体由单一的均匀材料组成 约束方式:为“0”位移固定约束,为保证无刚体位移,弹性体至少应有对三个自由度的独立约束 方程求解:针对半带宽刚度方程的Gauss消元法

输入文件:由手工生成节点信息文件NODE.IN,和单元信息文件ELEMENT.IN 结果文件:输出一般的结果文件DATA.OUT 程序的原理如框图:

弹性力学与有限元理论部分考试题2123

弹性力学与有限元(理论部分)考试题 姓名: (90分钟) 一、填空题(25分) 1.弹性力学的基本任务是:。(1分)2.弹性力学的基本假设是:,, ,,。(2.5分)3.应力分量包括:。(2分)应变分量包括:。(2分)位移分量包括:。(2分)4.平衡微分方程反映的是分量和分量的关系,有个方程。(1.5分)几何方程反映的是分量和分量的关系,有个方程。(1.5分)物理方程反映的是分量和分量的关系,有个方程。(1.5分)5.在对受力体进行有限元分析划分网格时,网格划分较密时,优点是:,缺点是:,反之亦然,因此划分适合的网格密度十分重要。(2分) 6.对于平面问题,三角形单元是最简单、最常用的单元,在平面应力问题中单元形状为,在平面应变问题中单元形状为。(2分) 7.在有限元分析中,有六个常用矩阵:D,B,S,k,K,N,它们分别叫做矩阵,矩阵,矩阵,矩阵,矩阵和函数。(6分) 8.刚度矩阵的半带宽B与有关,B=2(d+1)。(1分) 二、简答题(33分) 1.弹性力学与材料力学有何异同?(5分) 2.什么叫单元的位移模式?分别写出三节点三角形单元,四节点矩形单元,六节点三角形单元,八节点矩形单元的位移模式。(10分)

y 3.平面应力问题和平面应变问题的特点(应力、应变)各是什么?(6分) 4.圣维南原理?并举例说明如何应用之(5分) 5.轴对称问题的特点是什么?(3分) 6.如何理解等参元法,简述用等参元法进行空间问题有限元分析的过程(4分) 三.如下图为一个受力体,其划分的单元和节点编号如图1所示,求出半带宽,写出其整体刚度矩阵。(10分) 图1 受力体单元的划分和编号图 四.在单元e 中,三角形单元三个节点分别为i 、j 、m ,,请把图2和图3的力简化到各节点上。(12分) 图2 单元受力图 图3单元受力图 y

有限元分析报告大作业

有限元分析》大作业基本要求: 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 几何建模

有限元分析大作业试题

有限元分析习题及大作业试题 要求: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)

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

有限元大作业

风电主轴承有限元分析 XXX 摘要:基于有限元法在接触问题中的应用,对风电主轴承进行非线性分析。以轴承外圈的内表面和内圈的外表面为目标面,以滚子为接触面创建接触对分析滚子的接触应力情况。最大应力值出现在滚子边缘出,对最大承载滚子环向接触应力分析表明,有限元分析结果与理论计算结果相近,验证了利用有限元法分析风电主轴承应力状态的可行性。 关键词:风电主轴承;接触应力;有限元分析 0 引言 随着传统能源的日益枯竭以及环境污染问题愈发严重,风能作为一种清洁的的可再生能源近些年受到越来越多的关注。风力发电技术已广泛运用于世界各地。一些发达国家风力发电产业已得到了迅猛发展,技术日趋成熟,并开始走向产业化规模化发展阶段[1-3]。 风电主轴承是风力发电机重要的组成部分。其结构形式图下图1所示。据统计,如今安装的所有风力发电机中,采用主轴轴承支撑原理的占总数的75-80%[4],这种支撑是轴承内圈安装在旋转的主轴上,外圈固定在单独的轴承座上,相对于圆锥滚子轴承或圆柱滚子轴承来说,主轴轴承位置处轴产生变形,需要轴承具有一定的调心作用,所以都采用了调心滚子轴承。近年来由于计算机技术的飞速发展,轴承的受力分析计算已经普遍采用有限元分析的方法,能够准确合理地解决轴承复杂的非线性接触问题,为轴承的分析和计算提供了一种新的方法,成为未来的一个发展方向。在机械设备的设计过程中,对受力较大且复杂的零件进行受力分析,校核其整体和局部强度并进行合理的布局设计,是为了防止因应力过大而导致在实际工作中损坏或寿命降低[5]。本文主要运用ANSYS Workbench有限元软件对风电主轴承进行静力学计算,分析轴承内部结构参数对轴承载荷分布和最大接触应力的影响规律。 图1 风电主轴承结构及安装图 1 有限元分析过程 1.1 风电轴承有限元分析基本步骤 不同的物理性质和数学模型的问题,有限元法求解的基本步骤是相同的,只不过 具体公式推导和运算求解不尽相同。有限元分析求解问题的基本计算步骤[6]: 1.问题及求解域定义; 2.求解域离散化; 3.确定状态变量及控制方法; 4.单元推导;

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。

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

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

课程研究报告 科目:有限元分析技术教师:阎春平姓名:色学号: 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

弹性力学与有限元法分析及实例讲解

弹性力学与有限元法分析 弹性力学是固体力学的一个重要分支,是研究弹性固体在受外力作用、温度改变、边界约束或其他外界因素作用下而发生的应力、形变和位移状态的科学。有限单元法是力学、数学、物理学、计算方法、计算机技术等多种学科综合发展和结合的产物,是随着计算机技术的广泛应用而迅速发展起来的一种数值分析方法。有限元法的基本思想就是化整为零,分散分析,再集零为整。即用结构力学方法求解弹性力学问题,实质是将复杂的连续体划分为有限多个简单的单元体,单元体之间仅仅通过结点相连,实现化无限自由度问题为有限稀有度问题,将连续场函数的(偏)微分方程的求解问题转化为有限个参数的代数方程组的求解问题。 有限元方法经过近半个世纪的发展,目前已经成为各种工程问题特别是结构分析问题的标准分析方法,而有限元软件也已成为现代结构设计中不可缺少的工具。有限元软件是有限元理论通向实际工程应用的桥梁,它的应用极大地提高了力学学科解决自然科学和工程实际问题的能力,进一步促进了有限元方法的发展。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件,广泛用于机械制造、石油化工、航空航天、汽车交通、土木工程、造船、水利等一般工业及科学研究。 ANSYS软件的组成: (一)前处理模块 该模块为用户提供了一个强大的实体建模及网格划分工具,可以方便的构造有限元模型,软件提高了100种以上的单元类型,用来模拟工程中的各种结构和材料。包括: 1.实体建模:参数化建模,布尔运算及体素库,拖拉、旋转、拷贝、蒙皮、倒角等。 2.自动网格划分,自动进行单元形态、求解精度检查及修正。 3.在集合模型上加载:点加载、分布载荷、体载荷、函数载荷。 4.可扩展的标准梁截面形状库。 (二)分析计算模块 该模块包括结构分析(可进行线性分析、非线性分析和高度非线性分析)、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析,可模拟多种物理介质的相互作用,具有灵敏度分析及优化分析能力。 (三)后处理模块 将计算结果以彩色等值线、梯度、矢量、粒子流、立体切片、透明及半透明等图形方式显示出来,也可以用图表、曲线形式显示或输出。 由于现在只是对ANSYS工程软件有初步的了解和掌握,所以本次作业仅以(1)结构静力学分析为例,运用ANSYS软件对汽车连杆进行受力分析;(2)

有限元大作业matlab课程设计例子

有 限 元 大 作 业 程 序 设 计 学校:天津大学 院系:建筑工程与力学学院 专业:01级工程力学 姓名:刘秀 学号:\\\\\\\\\\\ 指导老师: 连续体平面问题的有限元程序分析 [题目]: 如图所示的正方形薄板四周受均匀载荷的作用,该结构在边界 上受正向分布压力, m kN p 1=,同时在沿对角线y 轴上受一对集中压 力,载荷为2KN ,若取板厚1=t ,泊松比0=v 。 [分析过程]: 由于连续平板的对称性, 只需要取其在第一象限的四分之一部分

参加分析,然后人为作出一些辅助线将平板“分割”成若干部分,再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性,四分之一部分的边界约束,载荷可等效如图所示。 [ 用和单元信息文件DATA.OUT。 位移模式:用用线性位移模式 载荷类型:节点载荷,非节点载荷应先换算为等效节点载荷 材料性质:弹性体由单一的均匀材料组成 约束方式:为“0”位移固定约束,为保证无刚体位移,弹性体至少应有对三个自由度的独立约束 方程求解:针对半带宽刚度方程的Gauss消元法 输入文件:由手工生成节点信息文件NODE.IN,和单元信息文件ELEMENT.IN 结果文件:输出一般的结果文件DATA.OUT 程序的原理如框图:

(1) ID : ID=2时为平面应变问题 (平面问题) ,LJK_ELE(I,1),LJK_ELE(I,2), X(I),Y(I)分别存放节点I 的x ,y 表示第I 个作用有节点载荷的节点x,y 方向的节点载荷数值 存放节点载荷向量,解方程后该矩 (2 READ_IN : 读入数据 BAND_K : 形成半带宽的整体刚度矩阵 FORM_KE : 计算单元刚度矩阵 FORM_P : 计算节点载荷 CAL_AREA :计算单元面积 DO_BC : 处理边界条件 CLA_DD : 计算单元弹性矩阵 SOLVE : 计算节点位移 CLA_BB : 计算单元位移……应变关系矩阵 CAL_STS :计算单元和节点应力 (3)文件管理: 源程序文件: chengxu.for 程序需读入的数据文件:

有限元分析及应用大作业

有限元分析及应用大作业 作业要求: 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

试题及其答案--弹性力学与有限元分析(DOC)

如下图所示三角形薄板,按三结点三角形单元划分后,对于与局部编码ijm 对应的整体编码,以下叙述正确的是( D )。 ① I 单元的整体编码为162 ② II 单元的整体编码为426 ③ II 单元的整体编码为246 ④ III 单元的整体编码为243 ⑤ IV 单元的整体编码为564 A. ①③ B. ②④ C. ①④ D. ③⑤ 一、填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、 形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相 适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规 定相适应。 4、物体受外力以后,其内部将发生内力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa ,50=y σMPa ,5010=xy τ MPa ,则主应力 =1σ150MPa ,=2σ0MPa ,=1α6135' 。 8、已知一点处的应力分量, 200=x σMPa ,0=y σMPa ,400-=xy τ MPa ,则主应力=1σ512 MPa ,=2σ-312 MPa ,=1α-37°57′。 9、已知一点处的应力分量,2000-=x σMPa ,1000=y σMPa ,400-=xy τ MPa ,则主应力 =1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三 套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、 应力边界条件和混合边界条件。 13、按应力求解平面问题时常采用逆解法和半逆解法。 14、有限单元法首先将连续体变换成为离散化结构,然后再用结构力学位移法进行求解。 其具体步骤分为单元分析和整体分析两部分。 15、每个单元的位移一般总是包含着两部分:一部分是由本单元的形变引起的,另一部

有限元分析报告大作业

基于ANSYS软件的有限元分析报告 机制1205班杜星宇U201210671 一、概述 本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变。然后与实际情况进行比较,证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据,并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想,对有限元分析方法的实际应用有一个大致的认识。 二、问题分析 已知:桌子几何尺寸如图所示,单位为mm。假设桌子的四只脚同地面完全固定,桌子上存放物品,物品产生的均匀分布压力作用在桌面,压力大小等于300Pa,其中弹性模量E=9.3GPa,泊松比μ=0.35,密度ρ=560kg/m3,分析桌子的变形和应力。

将桌脚固定在地面,然后在桌面施加均匀分布的压力,可以看作对进行平面应力分析,桌脚类似于梁单元。由于所分析的结构比较规整且为实体,所以可以将单元类型设为八节点六面体单元。 操作步骤如下: 1、定义工作文件名和工作标题 (1)定义工作文件名:执行Utility Menu/ File/Change Jobname,在弹出Change Jobname 对话框修改文件名为Table。选择New log and error files复选框。 (2)定义工作标题:Utility Menu/File/ Change Title,将弹出Change Title对话框修改工作标题名为The analysis of table。 (3)点击:Plot/Replot。 2、设置计算类型 (1)点击:Main Menu/Preferences,选择Structural,点击OK。

有限元法大作业

有限元法大作业 一平面刚架的程序 用Visual C++编制的平面刚架的源程序如下: ///////////////////////////////////////////////////////程序开始////////////////////////////////////////////////////////////////// #include"iostream.h" #include"math.h" #include"stdlib.h" #include"conio.h" //***************** //声明必要变量 //***************** #define PI 3.141592654 int NE; //单元数 int NJ; //节点数 int NZ; //支承数 int NPJ; //有节点载荷作用的节点数 int NPF; //非节点载荷数 int HZ; //载荷码 int E; //单元码 int fangchengshu; double F[303]; //各节点等效总载荷数值 int dym_jdm[100][2]; //单元码对应的节点码:dym_jdm[][0], dym_jdm[][1]分为前后节点总码 int zhichengweizhi[300]; //记录支持节点作用点的数组 int fjzhzuoyongdanyuan[100]; //非节点载荷作用单元 int fjzhleixing[100]; //非节点载荷类型:1-均布,2-垂直集中,3-平行集中,4-力偶,5-角度集中 double fjzhzhi[100]; //非节点载荷的值 double fjzhzuoyongdian[100]; //非节点载荷在各竿的作用点 double fjzhjiaodu[100]; //非节点载荷作用角度 int jdzhzuoyongdian[100]; //节点载荷作用的节点数组 double jiedianzaihe[101][3];//节点载荷值,其jiedianzaihe[][0]-- jiedianzaihe[][2]分别为U, V, M double zhengtigangdu[303][303]; //整体刚度数组 double changdu[100]; //各单元竿长数组 double jiaodu[100]; //各单元角度数组 double tanxingmoliang[100]; //各单元弹性模量数组 double J_moliang[100]; //各单元J模量数组 double mianji[100]; //各单元面积数组 double weiyi[303]; //记录各个节点位移的数组 double dy_weiyi[100][6]; //各个单元在局部坐标系中的位移数组dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的u1,v1,@1,u2,v2,@2 double dy_neili[100][6];//各个单元在局部坐标系中的固端内力dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的U1,V1,M1,U2,V2,M2 double gan_neili[100][6];//各个单元的竿端内力数组,gan_neili[i][6]表示第i+1单元的6内力. //*******************

弹性力学与有限元分析试题答案

最新弹性力学与有限元分析复习题及其答案 一、 填空题 1、弹性力学研究弹性体由于受外力作用、边界约束或温度改变等原因而发生的应力、 形变和位移。 2、在弹性力学中规定,线应变以伸长时为正,缩短时为负,与正应力的正负号规定相 适应。 3、在弹性力学中规定,切应变以直角变小时为正,变大时为负,与切应力的正负号规 定相适应。 4、物体受外力以后,其部将发生力,它的集度称为应力。与物体的形变和材料强度直接有关的,是应力在其作用截面的法线方向和切线方向的分量,也就是正应力和切 应力。应力及其分量的量纲是L -1MT -2。 5、弹性力学的基本假定为连续性、完全弹性、均匀性、各向同性。 6、平面问题分为平面应力问题和平面应变问题。 7、已知一点处的应力分量100=x σMPa ,50=y σMPa ,5010=xy τ MPa ,则主应力 =1σ150MPa ,=2σ0MPa ,=1α6135' 。 8、已知一点处的应力分量, 200=x σMPa ,0=y σMPa ,400-=xy τ MPa ,则主应力 =1σ512 MPa ,=2σ-312 MPa ,=1α-37°57′。 9、已知一点处的应力分量,2000-=x σMPa ,1000=y σMPa ,400-=xy τ MPa ,则主应 力=1σ1052 MPa ,=2σ-2052 MPa ,=1α-82°32′。 10、在弹性力学里分析问题,要考虑静力学、几何学和物理学三方面条件,分别建立三 套方程。 11、表示应力分量与体力分量之间关系的方程为平衡微分方程。 12、边界条件表示边界上位移与约束,或应力与面力之间的关系式。分为位移边界条件、 应力边界条件和混合边界条件。

相关主题
相关文档 最新文档