有限元程序设计报告材料报告材料

《有限元程序设计报告材料报告材料》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、实用标准文案有限元程序设计报告课程名称:有限元程序设计指导教师：张亮学校：重庆大学专业：工程力学01班姓名：苏世宏学号：201266992015年7月8日文档实用标准文案有限元程序设计报告一、前言有限元方法（theFiniteElementMethod）是起源与上个世纪50、60年代，基于弹性力学变分原理的一种近似计算方法，也是当今工程分析中获得最广泛应用的数值计算方法。由于它的通用性和有效性，受到工程技术界的高度重视。伴随着计算机科学和技术的快速发展，现已成为计算机辅助设计（CAD）和计算机辅助

2、制造（CAM）的重要组成部分。有限元程序系统通常包括前处理、有限元程序本体和后处理三部分。前处理包括几何实体模型的建立、材料参数的赋值、位移边界条件的定义、载荷的定义、分析问题类型的定义、单元类型的选择和网格的划分等。（分析问题类型如静力分析、动力特性分析、动力响应、温度场分析、电磁场分析、流体动力学分析等）有限元程序本体是有限元程序系统的核心部分，其功能是实现各种问题的计算。后处理则是将计算结果用图形、曲线和表格的形式表达。（通常包括结构的变形图、应力、应变分布云图等）本课程设计则是针对有限元程

3、序本体，参照教学程序（FEATP），编写简单的有限元程序以计算简单的平面应力、平面应变和轴对称问题，并将其结果与有限元商用软件（ANSYS）的计算结果，以及问题的理论值进行比较，从而验证程序以及问题模型建立的正确性。1．设计目的1）通过编写简单的有限元程序熟悉用有限元方法解决实际问题的基本步骤和过程，体会这种方法的处理手段。2）在VisualFortran中编写程序，熟悉并巩固Fortran语言的语法、算法，学习程序的调试方法，并体会其在执行某个具体算例时，文件的输入、输出以及程序的执行过程。2．

4、设计内容1）以教学程序（FEATP）为参照，编写程序，计算简单的平面应力（PlaneStress），平面应变(PlaneStrain)问题，验证程序的正确性。2）在具体的算例中，对同一问题，在程序和ANSYS中采用不同的单元和网格划分方式，将其结果与理论值进行对比，体会不同的单元和网格划分对问题解的影响，从而判断模型的正确性和合理性。3）总结在编写程序和算例中遇到的问题和解决方法，写出自己的心得体会。二、弹性力学平面问题有限元方法的基本公式文档实用标准文案平面问题1．三角形单元文档实用标准文案文档

5、实用标准文案(1)位移éuiù文档实用标准文案êúêviúéùéùêú文档实用标准文案uNiu=êú=ê0Nj0Nm0ujúêú文档实用标准文案ëvûë0Ni0Nj0Nmûêvjú文档实用标准文案êêú文档实用标准文案éaiùumúêëvmúû文档实用标准文案êiú=[ININjINm]êajú文档实用标准文案ij=[NNêëamúûmN]ae=Nae文档实用标准文案其中，形函数ai,bi,ci是取决于节点坐标的常数应变矩阵B的分块子矩阵是文档实用标准文案3节点单元的应变矩阵

6、是1文档实用标准文案ú文档实用标准文案(2)形成单元的刚度矩阵和等效节点载荷列阵òKe=BTDBdV=BTDBtAVePe=Pe+Pe+Pe+Pefss0e0文档实用标准文案BT(3)集成结构的刚度矩阵和等效节点载荷列阵文档实用标准文案K=åKe=åòVDBdV文档实用标准文案eP=Pfe0+PsPsee+P0+PF文档实用标准文案åfSs0e0F文档实用标准文案=(Pe+PePee+Pe)+P文档实用标准文案(4)引入强制边界条件（消除K的奇异性）(5)求解有限元解方程，得到节点位

7、移aKa=P(6)计算单元应变和应力éexù文档实用标准文案ê[yme=êeêgúeú=Lu=LNa=LNiNjúN]ae文档实用标准文案mëi=[BxyûBjB]ae=Bae文档实用标准文案úésxù文档实用标准文案s=ê文档实用标准文案êsyêtú=De=DBae=Saeú文档实用标准文案ëxyûijmijm应力矩阵S=DB=D[BBB]=[SSS]文档实用标准文案S的分块子矩阵为éêbiùn0ciú文档实用标准文案Si=DBi=E0ê2ên0bicú（i,j,m）iú文档

8、实用标准文案2(1-n0)Aê1-n01-n0ú文档实用标准文案ëê2ci2biúû文档实用标准文案对于平面应力问题E0=E,n0=n文档实用标准文案文档实用标准文案对于平面应变问题3．四边形单元用相同的推导方法可以得到四边形单元位移、应变、应力的有限元表达格式，它具有和三角形单元相似的表达形式，这里就不一一列举了。三、有限元程序设计1．程序功能本程序是在教学程序（FEATP）的基础上修改，删减而成的，主要应用于解决各项同性的弹性力学二维问题，平面应力问题，平面应变问题和轴对称