除尘器箱体立柱在横向荷载作用下的抗剪强度计算方法与流程

本发明涉及一种除尘器箱体(亦可称壳体)立柱在横向荷载作用下的抗剪强度计算方法，属于结构技术领域。

背景技术

除尘器是广泛应用于火电、冶金、化工和建材等行业中以消除烟尘的主力环保装备，箱体是其中最重要的工艺部件。当箱体围护结构采用带加劲肋的平直钢板墙板—h形(或工字形)截面立柱结构体系时，h形截面立柱一侧翼缘与钢板墙板连续焊接连接，形成受力整体。箱体墙板主要直接承受由于内外温差而形成的空气负压(由外向内作用)和风荷载等横向荷载，由于立柱是墙板两侧的支承边界，这部分荷载会传递到立柱，使得立柱承担横向荷载。

横向荷载时箱体立柱抗剪强度计算方法与跨度、墙板宽度、加劲肋间距以及施加到墙板上的均布面荷载有关；所以对于除尘器箱体支承立柱进行结构设计时，需要有准确可靠的箱体墙板上直接作用横向荷载时箱体立柱抗剪强度计算方法。目前，工程设计中常用的立柱抗剪强度简化计算方法为：偏于简单和保守地将墙板视作均匀同性的单向受力板，其所受横向荷载均匀分配到相邻两侧立柱上，将立柱视作一有跨间支承的受均布线荷载的独立工作等截面多跨连续梁，其最大剪力vs可根据线弹性结构力学方法得到。该简化计算方法传力机制不够准确，除尘器箱体墙板并非单向传力板，应考虑墙板上加劲肋的传力作用，因此这种工程简化计算方法会造成计算剪力过大，对与计算得到的剪力过大会导致在设计时会使墙板宽度过大，进而会造成浪费。

鉴于以往计算方法的准确度不够，本发明对除尘器箱体立柱在横向荷载作用下的抗剪强度提出了一种更为准确可靠和考虑全面的计算方法。

技术实现要素：

本发明目的是针对现有除尘器箱体立柱在横向荷载作用下抗剪强度的工程简化计算方法存在的不足，提出一种受力模型更合理、结果准确度更高的抗剪强度计算方法，可以一定程度上优化立柱截面设计。本发明在推导形成过程中，依据实际工程的除尘器构造，适用范围涵盖了小、中、大型除尘器的几何尺寸。

一种除尘器箱体或壳体立柱在横向荷载作用下的抗剪强度计算方法；其步骤为：第一步，不考虑墙板与立柱的协同工作，将上、下相邻加劲肋和左、右两侧立柱围成的墙板区格作为一块四边简支受横向均布荷载作用的独立工作弹性板，求解边界横向反力；第二步，将各区格墙板边界横向反力反向作用到立柱，立柱作为一独立工作的等截面多跨连续梁求解其截面剪力；第三步，不考虑墙板与立柱协同受力的作用，计算立柱截面上最大剪应力，完成立柱抗剪强度验算。

在一种实施方式中，所述的除尘器箱体或壳体其结构为：除尘器箱体墙板为带有加劲肋的钢板，墙板上作用横向均布荷载，所述的横向是垂直于墙板方向，墙板与立柱一侧翼缘连续焊接连接，立柱为轧制h型钢或轧制i型钢或焊接h形截面，除尘器箱体内部对于立柱布置等间距的垂直墙板方向支撑。

在一种实施方式中，所述的上、下相邻加劲肋和左、右两侧立柱围成的墙板区格作为一块四边简支受横向均布荷载作用的独立工作弹性板，该墙板区格向左侧或者右侧立柱传递非均匀分布横向荷载qc；该墙板区格向上侧或下侧加劲肋边界传递横向荷载的合力为fs，该墙板区格向板件四个角点传递的横向集中力为-fc，该墙板区格通过上侧或者下侧的角钢加劲肋向左侧或者右侧的立柱传递集中力为f＝fs/2-fc。

在一种实施方式中，除尘器箱体墙板受横向荷载作用时立柱最大剪力发生在立柱顶部第一道跨间支撑处截面，中间立柱最大剪力vtm可以按照下式计算，边缘立柱最大剪力vte＝vtm/2：

vtm＝β[(n×a)·pb](1)

式中，n为一跨立柱范围内的墙板区格数量；

a为墙板上加劲肋的间距，单位为：mm；

b为相邻立柱间墙板的宽度，单位为：mm；

p为横向均布荷载值(由负压和风荷载等组成)，单位为：n/mm²；

β为立柱截面最大剪力的计算系数，见表1～表5：

表1两跨立柱截面最大剪力计算系数

表2三跨立柱截面最大剪力计算系数

表3四跨立柱截面最大剪力计算系数

表4五跨立柱截面最大剪力计算系数

表5六跨立柱截面最大剪力计算系数

在一种实施方式中，利用中间立柱控制截面最大弯矩mt，对于一跨范围内对应n个墙板区格的柱段，受到两侧墙板向立柱直接传递的非均匀分布荷载2qc和在其加劲肋与立柱连接点位置向立柱传递的集中力4f作用，推导得到中间立柱在横向荷载作用下顶部第一道跨间支承位置最大剪力理论计算公式为：

在一种实施方式中，中间立柱截面最大剪力表达式改写为：

式中，l为立柱一跨长度，单位为：mm；

i为中间过程系数(i＝0,1,2,3,4,5,6……)；

x为分布荷载计算微元段到计算截面的距离，公式(3)能够方便编制计算机程序对剪力进行数值计算。

在一种实施方式中，本发明在计算得到立柱截面最大剪力后，不考虑墙板的抗剪作用，最大剪应力按照h形截面(或i形截面)立柱独立承担计算，利用立柱控制截面最大剪力vt，当立柱为中间立柱时所述的vt为vtm，当立柱为边缘立柱时所述的vt为vte，立柱截面最大剪应力τt可以利用以下公式求得；ζ为立柱剪应力修正系数，取值为1.1。

式中，smax为立柱h形截面(或i形截面)的最大面积矩，单位为：mm³；

iy为立柱h形截面(或i形截面)的惯性矩，单位为：mm⁴；

t1为h形截面(或i形截面)立柱的腹板厚度，单位为：mm。

在一种实施方式中，利用边缘立柱控制截面最大弯矩mt/2，对于一跨范围内对应n个墙板区格的边缘立柱柱段，受到仅一侧墙板向立柱直接传递的非均匀分布荷载qc和在其加劲肋与立柱连接点位置向立柱传递的集中力2f作用，边缘立柱在横向荷载作用下顶部第一道跨间支承位置最大剪力理论计算公式为：

其中，l为一跨立柱的长度，单位为：mm；x为分布荷载计算微元段到计算截面的距离。

在一种实施方式中，边缘立柱在横向荷载作用下顶部第一道跨间支承位置最大剪力计算公式为：

其中，i为中间过程系数，其中i＝0,1,2,3,4,5,6……，公式(19)能够方便编制计算机程序对剪力进行数值计算。

借助以上技术方案，本发明的有益效果为：

1、准确性高，考虑因素全面：考虑墙板上加劲肋的传力作用，传力机制更加准确；较目前通行的工程简化计算方法更为准确可靠。

2、使用方便：根据本发明提出的计算公式和计算表格可直接求得除尘器箱体墙板受横向荷载作用时中间立柱最大剪力值，取立柱h形截面为受剪截面，依据公式(4)可方便准确地进行截面抗剪强度计算。

附图说明

图1为本发明的结构模型与位移坐标系。

图2为本发明中四边简支板边界反力分布及坐标系图。

图3为本发明中横向荷载传递机制示意图。

图4为本发明中立柱受到一侧墙板传递的横向荷载作用时受力简图。

图5为本发明中阶数上限m’取不同值时各项荷载结果比较。

图6为本发明中中间立柱最大剪力计算受力简图。

具体实施方式

下面结合附图并通过具体实施例，对本发明中的技术方案进行完整详细的描述，以进一步说明本发明的技术方案特征及其形成过程。可以理解的是，此处所描述的实施例仅仅用于解释本发明，而非对本发明的限定。

计算墙板所受横向荷载分配传递至立柱后引起的立柱最大剪力。

除尘器箱体墙板-立柱结构体系如附图1所示。除尘器箱体墙板直接承受横向均布荷载p作用，横向荷载会分配传递到墙板两侧立柱，使立柱受到剪力。本发明第一步计算横向荷载向立柱的传递，不考虑墙板与立柱的协同工作，将上、下相邻加劲肋和左、右两侧立柱围成的墙板区格作为一块四边简支受横向均布荷载作用的独立工作弹性板，求解边界横向反力。四边简支板宽度为一跨墙板宽度b，其高度为墙板加劲肋间距a，墙板厚度为t，一个墙板区格受力如附图2所示。负压和风荷载组成的横向均布荷载p直接作用到箱体每一个墙板区格上，继而向左、右两边的立柱和上、下两端的加劲肋传递。墙板区格左、右两侧立柱提供的横向分布边界反力为vy，将其反向作用到立柱上作为立柱直接承受的横向线荷载qc。墙板区格上、下端加劲肋提供的横向边界反力为vx，其合力大小为fs，加劲肋与两边立柱通过连接板连接，将一侧加劲肋提供边界反力的合力fs/2反向作用到加劲肋与立柱连接点，作为各墙板区格上、下端边界传递的集中横向荷载。此外在各墙板区格的四个角点上还有集中反力r，也将其反向作用到加劲肋与立柱连接点，作为墙板区格角点向立柱传递的集中横向荷载fc。横向荷载传递机制如附图3所示。立柱沿高度方向设有等间距的横向支撑，为其提供垂直墙板方向的约束，立柱跨度即为支撑间距l，每一跨对应有n个墙板区格。以三跨中间立柱为例，受到单侧墙板传递的横向荷载时受力简图如附图4所示。

一块墙板区格的边界横向反力采用对边简支矩形板的单三角级数列维解。根据常规除尘器结构几何尺寸，墙板宽b大于加劲肋间距a。墙板区格上、下端加劲肋提供的横向边界反力vx按下式计算：

式中，m为级数的阶数，取1,3,5…的奇数；系数αm、am’和bm’分别按下式计算：

则对应一块墙板区格的一端加劲肋传递到立柱的横向集中力为：

一块墙板区格(0≤x≤a)传递到一侧立柱上的横向分布荷载qc在数值上等于墙板区格左侧或者右侧边界横向反力vy，按下式计算：

各墙板区格角点向立柱传递的集中横向荷载-fc在数值上等于墙板区格四个角点上存在的集中反力-r，按下式计算：

上述解答均为无穷级数形式，在实际工程计算中，m取值总有一个上限。若m上限值取太大，则计算效率低下；若m上限值取太小，则结果不够精确。本发明采用matlab语言编写计算程序，通过对m取不同上限值时各项计算结果的数值计算，比较分析得到m上限值的影响，继而确定合理的m上限值。

取算例墙板区格a＝1000mm，均布荷载p＝9000pa，计算在m上限取值m’变化时各项荷载结果分别如附图5中(a)、(b)、(c)所示。

图5(a)中m上限取值不同情况的曲线几乎重合，表明改变解答中级数的求和阶数对于墙板区格上、下侧边界反力的合力fs大小几乎没有影响。图5(b)表明，墙板区格宽高比变化时，m上限值取较小时，墙板区格立柱侧边界最大横向反力vy,max差异较大；当m上限值取大于等于21时，数值计算结果趋于相等。图5(c)表明，当m上限值取大于等于21时，墙板区格角点处集中反力r的数值计算结果也趋于相等。以较为常规的除尘器箱体墙板区格高宽比b/a＝4的情况为例，当m上限取值m’＝21增大至m’＝201时，墙板区格上、下侧边界反力的合力fs由16.413kn改变为16.412kn，降低0.01％；墙板区格立柱侧边界最大横向反力vy,max由4.515kn/m改变为4.51kn/m，降低0.11％；墙板区格角点处集中反力r由0.8543kn改变为0.8547kn，提高0.05％。综合考虑数值计算效率和精度，本发明在后续计算对于各项解答均取m上限值为21。

立柱各处截面内力求解时可以将其视为一等截面多跨连续梁。计算表明，立柱上第一道跨间支承位置有最大剪力值，因此对于等截面立柱而言，强度设计时只需要验算此处截面的抗剪强度，后续分析主要针对立柱第一道跨间支承位置的剪力。

本发明需要首先计算中间立柱第一道跨间支承截面弯矩mt，在此基础上推导得到其最大剪力vt；边缘立柱最大内力值按照中间立柱值减半即可。mt推导计算过程如下：

考虑对于各跨立柱，其几何尺寸以及所受荷载的形式与大小都相同，且支承处无侧移，因此采用力矩分配法求解立柱控制截面最大弯矩。所需解决的是一跨段内立柱在荷载作用下的固端弯矩，包括远端为固定和远端为铰支两种约束情况。

取一跨度为l的除尘器箱体中间立柱，假设其两端固定，跨间对应n个墙板区格，受到一侧墙板向立柱直接传递的分布荷载qc和一侧墙板区格在其加劲肋边界与立柱连接点位置向立柱传递的集中力2f作用。根据平衡关系和位移协调条件，可以推导得到该跨立柱在横向荷载作用下，其固端弯矩如下式所示：

f＝fs/2-fc(13)

式中，k为所计算的微元段至a端所包括的墙板区格数目。

对于一跨一端固定另一端铰支的立柱，可推导得到其固端弯矩

采用matlab语言对固端弯矩进行数值计算时，由于式(12)和式(14)中的表达式带有绝对值项，无法进行直接积分求解，因此将其进行变换，两端固定情况的立柱固端弯矩表达式改写为：

式中，i为中间过程系数。

一端固定另一端铰支情况的固端弯矩表达式改写为：

立柱各跨段线刚度相等，根据总跨数和远端约束情况即可确定各跨间支承位置截面两侧的弯矩分配系数，继而可以求得立柱各跨间支承位置截面的弯矩。计算表明，立柱上第一道跨间支承位置有最大弯矩值。上述计算得到的弯矩只考虑了立柱单侧墙板传递的横向荷载，实际上箱体中间立柱两侧都有墙板，要考虑两侧墙板传递来的荷载，将此结果乘以2，即可求解出中间立柱截面最大弯矩值mt。

对于跨间对应n个墙板区格的柱段，受到两侧墙板向立柱直接传递的非均匀分布荷载2qc和在其加劲肋与立柱连接点位置向立柱传递的集中力4f作用，受力简图如附图6所示。根据平衡关系和位移协调条件，可以推导得到中间立柱在横向荷载作用下顶部第一道跨间支承位置最大剪力vtm如下式所示：

式中，l为立柱一跨长度，单位为：mm；

采用matlab语言对立柱控制截面剪力进行数值计算时，由于式(2)中的表达式带有绝对值项，无法进行直接积分求解，因此将其进行变换，中间立柱截面最大剪力表达式改写为：

式中，i为中间过程系数。

下述实施例考虑除尘器箱体立柱的跨数、每跨立柱范围内对应墙板区格数目以及墙板区格宽高比后总结归纳制定立柱控制截面的最大剪力计算方法。

实施例1：

除尘器箱体立柱为两跨立柱；立柱每跨间的墙板区格数量n为1；墙板上加劲肋的间距a为600mm，一跨墙板的宽度b为600mm，即b/a＝1；立柱跨度l＝n×a＝600mm；负压和风荷载等横向均布荷载的值p为0.008n/mm²。根据前述计算公式，用matlab语言编程计算出的中间立柱截面最大剪力vtm如表1所示，用上述方法计算vtm值时与立柱截面尺寸无关。

实施例2～实施例54：

实施例2～实施例54相对于实施例1仅改变立柱每跨间的墙板区格数量n和一跨墙板的宽度b，继而改变了b/a的大小，具体构造参数及用matlab语言编程计算出的中间立柱截面最大剪力vtm如表1所示。

实施例55：

除尘器箱体立柱为两跨立柱；立柱每跨间的墙板区格数量n为1；墙板上加劲肋的间距a为1000mm，一跨墙板的宽度b为1000mm，即b/a＝1；立柱跨度l＝n×a＝1000mm；负压和风荷载等横向均布荷载的值p为0.009n/mm²。用matlab语言编程计算出的中间立柱截面最大剪力vtm如表1所示。

实施例56～实施例108：

实施例56～实施例108相对于实施例55仅改变立柱每跨间的墙板区格数量n和一跨墙板的宽度b，继而改变了b/a的大小，具体构造参数及用matlab语言编程计算出的中间立柱截面最大剪力vtm如表1所示。

表1两跨除尘器箱体立柱结构在横向荷载作用下最大剪力计算值

通过分析比较两跨结构实施例中间立柱截面最大剪力vtm，发现中间立柱截面最大剪力主要与其跨度、墙板宽度、加劲肋间距以及施加到墙板上的均布面荷载有关，对结果的统计归纳表明，vtm/[(n×a)pb]之值仅与n和b/a相关，因此根据立柱跨数、n、a、b和荷载p即可得到最大剪力vtm。

实施例109～实施例216：

实施例109～实施例216相对于实施例1～实施例108仅改变立柱为三跨立柱，其余参数均不改变。具体构造参数及用matlab语言编程计算出的中间立柱截面最大剪力vtm如表2所示。

表2三跨除尘器箱体立柱结构在横向荷载作用下最大剪力计算值

实施例217～实施例324：

实施例217～实施例324相对于实施例1～实施例108仅改变立柱为四跨立柱，其余参数均不改变。具体构造参数及用matlab语言编程计算出的中间立柱截面最大剪力vtm如表3所示。

表3四跨除尘器箱体立柱结构在横向荷载作用下最大剪力计算值

实施例325～实施例432：

实施例325～实施例432相对于实施例1～实施例108仅改变立柱为五跨立柱，其余参数均不改变。具体构造参数及用matlab语言编程计算出的中间立柱截面最大剪力vtm如表4所示。

表4五跨除尘器箱体立柱结构在横向荷载作用下最大剪力计算值

实施例433～实施例540：

实施例433～实施例540相对于实施例1～实施例108仅改变立柱为六跨立柱，其余参数均不改变。具体构造参数及用matlab语言编程计算出的中间立柱截面最大剪力vtm如表5所示。

表5六跨除尘器箱体立柱结构在横向荷载作用下最大剪力计算值

本发明研发过程中计算了以上实施例中间立柱截面最大剪力vtm，对结果的统计归纳表明，采用本发明计算方法得到的立柱截面最大剪力，其主要影响参数有立柱的跨数、跨度l＝n×a、墙板宽度b、墙板加劲肋间距a以及施加到墙板上的均布面荷载p。由于立柱承受荷载作用形式是一定的，为便于工程计算应用，本发明引入除尘器箱体墙板受横向荷载作用时中间立柱截面最大剪力的计算系数β，利用公式(1)，可求得中间立柱截面最大剪力的理论计算值。根据除尘器箱体立柱的跨数、每跨立柱范围内对应墙板区格数目以及墙板区格宽高比编制成计算系数β表格，如表6-表10所示。

表6两跨立柱截面最大剪力计算系数

表7三跨立柱截面最大剪力计算系数

表8四跨立柱截面最大剪力计算系数

表9五跨立柱截面最大剪力计算系数

表10六跨立柱截面最大剪力计算系数

下述实施例验证了本发明立柱截面最大剪力计算方法的准确可靠性。

该实施例中使用有限元计算分析方法对本发明中的计算方法进行验证。有限元方法是一种利用计算机高效计算手段以用于求解复杂固体力学问题的数值分析方法，其准确性较高，但是建模、网格划分、加载和求解过程非常复杂，需要耗费大量时间精力，且操作难度较大，不便于复杂结构问题的工程普遍应用。因此，工程上需要准确可靠而简单便利的计算方法。

依据实际工程的除尘器构造，构建了5例箱体墙板-立柱结构体系计算模型，涵盖了小、中、大型除尘器的几何尺寸，可以使研究结果具有普遍参考意义。有限元方法可以准确反映除尘器箱体墙板-立柱结构体系中的内力分布情况，将有限元计算得到的立柱截面最大剪力值vfem认定为真值，与本发明理论计算方法求出的立柱控制截面最大剪力vt进行比较，验证本发明提出立柱截面最大剪力计算方法的准确性。有限元计算分析过程说明如下：

(1)定义单元：所有结构部件均采用shell181单元模拟。

(2)定义材料：由于横向荷载对结构体系的作用效应较小，变形和应力水平均较低，因此进行线弹性计算。制作除尘器一般采用q235钢材，其屈服强度fy＝235mpa，弹性模量e＝2.06×105mpa，泊松比ν＝0.3。

(3)施加约束情况：除尘器箱体墙板顶端与箱体加劲顶板连接，因此在墙板顶部边界施加垂直墙板方向(z向)的平动约束。墙板底端与灰斗加劲壁板连接，因此在墙板底端边界施加垂直墙板方向的平动约束。立柱受等间距布置的横向支撑(垂直于墙板方向)约束，在立柱与横向支撑连接处施加垂直墙板方向的平动约束。在中间立柱柱底施加三个方向的平动约束。由于箱体内烟气往往是高温，为了释放温度变形，两侧边缘立柱底部仅施加沿墙板高度方向(x向)和垂直于墙板方向的约束，以实现结构在墙板平面内(y向)可以伸缩变形。

(4)施加荷载情况：除尘器箱体墙板在运行过程中受到风荷载以及内外压差(负压)产生垂直于墙板的横向均布荷载，在墙板上施加横向均布荷载0.009mpa。

实施例541：

除尘器箱体墙板厚度t为6mm，墙板宽度b为3500mm，角钢加劲肋间距a为1170mm，立柱每跨内加劲肋区格数量n为3，立柱横向支撑间距为l为3510mm，立柱为三跨立柱，立柱总高度h为11990mm，立柱截面为h200mm×150mm×6mm×9mm，立柱截面面积a为3846mm²，立柱截面惯性矩iy为3.1×10⁷mm⁴，立柱截面模量wh为2.9×10⁵mm³。有限元方法计算出中间立柱截面最大剪力值vfem与本发明提出的剪力计算方法计算得到中间立柱截面最大剪力值vtm分别如表11所示。

实施例542：

除尘器箱体墙板厚度t为6mm，墙板宽度b为3850mm，角钢加劲肋间距a为1126mm，立柱每跨内加劲肋区格数量n为4，立柱横向支撑间距为l为4504mm，立柱为三跨立柱，立柱总高度h为14972mm，立柱截面为h250mm×175mm×7mm×11mm，立柱截面面积a为5523mm²，立柱截面惯性矩iy为6.8×10⁷mm⁴，立柱截面模量wh为5.2×10⁵mm³。有限元方法计算出中间立柱截面最大剪力值vfem与本发明提出的剪力计算方法计算得到中间立柱截面最大剪力值vtm分别如表11所示。

实施例543：

除尘器箱体墙板厚度t为7mm，墙板宽度b为4200mm，角钢加劲肋间距a为1000mm，立柱每跨内加劲肋区格数量n为5，立柱横向支撑间距为l为5000mm，立柱为三跨立柱，立柱总高度h为16460mm，立柱截面为h250mm×250mm×9mm×14mm，立柱截面面积a为9124mm²，立柱截面惯性矩iy为1.2×10⁸mm⁴，立柱截面模量wh为9.0×10⁵mm³。有限元方法计算出中间立柱截面最大剪力值vfem与本发明提出的剪力计算方法计算得到中间立柱截面最大剪力值vtm分别如表11所示。

实施例544：

除尘器箱体墙板厚度t为7mm，墙板宽度b为4030mm，角钢加劲肋间距a为1040mm，立柱每跨内加劲肋区格数量n为5，立柱横向支撑间距为l为5200mm，立柱为三跨立柱，立柱总高度h为13940mm，立柱截面为h294mm×200mm×8mm×12mm，立柱截面面积a为7056mm²，立柱截面惯性矩iy为1.2×10⁸mm⁴，立柱截面模量wh为7.8×10⁵mm³。有限元方法计算出中间立柱截面最大剪力值vfem与本发明提出的剪力计算方法计算得到中间立柱截面最大剪力值vtm分别如表11所示。

实施例545：

除尘器箱体墙板厚度t为7mm，墙板宽度b为4550mm，角钢加劲肋间距a为1000mm，立柱每跨内加劲肋区格数量n为5，立柱横向支撑间距为l为5000mm，立柱为四跨立柱，立柱总高度h为21460mm，立柱截面为h300mm×300mm×10mm×15mm，立柱截面面积a为11850mm²，立柱截面惯性矩iy为2.2×10⁸mm⁴，立柱截面模量wh为1.4×10⁶mm³。有限元方法计算出中间立柱截面最大剪力值vfem与本发明提出的剪力计算方法计算得到中间立柱截面最大剪力值vtm分别如表11所示。

表11有限元方法和本发明计算方法得到立柱截面最大剪力对比

考察比较实施例组541、542、543、544、545，本发明提出的计算方法得到的剪力值和有限元方法得到的剪力值两者相对误差在2％以内，是准确可靠的。有限元方法读取的是立柱截面的剪力值，本发明方法不考虑墙板-立柱的协同工作效应计算出的剪力值和有限元方法得到的剪力值差异很小，表明墙板在抵抗横向荷载引起的剪力时贡献不大，且立柱抗剪主要是由h型钢腹板来承担，翼缘的抗剪作用较弱，可以视作立柱一侧翼缘延伸的墙板对立柱的抗剪性能影响不大，因此立柱内力理论计算方法的剪力解答是准确可靠的。

实施例组541、542、543、544、545在计算得到立柱截面最大剪力后，最大剪应力按照h形截面立柱独立承担计算，不考虑墙板的抗剪作用，根据材料力学基本理论，立柱截面最大剪应力按照式(17)计算，材料力学基本理论计算方法算出中间立柱截面最大剪应力值τ与有限元方法计算出中间立柱截面最大剪应力值τfem分别如表12所示。

式中，vt对于中间立柱取vtm，对于边缘立柱取vte。

smax为立柱h形截面的最大面积矩，单位为：mm³；

iy为立柱h形截面的惯性矩，单位为：mm⁴；

t1为h形截面立柱的腹板厚度，单位为：mm。

表12有限元方法和材料力学计算方法得到截面最大剪应力对比

由表12可以发现，材料力学基本理论计算方法与有限元方法计算出中间立柱截面最大剪应力值有一定差异，一是因为实际承担剪力的是立柱与墙板的组合结构，并非单独立柱截面；二是因为在控制截面位置，有连接板等构造，使得应力分布相对复杂。因此本发明以有限元计算剪应力为真值，在式(17)基础上提出立柱剪应力修正系数ζ，取ζ＝1.1，立柱截面最大剪应力按下式计算：

依据本发明提出的立柱截面最大剪应力计算公式(4)算出立柱截面最大剪应力值τt与有限元方法计算出中间立柱截面最大剪应力值τfem对比如表13所示。

表13有限元方法和本发明计算方法得到截面最大剪应力对比

考察比较实施例组541、542、543、544、545用有限元方法计算出的中间立柱截面最大剪应力值与利用本发明提出的剪应力计算公式计算得到的中间立柱截面最大剪应力值，两者相对误差在6％以内，因此本发明提出的立柱截面最大剪应力计算公式(4)是准确可靠的。

虽然本发明已以较佳实施例公开如上，但其并非用以限定本发明，任何熟悉此技术的人，在不脱离本发明的精神和范围内，都可做各种的改动与修饰，因此本发明的保护范围应该以权利要求书所界定的为准。