有限元分析题及大作业题答案
有限单元法考试题及答案
有限单元法考试题及答案一、单项选择题（每题2分，共10分）1. 有限元法中，单元刚度矩阵的计算是基于（）。
A. 位移法B. 势能原理C. 能量守恒定律D. 牛顿第二定律答案：B2. 在有限元分析中，以下哪项不是网格划分时需要考虑的因素？（）A. 网格数量B. 网格形状C. 材料属性D. 边界条件答案：C3. 有限元分析中，以下哪项不是结构分析的基本步骤？（）A. 离散化B. 求解C. 后处理D. 优化设计答案：D4. 在有限元分析中，以下哪种类型的单元不适用于平面应力问题？（）A. 三角形单元B. 四边形单元C. 六面体单元D. 楔形单元答案：C5. 有限元分析中，以下哪种边界条件不属于几何边界条件？（）A. 固定支座B. 压力C. 温度D. 位移答案：C二、多项选择题（每题3分，共15分）6. 有限元法中，以下哪些因素会影响单元的精度？（）A. 单元形状B. 单元数量C. 材料属性D. 网格划分答案：ABD7. 在有限元分析中，以下哪些是常见的数值积分方法？（）A. 一阶积分B. 二阶积分C. 高斯积分D. 牛顿-莱布尼茨积分答案：ABC8. 有限元分析中，以下哪些是常见的单元类型？（）A. 线性单元B. 二次单元C. 三次单元D. 非线性单元答案：ABCD9. 在有限元分析中，以下哪些是常见的后处理技术？（）A. 应力云图B. 位移云图C. 模态分析D. 热分析答案：ABC10. 有限元分析中，以下哪些是常见的非线性问题？（）A. 几何非线性B. 材料非线性C. 接触非线性D. 热应力问题答案：ABCD三、填空题（每题2分，共20分）11. 有限元法中，单元刚度矩阵的计算通常基于___________原理。
答案：势能12. 在有限元分析中，网格划分的目的是将连续的___________离散化为有限数量的单元。
答案：域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个单元。
有限元大作业
250250试题 5：图示为带方孔（边长为 80mm ）的悬臂梁，其上受部分均布载荷（P=10KN/m ）作用，试采用一种平面单元，对图示两种结构进行有限元分析，并就方孔的布置（即方位）进行分析比较，如将方孔设计为圆孔，结果有何变化？（板厚为 1mm ，材料为钢）。
3001KN9003001KN图6-1一、几何建模与分析由图6-1及问题描述可知，板的长宽尺寸远远大于厚度，研究结构为一很薄的等厚度薄板，满足平面应力的几何条件；作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面上无外力作用，满足平面应力的载荷条件。
故该问题属于平面应力问题，薄板所受的载荷为面载荷，分布情况及方向如图6-1所示，建立几何模型，进行求解。
薄板的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3。
二、有限元分析及其计算结果选取PLANE182作为分析的单元，来分析薄板的位移和应力，由于此问题是平面应力问题，并在K3选择str w/thk ，设置THK 为1。
1）方孔竖直制，划分方式采用自由方式，划分后网格的模型如图6-2所示。
计算得到的位移和应力分布如图6-3所示。
图6-2 方孔竖直的网格划分图6-3 位移及应力分布云图2）方孔正直制，划分方式采用自由方式，划分后网格的模型如图6-4所示。
计算得到的位移和应力分布如图6-5所示。
图6-4 方孔正直的网格划分图6-5 位移及应力分布云图3）圆孔按图6-1所示模型进行建模。
并用PLANE182单元进行划分网格，网格大小采用全局网格控制，划分方式采用自由方式，划分后网格的模型如图6-6所示。
计算得到的位移和应力分布如图6-7所示。
图6-4 方孔正直的网格划分图6-5 位移及应力分布云图根据以上的模型分析的位移和应力图，可以得出方孔竖直、方孔正直、圆孔的最大最小位移应力的分布如表6-1所示。
三、比较与分析1）方孔竖直与方孔正直的比较，发现方孔正直的位移变形较小，应力相差不大2）圆孔与方孔比较，发现圆孔的位移变性最小，应力也最小，故可以得出圆孔的布置结构对整体布置的效果最好。
有限元实例分析
作业一：有限元分析实例实例：请对一个盘轴配合机构进行接触分析。
轴为一等直径空心轴，盘为等厚度圆盘，其结构及尺寸如图所示。
盘和轴为一种材料，材料参数为：弹性模量Ex=2.5E5，泊松比NUXY=0.3，摩擦系数MU=0.25，试采用有限元计算方法分析轴和盘在过盈配合时的应力应变分布以及将轴从盘心拔出时轴和盘的接触情况。
问题分析说明（1）本题主要分析装配过程中结构的静态响应，所以分析步选择通用静态分析步。
由于为过盈配合，属于大变形，故应考虑几何非线性的影响。
（2）模型具有轴对称性，所以可以采取轴对称模型来进行分析，先建立二维模型计算，再转换为三维模型计算，这样可以节省计算时间。
分析过程由两个载荷步组成, 第一个载荷步为过盈分析, 求解过盈安装时的情况。
第二个载荷步为将轴从盘心拔出时的接触分析, 分析在这个过程中盘心面和轴的外表面之间的接触应力。
它们都属于大变形问题, 属于非线性问题。
在分析时需要定义一些非线性选项来帮助问题的收敛。
（3）接触面之间有很大的相对滑动，所以模型要使用有限滑移。
模型建立的分析说明（1）进定义单元类型此项实例分析的问题中涉及到大变形, 故选用So li d185 单元类型来建立本实例入部件模块，的模型。
盘轴接触问题属于面面接触, 目标面和接触面都是柔性的,将使用接触单元T ARGET 170 和CO NTAT17 4来模拟接触面。
分别创建名为为part1、part2的部件。
（2）定义材料属性，在线性各向同性材料属性对话框中的EX (弹性模量) 文本框中输入 2 . 5E5,PRX Y (泊松比) 文本框中输入0 . 3,并将定义的材料属性赋予给part1和part2。
如下图所示。
（3）进入装配模块，创建两者间的装配关系。
（4）进入分析步模块定义名为step1和step2的两个分析步。
（5）进入相互作用模块，创建相互作用属性，设置摩擦系数；然后定义接触关系。
如下图所示。
（6）进入载荷模块，创建边界条件，依次定义名为BC -2（类型为：完全固定）、BC -3（类型为：位移/转角，约束U1、UR3），分析步均为Initial 。
有限元分析大作业
有限元大作业一题目要求：图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号节点的支反力如下图所示：。
有限元分析与应用大作业
有限元分析及应用大作业课程名称: 有限元分析及应用班级:姓名:试题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。
图2-2 数学模型二、第一问三节点常应变单元（2个和200个单元）三节点单元类型为PLANE42，设置好单元类型后，实常数设置板厚为0.3M 。
采用2个单元的网格划分后的结果如图2-3,200个单元的网格划分图如图2-6所示。
约束的施加方式和载荷分布如图2-2中所示。
约束右边线上节点全部自由度。
计算得到的位移云图分别如图2-4、7所示，应力云图如图2-5、8所示。
P=33.33MPa图2-3 2个三角形单元的网格划分图图2-5 2个三角形单元的应力云图图2-7 200个三角形单元的位移云图三、第二问四节点矩形单元的计算四节点单元类型为PLANE42，设置好单元类型后，实常数设置板厚为0.3M。
采用1个单元的网格划分后的结果如图2-9,50个单元的网格划分图如图2-12所示。
约束的施加方式和载荷分布如图2-2中所示。
约束右边线上节点全部自由度。
计算得到的位移云图分别如图2-10、11所示，应力云图如图2-13、14所示。
图2-9 1个四边形单元的网格划分图图2-11 1个四边形单元的应力云图图2-12 50个四边形单元的网格划分图图2-13 50个四边形单元的位移云图图2-14 50个四边形单元的应力云图四、第三问八节点等参单元的计算四节点单元类型为PLANE82，设置好单元类型后，实常数设置板厚为0.3M。
有限元习题及答案
有限元习题及答案有限元习题及答案有限元方法是一种常用的数值计算方法，用于求解各种工程和科学问题。
在学习有限元方法的过程中，练习习题是非常重要的，可以帮助学生巩固所学的知识，并提高解决实际问题的能力。
本文将介绍一些有限元习题及其答案，希望对学习有限元方法的同学有所帮助。
习题一：一维热传导问题考虑一个长度为L的一维杆，其两端固定，杆上的温度满足以下热传导方程：∂²T/∂x² = 0，其中T为温度，x为位置。
已知杆的两端温度分别为T1和T2，求解杆上的温度分布。
解答一：根据热传导方程，可以得到温度分布的一般解为T(x) = Ax + B，其中A和B为常数。
根据边界条件，可以得到方程组：T(0) = B = T1T(L) = AL + B = T2解方程组可得A = (T2 - T1) / L，B = T1。
因此，温度分布为T(x) = ((T2 - T1) / L) * x + T1。
习题二：二维弹性问题考虑一个矩形薄板，其长为L，宽为W，材料的弹性模量为E，泊松比为ν。
已知薄板的边界上施加了一定的边界条件，求解薄板上的位移场。
解答二：对于二维弹性问题，可以使用平面应力假设，即假设薄板内部的应力只有两个分量σx和σy，并且与z轴无关。
根据平面应力假设和胡克定律，可以得到位移场的偏微分方程：∂²u/∂x² + ν * (∂²u/∂y²) + (1 - ν) * (∂²v/∂x∂y) = 0∂²v/∂y² + ν * (∂²v/∂x²) + (1 - ν) * (∂²u/∂x∂y) = 0其中u和v分别为位移场在x和y方向上的分量。
边界条件根据具体情况给定。
通过数值方法，如有限元方法，可以求解位移场的近似解。
习题三：三维流体力学问题考虑一个三维流体力学问题，流体在一个封闭容器内流动，容器的形状为一个长方体，已知流体的速度场和压力场的初始条件，求解流体的运动状态。
电机组件的有限元分析考核试卷
B.力边界
C.位移边界
D.材料属性
7.在电磁场分析中，哪些物理量通常作为输入参数？()
A.磁势
B.电势
C.电流密度
D.材料电阻
8.哪些方法可以用来提高电机组件有限元分析的收敛性？()
A.优化网格划分
B.调整边界条件
C.选择合适的求解器
D.增加计算机硬件配置
9.以下哪些是电机组件动态分析的主要内容？()
A.检查网格质量
B.检查边界条件设置
C.对比理论值和计算结果
D.优化计算模型
19.以下哪种方法适用于电机组件的疲劳分析？()
A.静态分析
B.动态分析
C.稳态热分析
D.疲劳分析
20.在电机组件的有限元分析中，以下哪个参数代表材料的磁性能？()
A.磁导率
B.电导率
C.比热容
D.电阻率
二、多选题（本题共20小题，每小题1.5分，共30分，在每小题给出的四个选项中，至少有一项是符合题目要求的）
3.热传递机制包括传导、对流和辐射。优化包括改善材料布局、散热设计和边界条件设置。
4.耦合效应处理包括顺序耦合和同时耦合。解决方法有直接耦合法和间接耦合法。
1. ×
2. ×
3. √
4. ×
5. √
6. √
7. ×
8. √
9. ×
10. ×
五、主观题（参考）
1.基本流程包括建模、网格划分、材料属性定义、边界条件设定、加载求解和结果分析。挑战包括模型简化、网格质量、收敛性和计算资源需求。
2.选择求解器需考虑问题类型、求解精度和计算效率。策略包括预处理、迭代方法和后处理优化。
5.在有限元分析中，所有的边界条件都必须在加载步骤中定义。（）
汽车结构有限元分析试题及答案(精华)
一、20分)(×) 1. 节点的位置依赖于形态，而并不依赖于载荷的位置( √ ) 2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元(×) 3. 不能把梁单元、壳单元和实体单元混合在一起作成模型( √ ) 4. 四边形的平面单元尽可能作成接近正方形形状的单元(×) 5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化处理的话会得到一样的答案(×) 6. 用有限元法不可以对运动的物体的结构进行静力分析( √ ) 7. 一般应力变化大的地方单元尺寸要划的小才好(×) 8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度( √ ) 9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小( √ ) 10 一维变带宽存储通常比二维等带宽存储更节省存储量。
二、填空(20 分)1．平面应力问题与薄板弯曲问题的弹性体几何形状都是薄板，但前者受力特点是：平行于板面且沿厚度均布载荷作用，变形发生在板面内；后者受力特点是：垂直于板面的力的作用，板将变成有弯有扭的曲面。
2 ．平面应力问题与平面应变问题都具有三个独立的应力分量：σx，σy，τxy ，三个独立的应变分量：εx，εy，γxy，但对应的弹性体几何形状前者为薄板，后者为长柱体。
3．位移模式需反映刚体位移，反映常变形，满足单元边界上位移连续。
4 ．单元刚度矩阵的特点有：对称性，奇异性，还可按节点分块。
5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为二维问题处理。
6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。
等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。
7．有限单元法首先求出的解是节点位移，单元应力可由它求得，其计算公式为} = [D][B]6}e 。
有限元作业试题及答案.doc
2
答：一般选用三角形或四边形单元，在满足一定精度情况,
有限元划分网格的基本原则是:
1、拓朴正确性原则。即单元间是靠单元顶点、或单元边、或单元面连接
2、几何保形原则。即网格划分后，单元的集合为原结构近似
3、特性一致原则。即材料相同，厚度相同
4、单元形状优良原则。单元边、角相差尽可能小
c j二elcm= —a
Ni = l/a2 • a x = x/a
同理可得：Nj二y/a
有限元方法及应用试题
1
答：单元离散（划分、剖分）一单元分析一整体分析
有限元分析的主要步骤主要有：
A结构的离散化
B单元分析。选择位移函数、根据几何方程建立应变与位移的关系、根据物理方程建立应力
与位移的关系、根据虚功原理建立节点力与节点位移的关系（单元刚度方程）
C等效节点载荷计算
D整体分析，建立整体刚度方程
7、图示三角形ijni为等边三角形单元，边长为1,单位面积材料密度位P,集 中力F垂直作用于nij边的中点，集度为q的均布载荷垂直作用于im边。写出三 角形单元的节点载荷向量。
q:移到m, i点F:移到m, j点重力：移到m, I, j点
要证{8}=0
只需证，Nm = 0
Nm= 1/2A (am+bmx +cmy)
(d)平面三角形单元，29个节点，38个自由度
4、什么是等参数单元？。
如果坐标变换和位移插值采用相同的节点，并且单元的形状变换函数与位移插值的形函
数一样，则称这种变换为等参变换，这样的单元称为等参单元。
5பைடு நூலகம்
v(x, y)=
答：不能取这样的位移模式，因为在平面三节点三角形单元中，位移模式应该是呈线性的。
有限元受力分析大作业
１. 定义分析类型
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）
四. 加载和求解
１. 定义分析类型
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分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算;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个单元。
有限元分析与应用大作业
有限元分析及应用大作业课程名称: 有限元分析及应用班级:姓名:试题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。
中南大学有限元习题与答案(Word最新版)
中南大学有限元习题与答案通过整理的中南大学有限元习题与答案相关文档，希望对大家有所帮助，谢谢观看！中南大学有限元习题与答案习题 2.1 解释如下的概念：应力、应变，几何方程、物理方程、虚位移原理。
解应力是某截面上的应力在该处的集度。
应变是指单元体在某一个方向上有一个ΔU的伸长量，其相对变化量就是应变。
表示在x轴的方向上的正应变，其包括正应变和剪应变。
几何方程是表示弹性体内节点的应变分量与位移分量之间的关系，其完整表示如下：物理方程：表示应力和应变关系的方程某一点应力分量与应变分量之间的关系如下：虚位移原理：在弹性有一虚位移情况下，由于作用在每个质点上的力系，在相应的虚位移上虚功总和为零，即为：若弹性体在已知的面力和体力的作用下处于平衡状态，那么使弹性体产生虚位移，所有作用在弹性体上的体力在虚位移上所做的工就等于弹性体所具有的虚位能。
2.2说明弹性体力学中的几个基本假设。
连续性假设：就是假定整个物体的体积都被组成该物体的介质所填满，不存在任何间隙。
完全弹性假设：就是假定物体服从虎克定律。
各向同性假设：就是假定整个物体是由同意材料组成的。
小变形和小位移假设：就是指物体各点的位移都远远小于物体原来的尺寸，并且其应变和转角都小于1。
2.3简述线应变与剪应变的几何含义。
线应变：应变和刚体转动与位移导数的关系，剪应变表示单元体棱边之间夹角的变化。
2.4 推到平面应变平衡微分方程。
解：对于单元体而言其平衡方程：在平面中有代入上式的2.5 如题图2.1所示，被三个表面隔离出来平面应力状态中的一点，求和的值。
解：x方向上：联立二式得：2.6相对于xyz坐标系，一点的应力如下某表面的外法线方向余弦值为，，求该表面的法相和切向应力。
解：该平面的正应力全应力该平面的切应力2.7一点的应力如下MP 求主应力和每一个主应力方向的方向余弦；球该店的最大剪应力。
解：设主平面方向余弦为，由题知将代入得即，。
最大剪应力（1）当时代入式（2.21）（2）当时代入式（2.21）且2.8已知一点P的位移场为，求该点p(1,0,2)的应变分量。
有限元分析题及大作业题答案
有限元分析及应用作业报告试题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所示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
有限元分析及应用作业报告
试题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 →OK
3、生成几何模型：
a.创建关键点点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入16个点的坐标→OK
b、将这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 →拾取模型左部的竖直边→OK
B、给模型施加载荷
ANSYS Main Menu: Solution →Define Loads →Apply →Structural →force→on keypoints→
拾取上面左端关键点→700N/mm→ok
C、分析计算：ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK
6)结果显示：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所示。
3、结果分析
由以上结果分析可知：
1）扳手柄根部的转折点处的应力比其它部位大，不利于结构承受载荷；扳手卡紧螺母部位的转折角应力也较大，同时扳手头部的应力也较大。
2）扳手柄的应变最大，扳手头部的应变也较大。
3）正确性分析评判。
由以上应变和应力分析结果可以看出，如果要是使最大应力最小，一般选用圆弧过渡，折线次之，方孔方角最大。
在选择圆弧时，尽量使圆弧半径大些，这样可以最大限度的减小应力。
同时，圆弧对零件内部的应力分布改善效果也较为明显。
对于应变的处理，应选用刚度强的材料，以减少应变。
五、多方案计算比较
1）不同的网格划分对最大应力和最大位移的位置及大小影响较小，但是对零件内部的
应力分布影响较大。
2）结点规模增减对精度也有一定的的影响。
3）单元改变对精度的影响较大。
六、建议与体会
通过以上分析情况可以看出，如果要使分析结果较为精确，首先是单元的类型选择要恰当，不同的单元类型会造成结果的不同，节点较多可以保证计算精度较高；由计算结果可知，划分网格时，单元数目也不能太少，单元数目的增加也可以提高计算的精度；但是对于实际工程而言，采用较多节点的单元反而会影响计算的工作量，这是不经济不必要的。
因此在保证网格划分大小适当和均匀的前提下，使应力集中处划的密集些，这样也能得到较为精确的结果。