第2期 马明军等：岩体动态响应分析程序的微机化及其应用 .109. 数组也会遇到同样的问题，而且大量数组的虚实结合还会降低运行速度， 2 DYNPAK程序的微机化改造 二维非线性动态有限元程序DYNPAK是用FORTRAN语言写成的，原始程序由Owen 和Hinton于1980年给出).该程序可模拟材料的物理非线性和几何非线性，以及对地震载 荷或脉冲载荷的响应·微机化改造的主要内容是引进前面讨论的内存空间数据处理方 法，以便：(1)该程序能在普通微机上运行；(2)扩大程序的解题能力；(3)增强程序的 可读性，便于修改和扩充，此外，还对该程序的输入系统做了改进，将传统的固定格式输 入改为可标识的自由格式输人，大大地方便了输入数据的准备，下面是改造后的部分主程序， PROGRAM DYNPAK IMPLICIT REAL (A-H,O-Z) COMMON MTOT,NP,IA (100 000) COMMON/CONTRL/NPOIN,NELEM,NDOFN,NMAT,NVFIX,... COMMON/IOFILE/IPR,NF1,NF2,NF3,. … CALL ARRDFN (KCORD,KIFPR,KMATN,. CALL INPUTD (KCORD,KIFPR,KMATN,.. …… CALL LUMASS (KCORD,KMATN,.. … STOP END SUBROUTINE ARRDFN(KCORD,KIFPR,KMATN,. IMPLICIT REAL (A-H,O-Z) COMMON MTOT,NP,IA (1) COMMON/CONTRL/NPOIN,NELEM,NDOFN,NMAT,NVFIX,.. CALL DEFINE (CORD',KCORD,NDOFN,NPOIN) CALL DEFINE (IFPR,KIFPR,NDOFN,NPOIN CALL DEFINE (MATN',KMATN,1,NELEM ”+ RETURN END 改造后的DYNPAK程序可在普通IBM-PC/XT/AT/286/386微机上运行，由于保留 了原始程序的核心部分（语法按MS-FORTRAN77修改）和基本功能，示例计算的结果与 Owen和Hinton给出的结果吻合很好，限于篇幅，在此不一一给出.第 期 马 明军等 岩体动态响应分析程序 的微机化及其应用 数组也会遇到 同样 的 问题 ， 而 且大量 数组 的虚 实结合还 会 降低 运行 速 度 程序 的微机化 改造 二 维非 线性 动态有 限元程序 是 用 人 语言 写 成 的 ， 原 始 程 序 由 和 于 年给 出 「】 该程 序可模拟 材料 的物理 非 线性 和几何 非 线性 ， 以 及 对 地 震 载 荷 或脉 冲 载 荷 的 响应 微 机 化 改 造 的 主 要 内容 是 引 进 前 面 讨 论 的 内 存 空 间 数 据 处 理 方 法 ， 以便 该程序能在普通微 机上 运行 扩大程序 的解题 能 力 增 强 程 序 的 可读性 ， 便于修改和 扩充 此外 ， 还 对该 程 序 的输 入 系 统 做 了 改 进 ， 将 传 统 的 固 定 格 式输 人改 为可 标识 的 自由格式输人 ， 大大地方便 了输人数据 的准备 下面是 改造后 的部分主程序 仃 一 ， 一 ， ， 丁 刊田 呵 ， ， ， ， 』 ， … 刀〕 正 ， ， 于 ， ， ， … ， ， ， … ， ， ， … ， ， … 叮 ， ， ， … 吐 】 ， ， 一 均 ， ， 陌 ， ， ， ， ， … ’ ， ， ， ， ， ， ， 卜丁 卜 ’ ‘ ， ， ， 改造后 的 程 序 可 在 普 通 一 叉丁 戌 ， 微 机 上 运 行 由于 保 留 了原始程序 的核心部分 语法按 一 修 改 和基本功 能 ， 示 例计 算 的结 果 与 叭七 和 给 出的结果 吻合很 好 ， 限于 篇 幅 ， 在此不 一 一 给 出