[Fortran] 纯文本查看 复制代码
!*****节点应变的有限元计算，该程序参考曾攀的《有限元分析及应用》中的有限元程序编得，
!******理论基础可见于曾攀的《有限元分析及应用》与王勖成的《应用有限单元法》等。
Subroutine disp(x, y, mel, ak, stsmel, ijku)
!*****MEL 单元节点编号数组，X，Y 节点坐标数组，AK 整体刚度矩阵，JR 节点约束信息，P 节点载荷
!*****IJKU STSMEL 单元应变，STSNP 节点应变
Dimension mel(5000, 3), x(5000), y(5000), u(5000, 3), p(5000, 3), ijku(5000, 3), ren(5000), ak(5000, 5000), stsmel(5000, 3), stsnp(5000, 3), jr(5000, 3)
Common /ca/np, ne, nr, nl, band, ndof
Common /cb/c11, c12, c44, q11, q12, q44
Common /cc/n, m, time, h, hx, hy
!*****计算节点载荷
Call form(nload, mel, x, y, resul, ae)
!*****处理边界条件
Call dobc(ijku, ak, resul)
!*****计算节点位移
Call slove(ak, resul)
!*****计算单元与节点应变
Call calsts(mel, x, y, resul, stsnd)
Return
End Subroutine disp