清华曾攀ansys二次开发资料.pdfVIP

自主程序开发使用说明(fortran,C,ANSYS 平台衔接) 清华大学机械工程系 曾攀 《有限元分析及应用》 自主程序开发使用说明 (fortran, C, 与 ANSYS 平台衔接) 曾 攀 (清华大学机械工程系) 1 连续体平面问题的有限元分析程序(Fortran) 源程序名称：FEM2D.FOR Fortran 源程序、相关接口程序和数据文件见本配书盘中目录 /有限元分析源程序(f,c,ANSYS 衔接)/fortran 源程序及与 ANSYS 衔接(FEM2D)/ 下面将介绍 FEM2D.FOR 程序的实现原理。该程序可以接受由 ANSYS 前处理所输出 的节点信息文件 NODE_ANSYS.IN 和单元信息文件 ELEMENT_ANSYS.IN，经过计算分析 后，输出一个一般性的结果文件 DATA.OUT 和一个专供 ANSYS 进行后处理的结果数据文 件 FOR_POST.DAT 。该程序所使用的需手工生成的模型基本信息文件为 BASIC.IN ，与 ANSYS 后处理衔接的接口程序程序为 USER_POST.LOG 。 1.1 程序原理及实现 该程序的特点如下： 问题类型：可用于计算弹性力学平面应力问题和平面应变问题 单元类型：采用常应变(CST)三角形单元 位移模式：采用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为 “0 ”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独 立约束 方程求解：针对半带宽刚度方程的 Gauss 消元法 节点信息：可以读入由 ANSYS 前处理导出的节点信息文件 NODE_ANSYS.IN ，或手工生 成 单元信息：可以读入由 ANSYS 前处理导出的单元信息文件 ELEMENT_ANSYS.IN ，或手 工生成 结果文件：输出一般的结果文件 DATA.OUT ，还输出供 ANSYS 进行后处理 的文件 FOR_POST.DAT 该程序的原理如框图 1 所示。 1 自主程序开发使用说明(fortran,C,ANSYS 平台衔接) 清华大学机械工程系 曾攀 开 始 输入数据(子程序READ_IN) BASIC.IN （基本信息文件） NODE_ANSYS.IN （节点信息文件，可由ANSYS前处理得到） ELEMENT.ANSYS （单元信息文件，可由ANSYS前处理得到） 形成单元刚度矩阵(子程序FORM_KE) 以半带存储方式形成整体刚度矩阵(BAND_K) 形成节点载荷向量(子程序FORM_P) 处理边界条件(子程序DO_BC) 求解方程获得节点位移(子程序SOLVE) 计算单元及节点应力(子程序CAL_STS) 输出文件 （FOR_POST.DAT,供ANSYS后处理用） （DATA.OUT，一般的结果文件）