爆破拆除钢结构建筑计算机仿真算法及程序实现.pdf

第 “ 卷第 期 ““ 年 月 爆破 “ ,-;.06 “ , 8A;.5-2.3, 2 ,。 8; 数据录入系统中， 在窗体上添加 个 320 组件、 ’ 个 .50320 组件， 在其中依次输入最小边长 系数， 节理线及边界线总数， 材料线数目， 交叉单元 数， 约束点数， 荷载点数， 观测点数， 空洞点数； 在窗 体上添加 个 A0B/320 组件， 在 A0B/320 框中依序 按格式逐行输入节理线的第一端点 、 “ 坐标、 第二 端点 、 “ 坐标、 节理材料号； 材料线的第一端点 、 “ 坐标、 第二端点 、 “ 坐标、 块体材料号； 约束线的 第一端点 、 “ 坐标、 第二端点 、 “ 坐标、 荷载点的 、 “ 坐标、 观测点的 、 “ 坐标、 空洞点的 、 “ 坐 标。 为 “另存为” 菜单添加事件， 在其 C8-0BD 事件的 处理过程中完成如下代码 G,0HE6I“ *,J1,*8/BD（.1216EKL M1B） ； J,6 0 N 01O16； K1O0 0P *,J120,-EO“141BF1 /1 Q0/ I1IE“-01* 2E K1O0 B-1,6； ,5512 （120“14） ； ,5512 （120“14） ； ,5512 （*50120“14） ； ,5512 （*501209“14） ； ,5512 （*50120“14） ； ,5512 （*50120R“14） ； ,5512 （*50120’“14） ； ,5512 （*50120S“14） ； PE6 0 N E 60B/120“-01*“BEF 2E ,5512 （60B/120“-01* ［0］ ） ； *,J1EP0-1 （*,J120,-EO“P0-1,I1） ； 12； 12； “9H; 程序 笔者采用 1-5/0 ’“ 开发了 H; 数据录入 系统， 为 H 程序提供原始数据文本文件 H4“。 H; 数据录入系统中， 在窗体上添加 个 .50320 组件、 9 个 320 组件， 在其中依次输入计 算方式控制参数， 计算时间步数， 块体材料号数， 节 理材料号数， 假设的最大变形率， 每时步中时间间隔 上限， 接触弹簧刚度； 在窗体上添加一个 A0B/320 组件， 在 A0B/320 框中依序按格式逐行输入 R （） 行， 时间 、 在 时的 、 “ 方向的位移功荷载、 单位 面积质量、 单位面积在 、 “ 方向的重量、 弹模、 泊松 比、 、 “ 方向的初始应力、 垂直于 方向平面的初 始剪应力、 、 “ 方向的初始速度、 初始角速度、 摩擦 角、 粘滞力、 抗拉强度、 超松弛系数。 并为 “另存为” 菜单添加事件 （与 8; 类似） 。 “8 程序 8 程度是 H 程序的一个预处理程序， 它通过 输入的节理线把要计算的节理域分割成一个个独立 的块体。其原始数据文本文件由 8; 数据录 入系统生成。 “RH 程序 H 程序是用来进行静力分析或动力分析的。 其物理信息数据来自由 H; 数据录入系统生 成的文本文件。在每一个时间步中， 都要搜索一遍 块体与块体之间的接触关系， 然后计算出刚度子矩 阵、 荷载子矩阵、 质量子矩阵和接触子矩阵并加入到 平衡方程的力矩阵或系数矩阵 ［R， ’］ ， 每一个计算时 段， 需要进行开L合迭代， 以保证块体之间没有嵌入， 它是指块体之间刚性弹簧的加设和移开过程。计算 中需要检查相邻的顶点与顶点， 或顶点与边的接触 关系， 以确定如何放置刚性弹簧来保证块体与块体 之间没有相互嵌入。若迭代计算中出现嵌入， 则要 加一个刚性弹簧将其顶回去； 若原有的接触在法向 弹簧上出现拉力， 则要移开刚性弹簧， 然后相邻块体 将分离。在每一个计算时段末， 新生成的变形块体 爆破9 年 月 万方数据 的几何形状以及块体的应力、 速度和接触力信息会 被留作下一步的起始状态。在每一个时间步的开 合迭代中， 平衡方程既可采用直接方法求解， 又可采 用 “ 迭代方法求解。 ］ 曹金河， 于亚伦7第 A 卷第 8 期胡文军爆破拆除钢结构建筑计算机仿真算法及程度实现 万方数据