急！！！求救各位大神

wutianwen

subroutine Get_dN_dxi(dN_dxi,xi,eta) ! 得到形函数对局部坐标系的导数
implicit none
real*8 :: dN_dxi(:,, xi,eta
real*8:: x1, x2, x3, x4, x5 ,x6 ,x7 ,x8
x1=.25*(1.-eta);
x2=.25*(1.+eta);
x3=.25*(1.-xi);
x4=.25*(1.+xi)

dN_dxi(1,1)=x1*(2.*xi+eta)
dN_dxi(1,2)=-8.*x1*x2
dN_dxi(1,3)=x2*(2.*xi-eta)
dN_dxi(1,4)=-4.*x2*xi
dN_dxi(1,5)=x2*(2.*xi+eta)
dN_dxi(1,6)=8.*x2*x1
dN_dxi(1,7)=x1*(2.*xi-eta)
dN_dxi(1,8)=-4.*x1*xi
dN_dxi(2,1)=x3*(xi+2.*eta)
dN_dxi(2,2)=-4.*x3*eta
dN_dxi(2,3)=x3*(2.*eta-xi)
dN_dxi(2,4)=8.*x3*x4
dN_dxi(2,5)=x4*(xi+2.*eta)
dN_dxi(2,6)=-4.*x4*eta
dN_dxi(2,7)=x4*(2.*eta-xi)
dN_dxi(2,8)=-8.*x3*x4
return
end subroutine Get_dN_dxi
各位大神看看这里面eta的一串语句表示什么意思?这是一段八节点单元形函数对坐标轴求导的fortran程序

penultimate

把偏导的公式写出来，对一下程序就知道了，就是对两个变量的偏导在xi，eta处的值

wutianwen

penultimate 发表于 2017-5-13 11:26
把偏导的公式写出来，对一下程序就知道了，就是对两个变量的偏导在xi，eta处的值 ...

我完全是个新手，能否解释得再详细一点

cadence

real*8 :: dN_dxi(:  ,   :    )    ,    xi,eta，这个是可变数组的写法
也可以写成
real*8 :: dN_dxi(2,8) , xi,eta
我一般写固定数组形式，免得出错。
另外，我发现用老的FORTRAN77格式写的程序，最保险
好查错。

wutianwen

cadence 发表于 2017-5-13 23:19
real*8 :: dN_dxi(:  ,   :    )    ,    xi,eta，这个是可变数组的写法
也可以写成
real*8 :: dN_dxi(2,8) ...

谢谢老哥，那xi+eta和eta-xi是什么意思呢

cadence

看一下有限元的形函数部分表达式