uke1.f90

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

!***************************************************************************!
!                       subroutine UKE1(x,y,ke)                             !
!   功能   三角单元、线性插值时计算单元系数矩阵ke。                         !
!   参数说明                                                                !
!    x----3个元的一维实数组,输入参数,存放单元节点的x坐标,存放顺序xi,xj,xm.  !
!    y----3个元的一维实数组,输入参数,存放单元节点的y坐标,存放顺序yi,yj,ym.  !
!    ke---3*3的二维实数组,输出参数,存放单元系数矩阵.                        !
!***************************************************************************!

subroutine UKE1(x,y,ke)

  dimension x(3),y(3),a(3),b(3)
  real ke(3,3)

  a(1)=y(2)-y(3)
  a(2)=y(3)-y(1)
  a(3)=y(1)-y(2)
  b(1)=x(3)-x(2)
  b(2)=x(1)-x(3)
  b(3)=x(2)-x(1)
  s=2.*(a(1)*b(2)-a(2)*b(1))
  
    do i=1,3
	  do j=1,i
	    ke(i,j)=(a(i)*a(j)+b(i)*b(j))/s
	    return
	  end do
	end do
end subroutine UKE1