本发明提供一种正交异性钢桥面板的计算方法，包括将纵肋保持高度不变，按刚度与质量等效为矩形肋，视为Euler‑Bernoulli梁单元进行分析，并在分析中考虑其偏心效应；将横隔板按实际的弹性约束刚度简化为弹性支承，母板沿横向支承在两侧腹板上，腹板亦按实际的弹性约束刚度简化为弹性支承；母板运用Kirchhoff薄板理论分析。分析中，考虑横截面的纵向焊接残余应力效应，由于正交异性钢桥面板的母板整体受压而纵肋整体受拉，因此，可视等效截面的母板均布受压而纵肋均布受拉。在此基础上，基于能量原理建立正交异性钢桥面板的刚度方程，并可采用MATLAB软件方便地计算结构刚度矩阵。
1.一种正交异性钢桥面板的计算方法，其特征在于，包括：步骤1.搭建理论模型
正交异性钢桥面板，由母板、x方向的纵肋、y方向的横隔板焊接而成，将x方向的纵肋保持高度不变，按刚度与质量等效为矩形肋，视为Euler‑Bemoulli梁单元进行分析，其偏心效应通过中性面的位置进行考虑：将y方向的横隔板按实际的弹性约束刚度简化为弹性支承，母板沿横向弹性支承在两侧腹板上，两侧腹板亦按实际的弹性约束刚度简化为弹性支承；
1)母板x与y方向的长度分别为a与b，x、y、z方向的位移分别用u、v、w表示；
2)正交异性钢桥面板的密度、弹性模量、剪切模量、泊松比分别为ρ、E、G、μ；
3)x方向的纵肋、y方向的横隔板的数量分别为Nx、Ny，x方向的纵肋、y方向的横隔板的间距分别为b1、a1；
4)第i条x方向纵肋及第i条y方向横隔板坐标分别为yi与xi；
5)母板的厚度为tp，x方向的纵肋的等效厚度为tr、x方向的纵肋的高度为hs、截面积为Asx；
7)正交异性钢桥面板垂直于x轴截面的中性面距母板顶缘hx；
8)母板的弯曲刚度为 x方向的纵肋的拉压刚度、弯曲刚度、扭转刚度分别为EAsx、EIsx、GJsx；
9)y方向的横隔板及两侧腹板作为弹性支承考虑，其弹性约束刚度由各自的几何、材料参数计算得到，y方向的横隔板的竖向拉压刚度及绕y轴的转动刚度分别用kda、kdy表示，两侧腹板的竖向拉压刚度及绕x轴的转动刚度分别用kwa、kwx表示，y方向的横隔板及两侧腹板的其余刚度对结构的影响不予考虑；
步骤2.对母板运用Kirchhoff薄板理论分析考虑母板横截面的纵向焊接残余应力效应，由于正交异性钢桥面板的母板整体受压而x方向的纵肋整体受拉，因此，视等效截面的母板均布受压而x方向的纵肋均布受拉，具体分析如下：
1)令母板的均布纵向初始应力为σpx0、x方向的纵肋的均布纵向初始应力为σsx0，由初始应力自平衡条件得到：σpx0b1tp+σsx0hstr＝0 (1)
2)沿着正交异性钢桥面板的纵向，考虑母板与x方向的纵肋中的纵向焊接残余应力不变；
3)正交异性钢桥面板的横向焊接残余应力效应较小，对结构的影响不予考虑；
步骤3.基于能量原理建立正交异性钢桥面板的方程具体过程如下，根据Hooke定律，母板与加劲肋的应力与应变关系如下：σsx＝Eεsx (2‑4)式(2‑1)至(2‑4)中：σpxx与σpyy分别为母板内x方向、y方向的正应力；τpxy为母板面内剪应力，σsx与σsy分别为x方向、y方向加劲肋的轴向应力；εpxx与εpyy分别为母板内x方向、y方向的正应变，γpxy为母板面内剪应力，εsx与εsy分别为x方向、y方向加劲肋的轴向应变；
由于母板与纵肋的纵向初始应力分别为σpx0、σsx0，因此，根据Hooke定律，得母板与加劲肋的初始应变为：式(3‑1)至(3‑3)中：εpx0为母板的纵向初始应变；σpy0是母板的横向初始应力；εpy0为母板的横向初始应变；εsx0是纵肋的纵向初始应变；
根据针对弹性支承加劲板，经过简单的几何运算，再结合式(2‑1)、(2‑4)、(3‑1)、(3‑3)得到纵肋中性轴及母板中性面的应变分别为：母板的应变能由弯曲应变能及中性面的拉压应变能组成：纵肋的应变能为：
y方向横隔板与两侧腹板作为弹性支承，取其中心位移计算各自的等效应变能：综上，正交异性钢桥面板的总应变能为：
U＝Up+Usx+Udw (8)对总应变能U进行变分运算，即得刚度矩阵：P sx dw
[K]＝[K]+[K ]+[K ] (9)p sx dw
式中：[K]为刚度矩阵，[K]、[K ]、[K ]分别为母板、x方向的纵肋、y方向的横隔板与两侧腹板对应的刚度矩阵，根据母板的横向位移w，以及正交异性钢桥面板的详细结构与材料参数，得到刚度矩阵；
式中：wmn为广义坐标，对应于第(m，n)阶屈曲模态的幅值， 和φn(y)分别为沿着x方向与y方向的位移函数，表示如下：φn(y)＝Ancos(any)+Bnsin(any)+Cncosh(any)+Dnsinh(any) (12)式中：常数Am、常数Bm、常数Cm、常数Dm、常数am，根据x方向的纵肋的边界条件确定，常数An、常数Bn、常数Cn、常数Dn、常数an，根据y方向的横隔板的边界条件确定；
对于x＝0及x＝a边，根据边界条件确定各参数，如：两边简支：
对于y＝0及y＝b边，两边由两侧腹板支承，力的边界条件为：y＝0边：支承力Qy＝kwaφn(n)、力矩y＝b边：支承力Qy＝kwaφn(b)、力矩刚度矩阵的形式根据应变能及母板的横向位移确定，母板的刚度矩阵为：其中，元素 的表达形式为：
y方向的横隔板与两侧腹板的刚度矩阵为：其中，元素 的表达形式为：
步骤4.运用MATLAB软件计算刚度矩阵，进而对正交异性钢桥面板进行详细地局部计算。
一种正交异性钢桥面板的计算方法
技术领域
[0001] 本发明属于桥梁结构分析技术领域，尤其涉及一种正交异性钢桥面板的计算方法。
背景技术
[0002] 近年来，正交异性钢桥面板因具有自重小、整体性好、承载力高、加工制造方便等众多优点而广泛应用于大跨度桥梁中。目前，正交异性钢桥面板的制造已基本采用焊接，在焊接过程中将产生焊接残余应力，焊接残余应力的存在必然使得正交异性钢桥面板的力学特性发生变化，如果仍然按照理想状态进行结构分析，结果将偏离实际。因此，实际分析中需要考虑焊接残余应力的影响。目前，正交异性钢桥面板的结构分析方法有等效正交异性板法、等效格子梁法、有限单元法等，其中尤以有限单元法的应用最为广泛。尽管有限单元法具有计算方便、高效等优点，但是一方面，精细分析时单元划分的数量非常多，计算工作量大，从而使得计算分析时间增加；另一方面，作为数值方法，有限单元法很难直观地反映正交异性钢桥面板的受力机理与受力特征。
发明内容
[0003] 本发明的目的在于解决上述现有技术存在的缺陷，提供一种正交异性钢桥面板的计算方法。
[0004] 本发明采用如下技术方案：
[0005] 一种正交异性钢桥面板的计算方法，包括：
[0006] 步骤1.搭建理论模型
[0007] 正交异性钢桥面板，由母板、x方向的纵肋、y方向的横隔板焊接而成，将纵肋保持高度不变，按刚度与质量等效为矩形肋，视为Euler‑Bemoulli梁单元进行分析，其偏心效应通过中性面的位置进行考虑：
[0008] 将横隔板按实际的弹性约束刚度简化为弹性支承，母板沿横向支承在两侧腹板上，腹板亦按实际的弹性约束刚度简化为弹性支承；
[0009] 参数解释如下：
[0010] 1)母板x与y方向的长度分别为a与b，x、y、z方向的位移分别用u、v、w 表示；
[0011] 2)正交异性钢桥面板的密度、弹性模量、剪切模量、泊松比分别为ρ、E、 G、μ；
[0012] 3)纵肋、横隔板的数量分别为Nx、Ny，纵肋、横隔板的间距分别为b1、a1；
[0013] 2)第i条纵肋及第i条横隔板坐标分别为yi与xi；
[0014] 5)母板的厚度为tp，纵肋的等效厚度为tr、纵肋的高度为hs、截面积为Asx；
[0015] 6)纵肋横截面的中性轴距母板顶缘hsx；
[0016] 7)正交异性钢桥面板垂直于x轴截面的中性面距母板顶缘hx；
[0017] 8)母板的弯曲刚度为 纵肋的拉压刚度、弯曲刚度、扭转刚度分别为EAsx、EIsx、GJsx；
[0018] 9)横隔板及两侧腹板作为弹性支承考虑，其弹性约束刚度由各自的几何、材料参数计算得到，横隔板的竖向拉压刚度及绕y轴的转动刚度分别用kda、kdy表示，腹板的竖向拉压刚度及绕x轴的转动刚度分别用kwa、kwx表示，横隔板及腹板的其余刚度对结构的影响不予考虑；
[0019] 步骤2.对母板运用Kirchhoff薄板理论分析
[0020] 考虑横截面的纵向焊接残余应力效应，由于正交异性钢桥面板的母板整体受压而纵肋整体受拉，因此，视等效截面的母板均布受压而纵肋均布受拉，具体分析如下：
[0021] 1)令母板的均布纵向初始应力为σpx0、纵肋的均布纵向初始应力为σsx0，由初始应力自平衡条件得到：
[0022] σPx0b1tp+σsx0hstr＝0 (1)
[0023] 2)沿着正交异性钢桥面板的纵向，考虑母板与纵肋中的纵向焊接残余应力不变；
[0024] 3)正交异性钢桥面板的横向焊接残余应力效应较小，对结构的影响不予考虑；
[0025] (4)基于能量原理建立正交异性钢桥面板的方程，具体过程如下：
[0026] 根据Hooke定律，母板与加劲肋的应力与应变关系如下：
[0027]
[0028]
[0029]
[0030] σsx＝Eεsx (2‑4)
[0031] 式(2‑1)至(2‑4)中：σpxx与σpyy分别为母板内x、y方向的正应力；τpxy为母板面内剪应力，σsx与σsy分别为x、y方向加劲肋的轴向应力；εpxx与εpyy分别为母板内x、y方向的正应变，γpxy为母板面内剪应力，εsx与εsy分别为x、y方向加劲肋的轴向应变；
[0032] 由于母板与纵肋的纵向初始应力分别为σpx0、σsx0，因此，根据Hooke定律，得母板与加劲肋的初始应变为：
[0033]
[0034]
[0035]
[0036] 式(3‑1)至(3‑3)中：εpx0为母板的纵向初始应变；σpy0是母板的横向初始应力；εpy0为母板的横向初始应变；εsx0是纵肋的纵向初始应变；
[0037] 根据图2(a)‑图2(c)所示的加劲板，经过简单的几何运算，再结合式 (2‑1)、(2‑4)、(3‑1)、(3‑3)得到纵肋中性轴及母板中性面的应变分别为：
[0038]
[0039]
[0040]
[0041]
[0042] 母板的应变能由弯曲应变能及中性面的拉压应变能组成：
[0043]
[0044] 纵肋的应变能为：
[0045]
[0046] 横隔板与腹板作为弹性支承，取其中心位移计算各自的等效应变能：
[0047]
[0048] 综上，正交异性钢桥面板的总应变能为：
[0049] U＝Up+Usx+Udw (8)
[0050] 对总应变能U进行变分运算，即得刚度矩阵：
[0051] [K]＝[KP]+[Ksx]+[Kdw] (9)
[0052] 式中：[K]为刚度矩阵，[Kp]、[Ksx]、[Kdw]分别为母板、纵肋、横隔板与腹板对应的刚度矩阵，根据母板的横向位移w，以及正交异性钢桥面板的详细结构与材料参数，得到刚度矩阵；
[0053] 母板的横向位移表示如下：
[0054]
[0055] 式中：wmn为广义坐标，对应于第(m,n)阶屈曲模态的幅值， 和φn(y)分别为沿着x与y方向的位移函数，表示如下：
[0056]
[0057] φn(y)＝An cos(any)+Bn sin(any)+Cn cosh(any)+Dn sinh(any) (12)[0058] 式中：常数(Am、Bm、Cm、Dm、am)根据x方向板的边界条件确定，常数(An、 Bn、Cn、Dn、an)根据y方向板的边界条件确定；
[0059] 对于x＝0及x＝a边，根据边界条件确定各参数，如：
[0060] 两边简支：
[0061] 两边固支：
[0062] 对于y＝0及y＝b边，两边由腹板支承，力的边界条件为：
[0063] y＝0边：支承力Qy＝kwaφn(n)、力矩
[0064] y＝b边：支承力Qy＝kwaφn(b)、力矩
[0065] 刚度矩阵的形式根据应变能及母板的横向位移确定，母板的刚度矩阵为：
[0066]
[0067] 其中，元素 的表达形式为：
[0068]
[0069]
[0070] 纵肋的刚度矩阵为：
[0071]
[0072] 其中，元素 的表达形式为：
[0073]
[0074] 横隔板与腹板的刚度矩阵为：
[0075]
[0076] 其中，元素 的表达形式为：
[0077]
[0078] 步骤3.运用MATLAB软件方便地计算刚度矩阵，进而对正交异性钢桥面板进行详细地局部分析。(计算出正交异性钢桥面板的刚度矩阵之后，才能进行下一步的计算，即对正交异性钢桥面板进行受力计算)。
[0079] 进一步的，步骤3包括程序段：
[0080] 1)定义结构参数
[0081] 定义变量：syms E G μ ρ a b a1 b1 tp tr hs Asx hsx hx kda kdy kwa kwx σpx0 σsx0[0082] 节点坐标：xyz＝[；]
[0083] 材料参数：mat＝[E；μ；ρ]
[0084] 1)计算刚度矩阵
[0085] 2)定义函数：funetionK_total＝Node_total(Ke)
[0086] 计算单元数量及整体刚度矩阵：
[0087] J＝length(Ke)；
[0088] K_total(1:1+J,1:1+J)＝0
[0089] Ke_J(1:2,1:2)＝0；
[0090] forI＝1:J
[0091] Ke_J＝Ke(I)
[0092] forN1＝1:2
[0093] forN2＝1:2
[0094] K_total(N1+I‑1,nN2+I‑1)＝K_total(N1+I‑1,N2+I‑1)+Ke_J{1,1}(N1,N2)；
[0095] end
[0096] K_total。
[0097] 本发明的有益效果：
[0098] 本发明能够考虑纵肋的扭转刚度及偏心效应，并且可以简单地将横隔板的约束作用体现出来，同时亦考虑焊接残余应力效应。本发明不仅物理意义明确，考虑因素全面，而且计算分析效率高。
附图说明
[0099] 图1(a)‑图1(c)为正交异性钢桥面板示意图；
[0100] 图2(a)‑图2(c)为弹性支承加劲板计算示意图；
[0101] 图3为等效截面的纵向应力分布图；
[0102] 图4为梯形加劲板结构示意图；
[0103] 图5为本发明的步骤流程图。
具体实施方式
[0104] 为使本发明的目的、技术方案和优点更加清楚，下面本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0105] 注：“x方向的纵肋”和“y方向的横隔板”是相对于图1而言的，即x 方向的加劲肋称为“纵肋”，y方向的加劲肋称为“横隔板”，以下可以将“x方向的纵肋”简称为“纵肋”，将“y方向的横隔板”简称为“横隔板”。
[0106] 弹性支承”是边界支承条件的一种。
[0107] “正交异性钢桥面板”的两侧是腹板。
[0108] “纵肋”与“横隔板”均是“加劲肋”的一种，即“加劲肋”包括“纵肋”与“横隔板”。
[0109] 如图1(a)‑图1(c)、图5所示，本发明的一种正交异性钢桥面板的计算方法，包括：
[0110] (1)正交异性钢桥面板，由母板、x方向的纵肋、y方向的横隔板焊接而成，将纵肋保持高度不变，按刚度与质量等效为矩形肋，视为Euler‑Bemoulli 梁单元进行分析，其偏心效应通过中性面的位置进行考虑：
[0111] 将横隔板按实际的弹性约束刚度简化为弹性支承，母板沿横向支承在两侧腹板上，腹板亦按实际的弹性约束刚度简化为弹性支承。
[0112] 参数解释如下：
[0113] 1)母板x与y方向的长度分别为a与b，x、y、z方向的位移分别用u、 v、w表示；
[0114] 2)正交异性钢桥面板的密度、弹性模量、剪切模量、泊松比分别为ρ、 E、G、μ；
[0115] 3)纵肋、横隔板的数量分别为Nx、Ny，纵肋、横隔板的间距分别为b1、 a1；其中b3表示梯形肋的底板宽度。
[0116] 4)第i条纵肋及第i条横隔板坐标分别为yi与xi；
[0117] 5)母板的厚度为tp，纵肋的等效厚度为tr、纵肋的高度为hs、截面积为 Asx；
[0118] 6)纵肋横截面的中性轴距母板顶缘hsx；
[0119] 7)正交异性钢桥面板垂直于x轴截面的中性面距母板顶缘hx；
[0120] 8)母板的弯曲刚度为 纵肋的拉压刚度、弯曲刚度、扭转刚度分别为EAsx、EIsx、GJsx；
[0121] 9)横隔板及两侧腹板作为弹性支承考虑，其弹性约束刚度可由各自的几何、材料参数计算得到，横隔板的竖向拉压刚度及绕y轴的转动刚度分别用kda、 kdy表示，腹板的竖向拉压刚度及绕x轴的转动刚度分别用kwa、kwx表示，横隔板及腹板的其余刚度对结构的影响可以不予考虑。
[0122] (2)母板运用Kirchhoff薄板理论分析。
[0123] (3)分析中，考虑横截面的纵向焊接残余应力效应，由于正交异性钢桥面板的母板整体受压而纵肋整体受拉，因此，可视等效截面的母板均布受压而纵肋均布受拉，具体分析如下：
[0124] 1).如图3所示，令母板的均布纵向初始应力为σpx0、纵肋的均布纵向初始应力为σsx0，由初始应力自平衡条件得到：
[0125] σpx0b1tp+σsx0hstr＝0 (1)
[0126] 2).沿着正交异性钢桥面板的纵向，考虑母板与纵肋中的纵向焊接残余应力不变；
[0127] 3).正交异性钢桥面板的横向焊接残余应力效应较小，对结构的影响不予考虑。
[0128] (4).基于能量原理建立正交异性钢桥面板的方程，具体过程如下：
[0129] 根据Hooke定律，母板与加劲肋的应力与应变关系如下：
[0130]
[0131]
[0132]
[0133] σsx＝Eεsx (2‑4)[0134] 式(2‑1)至(2‑4)中：σpxx与σpyy分别为母板内x、y方向的正应力；τpxy为母板面内剪应力，σsx与σsy分别为x、y方向加劲肋的轴向应力；εpxx与εpyy分别为母板内x、y方向的正应变，γpxy为母板面内剪应力，εsx与εsy分别为x、y方向加劲肋的轴向应变。
[0135] 由于母板与纵肋的纵向初始应力分别为σpx0、σsx0，因此，根据Hooke定律，可得母板与加劲肋的初始应变为：
[0136]
[0137]
[0138]
[0139] 式(3‑1)至(3‑3)中：εpx0为母板的纵向初始应变；σpx0为母板的纵向初始应力；σpy0是母板的横向初始应力；εpy0为母板的横向初始应变；εsx0是纵肋的纵向初始应变；σsx0为纵肋的纵向初始应力。
[0140] 根据图2(a)‑图2(c)所示的加劲板，经过简单的几何运算，再结合式(2‑1)、 (2‑4)、(3‑1)、(3‑3)可得到纵肋中性轴及母板中性面的应变分别为：
[0141]
[0142]
[0143]
[0144]
[0145] 母板的应变能由弯曲应变能及中性面的拉压应变能组成：
[0146]
[0147] 纵肋的应变能为：
[0148]
[0149] 横隔板与腹板作为弹性支承，取其中心位移计算各自的等效应变能：
[0150]
[0151] 综上，正交异性钢桥面板的总应变能为：
[0152] U＝Up+Usx+Udw (8)
[0153] 对总应变能U进行变分运算，即可得刚度矩阵：
[0154] [K]＝[KP]+[Ksx]+[Kdw] (9)
[0155] 式中：[K]为刚度矩阵，[Kp]、[Ksx]、[Kdw]分别为母板、纵肋、横隔板与腹板对应的刚度矩阵，根据母板的横向位移w，以及正交异性钢桥面板的详细结构与材料参数，可得到刚度矩阵；Up为母板的应变能；Usx为纵肋的应变能；Udw为等效应变能；U为总应变能。
[0156] 母板的横向位移可以表示如下：
[0157]
[0158] 式中：w(x，y)表示母板的横向位移，wmn为广义坐标，对应于第(m,n)阶屈曲模态的幅值， 和φn(y)分别为沿着x与y方向的位移函数，表示如下：
[0159]
[0160] φn(y)＝An cos(any)+Bn sin(any)+Cn cosh(any)+Dn sinh(any) (12)[0161] 式中：常数(Am、Bm、Cm、Dm、am)根据x方向板的边界条件确定，常数(An、 Bn、Cn、Dn、an)根据y方向板的边界条件确定。
[0162] 对于x＝0及x＝a边，根据边界条件确定各参数，如：
[0163] 两边简支：
[0164] 两边固支：
[0165] 对于y＝0及y＝b边，两边由腹板支承，力的边界条件为：
[0166] y＝0边：支承力Qy＝kwaφn(n)、力矩
[0167] y＝b边：支承力Qy＝kwaφn(b)、力矩
[0168] 刚度矩阵的形式根据应变能及母板的横向位移确定，母板的刚度矩阵为：
[0169]
[0170] 其中，元素 的表达形式为：
[0171]
[0172]
[0173] 纵肋的刚度矩阵为：
[0174]
[0175] 其中，元素 的表达形式为：
[0176]
[0177] 横隔板与腹板的刚度矩阵为：
[0178]
[0179] 其中，元素 的表达形式为：
[0180]
[0181] (5)运用MATLAB软件可以方便地计算刚度矩阵，进而对正交异性钢桥面板进行详细地局部分析，以下简要列出主要程序段。
[0182] 1)定义结构参数
[0183] 定义变量：syms E G μ ρ a b a1 b1 tp tr hs Asx hsx hx kda kdy kwa kwx σpx0 σsx0[0184] 节点坐标：xyz＝[；]
[0185] 材料参数：mat＝[E；μ；ρ]
[0186] 3)计算刚度矩阵
[0187] 4)定义函数：funetionK_total＝Node_total(Ke)
[0188] 计算单元数量及整体刚度矩阵：
[0189] J＝length(Ke)；
[0190] K_total(1:1+J,1:1+J)＝0
[0191] Ke_J(1:2,1:2)＝0；
[0192] forI＝1:J
[0193] Ke_J＝Ke(I)
[0194] forN1＝1:2
[0195] forN2＝1:2
[0196] K_total(N1+I‑1,nN2+I‑1)＝K_total(N1+I‑1,N2+I‑1)+Ke_J{1,1}(N1,N2)；
[0197] end
[0198] K_total；
[0199] 实施例
[0200] 以一座单索面钢箱梁三塔斜拉桥为分析对象，主梁横截面为单箱三室钢箱梁，顶板宽38.8m、厚18mm，顶板纵肋为梯形肋，肋厚8mm、高300mm，梯形肋的上下口宽分别为300mm、180mm，梯形肋的横向间距为600mm。桥面板结构如图4所示。相应的材料及几何参数
3
为：ρ＝7850kg/m 、E＝206GPa、μ＝0.3、a＝6.4m、a1＝3.2m、b＝2.4m、b1＝600mm、b2＝300mm、b3＝180mm、hs＝280mm、tr＝8mm、tp＝18mm。
[0201] 计算考虑：考虑母板中的初始压应力为σpx0＝50MPa，分析母板中心承受横向集中力1KN作用下，沿着母板中心x方向6个点(x＝0,x＝0.64m,x＝1.28m,x＝1.92m, x＝2.56m,x＝3.2m)的竖向位移，并将其与ANSYS有限元计算结果对比验证，计算结果如下表1所示。
[0202] 表1母板中心各点竖向位移(mm)
[0203]位置 本发明方法 ANSYS计算值
x＝0 8.2 8.0
x＝0.64m 6.8 6.6
x＝1.28m 5.1 4.9
x＝1.92m 4.6 4.5
x＝2.56m 3.0 2.9
x＝3.2m 1.3 1.3
[0204] 从表1可以看出，本发明方法与ANSYS计算值非常接近，说明本发明方法比较准确；另外，相比于ANSYS等有限元计算方法，本发明无需建立精细有限元模型，能够节省运算时间；相比于其他常用的解析法及半解析法而言，则考虑的结构因素更为全面，从而使得计算结果更符合实际。
[0205] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
