有限元法分析与建模课程设计报告学院：机电学院专业：机械制造及其自动化指导教师：****学生：* ***学号：***********2015-12-31摘要本文通过ANSYS10.0建立了标准光盘的离心力分析模型，采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析，同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。
力求较为真实地反映光盘在光驱中实际应力和应变分布情况，为人们进行合理的标准光盘结构设计和制造工艺提供理论依据。
关键词：ANSYS10.0；光盘；应力；应变。
目录第一章引言 (3)1.1 引言 (3)第二章问题描述 (4)2.1有限元法及其基本思想 (4)2.2 问题描述 (4)第三章力学模型的建立和求解 (5)3.1设定分析作业名和标题 (5)3.2定义单元类型 (6)3.3定义实常数 (9)3.4定义材料属性 (12)3.5建立盘面模型 (14)3.6对盘面划分网格 (22)3.7施加位移边界 (27)3.8施加转速惯性载荷并求解 (30)第四章结果分析 (32)4.1 旋转结果坐标系 (32)4.2查看变形 (33)4.3查看应力 (35)总结 (38)参考文献 (39)第一章引言1.1 引言光盘业是我国信息化建设中发展迅速的产业之一，认真研究光盘产业的规律和发展趋势，是一件非常迫切的工作。
光盘产业发展的整体性强，宏观调控要求高，因此，对于光盘产业的总体部署、合理布局和有序发展等问题，包括节目制作、软件开发、硬件制造、节目生产、技术标准等。
在高速光盘驱动器中，光盘片会产生应力和应变，在用ANSYS分析时，要施加盘片高速旋转引起的惯性载荷，即可以施加角速度。
需要注意的是，利用ANSYS施加边界条件时，要将内孔边缘节点的周向位移固定，为施加周向位移，而且还需要将节点坐标系旋转到柱坐标系下。
本文通过ANSYS10.0建立了标准光盘的离心力分析模型，采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析，同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。
第二章 问题描述2.1有限元法及其基本思想有限元法的基本思想是将连续的求解区域离散为一组有限个、且按一定方式相互联结在一起的单元组合体。
由于单元能按不同的联结方式进行组合，且单元本身可以有不同形状，因此可以模型化几何形状复杂的求解区域。
有限元法作为数值分析方法的一个重要特点是利用在每一个单元内假设的近似函数，分片地表示全求解域上待求的未知场函数。
单元内的近似函数通常由未知场函数或其导数在单元的各个节点的数值和其插值函数表达。
这样，一个问题的有限元分析中，未知场函数或其导数在各个节点上的数值就成为新的未知量（即自由度），从而使一个连续的无限自由度问题变成离散的有限自由度问题。
一经求解出这些未知量，就可通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上的近似解。
显然，随着单元数目的增加，即单元尺寸的缩小，或者随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断改进。
如果单元是满足收敛要求的，近似解最后将收敛于精确解。
2.2 问题描述标准光盘，置于52倍速的光驱中处于最大读取速度（约为1000min r ），计算其应力分布。
标准光盘参数：外径：120mm 内孔径：15mm 厚度：1.1mm 弹性模量：a MP 4106.1⨯ 密度：33102.2m Kg ⨯第三章力学模型的建立和求解3.1设定分析作业名和标题（1）从常用菜单中选择File--Change jobname命令，将弹出Change jobname对话框，在Enter new jobname文本框中输入文字“CH01”为数据库文件名，单击OK完成文件名的修改，如图所示。
（2）从实用菜单中选择File--Change Title命令，将弹出Change Title对话框，在Enter new title文本框中输入文字“CD”为标题名，单击OK完成文件名的修改，如图所示。
（3）选取菜单路径Utility Menu/Plot /Replot，指定的标题“CD”将显示在图形窗口的左下角，如图所示。
（4）选取菜单路径Main Menu/Preference，将弹出Preference of GUI Filtering（菜单过滤参数选择）对话框，选中Structural复选框，单击OK按钮确定，如图所示。
3.2定义单元类型在进行有限元分析时，首先应根据分析问题的几何结构、分析类型和所分析的问题精度要求等，选定适合具体分析的单元类型。
本课题用选4节点四边形板单元PLANE42。
（1）从主菜单中选择Main Menu/Preprocessor/Element/Add/Edit /Delete命令，弹出Element Types（单元类型）对话框，如图所示。
（2）单击Add...按钮，将弹出Library of Element Types（单元类型库）对话框，如图所示。
（3）在左边的列表框中选择“Solid”选项，选择实体单元类型。
在右边的列表框中选择“Quad 4node 42”选项，选择4节点四边形板单元PLANE42，单击OK按钮，将PLANE42单元添加，并关闭单元类型库对话框，同时返回到第一步弹出的单元类型对话框，如图所示。
（4）单击Options…弹出如图所示的PLANE42 element type options（单元选项设置）对话框，对PLANE单元进行设置，使其可用于计算平面应力问题，如图所示。
（5）在Element behavior（单元行为方式）下拉列表框中选择Plane-stress w/thk(带有厚度的平面应力)选项。
单击OK按钮，接受选项，关闭单元选项设置对话框，返回到单元类型对话框，完成单元类型定义，如图所示。
（6）单击Close按钮，关闭单元类型对话框，结束单元类型的添加，如图所示。
3.3定义实常数本课题中选用带有厚度的平面应力行为方式的PLANE42单元，需要设置其厚度实常数。
（1）从主菜单中选择Main Menu/Preprocessor/Real Constants/Add/Edit/Delete命令，弹出Real Constants(实常数)对话框，如图所示。
（2）单击add…按钮，弹出如图所示Element Type for Real Constants （实常数单元类型）对话框，要求选择欲定义实常数的单元类型。
本课程设计中只定义了一种单元类型，在已知的单元类型列表中选择“Type 1 PLANE42”，将要为PLANE42单元类型定义实常数。
单击OK弹出Real Constant Set（实常数集）对话框，如图所示。
（3)在Thickness（厚度）文本框中输入“1.1”，单击OK，关闭实常数输入对话框，返回到实常数对话框，如图所示已经定义了一组实常数，如图所示。
(4)单击Close按钮，关闭实常数设置对话框，如图所示。
3.4定义材料属性考虑惯性力的静力分析中必须定义材料的弹性模量和密度。
具体步骤如下：（1）从主菜单中选择Main Menu/Preprocessor/Material Props/Material Models命令，将弹出Define Material Model Behavior(定义材料模型)窗口，如图所示。
（2）依次双击Structural/Linear/Elastic/Isotropic,将弹出1号材料的弹性模量EX和泊松比PRXY的定义对话框，如图所示。
（3）在对话框EX文本框中输入弹性模量EX为1.6e4，在PRXY文本框中输入泊松比0.3，单击OK按钮，关闭对话框，并返回到定义材料属性对话框，在定义材料属性会话框的左边一栏出现刚刚定义的参与号为1的材料属性，如图所示。
（4）依次双击Structural/Density弹出定义密度对话框，在DENS文本框中输入密度数值“2.2e-9”，单位为。
单击OK按钮，退出定义材料模型属性窗口，完成对材料模型属性的定义，如图所示。
3.5建立盘面模型本课程设计的模型比较简单，仅有一个平面。
需要注意的是使用PLANE系列单元时，要求模型必须位于全局XY平面内。
默认的工作平面即位于全局XY 平面内，因此可以直接在默认的工作平面内创建圆环面，为了对圆环面划分有限元网格，还需要将圆环面分成两半。
（1）从主菜单中选择Main Menu/Preprocessor/Modeling/Create/Areas/Circle/ Annulus命令，弹出Annular Circular Area(圆环面)对话框，如图所示。
（2）在WP X文本框中输入圆环面圆心的X坐标（相对于工作平面）“0”；在WP Y文本框中输入圆环面圆心的Y坐标（相对于工作平面）“0”；在Rad-1文本框中输入内孔半径“7.5”；在Rad-2文本框中输入内孔半径“60”。
单击OK按钮确认，关闭对话框。
同时ANSYS在当前工作平面的原点位置创建内径为15mm，外径为120mm的圆环面，如图所示。
（3）从实用菜单中选择Main Menu/PlotCtrls/Numbering命令，弹出Plot Numbering Controls（显示图元编号控制）对话框，如图所示。
选中Keypoint numbers(关键点编号)复选框。
单击OK按钮，使之生效，并对图形窗口进行重绘，如图所示。
图为显示图元编号控制对话框图为创建的圆环面（4）从主菜单中选择Main Menu/Modeling/Create/Lines/Lines /Straight Line命令，弹出关键点选择对话框，要求选择要创建的直线的两个端点。
用鼠标在图形窗口中选择选取关键点1、3，创建出端点为关键点1,3的直线，如图所示。
（5）单击OK按钮，关闭选择对话框，如图所示。
（6）从实用菜单中选择Utility Menu/Plot/Multi-Plots命令，在图形窗口显示所有图元，如图所示。
图为创建的盘面和切割盘面的线（7）从主菜单中选择Main Menu/Preprocessor/Modeling/Operate/Booleans/Divide/Area by Line命令，弹出面选择对话框，要求选择欲被分割的面，如图所示。
(8) 选择创建的圆环面，并单击选择对话框的OK按钮，又弹出选择对话框，要强求选择分割面的线，如图所示。
（9）选择在第（4）步中创建的线，并单击选择对话框的OK按钮，如图所示。
（10）ANSYS将进行布尔运算，将选定的面分割开来成为两个面，如图所示。
如图是被分割后的两个面3.6对盘面划分网格（1）从主菜单中选择Main Menu/Preprocessor/Meshing/MeshTool 命令，弹出Mesh Tool，单击Lines域Set按钮，点击Pick All按钮，弹出Element Sizes on Picked Lines对话框，在NDIV No. of element divisions 文本框中输入“10”单击OK按钮，如图所示。