自重的等效节点荷载问题
本帖最后由 ccdier 于 2012-8-13 15:31 编辑为何计算得到的重力的等效节点荷载，各节点的荷载值不一致呢？？
以8节点实体（六面体）单元为例，对某一单元，代码如下：
DO IGAUS=1,NGAUS
DO JGAUS=1,NGAUS
DO KGAUS=1,NGAUS
EXISP=POSGP(IGAUS)
ETASP=POSGP(JGAUS)
EYBSP=POSGP(KGAUS)
CALL SFR2(EXISP,ETASP,EYBSP)
CALL JACOB2(NNNODE,DJACB)
DVOLU=DJACB*WEIGP(IGAUS)*WEIGP(JGAUS)*WEIGP(KGAUS)
DO IEVAB=1,NDOFN*NNNODE
FBEMATX(IEVAB,1)=0.0
END DO
DO INODE=1,NNNODE ！NNNODE=8
IEVAB=(INODE-1)*NDOFN+1 !NDOFN=3
JEVAB=(INODE-1)*NDOFN+2
KEVAB=(INODE-1)*NDOFN+3
FBEMATX(IEVAB,1)=FBEMATX(IEVAB,1)+SHAPEN(INODE)*FBMATX(1)*DVOLU!FBMATX=T重度
FBEMATX(JEVAB,1)=FBEMATX(JEVAB,1)+SHAPEN(INODE)*FBMATX(2)*DVOLU
FBEMATX(KEVAB,1)=FBEMATX(KEVAB,1)+SHAPEN(INODE)*FBMATX(3)*DVOLU
END DO
END DO
END DO
END DO
其中，单元的法向轴线和z轴是平行的。。这样理论上不是各节点力应该相同吗？
另外，对于温度变化所产生的温度变形，转换成各单元节点等效荷载时，各节点力也是不一致的。。其中，温度变形只考虑X,Y,Z方向的变形，即ε=[αΔT,αΔT,αΔT,0,0,0]
计算公式为：
f=∫BT D ε dV
页:
[1]