模块分割法在断裂力学分析微机程序中的实现

<正> 为合理使用材料,针对微机普遍使用的现状,我们编制了二维弹塑性断裂力学的微机程序。断裂力学的计算基本是在大型计算机上进行的,因为它的内存大,计算速度快且容易编程。在结构分析中的本构方程为: 〔K〕{δ}={F} 其总刚只能容纳在大型机的内存中。一般二维问题总刚数为(2 n)~2次,三维的为(3n)~2次,例如节点1 200个则总刚数为576×10~4。这样大的数字在微机的内存中是无法存储的,即便使用大家熟知的变带宽、非零一维数组存储也需要很大的内存,在运算中也要解很多次的联立方程。为此,我们在二维弹塑性微机程序中使用了模块分割法。