【求助】单刚矩阵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
谢谢指点!! 这个只有自己慢慢看了,不过按照我的记忆应该是材料参数什么的。不过就我们用户而言,这些是不用修改的。
页:
[1]