主程序及子程序功能.txt

这是一个用fortran编写的“第一类边界条件的有限元算法”

!***************************************************************************!                   
!                               主程序功能：                                !
!                                                                           !
!         求解第一类边界条件、三角元剖分、线性插值的位场延拓                !
!***************************************************************************!
          
!***************************************************************************!                             
!                            子程序功能说明：                               !
  
1、 SUBROUTINE UKE1(x,y,ke)
 
!   功能     三角单元、线性插值时计算单元系数矩阵ke。                       !
!   参数说明                                                                !
!     x----3个元的一维实数组,输入参数,存放单元节点的x坐标,存放顺序xi,xj,xm. !
!     y----3个元的一维实数组,输入参数,存放单元节点的y坐标,存放顺序yi,yj,ym. !
!     ke---3*3的二维实数组,输出参数,存放单元系数矩阵.                       !

2、 SUBROUTINE MBW(ne,i3,iw)

!   功能    计算总体系数矩阵的半带宽。                                      !
!   参数说明                                                                !
!     ne---整型变量,输入参数,单元总数。                                     !
!     i3---3*ne的二维整数组,输入参数,存放单元节点编号。                     !
!     iw---整型变量,输出参数,半带宽。                                       !

3、 SUBROUTINE UK1(nd,ne,iw,i3,xy,sk)

!   功能     三角单元,线性插值时用定带宽储存的方法集成总体矩阵              !
!   参数说明                                                                !
!     nd---整型变量,输入参数,节点总数。                                     !
!     ne---整型变量,输入参数,单元总数。                                     !
!     iw---整型变量,输入参数,半带宽。                                       !
!     i3---3*ne的二维整数组,输入参数,存放单元节点编号,存放顺序:i1,j1,m1,    !
!          ...,ine,jne,mne。                                                !
!     xy---2*nd的二维整数组,输入参数,存放节点的xy坐标,存放顺序:x1,y1,       !
!          ...,xnd,ynd。                                                    !
!     sk---nd*iw的二维整数组,输出参数,存放定带宽储存的总体系数矩阵          !

4、  SUBROUTINE UB1(nd1,nb1,u1,nd,iw,sk,u)
!    功能    在定带宽储存的总体系数矩阵和右侧列向量上加上第一类边界条件。   !
!    参数说明                                                               !
!      nd1---整型变量,输入参数,第一类边界条件的节点数。                     !
!      nb1---nd1个元的一维整数组，输入参数，存放第一类边界条件的节点号。    !
!      u1----nd1个元的一维实数组，输入参数，存放第一类边界节点的场值。      !
!      nd----整型变量,输入参数,节点总数。                                   !
!      iw----整型变量,输入参数,半带宽。                                     !
!      sk----nd*iw的二维整数组,输入、输出参数,定带宽存放总体系数矩阵，输出  !
!                时,第一类边界条件已带入。                                  !
!      u-----nd个元的一维实数组，输出参数，存放加入第一类边界条件后的右端列 !
!                向量。                                                     ! 

5、    SUBROUTINE LDLT(a,n,iw,p,ie)
!      功能    对称带型线性方程组的系数矩阵的下三角部分被定带宽地储存在矩型 !
!              数组a中，利用a数组来解方程组。                               ! 
!      参数说明                                                             ! 
!        a-----n*iw的二维实数组,输入参数,存放对称带型线性方程组的系数矩阵的 !
!              下三角部分。                                                 !
!        n-----整型变量,输入参数,方程组的阶数。                             !
!        iw----整型变量,输入参数,对称带型线性方程组的半带宽。               !
!        p-----n个元素的一维实数组，输入、输出参数。开始存放方程组的右侧列  !
!              向量，工作结束时，存放解向量。                               !
!        ie----整型变量，输出参数，标志。ie=0时表示程序正常结束；           !
!              ie=1时，表示系数矩阵奇异。                                   !
!***************************************************************************!