1、有限元方法理论及其应用考试题目及要求南京理工大学机械工程学院研究生研究型课程考试题目及要求课程名称： 有限元方法理论及应用 考试形式：专题研究报告 论文 大作业 综合考试考试题目：“有限元方法理论及应用”理论研讨及上机实验试题及要求：一、课程论文：弹性力学有限元位移法原理（30分）撰写一篇论文，对有限元位移法的原理作一般性概括和论述。要求论文论及但不限于下列内容：1）弹性力学有限元位移法的基本思想和数学、力学基础；2）有限元法求解的原理和过程，推导计算列式；对基本概念和矩阵符号进行解释和讨论；3）等参单元的概念、原理和应用。二、分析与计算（40分）1、图示两个结构和单元相似，单元方位相同的平面
2、应力有限元模型，两模型的单元厚度和材料相同。两个模型右端单元边上受均匀剪切面力。对于下列2种情况，试根据有限元法和力学有关知识来分析论证两个模型求解后对应节点（节点1）的位移值和对应单元的应力值之间的关系：1）两个模型面力的合力相等；2）两个模型面力值相等。（10分） (4)对上图（a）,（b）求解刚度矩阵：单元编号ijm(1)124(2)314(3)346(4)536对于（a）(b)刚度矩阵相等结构总的刚度矩阵的组集： (5)外部载荷与约束力：对于第一种情况；（a） (b) 对于第二钟情况：(a) (b) （6）位移矩阵： 有约束条件可知：（7）根据最小势能原理： 进行求解（8）位移和应力值
3、的关系： 对于第一种情况： 节点1的位移： 单元（1）的应力值： 对于第二种情况： 节点1的位移： 单元（1）的应力值：2、证明3节点三角形单元满足协调性条件（相邻单元之间位移连续）。（10分）证明：假设任意两个相邻的三角形单元如图所示：这里采用x,y的一次多项式作为位移插值函数：将广义坐标换为单元节点自由度的二维插值： （1）相邻之间单元连续，先说明插值函数连续； 有方程可知道u,v函数在平面内是连续的。 （2）单元1(i-j-m)和单元2（i-j-n）在i-j边界处连续； （a）先证明u函数在单元间连续 把，代入上面的方程： 有上式可以看出在边界处的位移变化与m,n点的坐标值无关，只是与
4、i,j的坐标值和位移值有关，所以单元1与单元2在i-j边界处u值相等 (b)同理可证得：单元1与单元2在i-j边界处v值相等 (3)结论：由于插值函数连续，单元1与单元2在边界处u,v的值相等，所以三节点三角形单元相邻单元之间位移连续。3、对4节点四边形平面等参元，试验证等参变换能把平面上的正方形母单元映射成为平面上4节点任意四边形单元。（10分）证明：如上图所示，坐标系中的正方形单元1-2-3-4（图b）通过映射关系：可以得到坐标系中的矩形单元i-j-k-l(图a),并保证四个顶点间的映射关系为： 1i 2j 3k 4l 更一般地，如果假设坐标系xoy中的坐标x,y与原坐标系中的映射关系为：
5、 则可以实现坐标系中的正方形单元1-2-3-4（图b）向坐标系xoy中任意直边四边形i-j-k-l(图c)的映射。若进一步假定在两个不同坐标系中，四边形顶点的对应关系为： 1i 2j 3k 4l同时也可以将上式改写成插值函数形式：上式中中的 就是坐标系中的正方形单元顶点1，2，3，4上的拉格朗日插值基函数：所以平面上的正方形母单元经过等参变换能映射成为平面上4节点任意四边形单元4、图示一个一维直杆问题，杆的截面积为A,弹性模量为E。杆受线性变化的轴向线分布力。试构造一种三次杆单元求解该问题，单元有4个节点，节点间隔均匀，形函数可以由形函数性质直接构造或采用拉格朗日插值多项式。整个杆用1个单元离
6、散化。解出节点位移后，由单元有关方程导出单元上位移和应力的函数表达式，并将有限元解与精确解作比较。（10分）解：（1）对于三次单元，首先假设三次函数作为插值函数：（2）等分的四个节点分别为1，2，3，4整理为形函数形式： 采用拉格朗日插值多项式求解形函数：（3）求出形函数：(4)求解刚度矩阵： 经过计算得：（5）求解载荷矩阵： 经过计算得： (6)用最小势能原理得： 其中 由于可以划去第一个方程解出其余三个方程组得：(7)求解插值函数表达式： 所以： 精确解为： （8）单元应力： 经计算得： 精确解为：一 实验题目： 一个200mm200mm平板，中心有一个直径5mm圆孔，左右两边受面内均匀拉
7、伸载荷1MPa。建立平面应力问题有限元模型，分别采用3节点三角形单元和8节点四边形等参元计算孔边应力集中。二 实验目的： 通过采用3节点三角形单元和8节点四边形等参元计算孔边应力集中，对两种单元的求解精度进行比较。经过简单力学分析，可以知道本实验问题属于平面应力问题，基于结构和载荷的对称性，可以只取模型的1/4进行分析。以此来掌握平面问题的有限元分析方法和对称性问题建模的方法。三 建模概述:3.1 定义工作文件和工作标题3.1.1 定义工作文件名 执行Utility Menu-FileChange Jobname3-1,单击OK按钮。3.1.2 定义工作标题 执行Utility Menu-Fi
8、leChange Tilemy work，单击OK按钮。3.1.3 更改当前工作目录 执行Utility Menu-FileChange the working directoryE/STUDY/ANSYS/kaoshi。3.2 定义单元类型、实常数和材料属性3.2.1 设置计算类型 执行Main MenuPreferencesselect StructuralOK。3.2.2 选择单元类型 先：Main MenuPreprocessorElement TypeAdd/Edit/Delete Addselect SolidQuad 8node 182OK。（三节点三角形单元） 后：Main M
9、enuPreprocessorElement TypeAdd/Edit/Delete Addselect SolidQuad 8node 183OK。（八节点四边形单元）如图下图所示：3.2.3 定义实常数 执行Main MenuPreprocessorReal ConstantsAdd/Edit/DeleteAddOKClose。3.2.4 设置材料属性 执行Main MenuPreprocessorMaterial PropsMaterial ModelsStructuralLinearElasticIsotropic输入实常数（在EX框中输入210000，在PRXY框中输入0.3）OK。
10、3.3 创建几何模型3.3.1 生成矩形平面 执行Main MenuPreprocessorModelingCreateAreasRectangleBy 2 Corners输入尺寸OK。3.3.2 生成圆形 执行Main MenuPreprocessorModelingCreateAreasCircleSolid Circle输入圆点坐标及半径OK。3.3.3 布尔运算得到几何模型 执行Main MenuPreprocessorModelingOperateBooleansSubtractArea选择矩形面OK选择圆形OK。如图3-2所示：图3-2 3.4 生成有限元网格 执行Preproce
11、ssorMeshingMesh ToolMesh AreasTetFreeMesh拾取几何模型OK。如图3-3所示的模型。图3-3 初次划分网格后的单元由于板中间有孔存在，应力集中严重，所以应该将孔的边缘网格细化，执行PreprocessorMeshingMesh ToolMesh AreasRefine选中孔周围的网格（如图3-4）OK将划分网格精度定位2OK。得到如图3-5所示。图3-4孔周围要细化的单元图3-5细化后的单元3.5 加载并求解3.5.1 施加约束条件 执行Main MenuSolutionDefine loadsApplyStructuralDisplacementOn L
12、ines，弹出一个拾取框，拾取左边缘线，单击OK按钮，弹出Apply U.ROT on Line对话框，选择“UX”选项，单击OK按钮。 执行Main MenuSolutionDefine loadsApplyStructuralDisplacementOn Lines，弹出一个拾取框，拾取下边缘线，单击OK按钮，弹出Apply U.ROT on Line对话框，选择“UY”选项，单击OK按钮。3.5.2 施加载荷 执行Main MenuSolutionApplyStructuralPressureOn Lines，弹出一个拾取框，拾取右边缘线，单击OK按钮，弹出对话框，所示输入数据-1，单击
13、OK按钮.如图3-6所示图3-63.5.3 求解 执行Main MenuSolutionSolveCurrent LS，弹出一个提示框。浏览后执行file-close，单击OK按钮开始求解运算。出现一个【Solution is done】对话框是单击close按钮完成求解运算。四 计算结果分析与结论：（1）两种情况变形图。（2）两种情况X方向应力分量x应力云图：(3)分别绘制x沿Y轴和X轴的应力大小： （沿Y轴方向的应力大小） （沿X轴方向的应力大小） （4） 对结果的分析： （1）模型x应力分布：x应力集中分布于中心圆孔与x、y轴相交的地方，且与x轴相交处应力为负（应力值最小），与y轴相交处
14、应力为正（应力值最大）；沿圆周向周围，x有最大值迅速减小至最小值；沿y方向的x应力随着Y值的增加迅速下降，沿x方向的x应力随着X的增加先是急速增加后又有所下降最后有缓慢增加。 （2）应力最大值2.95539，最小值-0.0668。误差来源：有限元分析方法是将结构离散化，网格划分得越稀疏，计算出的结果就越 偏离理论值，分的越密集，结果越接近与真实值，三节点三角形与八节点四边形比起来离散的比较稀疏，计算的结果偏离真实值较大。五 实验体会与总结： 通过本次实验，对理论课所学有限元基本方法有了一个更加直观、深入的理解。通过对ANSYS软件处理平面孔的应力集中问题，了解了这款软件的基本应用和它对有限元的
15、一些很好的应用。试验中，遇到诸多问题，仔细思考，加之对讲义的理解，确实很有收获。更增加了对有限元的认识，和对其功能之强大有了更深的理解。一 实验题目： 一个空心球的外半径，内半径。内壁受均匀压力。试用有限元法计算该空心球体的应力分布情况。要求分别应用轴对称二次等参单元建立轴对称模型、应用二次六面体等参单元建立三维模型求解。二 实验目的： 掌握用轴对称二次等参单元建立轴对称模型，并且通过与二次六面体等参单元建立三维模型进新行比较，熟悉轴对称模型的优点。三 建模概述: 为了比较，分别采用三维模型和轴对称模型进行分析。考虑到到三维实体模型的对称性，只对1/8球体划分网格，然后在对称面上施加对称边界条
16、件。其中： 三维实体模型（1/8空心球体）使用Solid186(二次六面体单元)单元建立实体模型及划分网格。 轴对称模型（1/4圆环）使用Plane183（轴对称类型）单元建立实体模型及划分网格。两种模型划分网格如下：然后施加载荷和约束：三维实体模型 载荷为内表面上压力10Mpa；约束条件为三个对称截面的对称条件轴对称模型 载荷为内表面（圆弧段）上压力10Kpa；约束条件是两个对称轴的法向位移为0。四 计算结果分析与结论：1计算结果：（1）径向位移（2）等效应力2 对结果的分析： 对于两种模型的结果进行比较，径向位移的变化完全一样；等效应力的最小值分别为：15.7085和15.7231，对称模
17、型与三维模型比较增加了0.93%，等效应力的最大值分别为：30.7298和30.6975，对称模型与三维模型比较减少了1.1%，表明两种模型的结果基本一致，五 实验体会与总结： 通过本次实验对与二次六面体三维实体模型和轴对称模型有了更深一步的认识，学会了关于对称图形的处理方法，并对轴对称模型有了深一步的理解，轴对称模型与三维模型的解基本一致，但是轴对称模型建模容易，求解迅速，是一种不错的方法。一 实验题目： 一个矩形平板，长1000mm,宽100mm，厚度10mm。材料的E=200GPa， 。板的一侧短边固支，其它三边自由。在相同的较粗网格（厚度方向1层单元）下，分别用8节点六面体全积分等参元
18、和二次六面体等参元计算其前六阶自由振动频率和振型，计算结果列在表中，对计算结果作对比和分析。二 实验目的： 学习用有限元分析结构的自振频率和振型并且比较八节点六面体网格和二次六面体网格的区别。三 建模概述:1建立模型 分别用Solide185建立八节点六面体模型和Solide186建立二次六面体模型。并进行网格划分。如下图：2设置约束条件 执行Main MenuPreprocessorLoadsDefine loadsApplyStructuralDisplacementOn Areas3设置六阶自由振动模型：四 计算结果分析与结论：1.计算结果：（1）两种模型的频率： 八节点六面体单元自由振动频率 二次六面体单元自由振动频率（2）两种的振型 1 八节点六面体单元2二次六面体振型2.计算结果分析：阶数八节点六面体模型频率二次六面体模型频率10.235940.2601521.47821.629332.56982.567644.14194.562054.93614.911768.12808.9427 由两种模型的自由振动频率和振型可以看出两种模型的结果基本上相同五 实验体会与总结 通过本次实验，学习了经过有限元对结构自振频率和阵型的进一步理解，分别用一次单元和二次单元分析，发现结果差别并不是很大，本次实验采用的网格较粗，与真实值可能有些冲突。
