鼓或带。 程序根据轧辊的辊身和辊颈在z向和r向的单元数自动形成图3的计算网格和结点坐标， 减少了准备数据的工作量，每个展开项的{8m}和{on}需存于外部设备中，在计算总的位移 和应力时逐项调用，程序的主要输出结果是：子午而上冬点在n=0~1各展开项时的位移 {òn}和应力{·n},以及在给定不同0角下，子午面上各结点的总位移{δ}和总应力{g}。 四、算例 为了考核基木原理的适用情况，程序设计的可靠性和计第精度，通过：个算例加以考证。 1.厚壁筒受内压 严壁例筒可看作为一个没有辊颈部份的中 空轧辊，它在内压作用下的径向位移和应力有 精确解答，可按公式计算出来，并与S-NSL 序训第的结果相比较，这时单元刘分如图7、 径向（父分6个单元，约束尔件使厚壁圆筒处于 平应变东件下，筒内壁作用有面载荷q,=p =2公斤/毫米2，q。=q2=0,按(17)式把 q,胫开时只有n=0一项。理论值和用S-NSL 程序计算结果的比较见专1。 图7 表1 r(毫米) 150 175 200 225 250 275 300 径向位移 理论值 272.381 241.134218.730202.222 189.842 180.462 173.333 1×104毫米 有限元值 271.148240.119217.856201.445189.134179.806172.717 周向应力 理论值 3.3333 2.6258 2.1667:1.8518 1.6266 1.4600 1.3333 0 公斤/毫米2 有限元值3.17092.5875 2.1422 1.8350 1.6143 1.4506 1.3439 径向应力 理论值 -2.0000-1.2925,-0.8333-0.5185-0.2933-0.1267 0 公斤/毫米：： 有限元值 -1.9909 -1.3507-0.8668-0.5392-0.3069-0.1361-0.0268 2.四辊轧机的工作辊 四辊轧机的工作辊D (图8)？到轧件支撑批 的力。若研究一种简单督 况：①暂不为虚工作辊受到 的扭摩擦力。②轧件 (a) 支撑辊对1工作辊的压力在0 方向的分布新人小、两者相 同，在z方向都是均匀分布 图8 (b) 的。 191鼓 或千告 。 程 序 根 据 轧辊 的辊 身和辊颈 在 向和 向的单元 数 自动形成图 的计算网格 和 结点 坐标 ， 减 少 了准备数据 的工作 量 ， 每个展 开项 的 乙 。 和 。 需存于 外 部设 备 中 ， 在 汁算 憋的 位移 和应 力 时逐 项调 用 ， 程 序 的 主 要 输 出结果 是 子午 面 上 各 汽在 二 。 各展 开项 时 的 位移 各 。 和应 力笼口 ， 以 及在给定 不 同 角下 ， 子 午面 上 各结 点 的 总位 移 乙 和 总应 力 。 四 、 算例 为了考 伙垫术原理 的适 用情 况 ， 程 序 没计 的可 靠性 和 计算 精度 ， 通 过 三个 算例 加 以 考证 。 厚 壁 简受 内压 厚 壁 圆 筒可 看作 为一 个没 有辊颈 部份 的 中 空轧辊 ， 它 在 内压作用 下 的径 向位 移 和应 力有 书 ， 确解 答 ， 可按 公 式计算 出来 ， 并 与 一 程 序 计算 的 结 果 相 比 较 ， 这 时单元 划 分如 图 、 径 向 仪分 个单元 ， 约 束 条件 使厚 壁 圆筒处于 平 而应 变 条件下 ， 筒 内壁作用 有面 载荷 二 公 斤 毫 米 “ ， ， 按 式 把 ， 展 开时只 有 二 。 一项 。 理 论 值和 用 一 程序 计算结果 的 比 较见 岌 。 表 图 毫 米 理 论 值 有 限元 值 。 … 。 。 一 … 云百丁丽 一 丽万三诬下爪…下 … … · ， 。 · 。 … ， 径 向位 移 ， 一 ‘ 毫 米 周 向应 力 … 公 拼飞米 径 向应 力 理 论 值 … 有限 元 位 … … 。 理 仑值 一 · 。 … 一 ， ‘ 一 ‘ 一 ‘ 一 ” 。 “ 门甘︵ 一 口 公斤 毫米 、 有 限元 值 一 一 一 一 一 一 四 辊 轧 机 的工 作 辊 四辊 轧机 的 工 作 辊 图 受到 轧件 和 支撑辊 的 仄 力 。 若 研 下己一 种 简单情 况 ① 哲 不 考 虑 工作 辊 受 到 的 扭 知 和 摩 擦 力 。 ②轧件 和 支撑辊 对 一 几作辊 的 压 力在 方 向 的 分 布和 人 小 两者 相 同 ， 在 方 向都是 均 匀分布 的 。 图 工