C语言与有限元程序设计

合集下载

有限元基础编程百科全书

有限元基础编程百科全书
有限元基础编程是指使用有限元方法（FEM）进行工程分析和设计的计算机编程技术。

有限元方法是一种数值分析技术，用于解决复杂的工程和物理问题，例如结构分析、热传导、流体力学等。

下面我将从多个角度全面介绍有限元基础编程的百科全书。

首先，有限元基础编程百科全书应该包括对有限元方法的基本原理和数学基础的详细讲解。

这包括有限元离散化过程、单元和节点的概念、刚度矩阵和质量矩阵的推导，以及有限元解的数值求解技术等内容。

此外，对于常见的工程问题，如静力学、动力学、热传导和流体力学等，百科全书还应该包括有限元方法在这些领域的应用原理和算法。

其次，有限元基础编程百科全书还应该涵盖有限元程序的编写和实现。

这包括使用常见的有限元软件（如Abaqus、Ansys、Nastran等）进行编程的基本步骤和技巧，以及各种编程语言（如Fortran、C++、Python等）在有限元分析中的应用。

此外，还应该包括有限元程序的优化和并行计算技术，以提高计算效率和精度。

此外，有限元基础编程百科全书还应该介绍有限元分析在工程
实践中的应用。

这包括结构分析、材料力学、振动和声学分析、热传导和传热分析、流体力学和多物理场耦合分析等方面的工程案例和实际应用。

这些案例不仅可以帮助读者理解有限元方法的具体应用，还可以帮助他们将理论知识转化为实际工程问题的解决方案。

综上所述，有限元基础编程百科全书应该全面介绍有限元方法的理论基础、编程实现和工程应用，从而帮助读者全面深入地理解和掌握有限元分析技术。

希望这些信息对你有所帮助。

有限元教材PFF程序使用说明C版

1.8 平面刚架程序设计1.8.1 程序说明及总流程图依据前述平面刚架矩阵位移法原理，本节使用Visual C++ 6.0编译器编制了平面刚架先处理法计算程序PFF，该程序未使用C++面向对象的编程方法，而仍沿用兼容C的结构化程序设计方法，以便和本科生先修课程C语言相容。

本程序不建议使用TC++3.0及更低版本的编译器编译。

一．PFF的主要功能和特点(1)输入单元编号、结点编号和结点位移分量统一编码，单元材料信息；(2)先处理法形成结构刚度矩阵和综合结点荷载列阵；(3)Gauss消元法解线性代数方程组；(4)计算并输出结点位移和单元杆端力；二．PFF的使用方法(1)输入输出数据文件PFF从原始数据文件中读取结构的离散化信息，经计算后，将结果输出到结果数据文件中。

运行编译好的PFF可执行程序，会提示“Please input primary data file name!”，这时输入事先准备好的原始数据文件名，回车后，提示“Please input calculation result file name!”，输入一个文件名，PFF会自动以该文件名生成结果数据文件。

原始数据文件最好同PFF可执行程序放置在同一文件夹中，结果数据文件不应是已有的文件，而必须由PFF自动生成。

两文件建议使用英文和数字命名。

在输入文件名时，应包含后缀。

C语言以0为数组各维的开始下标，称为0基数组；而一般习惯是以1为数组各维的开始下标，称为1基数组。

PFF源程序编制时，仍以C语言为标准，使用0基数组，即PFF 程序内部处理机制是0基的。

但是为方便读者，原始数据文件和输出数据文件则按一般习惯，使用1基数组对结点、单元、结点位移分量、结点荷载和非结点荷载进行编号，即PFF程序外部表现是1基的。

(2)原始数据文件的格式下面依行序说明原始数据文件输入数据的含义：①总体信息ne,nj,n,np,nf,ene——单元总数；nj——结点总数；n——结点位移未知量总数；np——结点荷载总数；nf——非结点荷载总数；e——弹性模量②结点坐标(x[i],y[i])(i=0;i<nj;i++)x[i]——结点i+1的x坐标；y[i]——结点i+1的y坐标。

本科专业认证《计算机与程序设计基础》(C语言)教学大纲

《计算机与程序设计基础》课程教学大纲英文名称：Fundamentals of Computer and Programming课程编号：0812010005学分：2.5总学时：40。

其中，讲授学时24，上机学时16。

适用专业：全校非计算机专业理工科本科生先修课程：无开课学院、系：计通学院计算机应用系一、课程目标《计算机与程序设计基础》课程是新生入校的第一门计算机课程，也是大学本科理工科专业学生必修的公共基础课程。

《计算机与程序设计基础》的主要教学目的是使学生较全面、系统地掌握计算机的基本知识，理解计算机操作系统基础知识、掌握操作系统的基本应用，了解办公自动化系统的概念、基本组成和关键技术，熟练掌握文字处理软件、电子表格软件、演示文稿制作软件的使用方法，掌握C语言的基本语法、编程技术和基本算法，掌握程序设计的基本思想和方法，具备利用计算机求解实际问题的基本能力，能灵活应用C语言进行程序设计，有一定的编程和调试程序的能力。

使学生具有利用计算机分析问题、解决问题的意识与能力，提高学生的计算机素养，为将来应用计算机知识和技能解决专业中的实际问题打下必要的基础。

按照国家提出的课程培养目标、以及开设了本课程的相关专业2017培养方案中的毕业要求，制定本课程学习目标如表1所示：表1 《计算机与程序设计基础》（C语言）课程目标-毕业要求关系表二、课程内容及学时分配本课程内容、建议学时以及知识单元与课程目标支撑关系如表2所示。

表2 《计算机与程序设计基础》（C语言）课程内容及学时分配三、教学方法根据所面向学生的不同特点和专业要求，针对不同的教学内容采用不同的组织方式。

对理论性较强的内容以课堂教学为主；对应用性强的部分精讲多练，加强课内实验和课外自主练习；对较复杂的操作性内容采用任务驱动的案例教学；对知识性方面的内容以学生自学为主，同时让学生充分利用网络资源、教育平台进行启发式、讨论式、研究式的自主学习，在教学中给学生提供多样化的探索空间，鼓励他们进行个性化发展。

有限元程序设计大作业

有限元程序设计大作业1.不同板宽的孔边应力集中问题姓名：胡宇学号：21201201282.摘要本文采用MATLAB和FOTRAN四节点平面单元，利用有限元数值解法对不同板宽的孔边应力集中问题进行了数值模拟研究。

对于不同的板宽，并且与解析系数（半板宽b/孔半径r），得到了不同的应力集中系数1解进行了比较，验证了有限元解的正确性，并且得出了解析解的适用范围。

3.引言通常情况下的有限元分析过程是运用可视化分析软件(如ANSYS、ABAQUS、SAP等)进行前处理和后处理，而中间的计算部分一般采用自己编制的程序来运算。

具有较强数值计算和处理能力的Fortran语言是传统有限元计算的首选语言。

随着有限元技术的逐步成熟,它被应用在越来越复杂的问题处理中。

MATLAB是由美国MATHWORKS公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。

它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。

4.MATLAB部分1，计算模型本程序采用MATLAB编程，编制平面四边形四节点等参元程序，用以求解近似平面结构问题。

本程序的研究对象为中央开有小孔的长方形板，选取的材料参数为：板厚h=1、材料强度E=1.0e11 Pa、泊松比mu=0.3。

此外，为方便网格的划分和计算，本文所取板的长度与宽度相等。

其孔半径为r=1，板宽为2b待定。

由于本程序的目的在于验证有限元解的正确性和确定解析解的适用范围，因此要求网格足够细密，以满足程序的精度要求。

同时为了减小计算量，我采取网格径向长度递增的网格划分方法。

此种方法特点是，靠近小孔部分的网格细密，在远离小孔的过程中，网格逐渐变得稀疏。

《有限元及程程序设计》课程设计说明书

河北工程大学《有限元及程序设计》课程设计说明书课程设计题目：平面钢架有限元程序功能修改设计副标题：1增加非节点荷载为集中力偶、纵向均布力以及线性分布力时固端反力的计算功能2 结果以文件形式输出指导教师：班级：姓名：学号：摘要有限元法是现代工程数值分析中应用广泛的一种方法，本文根据线性有限元理论对受有五种不同荷载的三杆钢架进行静力分析，将结构离散为三个单元四个节点。

分别建立每个单元的单元刚度矩阵和节点荷载列阵，根据刚度集成法建立了结构的总刚度矩阵和节点荷载列阵，得出结构的平衡方程，并用对角线元素置一法引入边界条件，用高斯消元法求解平衡方程。

最后编写C程序求解此问题，并通过与手算的比较验证了程序的准确性，通过增加一个杆件说明程序的通用性键词：有限元法杆平面刚架刚度矩阵对角线元素置一法 C语言A b s t r a c tFinite Element Method (FEM) makes an extensive use in the numerical analysis of modern construction. In this paper, we study Static state of the Plane frame which is loaded by five different kinds of loads based on linear finite element theory loads; the structure is divided into three units and four nodes. The stiffness matrix array node load of each unit is set up, According to Stiffness integration method we establish the total stiffness matrix and load node array of the structure which is aimed drawing the balance equation of the structure, and then introduce the boundary conditions by buy-one-diagonal elements .we can use Gaussian elimination method for solving equilibrium equations. Finally, we program procedures for the preparation of C to solve this problem, and through comparison with the hand count to verify the accuracy of the procedure, by adding a bar to descript generic property of the procedureKey words: Finite Element Method (FEM) Planar rigid frame Buy-one-diagonal elements Stiffness matrices C language目录设计题目说明-------------------------------------------------4 1.用有限元法进行手算----------------------------------------41.1化分单元，选取坐标系----------------------------------------------4 1.2求局部坐标系下各单元的单元刚度矩阵--------------------------------5 1.3求整体坐标系下各单元的单元刚度矩阵--------------------------------5 1.4求整体刚度矩阵----------------------------------------------------6 1.5求非节点荷载引起的等效节点荷载及节点荷载列阵----------------------6 1.6列整刚方程，求节点位移--------------------------------------------9 1.7求单元内力----- --------------------------------------------------101.7.1转换位移列阵- ---------------------------------------------------101.7.2求内力- --------------------------------------------------------111.8画内力图并列表--------------------------------------------132.程序设计与上机调试结果-------------------------------------142.1说明与结果- ------------------------------------------------------14 2.2程序设计中一些问题的描述- ----------------------------------------142.2.1数字描述--------------------------------------------------------142.2.2程序总框图------------------------------------------------------16主要结论-----------------------------------------------------16 设计心得体会---------------------------------------------------------16参考文献-----------------------------------------------------17附录C程序源代码及修改注释- --------------------------------17求图示平面刚架节点位移及各杆的内力错误!具体参数：面积弹性模量E=27/101.2mKN ⨯惯性矩I=421016667.4m-⨯1用有限元法进行手算1.1划分单元，标出单元号码及节点号码；选取整体坐标系O x y ，局部坐标系Oxy ，并标上单元的局部节点码i(1),j(2)，见下页图。

c++面向对象的有限元程序设计

《C++面向对象的有限元程序设计》一、引言在计算机科学和工程中，有限元方法是一种数值分析技术，广泛应用于工程设计和科学研究领域。

C++作为一种流行的编程语言，在有限元程序设计中也扮演了重要角色。

本文将从深度和广度两个方面对C++面向对象的有限元程序设计进行全面评估，并撰写一篇有价值的文章，以帮助读者更全面、深刻地理解这一主题。

二、C++面向对象的有限元程序设计的基本概念1. 有限元方法的基本原理有限元方法是一种数值计算方法，用于求解偏微分方程和积分方程。

通过将求解区域分割为有限个单元，建立单元之间的联系，将连续的问题转化为离散的代数问题，从而得到数值解。

在有限元程序设计中，需要考虑如何有效地表示和处理单元、节点、边界条件等信息。

2. 面向对象的程序设计思想面向对象的程序设计思想强调将现实世界中的问题抽象成对象，通过封装、继承和多态等机制构建模块化、可复用的代码结构。

在C++中，类和对象是面向对象程序设计的核心概念，有限元程序设计可以通过抽象出单元、节点、网格等对象来实现。

三、深入探讨C++面向对象的有限元程序设计1. C++语言特性在有限元程序设计中的应用在C++语言中，有丰富的特性可以用于实现面向对象的有限元程序设计。

类的封装可以用于表示单元和节点对象的属性和行为，继承可以用于构建具体单元类型的层次结构，多态可以实现对不同单元类型的统一处理。

2. 优化设计思路下的C++面向对象有限元程序设计针对大规模的有限元计算，优化的设计思路是必不可少的。

C++中提供了丰富的性能优化手段，如模板元编程、内联函数、移动语义等，可以在面向对象的有限元程序设计中发挥重要作用。

四、总结和回顾在本文中，我们对C++面向对象的有限元程序设计进行了全面评估，并撰写了一篇有价值的文章。

通过深入探讨原理、语言特性和优化设计思路，帮助读者更全面地理解了这一主题。

从我的个人观点看，C++面向对象的有限元程序设计是一个值得深入研究的领域，它不仅涉及到程序设计技术，还涉及到数值计算和工程应用等多个领域的知识。

有限元单元法程序设计

有限元单元法程序设计有限元单元法是一种用于工程结构分析和设计的计算方法，它将大型结构分解为许多小的离散单元，通过分析单元之间的相互作用来预测结构的力学行为。

有限元单元法程序设计是指针对特定工程问题，编写计算机程序来实现有限元分析的过程。

下面将介绍有限元单元法程序设计的基本流程和关键要点。

一、问题建模和网格划分有限元单元法程序设计的第一步是对工程结构进行合理的建模和网格划分。

建模的目的是将实际结构抽象为适用于有限元分析的数学模型，包括定义结构的几何特征、材料属性、边界条件等。

网格划分是将结构分解为许多小的单元，每个单元具有一定的形状和尺寸，以便于数值计算。

常用的单元形状包括三角形、四边形、四面体、六面体等，根据结构的特点选择合适的单元形状和尺寸。

二、单元刚度矩阵和载荷矩阵的求解在有限元单元法程序设计中，需要编写算法来求解每个单元的刚度矩阵和载荷矩阵。

单元刚度矩阵描述了单元内部的力学性能，包括刚度、弹性模量、泊松比等，它们通常通过数学公式或有限元理论推导得到。

载荷矩阵描述了单元受到的外部荷载，可以是均匀分布载荷、集中载荷或者边界条件引起的约束力。

通过合适的数值积分方法，可以计算得到每个单元的刚度矩阵和载荷矩阵。

三、组装全局刚度矩阵和载荷向量在有限元单元法程序设计中，需要将所有单元的刚度矩阵和载荷向量组装成整个结构的全局刚度矩阵和载荷向量。

这涉及到单元之间的连接关系以及边界条件的处理。

采用适当的组装算法，可以将各个单元的刚度矩阵和载荷向量叠加在一起，形成整个结构的刚度矩阵和载荷向量。

四、求解位移和应力有限元单元法程序设计的最后一步是求解结构的位移和应力。

通过斯蒂芬-泰勒算法或者其他迭代算法，可以得到整个结构的位移分布，然后根据位移场计算各个点的应变和应力。

这一过程涉及到对整个结构刚度矩阵的求解和对位移的后处理。

有限元单元法程序设计是一个复杂而又精密的工作，需要深入理解有限元原理、结构力学知识和数学方法。

有限元分析程序设计

结构有限元分析程序设计绪论§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、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

Ｃ语言与有限元程序设计
叶又戚燕
（上海交通大攀国家棋其ＣＤ工程中心上海加０３）几００
（上海大学计算中心上海匆印力
摘要本文通过与ＦＲＲＮ语官的比说明ＯＴＡ７７较，利用Ｃ语言开发有限元程序的可行性。结果表明，语官不但完Ｃ全可以实现ＦＲＲＮ的语法功能，ＯＴＡ７７拥有更强大的实用编程环境，同时由于提供指针变ｆ、动态内存分配函数和结构变ｔ，使其编写的软件在维护性、可读性和内存利用效率方面具有明显的优势。文章最后还将Ｃ语言与近出最现的ＦＲＲＮ作了ＯＴＡ９０简单比较。
２语言与ＦＲＲＮ的主要语法结构比较Ａ７ＣＯＴ７
ｃ最初出现是用作ＵＩ操作系统的ＮＸ记述语言，发展之初一直作为“ 较低级语言” 编写与硬件关系密切的底层软件。由于ＵＩ成功和ＮＸ的广泛使用，语言开始流行起来，ｃ并成为一种普遍使用的程序设计语言。由于两者产生的背景不同，它们是存在差异的。主要表现在函数的调用和参数的传递，特别是多维数组的处理方面不如ＦＲＲＮ直观和易于理解。ＦＲＲＮＯＴＡＯＴＡ从一开始就用于科学计算，必须处理诸如矩阵一类的数学问题，多维数组是它的主要使用工具，语言一开而ｃ始主要是面向系统程序员的，在应用中很少会涉及多维数组。但这并不盒味着在处理多维数组时会遇到很多麻烦。３提供了利用指针数据类型处理多维数组及其文【〕在函数中传递的两种方法。利用这些方法可以方便地将ＦＲＲＮ的数组操作移植成ｃＯＴＡ语言
和数据类型，甚至允许用户利用基本数据类型自己定义数据类型。
ＦＲＲＮ的主要可执行语句与ｃ确的７ＯＴＡ７有明语句对应关系，并完成相应的功能，见表１，
袅１ＯＴＡ７主要句法与ｃ语官的对应关系ＦＲＲＮ７
ＣｎｒｇＦＲＲＮｔｓｒｗｔｓｔｉｙｇｎａｉｆｉａｐｉｗｈＴＡ７，ａｒｉｓｕａｌｕｎＣｕｅｎｅａｎｉＯｔ７ｈｐｅｅｅｈｉｉｓｌｇｇｎｉｐｖｅｂｔｉａｉｔ二ｄｅｐｅ．ｒｕｉｉｔ．ｔｎａｎｏｙｐｔａｉｅｔａｅｌｍｎＴｅｌｎｃｅｔＣｇｇｏｎｃｌｌｃｅｓｌｖｏｔｈｅｔｅｓｈｌｕｅｌｏｅｙｖｈｓｄａａｔｍｅｈｅｏｕｒ６ｌｒｍｇｗｒｌｅｒａｎ－ａｏａｎｄＦＲＲＮ，ａｏｐｖｅＯＴＡ７ｂｔｃｒｉ二ｐｆｐｃａｐｇｍｉ－－，７ｕｌａｏｄｓｎｄａｉｍｍｒａｏｔａｔｓｓｃｒｄａｓ咫ｐｇｍｉ阮ｙｍｅｏｃｍｆｃｎａｔｔｅａｔｅ阮ｎｃｙｌａｕｉｎｒｕｄｙ，ｎｏｄｕｔｐｒｒｏＣｏａｆｓ
ｉｒｍ：ｈ勺ｒｄｉ耐ｍｍａｉｉ．ｓ如ｃｐｓｎｔｎ，ａｌｎｎｎ山ｅ曲ｉａａｔｙｅｏｌｌＩｔｏａｏｙａｂｔｒｖａｉａ，ｍｒｎｙｉ
ｉｓｐｄｃｓ．ｓｌｉｕｅｉｙｓｄｍｓＦｉＥｍｎＰｇｍＯＴＡ７ＦＲＲＮｉｔｌｅｒｒ７ＴＡ９ｎｅｔａＦＲＲＮＯｅｏ０
Ｓｏ－１ａｔｗ－
．２１
万方数据
３动态内存分配技术
让我们先来研究一下利用ＦＲＲＮ编写的ＯＩＡ７７有限元程序所使用最广的一种内存动态管理方法。该方法在程序一开始就在无名公共块中定义一个大数组Ａ在程序不同，运行阶段存
放不同的内容。
ＣＲ１〔ＲＭＭ九ＮＡＬ
收期：一加。叶又，主要研究领城：维板稿日ｔ仍一＂博士，三料成形有限元棋拟系开发。数值统的
·２ · ０
万方数据
在各种机型和各种操作系统上都运行有ｃ语言的编译器，提供了丰富的库函数和实用程序，以改善用户一程序界面，ＯＵ效率，并简化编程。本文的主要目的是通过与ＦＲＲＮ的比较，ＯＴＡ７７说明利用ｃ语言开发有限元程序的可行性。ｃ语言不但包容了ＦＲＲＮ的全部语法功能，ＯＴＡ７７它提供的指针数据类型实现了真正意义上的内存动态分配。同合理有效地定义时，结构数据类型可大大增加程序可读性和维护效率。ＦＲＲＮ是ＦＲＲＮ０ＯＴＡ语言的最新标准版本。由ＯＴＡ９于采用了新的计算机软件科学理论，其功能实际上已超过ｃ语言。作者对两者作了概要的比从实用的角度提供了较，一些建议。
Ｎ６二褚＋ＮＵＰ－裸刃ＤＮ
ＣＬ１＂ＡＮ）ＡＮ）ＡＮ）ＡＮ）ＡＮ）ＡＮ）ＮＯ００ ቤተ መጻሕፍቲ ባይዱＬＮ（１，（２，（３，（４，（５，（６，ＤＦ · ）
ＲＵＲＥＴＮ
ＥＤＮ
ＩＰＮＳＰＮＮ，ＦＸＹＺＮＯ，．ＮＩＵＮ，ＰＭＯ，，ＤＦ．）ＭＰＤ，，．ＩＰａＲＡ－一，一）ＭＩＴＬＳＡＨ０ＺｄＥ（ＤＡ．ＳＮＳＰ１浏ＰＰ１，ＤＦＮＯ，ＩＭＩＮＮ（）Ｎ（）ＭＯ（ＤＦ１ＥＭＯＵ）
关扭词Ｃ语言有限元程序ＦＲＲＮＯＴＡ９ＯＴＡ７ＦＲＲＮ７０
ＣＮＵＧＡＤＮＴＦ．ｌＬＧＡＥＮＦＩＥＰｙＡＩｆＥＮｊＰＯＧＡ＇ＲＭＭＭＧＲ
ＹｏＱＹｒｅｕＹ，ｉａｎ
１引
言
ＦＲＲＮＯＴＡ是最早出现的高级编程语言之一。国际标准化组织先后推出了ＦＲＲＮＯＴＡＷ，ＦＲＲＮ以ＯＴＡ７及最近的ＦＲＲＮ三种标准版本，７ＯＴＡ９０至今已发展了三十余年，在科学计算领域有广泛的应用。当前有限元程序的抓代码几乎清一色是由ＦＲＲＮＯＴＡ写成。而Ｃ语言的发展仅有十几年，但却有一种后来居上的趋势，开始成为开发高质ｔ软件的主要编程语盲。目前
ＦＯＲＴＡＮ７Ｒ７
ｃ
二二二（》，。
ｆ．）｝Ｖ（：１．７
曰只份
｝一
ｅｄ．ｌ（）ｓ｛ｅ
ｉ．）ｆ．ｆ（．
日．｝｝
ＥＤＮＤ＇
ＤＩ，，Ｏ６ｃ＝ｎＤ二ｂＯＪａ ’
ＥＤＮＤＯ
Ｉ
。，，二＊。（．＜、十）二１．二｛
ＤＥＳＮ）Ｙ１，１ＩＮＯＸ１，（）（）ＭＩ（Ｚ
由于Ａ数组的大小必须事先给定，为满足解题规模，会尽可能开得足够大。当求解一个较小间题时，它同样占据着很大的内无疑造成计算机资源的浪费，存，这种情况对于目前多任
·２．２
万方数据
务多进程的操作环境更为重要。同时，求解更大规模的问题，为了还必须调整Ａ数组的大小，重新编译源程序。实际计算问题的大小是很难确定的，利用上面方法编写的程序实际上在生成执行代码时就决定需要内存的数龙。更一般的情况是，如果一个变量被指定为全局变量（如上面的Ａ数组）则它在整个程序运行期间都占，据存储单元。如果是自动变盆和形参，则在其所在的子程序执行期间，有的它占存储单元也是不释放的。这种方法处理问缺乏灵活性，题，往往会浪费许多内存，不是真正意义上的动态内存分配。人们应该设想能否找到这样一种方法，根据需要临时分配内存单元以存放有用的数据，数据不用时又可以释放存储单元。当此后这些存储单元又可用来分配给其它数据或程序使用，语言使用指针数据类型和标准的动态内存管理函ｃ数ｃｏ（可以非常方便地实现这一想法。ａｃ）Ｕ帕，ａｃｎｓ）记ｃｌ（，；ｌｏｉｚｅ分配ｎ个变Ａ的内存连续空间，每个变２字节大小是ｓ。由于此函ｉｚｅ数类型说明时是ｖｄｏ＊ｉ，故在使用时要强制转换为所用类型变Ａ的。指针·
ＮＮ－ＵＰ一一ＮＭＢＲＮＤＰＭＵＥＯＯＥＦＯ
ＮＯ－ＤＦ
ＮＯＳＤ－
ＮＭＵＭＲＤ！ＵＯＦＭＰＲＤＰＩＯｔＩ！Ｆ＂ＳＮＥＭＦＪ－＆ｍＥＯＯ
ＳＩＲＦＷｒＣＯＳＭＥＲＵＬＰＥＩＩＮＯＤＢＥＣＳＯＯＲ
气Ｎｒｏｆ（Ｄｉａ｝ｒＡ＆．ｉＤｅｗ｝Ｃｒ，ｈＪａｍｙＳｙ．０ｏｅ９ｗｍｗａｎ－ｉａ，ｈ２００Ｂｍ０３权Ｃ哪吨ｈｔ＇匆，－ｍｗｉ：，州２Ｍ） “ Ｃ．｝Ｓｈｍ１Ｖｍｔｏ
ｉ＊ｎＰ；ｔＩ
ｌｔｐ；ｆａ赞２ｏ
ｐ＝ｉ＊ｏ（ｌ；）１（ｔ）ｌｃ，ｎｃｎ２ａｐ二ｆａ）ｆ（４，；２（ｏ＊ｏｎ）ｌｔｃｃ２ａｌ
ｎ，的小根解的棋程运时值数Ｉ用ｎｘ个节内声共ｌ２大可据题规在序行喊。组Ｐ共占ｌ２宇存，ｎ占２字内。数ｌ２再要，们用空可单使下语用ｎｘ个节存当组ｐｐ不孺时它占的间简地用面句４，
的。
与ＦＲＲ７Ｎ一样，也是一种结构化的程序设计语言。结构化程序设计要求程序的ＯＴＡ７ｃ逻辑结构有顺序、选择和循环三种基本结构，提供了编写结构化程序所需要的语句。同时，ｃ使用ｃ也便于进行模块化程序设计，程序由ｃ众多的函数组成，数是进行模块化程序设计的函基本单位，ＯＴＡ与ＦＲＲＮ的子程序功能相似。另外，还提供了ＯＴＡ丰富ｃ较ＦＲＲＮ得多的运算符