《有限元法基础》考查要求及试卷分析12.docVIP

《有限元法基础》课程 考查考试要求 （2012年11月） 课程简介 该课程为“材料成型及控制工程专业”的专业前沿及选修课，同时也可作为材料学院其它专业化方向的选修课。该课程主要针对材料学科的发展与工程实际，以计算机辅助材料成形过程为背景，重点讲授有限元数值模拟方法的基础性知识，作为计算机虚拟和仿真材料成形过程的基础课，其主要内容如下：①有限元法的力学和数学基础；②有限元法的基本求解过程；③平面问题有限元法的基础知识；④轴对称问题有限元法的基础知识；⑤弹性问题有限元法的程序设计；⑥弹性问题有限元法的商用软件及工程应用。另外还每人安排10个机时的上机大作业，作为该课程的考查内容。 目的及要求 逐步了解有限元法的基础性知识，学习掌握弹性问题有限元法的求解过程，学会有限元程序的设计方法和步骤，了解常用的有限元分析软件，并能解决一些简单的工程问题。 通过课程学习和上机训练，使学生逐步了解有限元法的基础性知识，并能熟练掌握弹性问题有限元法的求解过程和程序设计步骤，熟悉有限元法在材料成型及控制工程专业领域内的应用，并能够分析解决一些简单的工程实际问题。 考试考查形式 重点考查学生在学习课程基础上的综合应用能力，其考查形式为上机考试。每位学生10个学时，在集中的时间内，能够针对某一确定的问题进行分析计算，并在线演示与讲解，回答老师的质询。 该考查课成绩的评定以上机成绩为主。上机考查的备选内容至少五类问题，按照难易程度分别限制最高得分为：60；70；80；90；100，学生根据自己的情况自由选择。每类问题只是给出问题的描述，并不限定具体的求解思路，相关材料的性能参数和求解问题的具体几何形状由学生根据自己掌握的知识自行确定。这样，一方面考查了每位学生解决实际问题的能力，又能避免同组同学的内容重复。 上机程序及说明 （一）程序说明 控制变量 NTYPE：问题类型判断变量； ILOAD：集中力判断数； IEDGE：分布面力判断数； L1：单元列数： L2：单元行数； NSTRE：应力分量数； NELEM：单元数； NPOIN：节点数； NREST：受约束自由度数； N：总自由度数。 总体变量 AA（I，J）：节点坐标数组； NOD（I，J）：单元节点数组； JA1（I）：受约束的自由度号； AREST（I）：受约束的位移值； B（I）：节点载荷列阵； GK（I，J）：整体刚度矩阵； A（I，J）：工作数组。 单元变量 KEM：循环单元号； DELTA：单元面积； T：平面问题时单元厚度； XY（6）：单元三个节点坐标； Z（6，6）：单元刚度矩阵； LEAD（6）：单元六个整体自由度号； A1（3）：存放ai，aj，am； BB（3）：存放bi，bj，bm； CC（3）：存放ci，cj，cm； D（I，J）：弹性矩阵； EB（I，J）：几何矩阵； EP（I）：单元应变列阵； ES（I）：单元应力列阵； S（I）：主应力列阵； EF（I）：单元节点力列阵； NLOAD：作用有集中力的自由度数； LOD（I）：作用有集中力的自由度号； POT（I）：集中力的具体数值； NEDGE：作用有分布面力的单元边数； NOP（I，J）：作用有分布面力的单元边的节点号； Q：分布面力集度； ALFA：分布面力与X轴的夹角。 材料常数 U：波松比； E：弹性模量。 主要子程序 INPUT子程序：生成原始数据； LOADS子程序：计算等效节点载荷； STIFF子程序：计算单元刚度矩阵； ASSEMB子程序：组装整体刚度矩阵； BOUND子程序：处理位移边界条件； GAUSS子程序：列主元高斯消去法解方程； STRES子程序：计算单元应变； XYABC子程序：计算单元的b，c，f，面积； DMATX子程序：计算弹性矩阵； BMATX子程序：计算几何矩阵； MAINS子程序：计算主应力、主方向及等效应力。 （二）上机程序 C MAIN PROGRAM OF PLANE OR AXIS-SYMMETRIC PROBLEMS DIMENSION AA(50,2),NOD(50,3),JA1(20),B(100),GK(100,100), Z(6,6),AREST(20) COMMON N,NELEM,NPOIN,NREST/C1/NTYPE,NSTRE/C2/U,E,T /C3/AA/C4/NOD/C5/GK OPEN(8,FILE= ) OPEN(9,FILE= ,STATUS=NEW) CALL INPUT(JA1,AREST) CALL LOADS(B) DO 10 KEM=1,NELEM CAL