中科大有限元作业答案(第一次到第五次)

清华大学有限元大作业

??????? ????? ???グ?? ? ??? ????φ??? ????φ??? ????φ??? ??? 1 Corresponding author: zzq_890709@https://www.360docs.net/doc/c09699717.html,

?? Contents 0 ????????????? (2) 0.1 ?????? (2) 0.2 ????????? (2) 0.3 ?????? (3) 1????????仈 (4) 1.1 ???? (4) 1.2 ?????????????? (5) 1.3 ???? (7) 1.4 ???????ウ (13) 2 ???????????仈 (14) 2.1 ???? (14) 2.2 ???? (15) 2.3 ???? (16) 3 й???????????仈 (19) 3.1 ??????? (19) 3.2 ?????????????? (19) 3.3 ???? (20) 3.4 ???????ウ (22) 4 ??????????仈 (23) 4.1 ???? (23) 4.2 ???? (24) 4.3 ???? (25) 5 ??????????????仈 (32) 5.1 ???? (32) 5.2 ???? (33) 6 ???仈 (39) 6.1 ???? (39) 6.2 ??傼? (41) 6.3 ???? (44) 6.4 ???? (45) 7 ?? (48)

0 ????????????? 0.1 ?????? ????????????????лй??僔? 1????——?????????????Gambit. 2??????——????Matlab??????? 3????——??Matlab??.dat?????Tecplot??? 0.2 ????????? ?0 Matlab????????? ?????bandwidth.m ????????? elem_B_plane.m ????????B elem_D_plane.m ??????D elem_MK_plane.m ????????K elem_stressIntp.m ??儈??????? input.m ?????? input_readload.m ???? integrate_gauss_1.m 1?儈??? integrate_gauss_2.m 2?儈??? integrate_hammer.m й??Hammer?? main_AssemSolve.m ???????????Ka = P main_CalStress.m ?????????? output2tec.m ??.dat???Tecplot?? plot_mesh.m ?Matlab????? plot_stress.m ?Matlab??????? shape_cubic_820.m 3?8?20??Serendipity??????shape_lagrange.m 1?Lagrange???? shape_quadrangle_9.m ???9???????? shape_quadrangle_48.m ???4?8??Serendipity??????shape_triangle_36.m й??3?6????????

CATIA有限元分析计算实例-完整版

CATIA有限元分析计算实例 CATIA有限元分析计算实例 11.1例题1 受扭矩作用的圆筒 11.1－1划分四面体网格的计算 （1）进入【零部件设计】工作台 启动CATIA软件。单击【开始】→【机械设计】→【零部件设计】选项，如图11－1所示，进入【零部件设计】工作台。 图11－1单击【开始】→【机械设计】→【零部件设计】选项 单击后弹出【新建零部件】对话框，如图11-2所示。在对话框内输入新的零件名称，在本例题中，使用默认的零件名称【Part1】。点击对话框内的【确定】按钮，关闭对话框，进入【零部件设计】工作台。 （2）进入【草图绘制器】工作台 在左边的模型树中单击选中【xy平面】, 如图11-3所示。单击【草图编辑器】工具栏内的【草图】按钮，如图11-4所示。这时进入【草图绘制器】工作台。 图11－2【新建零部件】对话框

图11－3单击选中【xy平面】 （3）绘制两个同心圆草图 点击【轮廓】工具栏内的【圆】按钮，如图11-5所示。在原点点击一点，作为圆草图的圆心位置，然后移动鼠标，绘制一个圆。用同样分方法再绘制一个同心圆，如图11-6所示。 图11－4【草图编辑器】工具栏 图11－5【轮廓】工具栏 下面标注圆的尺寸。点击【约束】工具栏内的【约束】按钮，如图11-7所示。点击选择圆，就标注出圆的直径尺寸。用同样分方法标注另外一个圆的直径，如图11-8所示。 图11－6两个同心圆草图 图11－7【约束】工具栏 双击一个尺寸线，弹出【约束定义】对话框，如图11－9所示。在【直径】数值栏内输入100mm，点击对话框内的【确定】按钮，关闭对话框，同时圆的直径尺寸被修改为100mm。用同样的方法修改第二个圆的直径尺寸为50mm。修改尺寸后的圆如图11-10所示。

《结构分析中的有限元法》2015-有限元习题-参考答案

本科有限元习题参考答案

2015年3月10日作业 1、简述力学课程中介绍的各种力学模型的简化条件、基本假设和适用范围（包括有拉压杆模型、弯曲梁模型、平面应力和平面应变模型、轴对称模型、板模型、壳模型等） 2、给出弹性力学问题中平衡方程、几何方程、物理方程的表达式及其意义。 （1）平衡方程：

zy yz xz zx yx xy z yz xz z y xy zy y x zx yx x f y x z f x z y f z y x ττττττττσττσττσ====+??+??+??=+??+??+??=+??+??+??,000， 物理意义：应力分量与体力分量之间的关系。 （2）几何方程： z u x w y w z v x v y u z w y v x u zx yz xy z y x ??+??=??+??=??+??=??=??=??=γγγεεε,,,, 物理意义：应变分量与位移分量之间的关系。 （3）物理方程： [] [] [] zx zx yz yz xy xy y x z z z x y y z y x x G G G E E E τγτγτγσσμσεσσμσεσσμσε1,1,1) (1 ) (1 )(1 ===+-=+-=+-= 物理意义：应变分量与应力分量之间的关系。 3、简述最小势能原理的主要内容和主要公式。 根据虚功原理得到：??=-Γ T Ω T T 0Td Γδu d Ω）F δu -σδε（，由 )(21εδσεδδεU T T =?? ? ??=则0)21((=Γ-Ω-=∏??ΩΓ）Td u d F u T T T p σεδδ 其中，??ΩΓ Γ-Ω-=∏Td u d F u T T T p )21 (σε即为系统的总势能，它是弹性体变 形势能和外力势能之和。上面变分为零式表明：在所有区域内满足几何关系，在边界上满足给定位移条件的可能位移中，真实位移使系统的总势能取驻值(可证

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

重庆大学有限元第一次作业

有限元分析技术课程大作业 科 目：有限元分析技术 教 师： 姓 名： 学 号： 专 业： 机械设计及理论 类 别： 学 术 上课时间： 2016 年 11 月至 2017 年 1 月 考 生 成 绩： 阅卷评语： 阅卷教师 (签名) 重庆大学研究生院

第一章 问题提出 1.1工程介绍 某露天大型玻璃平面舞台的钢结构如图1所示，每个分格（图2中每个最小的矩形即为一个分格）x 方向尺寸为1m ，y 方向尺寸为1m ；分格的列数（x 向分格）=学生序号的百位数值×10+十位数值+5，分格的行数（y 向分格）=学生序号的个位数值+4，如序号为041的同学分格的列数为9，行数为5，111号同学分格的列数为16，行数为5。 钢结构的主梁（图1中黄色标记单元）为高160宽100厚14的方钢管，其空间摆放形式如图3所示；次梁（图1中紫色标记单元）为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间（如不是正处于X 方向正中间，偏X 坐标小处布置）的次梁的两端，如图2中标记为UxyzRxyz 处。 玻璃采用四点支撑与钢结构连接（采用四点支撑表明垂直作用于玻璃平面的面载荷将传递作用于玻璃所在钢结构分格四周的节点处，表现为点载荷，如图4所示）；试对在垂直于玻璃平面方向的22 /KN m 的面载荷（包括玻璃自重、钢结构自重、活载荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析.（每分格面载荷对于每一支撑点的载荷可等效于0.5KN 的点载荷）。 1.2 作业内容 （1）屏幕截图显示该结构的平面布置结构，图形中应反映所使用软件的部分界面，如图1-2； （2）该结构每个支座的支座反力； （3）该结构节点的最大位移及其所在位置； （4）对该结构中最危险单元（杆件）进行强度校核。 图1-1

有限元法的基本思想及计算 步骤

有限元法的基本思想及计算步骤 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。 用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1）连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 2）单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个结点有两个位移分量u，v和两个结点力分量F x，F y。三个结点共六个结点位移分量可用列

现代设计方法(关于有限元)作业

《现代设计方法》作业关于有限元法的研究 学院：机械工程学院 专业：机械制造及其自动化

0.有限元法 有限元法分析起源于50年代初杆系结构矩阵的分析。随后，Clough于1960年第一次提出了“有限元法”的概念。其基本思想是利用结构离散化的概念，将连续介质体或复杂结构体划分成许多有限大小的子区域的集合体，每一个子区域称为单元（或元素），单元的集合称为网格，实际的连续介质体（或结构体）可以看成是这些单元在它们的节点上相互连接而组成的等效集合体；通过对每个单元力学特性的分析，再将各个单元的特性矩阵组集成可以建立整体结构的力学方程式，即力学计算模型；按照所选用计算程序的要求，输入所需的数据和信息，运用计算机进行求解。 当前，有限元方法/理论已经发展的相当成熟和完善，而计算机技术的不断革新，又在很大程度上推进了有限元法分析在工程技术领域的应用。然而，如此快速地推广和应用使得人们很容易忽视一个前提，即有限元分析软件提供的计算结果是否可靠、满足使用精度的前提，是合理地使用软件和专业的工程分析。有限元法分析一般包括四个步骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算结果的分析。每一个步骤在操作过程中都或多或少地引入了误差，这些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我们的认识和判断。 1.受内压空心圆筒的轴对称有限元分析 例图1.1所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径为120mm的刚性圆孔中，试求圆筒内径处的位移。结构的材料参数

为：200 =，0.3 E GPa μ=。 图1 结构图 对该问题进行有限元分析的过程如下。 （1）结构的离散化与编号 由于该圆筒为无限长，取出中间一段（20mm高），采用两个三角形轴对称单元，如图1.2所示。对该系统进行离散，单元编号及结点编号如图1.3所示，有关结点和单元的信息见表1.1。 图1.2 有限元模型

有限元分析与应用详细例题

《有限元分析与应用》详细例题 试题1：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比 较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 一．问题描述及数学建模 无限长的刚性地基上的三角形大坝受齐顶的水压作用可看作一个平面问题，简化为平面三角形受力问题，把无限长的地基看着平面三角形的底边受固定支座约束的作用，受力面的受力简化为受均布载荷的作用。 二．建模及计算过程 1. 分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算 下面简述三节点常应变单元有限元建模过程（其他类型的建模过程类似）： 1.1进入ANSYS 【开始】→【程序】→ANSYS 10.0→ANSYS Product Launcher →change the working directory →Job Name: shiti1→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 单元是三节点常应变单元，可以用4节点退化表示。 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4 node 42 →OK (back to Element Types window)→Options… →select K3: Plane Strain→OK→Close (the Element Type window) 1.4定义材料参数

有限元分析报告大作业

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

有限元分析软件比较分析

有限元分析软件 有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司，其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析： 1．在世界范围内的知名度：两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS，并且在五年前ANSYS 的界面是当时最好的界面之一，所以在中国，ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域：ANSYS 软件注重应用领域的拓展，目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究，致力于解决该领域的深层次实际问题。 3．性价比：ANSYS 软件由于价格政策灵活，具有多种销售方案，在解决常规的

有限元大作业

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

(完整word版)有限元分析软件的比较

有限元分析软件的比较（购买必看）－转贴 随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element A nalysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面： 增加设计功能，减少设计成本； 缩短设计和分析的循环周期； 增加产品和工程的可靠性； 采用优化设计，降低材料的消耗或成本； 在产品制造或工程施工前预先发现潜在的问题； 模拟各种试验方案，减少试验时间和经费； 进行机械事故分析，查找事故原因。 在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件

有限元法课后习题答案

1、有限元是近似求解一般连续场问题的数值方法 2、有限元法将连续的求解域离散为若干个子域,得到有限个单元,单元和单元之间用节点连接 3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个. 4、平面刚架结构在外力的作用下,横截面上的内力有轴力、剪力、弯矩. 5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 6、平面刚架有限元分析，节点位移有轴向位移、横向位移、转角。 7、在弹性和小变形下，节点力和节点位移关系是线性关系。 8、弹性力学问题的方程个数有15个，未知量个数有15个。 9、弹性力学平面问题方程个数有8，未知数8个。 10、几何方程是研究应变和位移之间关系的方程 11、物理方程是描述应力和应变关系的方程 12、平衡方程反映了应力和体力之间关系的 13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态 14、9形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 15、形函数是_三角形_单元内部坐标的_线性_函数,他反映了单元的_位移_状态 16、在进行节点编号时,同一单元的相邻节点的号码差尽量小. 17、三角形单元的位移模式为_线性位移模式_- 18、矩形单元的位移模式为__双线性位移模式_

19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 20、单元刚度矩阵描述了_节点力_和_节点位移之间的关系 21、矩形单元边界上位移是连续变化的 1.诉述有限元法的定义 答： 有限元法是近似求解一般连续场问题的数值方法 2.有限元法的基本思想是什么 答： 首先，将表示结构的连续离散为若干个子域，单元之间通过其边界上的节点连接成组合体。其次，用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3.有限元法的分类和基本步骤有哪些 答： 分类： 位移法、力法、混合法；步骤： 结构的离散化，单元分析，单元集成，引入约束条件，求解线性方程组，得出节点位移。 4.有限元法有哪些优缺点 答： 优点：

有限元作业

有限元作业

有限元分析大作业 学院： 班级： 姓名： 学号： 日期：

试题一（对应第二章） 如图所示，有一受轴向拉伸载荷2000P N =作用的变截面杆件，在0x =处，杆件截面积为2020A mm =，在180x L mm ==处，杆件截面积为201102 A mm =，杆件弹性模量为200GPa ，泊松比为0.3，试建立该杆件的有限元模型，并计算端部位移。（在划分网格时，沿长度方向取三个等长度杆单元） x P 0A 012A L 解：计算分析 000()(1)2(1)2(1)2x x x x A P x A A x A L P x A L P x E EA L σσσε===-=-==-

[ ]00 000022()[ln(2)]ln 2ln(2)12x x x x P dx PL PL u x dx L x L L x x EA EA EA L ε===--=--??- ??? ?? 0() 1.3860.1242mm PL u L EA == 数学建模：将其用二维模型进行降维处理，分为四个节点，三个等长度单元。 后处理

读出最大应力：1.750*10^2mpa 则计算得到右端部位移u（L）=0.12683 轴向位移随杆长变化图如下：

试题二（对应第三章） 一正方形平板,尺寸为40 mm×40 mm,厚度为2 mm,板中央有直径为d的圆孔如下图所示,板材弹性模量为200GPa,泊松比为0.3,在板的左端和右端分别施加20 MPa的拉力载荷.试建立该平板的有限元模型，并分别计算圆孔直 d=5,10,15,20和25mm时，平板开孔应力集中系数。

第三章平面问题的有限元法作业及答案

第三章 平面问题的有限元法作业 1. 图示一个等腰三角形单元及其节点编码情况，设μ=0，单元厚度为t 。求 1）形函数矩阵[]N ；2）应变矩阵[]B ；3）应力矩阵[]S 。 4 第1题图 第2题图 2. 如题图所示，结构为边长等于a 的正方形，已知其节点位移分别为：11(,)u v 、 22(,)u v 、33(,)u v 、44(,)u v 。试求A 、B 、C 三点的位移。其中A 为正方形形心，B 为三角形形心。 3.直角边边长为l 的三角形单元，如题图所示。试计算单元等效节点载荷列阵（单元厚度为t ，不计自重）。 第3题图 第4题图 4. 如题图所示，各单元均为直角边边长等于l 的直角三角形。试计算（1）单元等效节点载荷列阵；（2）整体等效节点载荷列阵。已知单元厚度为t ，不计自重。

5．下列3个有限元模型网格，哪种节点编号更合理？为什么？ 9 34 6 7912 11 34 6 12142 (a) (b) (c) 第5题图 6.将图示结构画出有限元模型；标出单元号和节点号；给出位移边界条件；并计算半带宽（结构厚度为t ）。 2a (a) (b) 无限长圆筒 (c) 第6题图 7. 结构如图所示，已知结构材料常数E 和 ，单元厚度为t 。利用结构的对称性，采用一个单元，分别计算节点位移和单元应力。 第7题图

答案： 1. 1）形函数 i x N a = ， j y N a = ， 1m x y N a a =-- 2）应变矩阵 []1000101 000101011011B a -????=-??--???? 3）应力矩阵 []100010100 01 0111 110022 2 2S a ? ???-? ?=-????- -? ?? ? 2. A 点的位移为 ()2312A u u u = + ， ()231 2A v v v =+ B 点的位移为 ()24313B u u u u = ++ ， ()2431 3B v v v v =++ C 点的位移为 ()1223C a u u u = + ， ()C 1223 a v v v =+ 3. 单元等效节点载荷列阵为 {}11 11 00003 663 T e i j i j R q q q q ?? =++?? ?? 4. （2）整体等效节点载荷向量为 {}111100006 322T R qlt P qlt P P qlt qlt ?? =-???? 7. （1） 减缩后的整体刚度方程 22 12 2 1222 22221110222021102(1)2 2102x x b b ab R b ab b P v Et ab a b ab ab R v b a μμμ μμμμμμ---??- - ??????????--?????? -??? ?=????---+ +? ???? ?????????-????+?? ? ? 节点位移

有限元答案

1.1有限单元法中“离散”的含义是什么？有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的？位移有限元法的标准化程式是怎样的？ （1）离散的含义即将结构离散化，即用假想的线或面将连续体分割成数目有限的单元，并在其上设定有限个节点；用这些单元组成的单元集合体代替原来的连续体，而场函数的节点值将成为问题的基本未知量。 （2）给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律，即通过插值以单元节点位移表示单元内任意点的位移。因节点位移个数是有限的，故无限自由度问题被转变成了有限自由度问题。 （3）有限元法的标准化程式：结构或区域离散，单元分析，整体分析，数值求解。 1.3单元刚度矩阵和整体刚度矩阵各有哪些性质？各自的物理意义是什么?两者有何区别? 单元刚度矩阵的性质：对称性、奇异性（单元刚度矩阵的行列式为零）。 整体刚度矩阵的性质：对称性、奇异性、稀疏性。 单元Kij物理意义Kij即单元节点位移向量中第j个自由度发生单位位移而其他位移分量为零时，在第j个自由度方向引起的节点力。 整体刚度矩阵K中每一列元素的物理意义是：要迫使结构的某节点位移自由度发生单位位移，而其他节点位移都保持为零的变形状态，在所有个节点上需要施加的节点荷载。 2.2什么叫应变能？什么叫外力势能？试叙述势能变分原理和最小势能原理，并回答下述问题：势能变分原理代表什么控制方程和边界条件？其中附加了哪些条件？ (1)在外力作用下，物体内部将产生应力ζ和应变ε，外力所做的功将以变形能的形式储存起来，这种能量称为应变能。 (2)外力势能就是外力功的负值。 (3)势能变分原理可叙述如下：在所有满足边界条件的协调位移中，那些满足静力平衡条件的位移使物体势能泛函取驻值，即势能的变分为零 δΠp=δUε+δV=0 此即变分方程。对于线性弹性体，势能取最小值，即 δ2ΠP=δ2Uε+δ2V≧0 此时的势能变分原理就是著名的最小势能原理。 势能变分原理代表平衡方程、本构方程和应力边界条件，其中附加了几何方程和位移边界条件。 2.3什么是强形式？什么是弱形式？两者有何区别？建立弱形式的关键步骤是什么？ 等效积分形式通过分部积分，称式 ∫ΩC T(v)D(u)dΩ+∫ΓE T(v)F(u)dΓ为微分方程的弱形式，相对而言，定解问题的微分方程称为强形式。 区别：弱形式得不到解析解。 建立弱形式的关键步骤：对场函数要求较低阶的连续性。 2.4为了使计算结果能够收敛于精确解，位移函数需要满足哪些条件？为什么？ 只要位移函数满足两个基本要求，即完备性和协调性，计算结果便收敛于精确解。 2.6为什么采用变分法求解通常只能得到近似解？变分法的应用常遇到什么困难？Ritz法收敛的条件是什么？ （1）在Ritz 法中，N决定了试探函数的基本形态，待定参数使得场函数具有一定的任意性。如果真实场函数包含在试探函数之内，则变分法得到的解答是精确的；如果试探函数取自完全的函数序列，则当项数不断增加时，近似解将趋近于精确解。然而，通常情况下试探函数不会将真实场函数完全包含在内，实际计算时也不可能取无穷多项。因此，试探函数只能是真实场函数的近似。可见，变分法就是在某个假定的范围内找出最佳解答，近似性就源于此。 （2）采用变分法近似求解，要求在整个求解区域内预先给出满足边界条件的场函数。通常情况下这是不可能的，因而变分法的应用受到了限制。 （3）Ritz 法的收敛条件是要求试探函数具有完备性和连续性，也就是说，如果试探函数满足完备性和连续性的要求，当试探函数的项数趋近于无穷时，则Ritz 法的近似解将趋近于数学微分方程的精确解。 3.1构造单元形函数有哪些基本原则？形函数是定义于单元内坐标的连续函数。单元位移函数通常采用多项式，其中的待定常数应该与单元节点自由度数相等。为满足完备性要求，位移函数中必须包括常函数和一次式，即完全一次多项式。多项式的选取应由低阶到高阶，尽量选择完全多项式以提高单元的精度。若由于项数限制而不能选取完全多项式时，也应使完全多项式具有坐标的对称性，并且一个坐标方向的次数不应超过完全多项式的次数。有时为了使位移函数保持一定阶次的完全多项式，可在单元内部配置节点。然而，这种节点的存在将增加有限元格式和计算上的复杂性，除非不得已才加以采用。形函数应保证用它定义的位移函数满足收敛要求，即满足完备性要求和协调性条件。 3.1构造单元形函数有哪些基本原则？试采用构造单元的几何方法，构造T10 单元的形函数，并对其收敛性进行讨论。 通常单元位移函数采用多项式，其中的待定常数由节点位移参数确定，因此其个数应与单元节点自由度数相等。根据实体结构的几何方程，单元的应变是位移的一次导数。为了反映单元刚体位移和常应变即满足完备性要求，位移函数中必须包含常数项和一次项，即完全一次多项式。 3.3何谓面积坐标？其特点是什么？为什么称其为自然坐标或局部坐标？ （1）三角形单元中，任一点P(x,y)与其3个角点相连形成3个子三角形,其位置可以用下述称为面积坐标的三个比值来确定： L1=A1/A L2=A2/A L3=A3/A 其中A1,A2,A3分别为P23,P31,P12的面积。 （2）面积坐标的特点： a T3单元的形函数Ni就是面积坐标Li b面积坐标与三角形在整体坐标系中的位置无关。 c三个节点的面积坐标分别为节点1（1, 0, 0）、节点2（0, 1, 0）、节点3（0, 0, 1），形心的面积坐标为（1/3, 1/3, 1/3）。 d单元边界方程为Li=0(i=1,2,3) e在平行于23边的一条直线上，所有点都有相同的面积坐标L1（L1对应的三角形具有相同的高和底边），而且L1就等于此直线至23边的距离与节点1至23边的距离之比值。

有限元作业

有限元作业答卷 一、问题解答 1、解：令2 21()()2()2dy p x q x y f x y dx ????=+-?? ??????? π 则可以得到()()y q x y f x π=-， ()y dy p x dx 'π=，()y d d dy p x dx dx dx '?? π= ??? 又有其Euler 方程公式为：0u u d dx ' ππ- = 综上得到原泛函问题的Euler 方程及其边界条件为： ()()(), [,] (1().1.1),()a b d dy p x q x y f x x a b dx dx y a y y b y ?????-+=∈ ?? ?==?? 2、（1）解：引入Sobolev 空间0V H ( )∞=Ω,任取V v ∈乘以方程两端积分： (((,))(,))(,)k x y u q x y u vdxdy f x y vdxdy Ω Ω -???+=???? 再利用格林公式得到： ((,)(,))(,) (,)u k x y u v q x y uv dxdy k x y vds f x y vdxdy n Ω Γ Ω ????+-=?????? 由边界条件得到： (,) ((,)(,))(,)(,) g x y k x y u v q x y uv dxdy f x y vdxdy k x y vds n Ω Ω Γ ????+=+?????? 令 A(,)((,)(,))(,) F()(,)(,)u v k x y u v q x y uv dxdy g x y v f x y vdxdy k x y vds n Ω ΩΓ ?=???+?????=+???????? 得变分方程 A(,)F(), (1.2.1)u v v v V =∈ 其解0u H ()∞∈Ω便为椭圆型方程第一边值问题的Galerkin 意义广义解。 （2）证：下面用Lax-Milgram 定理证明广义解的存在唯一性。 首先，由Hilbert 空间的Schwarz 不等式得到 1 (,)V f v f v f v v ≤≤?∈

CATIA有限元分析计算实例讲诉

CATIA有限元分析计算实例 11.1例题1 受扭矩作用的圆筒 11.1－1划分四面体网格的计算 （1）进入【零部件设计】工作台 启动CATIA软件。单击【开始】→【机械设计】→【零部件设计】选项，如图11－1所示，进入【零部件设计】工作台。 图11－1单击【开始】→【机械设计】→【零部件设计】选项 单击后弹出【新建零部件】对话框，如图11-2所示。在对话框内输入新的零件名称，在本例题中，使用默认的零件名称【Part1】。点击对话框内的【确定】按钮，关闭对话框，进入【零部件设计】工作台。 （2）进入【草图绘制器】工作台 在左边的模型树中单击选中【xy平面】, 如图11-3所示。单击【草图编辑器】工具栏内的【草图】按钮，如图11-4所示。这时进入【草图绘制器】工作台。

图11－2【新建零部件】对话框 图11－3单击选中【xy平面】 （3）绘制两个同心圆草图 点击【轮廓】工具栏内的【圆】按钮，如图11-5所示。在原点点击一点，作为圆草图的圆心位置，然后移动鼠标，绘制一个圆。用同样分方法再绘制一个同心圆，如图11-6所示。 图11－4【草图编辑器】工具栏 图11－5【轮廓】工具栏 下面标注圆的尺寸。点击【约束】工具栏内的【约束】按钮，如图11-7所示。点击选择圆，就标注出圆的直径尺寸。用同样分方法标注另外一个圆的直径，如图11-8所示。

图11－6两个同心圆草图 图11－7【约束】工具栏 双击一个尺寸线，弹出【约束定义】对话框，如图11－9所示。在【直径】数值栏内输入100mm，点击对话框内的【确定】按钮，关闭对话框，同时圆的直径尺寸被修改为100mm。用同样的方法修改第二个圆的直径尺寸为50mm。修改尺寸后的圆如图11-10所示。 图11－8标注直径尺寸的圆草图 图11－9【约束定义】对话框 （4）离开【草图绘制器】工作台 点击【工作台】工具栏内的【退出工作台】按钮，如图11-11所示。退出【草图绘制器】工作台，进入【零部件设计】工作台。 图11－10修改直径尺寸后的圆 图11－11【工作台】工具栏