等几何分析(IGA,Isogeometric Analysis)于2005年被提出,旨在将设计、分析和优化集成在一起。等几何分析利用准确描述几何形状的基函数(如NURBS, Non Uniform Rational B-spline)作为计算解空间的函数,统一的表达形式使得分析模型和几何模型可以无障碍地交互。随着工程实际中模型的日益复杂化,有限元分析中模型的网格生成耗费了大量的时间,而等几何分析将计算机辅助几何设计(CAD,Computer Aided Design)与计算机辅助工程分析(CAE,Computer Aided Engineering)进行了无缝结合,很好地解决了这一问题。等几何分析已成为当前工程分析发展的趋势,将会对CAD和CAE产生重大影响。现有的有限元环境很多基于Fortran语言,利用Fortran对等几何分析进行编程有希望实现等几何分析与现有有限元分析的集成。作者利用Fortran语言对NURBS等几何分析的理论以及实现方法进行了研究,并将基于Fortran的等几何分析用于不同的力学问题以提高等几何分析的计算效率和准确性。为了提高程序的运行效率,同时也为了使程序能够方便地被不同的算例调用,文中采用Fortran语言中广泛使用的模块化编程。为了使Fortran编程能对不同的算例兼容,作者最大化的将共用子程序,如高斯积分,NURBS基函数及其导数的计算,连接矩阵等,编入同一模块,对Fortran用于等几何分析的研究做了大量的基础性工作。另外,现有的等几何分析软件在解决大型稀疏矩阵问题时,计算效率低且不稳定,作者将一种大型稀疏矩阵的求解器GSS(Grus Sparse Solver)植入Fortran编程中,使得等几何分析的效率大大提高。文中提供可下载的相关算例的Fortran源程序。NURBS基函数通常不满足克罗内克函数的性质,即不具有插值性。控制点不在边界上时,难以直接施加位移边界条件。作者采用罚函数法处理边界条件,但罚因子取值的大小会影响计算结果。作者通过不同的算例探索合适的罚因子取值大小,结果表明,罚因子的取值在高于整体刚度矩阵中绝对值最大值2-3个数量级时,计算结果较为准确。此外,提高模型的网格质量可以减轻罚因子的影响,但高质量的网格势必会导致计算工作量的增加。作者完善了基于罚函数法处理边界条件的Kirchhoff-Love板壳理论。NURBS函数不仅用来描述单元形状和位移场,而且还提供了Kirchhoff-Love理论所需要的高阶连续函数。将基于Fortran的等几何分析用于不同的板壳算例模拟,证明了其准确性和快速收敛性,同时说明了等几何分析即使在粗糙的网格水平上,也能得到准确的计算结果。作者还结合扩展有限元与等几何分析对断裂力学问题进行了研究,证明了基于Fortran的扩展等几何分析模拟不连续问题的有效性,给出了选择强化控制点的方法并对裂纹不连续域和尖端位移场分别利用Heavisde方程和裂纹尖端方程进行了强化。通过与扩展有限元用于模拟相同的裂纹模型作比较,说明了扩展等几何分析仅需更少的单元便可获得准确的结果。在对不同算例位移和应力场的表达上,高阶函数的应用使得其光滑、连续。为了实现局部细化,作者利用多面片技术对带孔平板问题模型进行了分片处理,并在每一片模型上进行了位移和应力计算,得到了准确的结果。在模拟断裂力学问题时,利用线性节点值的插入对裂纹区域进行了局部细化,减少了计算误差。根据不同的算例应用不同程度的细化方法所获得的结果来看,提高网格的细化质量,可以使计算结果更加准确,位移和应力场表达的更加连续。综上,作者开发了基于Fortran的新的等几何分析工具,并将其用于不同的力学问题以验证其准确性和高效性。同时,完善了等几何分析中的相关理论,对等几何分析中出现的问题进行了分析,并提出合理的改进方法。然而由于时间和条件有限,还需对程序进行优化,也需要将Fortran编程应用于更多的等几何分析算法。