关于有限元程序中得到雅克比矩阵的子程序

zty · 发表于 2023-6-24 23:07:41

我无法得知二维情况下，循环中的公式是什么原理，或者说pn和co_xn是什么变量

[Fortran] 纯文本查看 复制代码

01subroutine get_jcob(detj,jcob,jinv,pn,co_xn)implicit none
02real(8)::co_xn(nnd_e,dimens),detj
03real(8)::jcob(dimens,dimens),jinv(dimens,dimens),pn(nnd_e,dimens)
04integer::i,j,k
05jcob=0.0;jinv=0.0
06if(dimens==2) then
07   do i=1,nnd_e
08      jinv(1,1)=jinv(1,1)+pn(i,2)*co_xn(i,2)
09      jinv(1,2)=jinv(1,2)-pn(i,1)*co_xn(i,2)
10      jinv(2,1)=jinv(2,1)-pn(i,2)*co_xn(i,1)
11      jinv(2,2)=jinv(2,2)+pn(i,1)*co_xn(i,1)
12   end do
13   detj=jinv(1,1)*jinv(2,2)-jinv(1,2)*jinv(2,1)

[子程序] 关于有限元程序中得到雅克比矩阵的子程序