格林积分在多边形截面特性计算的应用

在工程领域中，截面面积、形心、静矩、惯性矩等截面特性的计算十分广泛，目前工程设计或计算手册一般只给出简单截面的截面特性计算公式或常用截面的截面特性表，使用上有一定的局限性。

文中根据格林积分公式，推导任意多边形的截面特性计算公式及任意闭区域截面特性计算思路，并介绍该方法在理正结构工具箱软件中的应用。

1. 1. 基本原理

根据高等数学格林积分公式[1]可知，如图1所示，平面闭区域D的二重积分可通过沿闭区域D的边界曲线L上的积分表达。

图 1 平面闭区域D的二重积分

Fig. 1 Double integral of plane closed region D

1. 2. 平面闭区域截面特性积分式

根据平面闭区域截面特性的物理意义[2]

面积可表示为：

静矩可表示为：

形心可表示为：

惯性矩可表示为：

1. 3. 直边多边形截面特性数值积分

平面任意直边多边形，如图2所示，根据平面闭区域截面特性积分并结合格林积分公式，可将n边多边形截面特性积分表示为多边形每边计算，则n边多边形

图 2 平面直线边多边形

Fig. 2 Straight edge polygon

面积可表示为：

静矩可表示为：

形心可表示为：

惯性矩可表示为：

若要计算形心惯性矩可先将坐标原点移至形心点，再利用上式计算即可。

1. 4. 平面闭区域截面特性数值积分

对于曲线边闭区域截面特性计算，可根据计算精度要求，先将曲线边离散成多段直线段表示，用离散后的多边形近似求解。

矩形独立基础在荷载偏心比较大的情况下，基础与土体会脱空，地基应力需要按重分布进行计算，如图3所示。

图 3 基底应力重分布计算

Fig. 3 Calculation of stress redistribution of basement

基底应力重分布后，基底受力面已不是矩形，可能存在五边或三边的多边形状态，需要根据重分布后的基底受力面截面特性计算基底应力。图3为理正结构工具箱独立基础模块，利用文中方法计算重分布后基底压力分布情况。

文中方法可精确解出直边多边形截面特性，并能对曲边多边形截面特性离散求解近似解，适应各种截面形式的截面特性计算，尤其合适计算机编程实现，扩大了工程应用范围。

参考文献：

[1] 同济大学数学系. 高等数学（第六版）[M]. 北京：高等教育出版社，2007.4

[2] 孙训方，方孝淑，关来秦. 材料力学（第四版）[M]. 北京：高等教育出版社，2002.8