水平节，点力等丁零（无摩擦）。压头以外的 节点力为零。对于下部边界，因为对称，所以 开始 水平节点力为零，垂直位移为零。右边界为节 数据 点力为零。左边界，因为对称，所以垂直节点 自动划分元素，节点 力为零，水平位移为零。 编号，计算节点座标 6.解方。采用〔K)阵一维紧缩存 给定正头每次压下的位移 储，非用追赶法求解方程。 PP:=1 7.计算应变增量和加权系数。为得到较 好的m值，迭代2次。 YES 8.计算△{σ}，△{F}等。由计算所得的 No △{e},利用(13)，(9)和(5)式得 QQ:=1 A{o},{o},。并利用△{8}、〔K)和(6) 让饰元素刚度矩碎 式得A{F}。则第i次计算所得的总节点力为 {},={}:-1+△{F}:。由节点力就可以求 集合成结构刚度加阵 得压头总的压力，它除以压头面积，即得平均 运用边界条件 压力。 解方程，得△{8 9.输出结果。输出每个元素的m值，用 它来判别元素屈服与否。输出元素的应力，等 i计算A{c}.m 效应力，节点力，压头压力，平均压力和节 Q- No QQ:QQ+ 点位移等。 YES 10.压头每向下压一个位移，必须计算 计算△{o,{o}, 一次，PP丧示计算的次数。听PP=1时，可 可A{F},F) 把整个结构作为弹性体计算，得到每个元素的 E头正力半均团 等效应力g,确定·max,令0max/0,=L, 则把计算所得的值，除以L,即可得到被压材 料进入屈服时的位移，应力，应变节点力等。 No 全屈服> ppi=pp+1 11,片刚性区与边界相连，即达到全屈 YES 服状态时，此时的平均压力即为平面应变的变 停机 形抗力。 图2计算程序框图 四、压入问题计算结果 本文计流了已具方106潮件交形的工业纯铝的E入问题。其屈服极限σ。=6.17i公斤/ 毫米2，啪性模屉E=7201.6公斤/毫米°，泊松比v=0.33。它的应力一应变关系如图8中 曲线I所示。计算长度1=5mm。 取以下几种滑况进行计第。 1) 0=1毫米， h=2毫米， ,'h=1 2。 2) n=1毫米 h=1.5毫米， a/h=2 ● 3 3) a=1毫米 h=1毫米 ,a/h=1。 4) a=1毫米， h=0.5毫米，a/h=2。 173水 平节点 力件千 零 无 摩擦 。 压 头 以 外 的 竹点 力为零 。 对于 下 部 边 界 ， 因 为对 称 ， 所 以 水 平节点 力 为零 ， 垂 直 位 移 为零 。 右边 界 为节 点 力 与零 。 左 边 界 ， 因 为对 称 ， 所 以 垂 直 方点 力 为零 ， 水 平位 移 为零 。 解方 程 。 采 用 〔 〕矩 阵一 维 紧缩 存 储 ， 并 用 追 赶 法求解方 程 、卜算应 变增 。 和加 权 系数 。 为得 到较 好 的 故 ， 迭 代 次 。 计算 △ ， △ 等 。 由计算所得 的 。 卜 ， 利 用 ， 和 式 得 ， ， 。 并利用 △ 乙 、 〔 〕和 式得 州 卜 。 则 第 次 计 算 所 得 的 总节点 力 为 魂 ， 王 ‘ 一 ， △ ‘ 。 由节点 力 就可 以 求 得 仄 失总 的 力 ， 它除 以 压头 面积 ， 即 得 平均 压 力 。 输 出结果 。 输 出每个元 素 的 值 ， 用 它来 判别元 素屈 服 与 否 。 输 出元 素的应 力 ， 等 效应 力 ， 节点 力 ， 压 头 压 力 ， 平 均 压 力和节 点 位 移 等 。 压 少协 向下 压一 个位 移 ， 必 须 计 算 一 次 ， 友示 计算 的次 数 。 当 二 时 ， 可 把 整 个 结 构作 与弹 性 体 计算 ， 得 到 每个元 素 的 等效应 力 ， 确定 。 、 ， 令 。 、 、 ， 则 把 卜算所得 的 值 ， 除 以 ， 即可得 到被 压 材 料进 入 屈 服 时 的 位 移 ， 应 力 ， 应 变 和 节点 力 等 。 当 塑性 区 与 边 界相 连 ， 即 达 到 全 屈 服 状 态 时 ， 此 时 的 平均压 力即 为平 面应 变 的变 形 抗力 。 自动划分 元 素 ， 节点 编 号 ， 计 算节点座标 计算△ 魔 》 ， 笼 ， 下 八 卜 王 卜 图 计 算程 序框 图 四 、 压 入 问题 计 算结果 术 之计 牛 已 具 汀 丫 塑性 呀形 工 业 纯 铝 的 压 入 问 题 。 其屈 服 极 限 了 公斤 毫 米 “ ， 弹 性 模 鼠 落 公 斤 毫 米 巳 ， 泊 松 比 二 。 。 它 的应 力— 应 变关 系如 图 中 线 所 示 。 计算 长度 。 取 以 下 几 种 洁况 进子 一 八 一 算 。 一一一 ‘ ‘ 川 了 比且 幻幻勺 ‘ 二 毫 米 毫 米 ‘‘ 毫 米 二 毫米 毫米 毫米 ， 毫米 毫米 勺性八月 、、产产、产、声