有限元分析程序的设计
有限元分析实验报告

有限元分析实验报告有限元分析实验报告引言有限元分析是一种广泛应用于工程领域的数值计算方法,它可以通过将复杂的结构划分为许多小的有限元单元,通过计算每个单元的力学特性,来模拟和预测结构的行为。
本实验旨在通过有限元分析方法,对某一结构进行力学性能的分析和评估。
实验目的本实验的目的是通过有限元分析,对某一结构进行应力和变形的分析,了解该结构的强度和稳定性,为结构设计和优化提供参考。
实验原理有限元分析是一种基于弹性力学原理的数值计算方法。
它将结构划分为许多小的有限元单元,每个单元都有自己的力学特性和节点,通过计算每个单元的应力和变形,再将其组合起来得到整个结构的力学行为。
实验步骤1. 建立有限元模型:根据实际结构的几何形状和材料特性,使用有限元软件建立结构的有限元模型。
2. 网格划分:将结构划分为许多小的有限元单元,每个单元都有自己的节点和单元材料特性。
3. 材料参数设置:根据实际材料的力学特性,设置每个单元的材料参数,如弹性模量、泊松比等。
4. 载荷和边界条件设置:根据实际工况,设置结构的载荷和边界条件,如受力方向、大小等。
5. 求解有限元方程:根据有限元方法,求解结构的位移和应力。
6. 结果分析:根据求解结果,分析结构的应力分布、变形情况等。
实验结果与分析通过有限元分析,我们得到了结构的应力和变形情况。
根据分析结果,可以得出以下结论:1. 结构的应力分布:通过色彩图和云图等方式,我们可以清楚地看到结构中各个部位的应力分布情况。
通过对应力分布的分析,我们可以了解结构的强度分布情况,判断结构是否存在应力集中的问题。
2. 结构的变形情况:通过对结构的位移分析,我们可以了解结构在受力下的变形情况。
通过对变形情况的分析,可以判断结构的刚度和稳定性,并为结构的设计和优化提供参考。
实验结论通过有限元分析,我们对某一结构的应力和变形进行了分析和评估。
通过对应力分布和变形情况的分析,我们可以判断结构的强度和稳定性,并为结构的设计和优化提供参考。
有限元分析毕业设计

有限元分析毕业设计有限元分析毕业设计毕业设计是大学生在学业结束前的一项重要任务,也是对所学知识的综合应用和实践能力的考验。
在工程类专业中,有限元分析是一种常见的工程设计方法,被广泛应用于各个领域,如机械、土木、航空等。
本文将探讨有限元分析在毕业设计中的应用。
一、有限元分析的基本原理有限元分析是一种基于数值计算的工程设计方法,通过将复杂的结构划分为有限个简单的单元,利用数学方法求解各个单元的力学行为,最终得到整个结构的力学性能。
有限元分析的基本原理是将结构分割为有限个单元,每个单元都有一组未知的位移和应力,通过建立单元之间的关系,利用数值方法求解出这些未知量。
二、有限元分析在毕业设计中的应用1. 结构强度分析在毕业设计中,结构强度分析是一个重要的环节。
通过有限元分析,可以模拟结构在不同载荷下的受力情况,评估结构的强度和稳定性。
例如,在机械工程的毕业设计中,可以利用有限元分析来评估零件的强度,确定合适的材料和尺寸,从而提高产品的可靠性和安全性。
2. 热传导分析热传导分析是另一个常见的应用领域。
在毕业设计中,有时需要对材料或结构在不同温度下的热传导性能进行分析。
有限元分析可以模拟材料的热传导行为,预测温度分布和热流量。
例如,在建筑工程的毕业设计中,可以利用有限元分析来评估建筑物的保温性能,优化建筑材料的选择和结构设计。
3. 流体力学分析流体力学分析是有限元分析的另一个重要应用领域。
在毕业设计中,有时需要对流体在管道、泵站、水利工程等系统中的流动行为进行分析。
有限元分析可以模拟流体的流动特性,预测流速、压力分布和流量。
例如,在水利工程的毕业设计中,可以利用有限元分析来评估水流在河道中的流动情况,优化河道的设计和水利工程的规划。
三、有限元分析的优势和局限性有限元分析作为一种数值计算方法,具有一些明显的优势。
首先,它可以模拟复杂的结构和物理现象,提供准确的数值结果。
其次,有限元分析具有灵活性,可以根据不同的需求进行模型的建立和分析。
有限元单元法程序设计

有限元单元法程序设计有限元单元法是一种用于工程结构分析和设计的计算方法,它将大型结构分解为许多小的离散单元,通过分析单元之间的相互作用来预测结构的力学行为。
有限元单元法程序设计是指针对特定工程问题,编写计算机程序来实现有限元分析的过程。
下面将介绍有限元单元法程序设计的基本流程和关键要点。
一、问题建模和网格划分有限元单元法程序设计的第一步是对工程结构进行合理的建模和网格划分。
建模的目的是将实际结构抽象为适用于有限元分析的数学模型,包括定义结构的几何特征、材料属性、边界条件等。
网格划分是将结构分解为许多小的单元,每个单元具有一定的形状和尺寸,以便于数值计算。
常用的单元形状包括三角形、四边形、四面体、六面体等,根据结构的特点选择合适的单元形状和尺寸。
二、单元刚度矩阵和载荷矩阵的求解在有限元单元法程序设计中,需要编写算法来求解每个单元的刚度矩阵和载荷矩阵。
单元刚度矩阵描述了单元内部的力学性能,包括刚度、弹性模量、泊松比等,它们通常通过数学公式或有限元理论推导得到。
载荷矩阵描述了单元受到的外部荷载,可以是均匀分布载荷、集中载荷或者边界条件引起的约束力。
通过合适的数值积分方法,可以计算得到每个单元的刚度矩阵和载荷矩阵。
三、组装全局刚度矩阵和载荷向量在有限元单元法程序设计中,需要将所有单元的刚度矩阵和载荷向量组装成整个结构的全局刚度矩阵和载荷向量。
这涉及到单元之间的连接关系以及边界条件的处理。
采用适当的组装算法,可以将各个单元的刚度矩阵和载荷向量叠加在一起,形成整个结构的刚度矩阵和载荷向量。
四、求解位移和应力有限元单元法程序设计的最后一步是求解结构的位移和应力。
通过斯蒂芬-泰勒算法或者其他迭代算法,可以得到整个结构的位移分布,然后根据位移场计算各个点的应变和应力。
这一过程涉及到对整个结构刚度矩阵的求解和对位移的后处理。
有限元单元法程序设计是一个复杂而又精密的工作,需要深入理解有限元原理、结构力学知识和数学方法。
有限元分析程序设计

结构有限元分析程序设计绪论§0.1 开设“有限元程序设计”课程的意义和目的§0.2 课程特点§0.3 课程安排§0.4 课程要求§0.5 基本方法复习$0.1 意义和目的1.有限元数值分析技术本身要求工程设计研究人员掌握1). 有限元数值分析技术的完善标志着现代计算力学的真正成熟和实用化,已在各种力学中得到了广泛的应用。
比如:,已杨为工程结构分析中最得以收敛的技术手段,现代功用大致有:a). 现代结构论证。
对结构设计从内力,位移等方面进行优劣评定,从而进行结构优化设计。
b)可取代部份实验,局部实验+有限元分析,是现代工程设计研究方法的一大特点。
c)结构的各种功能分析(疲劳断裂,可靠性分析等)都以有限元分析工具作为核心的计算工具。
2). 有限元数值分析本身包括着理论+技术实现(本身功用所绝定的)有限元数值分析本身包括着泛函理论+分片插值函数+程序设计2. 有限元分析的技术实现(近十佘年的事)更依赖于计算机程序设计有限元分析的技术取得的巨大的成就,从某种意义上说,得益于计算机硬件技术的发展和程序设计技术的发展,这两者的依赖性在当代表现得更加突出。
(如可视化技术)3.从学习的角度,不仅要学习理论,而且要从程序设计设计角度对这些理论的技术实现有一个深入的了解,应当致力于掌握这些技术实现能力,从而开发它,发展它。
(理论本身还有待于进一步完美相应的程序设计必须去开发)4.程序设计不仅是实现有限元数值分析的工具和桥梁,而且在以下诸方面也有意义:1). 精通基本概念,深化理论认识;2). 锻炼实际工程分析,实际动手的能力;3). 获得以后工作中必备的工具。
(作业+老师给元素库)目的:通过讲述有限元程序设计的技术与技巧,便能达到自编自读的能力。
§0.2 课程特点总描述:理论+算法+数据结构(程序设计的意义)理论:有限元算法,构造,步骤,解的等外性,收敛性,稳定性,误差分析算法;指求解过程的技术方法,含两方面的含义;a. 有限元数值分析算法,b, 与数据结构有关的算法(总刚稀疏存贮,提取,节点优化编号等)数据结构:指各向量矩阵存贮管理与实现,辅助管理结构(指针,数据记录等)具体特点:理论性强:能量泛函理论+有限元构造算法+数据结构构造算法内容繁杂:理论方法+技术方法+技术技巧技巧性强:排序,管理结构(指针生成,整型运算等)§0.3 课程安排①. 单元刚度矩阵及元素设计(单元刚阵算法,杆梁平面分析,板弯非协调元等)②. 总刚的形式及程序设计(单刚提前准备,技术复杂)③. l边界条件及程序设计(等效荷载计算,位移边界条件置入,多工况的对称性)④. 总刚线性方程组求解(LDL T分解,分块算法,子结构算法,波前法)⑤.单元应力计算+应力处理与改善。
有限元分析在工程设计中的应用案例分析

有限元分析在工程设计中的应用案例分析有限元分析,简称FEA(Finite Element Analysis),是一种利用数值计算方法对复杂结构进行力学分析的技术。
它基于物理学原理,利用离散化方法将连续的结构在有限元上分解成多个互相联系但是局部地独立的单元,再通过数学算法进行求解,最终得到整个结构的力学行为。
因为它可以减少试错周期、降低开发成本和提高产品性能,所以有限元分析已经成为当今工程设计和生产领域一项非常重要的技术。
本文将介绍一些有限元分析在工程设计中的具体应用案例。
1.汽车发动机壳体优化汽车发动机壳体是承载引擎所有关键部件的重要结构,其制造复杂度很高。
为了减少开发过程中的试验成本和时间,一家风机厂专门利用有限元分析技术对汽车发动机壳体进行优化设计。
更改前发动机壳体在经过一定的较高频振动时会存在密封性能下降的现象,需要进行加强设计。
利用有限元分析技术,他们对发动机壳体进行了动力学分析,并计算了各部位的振动位移和应力分布,通过不断地修改控制点的位置和形状来提高振动阻尼性能和密封性能。
最终确定了优化方案,成功地减少了振动,提高了发动机壳体的防震性能和密封性能。
2.建筑物钢框架分析建筑物钢框架是建筑结构的重要组成部分,其承载能力和组装结构设计都需要严格控制。
如何选取更好的工艺和材料来设计出更安全可靠的钢框架结构,被许多建筑设计公司所思考。
有限元分析技术的应用可以帮助工程师确定结构的承载能力,最大应力极限和变形情况,进而实现结构的优化。
一家建筑设施的设计公司利用有限元分析技术来优化钢框架的结构,计算具体承载状况,最终确定钢框架结构的有效设计方案。
这一个优化设计方案进一步增强了建筑物钢框架的承载能力,提高了项目的整体优势性。
3.飞机负荷分析航空工业是重要的现代国家产业之一。
飞机设计、测试和生产都需要极高的准确性,而这需要大量的场地、人力和物资投入。
一家工程公司成功地利用有限元分析技术对飞机进行负荷分析并评估整体结构的强度和刚度。
课程设计ANSYS有限元分析(最完整)

有限元法分析与建模课程设计报告学院:机电学院专业:机械制造及其自动化指导教师:****学生:****学号:2012011****2015-12-31摘要本文通过ANSYS10.0建立了标准光盘的离心力分析模型,采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析,同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。
力求较为真实地反映光盘在光驱中实际应力和应变分布情况,为人们进行合理的标准光盘结构设计和制造工艺提供理论依据。
关键词:ANSYS10.0;光盘;应力;应变。
目录第一章引言31.1 引言3第二章问题描述52.1有限元法及其基本思想52.2 问题描述5第三章力学模型的建立和求解63.1设定分析作业名和标题63.2定义单元类型73.3定义实常数103.4定义材料属性133.5建立盘面模型153.6对盘面划分网格233.7施加位移边界283.8施加转速惯性载荷并求解31第四章结果分析334.1 旋转结果坐标系334.2查看变形344.3查看应力36总结39参考文献40第一章引言1.1 引言光盘业是我国信息化建设中发展迅速的产业之一,认真研究光盘产业的规律和发展趋势,是一件非常迫切的工作。
光盘产业发展的整体性强,宏观调控要求高,因此,对于光盘产业的总体部署、合理布局和有序发展等问题,包括节目制作、软件开发、硬件制造、节目生产、技术标准等。
在高速光盘驱动器中,光盘片会产生应力和应变,在用ANSYS分析时,要施加盘片高速旋转引起的惯性载荷,即可以施加角速度。
需要注意的是,利用ANSYS施加边界条件时,要将孔边缘节点的周向位移固定,为施加周向位移,而且还需要将节点坐标系旋转到柱坐标系下。
本文通过ANSYS10.0建立了标准光盘的离心力分析模型,采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析,同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。
有限元分析实验报告

有限元分析实验报告引言有限元分析是一种工程设计和分析的常用方法。
它通过将结构或物体分割为有限数量的单元,利用数值方法计算每个单元的行为,最终得出整体结构的行为。
本实验使用有限元分析方法来研究一个特定的结构或物体。
实验目的本实验的目的是使用有限元分析方法研究一个给定的结构或物体。
通过实验,我们将探索结构的强度、刚度和变形等性能,评估其设计的合理性,并提出改进的建议。
实验步骤实验的步骤如下:1.准备工作:收集和整理所需的材料和数据,包括结构的几何形状、材料特性和加载条件等。
确保所收集的数据准确无误。
2.建立有限元模型:将结构的几何形状转化为有限元模型。
根据结构的复杂程度和要求,选择合适的单元类型和网格密度。
使用有限元软件,如ANSYS、ABAQUS等,建立有限元模型。
3.定义边界条件:根据实际应用场景,定义结构的边界条件。
这些条件包括约束边界条件和加载边界条件。
约束边界条件用于限制结构的自由度,加载边界条件用于施加外部载荷。
4.分析结构的行为:使用有限元软件进行结构的强度、刚度和变形等分析。
根据加载和边界条件,计算结构在不同工况下的应力、位移和应变等结果。
5.结果分析和讨论:评估结构的性能,比较不同工况下的结果,分析结构的弱点和改进的空间。
提出改进的建议,并讨论其可能的影响和成本。
6.撰写实验报告:根据实验结果和讨论,撰写实验报告。
报告应包括实验目的、方法、结果和讨论等部分。
确保报告的结构清晰,表达准确。
结果与讨论根据实验的结果和讨论,我们得出以下结论:1.结构的强度:分析结果显示,结构在给定的加载条件下具有足够的强度,能够承受预期的载荷。
然而,在某些关键部位,应力集中现象可能会导致局部的应力超过材料的极限强度。
2.结构的刚度:结构的刚度是指结构在受力下的变形情况。
分析结果显示,结构在加载后会发生一定的变形,但变形量较小,不会对结构的正常功能产生明显的影响。
3.结构的优化:根据分析结果和讨论,我们提出了改进结构的建议。
有限元分析课程设计活塞

有限元分析课程设计活塞一、教学目标本课程的教学目标是使学生掌握有限元分析的基本原理和方法,能够熟练运用有限元分析软件进行工程问题的分析和计算。
具体分为以下三个部分:1.知识目标:使学生了解有限元分析的基本概念、原理和方法,掌握有限元分析的基本步骤和技巧,熟悉常见的有限元分析软件。
2.技能目标:通过实例教学,使学生能够熟练运用有限元分析软件进行简单的工程分析和计算,能够独立完成有限元分析的基本操作。
3.情感态度价值观目标:培养学生对工程问题的分析和解决能力,提高学生的科学素养和创新能力,使学生能够认识到有限元分析在工程实际中的应用价值。
二、教学内容本课程的教学内容主要包括有限元分析的基本概念、原理和方法,有限元分析的基本步骤和技巧,以及有限元分析软件的应用。
具体安排如下:1.有限元分析的基本概念:介绍有限元分析的定义、特点和发展历程。
2.有限元分析的基本原理:讲解有限元分析的基本原理,包括离散化方法、节点和元素的概念。
3.有限元分析的基本方法:介绍有限元分析的基本方法,包括静态分析、动态分析和温度分析等。
4.有限元分析的基本步骤:讲解有限元分析的基本步骤,包括模型的建立、网格的划分、加载和求解等。
5.有限元分析软件的应用:介绍常见的有限元分析软件,如ANSYS、ABAQUS等,并通过实例教学,使学生能够熟练运用有限元分析软件进行简单的工程分析和计算。
三、教学方法为了实现本课程的教学目标,我们将采用多种教学方法,包括讲授法、案例分析法、实验法和讨论法等。
具体安排如下:1.讲授法:通过课堂讲授,使学生掌握有限元分析的基本概念、原理和方法。
2.案例分析法:通过分析实际工程案例,使学生了解有限元分析在工程实际中的应用和技巧。
3.实验法:通过上机实验,使学生能够熟练运用有限元分析软件进行简单的工程分析和计算。
4.讨论法:通过分组讨论和课堂讨论,激发学生的学习兴趣,培养学生的科学素养和创新能力。
四、教学资源为了支持本课程的教学内容和教学方法的实施,我们将准备以下教学资源:1.教材:选择合适的有限元分析教材,作为学生学习的主要参考资料。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
结构有限元分析程序设计绪论§0.1 开设“有限元程序设计”课程的意义和目的§0.2 课程特点§0.3 课程安排§0.4 课程要求§0.5 基本方法复习$0.1 意义和目的1.有限元数值分析技术本身要求工程设计研究人员掌握1). 有限元数值分析技术的完善标志着现代计算力学的真正成熟和实用化,已在各种力学中得到了广泛的应用。
比如:,已杨为工程结构分析中最得以收敛的技术手段,现代功用大致有:a). 现代结构论证。
对结构设计从内力,位移等方面进行优劣评定,从而进行结构优化设计。
b)可取代部份实验,局部实验+有限元分析,是现代工程设计研究方法的一大特点。
c)结构的各种功能分析(疲劳断裂,可靠性分析等)都以有限元分析工具作为核心的计算工具。
2). 有限元数值分析本身包括着理论+技术实现(本身功用所绝定的)有限元数值分析本身包括着泛函理论+分片插值函数+程序设计2. 有限元分析的技术实现(近十佘年的事)更依赖于计算机程序设计有限元分析的技术取得的巨大的成就,从某种意义上说,得益于计算机硬件技术的发展和程序设计技术的发展,这两者的依赖性在当代表现得更加突出。
(如可视化技术)3.从学习的角度,不仅要学习理论,而且要从程序设计设计角度对这些理论的技术实现有一个深入的了解,应当致力于掌握这些技术实现能力,从而开发它,发展它。
(理论本身还有待于进一步完美相应的程序设计必须去开发)4.程序设计不仅是实现有限元数值分析的工具和桥梁,而且在以下诸方面也有意义:1). 精通基本概念,深化理论认识;2). 锻炼实际工程分析,实际动手的能力;3). 获得以后工作中必备的工具。
(作业+老师给元素库)目的:通过讲述有限元程序设计的技术与技巧,便能达到自编自读的能力。
§0.2 课程特点总描述:理论+算法+数据结构(程序设计的意义)理论:有限元算法,构造,步骤,解的等外性,收敛性,稳定性,误差分析算法;指求解过程的技术方法,含两方面的含义;a. 有限元数值分析算法,b, 与数据结构有关的算法(总刚稀疏存贮,提取,节点优化编号等)数据结构:指各向量矩阵存贮管理与实现,辅助管理结构(指针,数据记录等)具体特点:理论性强:能量泛函理论+有限元构造算法+数据结构构造算法内容繁杂:理论方法+技术方法+技术技巧技巧性强:排序,管理结构(指针生成,整型运算等)§0.3 课程安排①. 单元刚度矩阵及元素设计(单元刚阵算法,杆梁平面分析,板弯非协调元等)②. 总刚的形式及程序设计(单刚提前准备,技术复杂)③. l边界条件及程序设计(等效荷载计算,位移边界条件置入,多工况的对称性)④. 总刚线性方程组求解(LDL T分解,分块算法,子结构算法,波前法)⑤.单元应力计算+应力处理与改善。
⑥. 数据处理(数据分类,压缩存贮,节点排序方法)⑦, 变带宽存贮的优化理论,图论的理论,有限元的图结构,存贮管理复核)⑧. 有限元议程全稀疏管理与求解策略。
说明:仅线性部份,复材,接触,弹塑性等不包括,基本部份。
实践性作业安排:1.作业:总的结构管理程序+子功能模块的编程,一个题的计算实践2.送有限元元素库。
§0.4 课程要求1.先行要求2.作业要求(计算机编程得出正确结果)3.课程校核要求§0.5 基本方法复习0.5.1 结构化程序设计方法0.5.2 有限元分析方法回顾0.5.3 Fortran语言回顾0.5.4 结构化程序设计1.基本结构:构成一个问题从输入到求解输出的基本程序形式:Imput——→Process——→Output(输入) (处理)(输出)三种基本形式:a). 树形结构(顺序执行结构)积木式(Fortran): 每个设定的功能分析团体的一个模块,每个模块又称作整体结构的素材,主模块象积木一样堆积.语体不联系,但有通讯方法沟通模块间信息,各种模块有各自的特征语体,main progame……,subroutine sub….嵌套式(Pascall): 主模块与子模块相互嵌套,各模块的特征,语体相同procedure mainprocefure ssub1……..procedure sub2……….End sub1………end sub2…….End main函数式(c语言): 主要特点是功能模块作为库函数调用,需用时在库内调用,每一个函数有表征语句,这种语体接近自然思维,而且对系统资源的调配应用更完善.面向对象的程序设计: 实施过程的可视化+控制性3.结构化程序设计方法a). Top—Down(自上面下),系统性强,选择性强.b). Critical Component First (关键部份优先),先抓主要矛盾,分清重,缓,急.c). 独立调试,总体联调,(软件设计的社会化作业).4.程序设计要点a). 自觉有意识地设计一个良好的程序结构,做到:易读,易懂,易管理,易修改,易发展.b). 做到逻辑清晰,说明完整.c). 要有工艺设计概念有框图,有步骤.5). 结构化程序设计原则a). 尽可能通用性好(适应各种规模的复题,?的扩大依据程序设计指标而定)b).整体精炼,清晰;避免GOTO。
c).省机时,省存贮,计算精度高,(算法上下功夫,要理论分析加技巧)d).输入数据少,格式简单。
e).输出结果简明,忌讳打印过多(与具体调试过程不一样)。
f).易读易维护,易发展。
§0.5.2.有限无方法求解过程回顾一.力学模型的分级管理有限无程序对力学模型的数据按一级:结构级(有点广义,不仅指具体结构,也指模型题目的规模)` 二级:单元级3. A). (1). 节点位置(总体坐标系下的坐标).(2). 节点局部坐标(按节点的约束方向制定的特殊坐标系x ',y ',z ',v 如斜支撑) (3). 节点的性质(自由,固定,指定位移,从属其它节点). (4). 节点力:(F x ,F y , F z ,M x ,M y ,M z ) (5). 节点位移:(u , v , w, θx ,θy ,θz ) B. 单元描述参数(1). 材料特性参数不清 E, G, γ→[D ](2). 节点的几何刚度参数(即面积A,板厚H,梁抗弯模量I) (3). 单元的局部坐标. (用于应力分析等,如图形曲面) (4). 单元的节点编号(5). 单元的几何矩阵营 (节点变形与应力关系矩阵) (6). 单元刚度矩阵 [K](7). 单元的应力,应变向量,(有限元分析多用向量,而不用矩阵(张量))结构描述参数单元总数,节点总数,单元娄型总数,结构材料种娄数,节点自由度数(控制题目规模) 二. 基本公式系统1. 单元刚度计算公式2. 单元刚阵组合 [K]=ΣA T KA3. 单元节点荷载计算4. 节点荷载组装:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛333231232221131211231312332211σσσσσσσσσσσσσσσ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛33323123222113121123131233221121212121212121εεεεεεεεεεεεεεε⎰=cV T e dvB D B K ]][[][][⎰=eS T e dsP N P }{][}{()∑=eAP P5. 位移约束关系:6. 总刚方程解:7. 应变计算:8. 应力计算:9. 支撑反力计算:三. 有限元分析的模块组织.{}⎭⎬⎫⎩⎨⎧=-d d []{}{}{}[]{}P K P K e =⇒=δδ{}[]{}δεB ={}[]{}[][]{}δεσB D D ==iii i d K R -=四.结构分析的原始输入数据1.题目规模` 节点数目:NNP 单元数目:NE2.节点数据单元人坐标:XE(NNP,3)3.单元数据单元节点编号:ME(NE,3)、ME(NE,2)材料特性:E、N单元几何参数:I、RI(惯矩)4.荷载数据外荷载作用点,坐标及大小:PA(NNP,1)§0. 5.3 Fortran语言回顾1.子模块(子程序)subroutinea.特点:独立性强,只要输入输出接口,象一个黑匣子,与外界无关。
b.作用:完成一个独立的功能(求应力,矩阵分解,投放等)c.格式:subroutine function(ip1,ip2,rp1,rp2,io1,io2,ro1,ro2).(其中ip1,ip2,rp1,rp2,是输入形参,io1,io2.rp1,rp2是输出形参)2.数据传递形式1). COMMON 公共块语句传递,(公共块的内容不能作为形参)a.公共块分为无名公共块和有名公共块b.公共块的参数不能作为子程序的参数出现,c.公共块名一致,其内容在不同公共块中可以标志符不同(但其长度应一致)d.通用原始数据放入公共块(作为实参错误率大)e.尽可能不放数值,安息组一般可作成可调长度f. 格式Common/comm/…….Subroutine fun()Common/comm/……..2).形参————实参对应a.实参不能开辟存贮单元,子程序内定义语句中的形参数组由主程序定义,在子程序中仅形式定义(即仅说明是数组,因而大小无所谓)b.格式:Dimension RP(1000),RO(1000)…….Call sub1(RP,RO)………ENDDIMENTION IBANKSUBROUTINE SUB1(RP,RO,NE)DIMENTION RP(1),RO(1),SP(50)DIMENTION RP(NE,1),RO(NE,1) (形参的动态定义,实参不能)3). 数组长度自动调整方法。
PROGRAM MAININPLICIT REAL*8 (A—H,O---Z)CHARACTAR*20 TRCOMMON/COMM/….DIMENTION IBANK( ),RBANK( ),IP1( ),IP2( )IP1(1)=… IP(N)=…..IP2(1)=…… IP2(N)=…..CALL SUB1(IBANK(IP1(1)),IBANK(IP1(N)),RBANK(IP2((1)),….)…..ENDSUBROUTINE SUB1(II1,IO2,….RI1,….RO1….NE) DIMENTION II1(1) ,IO2(NE,1), RI1(1), RO1(1)。