现代设计方法作业.docVIP

《现代设计方法》作业 关于有限元法的研究 学院：机械工程学院 专业：机械制造及其自动化 0 ?有限元法 有限元法分析起源于50年代初杆系结构矩阵的分析。随后， Clough于1960年第一次提出了 “有限元法”的概念。其基木思想是 利用结构离散化的概念，将连续介质体或复杂结构体划分成许多有限 大小的子区域的集合体，每一个子区域称为单元（或元素），单元的 集合称为网格，实际的连续介质体（或结构体）可以看成是这些单元 在它们的节点上相互连接而组成的等效集合体；通过对每个单元力学 特性的分析，再将各个单元的特性矩阵组集成可以建立整体结构的力 学方程式，即力学计算模型；按照所选用计算程序的要求，输入所需 的数据和信息，运用计算机进行求解。 当前，有限元方法/理论C经发展的相当成熟和完善，而计算机技 术的不断革新，又在很大程度上推进了有限元法分析在工程技术领域 的应用。然而，如此快速地推广和应用使得人们很容易忽视一个前提, 即有限元分析软件提供的计算结果是否可靠、满足使用精度的前提， 是合理地使用软件和专业的工程分析。有限元法分析一般包括四个步 骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算 结果的分析。每一个步骤在操作过程中都或多或少地引入了误差，这 些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我 们的认识和判断。 1 ?受内压空心筒的轴对称有限元分析 1 ?受内压空心 筒的轴对称有限元分析 例 图1?1所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径 为120mm的刚性圆孔中，试求圆筒内径处的位移。结构的材料参数 为:E = 200GPg, |1 = 0.3 o 120 图1结构图 对该问题进行有限元分析的过程如下。 结构的离散化与编号 由于该圆筒为无限长，取出中间一段(20mm高)，采用两个三 角形轴对称单元，如图1?2所示。对该系统进行离散，单元编号及结 点编号如图1?3所示，有关结点和单兀的信息见表1.1 o 图1.2有限元模型  图1.3节点位移编号及单元编号表1?1单元编号及结点编号 单元编号 结 点 编 号 ① 1 2 3 ② 2 3 4 结构的结点位移列阵为 / =[冷1鸭冷2 %知3嗎 —4 订(1-1) 结构的结点外载列阵 F = [FrX 0 伦 0 0 0 0 0]r (1.2) 你和?为由内压作用而等效在结点1和结点2上的载荷, 其大小为 fjn = fjnX￡p = 2^x40x20x2 (1. 3) ■ — jUZo 2 约束的支反力矩阵 /? = [0 Rzl 0 Rz2 心心心 r （1.4） (1.5)其中心和心为结点1和结点2在Z方向的约束支反力，(心，RJ 和(氐，RJ为结点3和结点4在厂方向和Z方向的约束支反力。 总的结点载荷列阵 (1.5) P=F + R =[巴 Rz\ Fj 2 坊3 Rz3 Fr4 Z?z4】 各个单元的描述 (3)单元的弹性矩阵为Z)= 105x2.691」5 (3)单元的弹性矩阵为Z)= 105x 2.69 1」5 1」5 0 1.15 2.69 1」5 0 1.15 1」5 2.69 0 0 0 0 0.77 B7 DBlTirdrdz , (1.6) 计算各个单元的刚度矩阵KJ A” K⑴=1。7 -2.58 -2.34 1.45 -1.93 1」3 8.46 1.37 -7.89 1.93 -0.565 2.30 -0.24 0.16 -1.13 对称 7.89 -1.93 0 2.26 0 4.03 0.565 (1.7) 0 -2.22 1.69 -0.085 -1.69 0.64 1.29 -0.645 一 1.29 0 0)2 5」1 -0.645 一 1.29 0 对称 9.66 1.05 -9.02 — ◎ 2.61 0.24 — 人4 9.02 — 4 2.05 (1.8) 建立整体刚度方程 组装整体刚度矩阵并形成整体刚度方程，有 K- 5= P (1.9) (8x8) (8x1) (8x1) 其中 K = K⑴ + K ⑵ （1.10） （4）边界条件的处理及刚度方程求解 边界条件为％=0,叫=。，vv4 = 0 , |ir4 =0 ,叫=。，将其代入方程（1?9） 中，有 「4.03 -2.43「 _5026_ -2.43 4.35 _ 5026 107x 对该方程进行求解得;::囂需? 52） 2 ?误差分析 在上面的例子中，发现有限元法解决实际工程问题产生的误差主 要可能出现在以下几个方面: 对受内压空心圆筒的轴进行了模型简化。简化为理想的不受边界 条件以及z方向变形的轴，施加了 W]=0,W2=0,W3=0,W4=0的边界条件。 理想化误差是在有限元法分析开始之前引入的，因此我们不可能通过 改进有限