【求助】单刚矩阵aeq4.for 中的变量

lyn

能解释一下单刚矩阵aeq4.for中 子程序 shapc中的变量coefc 代表什么？
  
     call shapn(nrefc,ncoor,4,rporo(1,iporo),cporo,crtr,1,3,3)
      call shapc(nrefc,ncoor,6,coefd,coefc,crtr,2,5,5)
      poron=coef(1)
      po1=coef(2)
      pc1=coef(3)
      po=coef(4)
      pc=coef(5)
      sfs=coef(6)
      weigh=det*gaus(igaus)
      fpo=po
      fpc=pc
      fsfs=sfs
      fporo=poron
      call calcoef(fporo,fpo,fpc,fsfs,
     &       a1,b1,c1,a2,b2,c2,d2,a3,a4,pkx,pky)
      b31=+(coefc(4,1)*pkx)**2+(coefc(4,2)*pky)
     & **2
  
b31 的原计算式是
b31=+({po/x_i}*pk_i)**2    ,po 是耦合场的值
  
我发现子程序 shapc是在传耦合场值po 之前就调用并生成了coefc数组，说明coefc与po 值没有直接联系？
  
我的vde 文件中的stif 段
stif
$cv fpo=po
$cv fpc=pc
$cv fsfs=sfs
$cv fporo=poron
$c6 call calcoef(fporo,fpo,fpc,fsfs,
$c0      &       a1,b1,c1,a2,b2,c2,d2,a3,a4,pkx,pky)
$cv b31=+({po/x_i}*pk_i)**2
$c6 b3=sqrt(b31)
dist=+[poro;poro]*0.0
  
load=+[poro]*a3*b3
  
end
  
谢谢指点！！

zhoulon

这个只有自己慢慢看了，不过按照我的记忆应该是材料参数什么的。不过就我们用户而言，这些是不用修改的。