lyn 发表于 2004-7-8 06:10:21

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

能解释一下单刚矩阵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=+*0.0

load=+*a3*b3

end

谢谢指点!!

zhoulon 发表于 2008-10-6 10:53:40

这个只有自己慢慢看了,不过按照我的记忆应该是材料参数什么的。不过就我们用户而言,这些是不用修改的。
页: [1]
查看完整版本: 【求助】单刚矩阵aeq4.for 中的变量