其中符号 ,门，，门，一为单元上结点i,j在局部坐标中的坐标值。 2a,2b一为真实单元在z向和r向的边长（图4）。 r。一单元的平均半径。 即 B,=5:51 B、=n,n) Cg=5,+ξ， Cn=n:+门1 dn =24a(3+B,〔C.+C.(C,-B,) (11)(1,-1) 3■ Co=In Tota 1 ro-a C1=2-C。 (-1,1)(-1,1) G,入一为拉梅常数 图4 E Eμ G=2(1+4),λ=a+)1-2μ) 单元刚度矩阵（8)为12×12的方阵，并应注意到(9)中h11,h22,h33和n有关， h2,h21,h23,h3z和n有关。 在轧辊的计算中只有表面力的作用，暂不考虑热负荷的影响，因而单元的负荷列阵为 4f)=((N)T4f}Rdz (11) 若面力作用于单元的2~3边上，且为均匀分布（图5a),则(11）式经积分后得 {fn}°=〔0,0,0，q:n,qgn,qzn,q:n,qen,gzn,0,0,0)T ×Rb (12) 若2~3边上有线性分布的方载荷（图5b),则(11） (a) 式积分后 i=0,00号q+9号+ 2 3939n, 26 2 1 1 21 2 922n+3932,3921n+33m,3920a+3930, 892n+ 1 (b) 393zn,0,0,0)TRb (13) 其中R一2~3边的半径。 2b (12)和(13)式相当于把面载荷的合力(q1+q2)2b乘R 后按静力学原理分配到2,3点上去，1,4点不受影响。 图5 由〔K),{fn}以及(8)式可解出{δa},若位移，应力，载荷按(3)~(5)式版 开到n1项，则子午面上任一点的总位移为 188其中符号 乙 ‘ ， ‘ ， 毛，， ” ， － 为单元 上结点 ， 在局 部坐标中的坐标值 。 ， － 为真实单元 在 向和 向的边 长 图 。 。 － 单元 的 平均半径 。 以， 甲 万 乙 ‘ 自 一 、 息 ‘ 息 ， 。 二 息 ‘ 七， ， ， ， ， 月 ‘ 月 题 ‘ ， 见三毕 ‘ “ ‘ 蠢 ， 〔 。 一 ，〕 。 二 一 一 。 ， 入－ 为拉梅 常数 右 ” 一 ， 一 ，一 图 件 久 卜 件 一 卜 单元 刚度矩 阵 为 的方 阵 ， 并应 注意 到 中 ， ， ， ， 和 “ 有关 ， ， ，， ， 和 有关 。 在 轧 辊 的计算 中只 有表面 力的作用 ， 暂 不 考虑 热 负荷 的影 响 ， 因而单元 的 负荷 列 阵 为 、“ ’ 〔 〕 ’ “ · ， 若面 力作用 于 单元 的 边 上 ， 且 为均匀分布 图 ， 则 式 经积 分后得 引上下 韶曰四 王 〔 ， ， ， 。 ， 。 ， ， ， ， ， ， ， ， 〕 若 边 上有 线性 分布的 方 载荷 式积 分 后 图 ， 则 玉 〔。 ， 。 ， 。粤口 不 ， 万口 、 。 ， 口 了 副上丁 了 ， 叭 一 百 一 。 ， ， 〕 下 一 一万 公 其 中 － 、 边 的半径 。 和 式 相 当于 把 面 载荷的 合 力蚤 乘 后 按 静力学原 理分 配到 ， 点 上去 ， ， 点 不 受影 响 。 由〔 。 〕 ， 。 以 及 式可 解 出 ， 若位 移 ， 应 力 ， 开到 项 ， 则 子午面 上任一 点 的 总 位 移 为 图 载荷按 式展