一种基于群论的对称索杆结构找形方法与流程

本发明属于空间结构、对称索杆结构的形态分析、设计与开发等领域，涉及一种对称索杆结构找形方法。

背景技术：

在传统的结构中，结构是从几何和材料中获取刚度的，这种主要从几何和材料中获取刚度的结构即所谓的刚性结构。而索杆结构的主要从预应力中获取刚度，确定索杆结构的几何形状和预应力状态的分析就是所谓的找形分析，也称为初始形态分析，其中形为初始几何形状，态则是初始预应力态。

现有的找形方法大多针对自平衡的索杆体系，而在土木工程中应用的结构大多有约束条件。对于有支座边界约束条件的索杆结构，一些学者通过加虚拟杆件的方式将支座约束转化为杆件约束，从而将结构转化成自平衡的张拉整体结构进行分析，但是虚拟杆的设置方式不同对结果的影响比较大。

对于对称索杆结构，充分利用固有对称性可显著简化其找形分析过程。群集理论作为一种系统分析对称性的重要数学工具，结合对称群的对称操作及不可约表示建立结构的对称坐标系，将笛卡尔坐标系下的稀疏、带状分布的相关矩阵转换成对称坐标系下的分块、对角化矩阵，将显著提高张拉整体结构找形分析的计算效率。近年来，一些学者已成功将群集理论引入张拉整体结构的找形分析中，显著提高了张拉整体结构找形分析的计算效率。

技术实现要素：

技术问题：本发明提供一种针对任意有支座边界约束的索杆结构，能准确、迅速地计算其杆件力密度和对应节点坐标的一种对称索杆结构找形方法。

技术方案：本发明针对任意有支座边界约束的索杆结构，确定其所属对称群后，求出考虑所有节点和考虑自由节点的与结构x/y/z方向刚体平动相关联的不可约表示的对称子空间，根据这些不可约表示对应的力密度分块子矩阵的秩亏条件及零空间，直接求得结构的可行力密度及节点坐标。

本发明的基于群论的对称索杆结构找形方法，包括以下步骤：

步骤1判断结构所属对称群，根据对称性对结构杆件进行分组；

步骤2根据结构在各不同对称操作下节点之间的转换关系，确定结构在所属对称群下所有独立对称操作对应的所有节点转换的置换矩阵pg和自由节点转换的置换矩阵p＇g，其中g表示结构所属对称群所有独立对称操作的集合；

步骤3根据结构x/y/z方向刚体平动相关联的不可约表示和所述步骤2确定的置换矩阵pg，求得所有节点对应于x/y/z方向坐标的转换矩阵v^(μ)i，根据结构x/y/z方向刚体平动相关联的不可约表示和所述步骤2确定的置换矩阵p′g，求得自由节点对应于x/y/z方向坐标的转换矩阵v^′(μ)i；

其中，μ表示与结构x/y/z方向坐标对应的不可约表示，i表示取不可约表示μ的第i行，d＝c^t*diag(q)*c，c^t表示拓扑矩阵c的转置，q为以符号表示的各杆件力密度列向量，diag()表示矩阵的对角线元素为括号中的元素，其余元素为0；

步骤5根据三个x/y/z方向坐标对应的力密度分块子矩阵和不可约表示a1对应的力密度分块子矩阵的总秩亏不小于d+1分析各力密度分块子矩阵，求得各组杆件力密度之间的解析关系，其中d为待求解结构的维度；

步骤6求得分别与x/y/z方向坐标对应的三个力密度分块子矩阵的零空间，即为对称坐标系下节点在x/y/z方向的的坐标然后根据下式求出笛卡尔坐标系下结构的节点坐标模态：

步骤7根据设计要求中给出的固定节点，选取其中任意两个固定节点a、b，找出所述步骤6求出的节点的x/y/z方向的坐标中与选定的固定节点a、b相应的节点a⁰、b⁰的坐标，然后根据下式求出比例放大系数：

根据下式求出结构实际构形下的坐标x^0′，y^0′，z^0′：

x^0′＝coe*x⁰

y^0′＝coe*y⁰

z^0′＝coe*z⁰

进行整体平移，即根据下式求出结构的n个自由节点在笛卡尔坐标系下的真实坐标xj，yj，zj，其中n为结构中自由节点总数，j＝1，2，3，…，n：

其中，xf，yf，zf为设计要求中给出的任意一个固定节点的坐标，为对应的该固定节点实际构型下的坐标，f为整数且f∈[1，nf]，nf为结构固定节点总数。

进一步的，本发明方法中，步骤3中，所有节点对应于x/y/z方向坐标的转换矩阵v^(μ)i按照如下方式求得：首先将不同对称操作下与x/y/z方向刚体平动相关联的不可约表示和对应于所有节点转换的置换矩阵pg求乘积，将这些乘积求和得出所有节点对应于x/y/z方向坐标的三个投影算子矩阵，然后根据所述三个投影算子矩阵的列空间分别求出所有节点对应于x/y/z方向坐标的三个转换矩阵v^(μ)i。

进一步的，本发明方法中，步骤3中，自由节点对应于x/y/z方向坐标的转换矩阵v^′(μ)i按照如下方式求得：首先将不同对称操作下与x/y/z方向刚体平动相关联的不可约表示和对应于自由节点转换的置换矩阵p′g求乘积，将这些乘积求和得出自由节点对应于x/y/z方向坐标的三个投影算子矩阵，然后根据所述三个投影算子矩阵的列空间分别求出自由节点对应于x/y/z方向坐标的三个转换矩阵v^′(μ)i。

进一步的，本发明方法中，步骤5中，按照如下方式求出各组杆件力密度之间的解析关系：与x/y/z方向坐标对应的三个力密度分块子矩阵中，非方阵自动满足秩亏要求，无需进行计算，方阵则令三个力密度分块子矩阵行列式等于0，联立三个力密度分块子矩阵对应的组成方程组，解该方程组，即得到各组杆件力密度之间的解析关系，其中，方阵为行数与列数相等的矩阵，非方阵为行数与列数不相等的矩阵。

进一步的，本发明方法中，步骤6中，按照以下方式求出与x/y/z方向坐标对应的三个力密度分块子矩阵的零空间：将符合所述步骤5求得的各组杆件力密度之间的解析关系，并且使拉索受拉为正值、压杆受压为负值的数值分别代入与x/y/z方向坐标对应的三个力密度分块子矩阵，进行零空间求解，得到零空间。-

本发明在分析和研究现有将群论引入对称张拉整体结构找形分析的基础上，将应用于自平衡张拉整体结构的找形方法引入有支座边界约束条件的对称索杆结构的找形问题中，把结构节点分为自由节点和有约束节点，从所有节点和自由节点两个角度出发，结合与结构x/y/z方向刚体平动相关联的不可约表示，提出了一种新的对称索杆结构找形方法，大大简化了对称索杆结构的找形分析。

有益效果：本发明与现有技术相比，具有以下优点：

本发明的优点在于将有边界支座约束的预应力索杆结构，从所有节点和自由节点两个角度出发，充分利用与结构x/y/z方向刚体平动相关联的不可约表示求解力密度分块子矩阵，直接对维度较小的力密度分块子矩阵进行分析。现有的方法将笛卡尔坐标系下的稀疏、带状分布的相关矩阵转换成对称坐标系下的分块、对角化矩阵时需要计算所有不可约表示对应的对称子空间，求出所有投影算子矩阵及其列空间组成整体转换矩阵，将力密度矩阵转换成分块、对角化矩阵时不能区分与x/y/z方向相关联或无关分块子矩阵，需要进行整体力密度矩阵秩的分析，对于高对称结构和杆件数较多的结构，需要计算大量投影算子矩阵和坐标转换矩阵，且进行力密度矩阵秩的分析时矩阵维度较大，计算量较大，因而计算效率较低。且现有找形方法主要针对自平衡的张拉整体结构，对有支座约束的索杆结构不适用。本发明所述方法首先将结构节点分为自由节点和约束节点，从结构所有节点和自由节点两个角度出发，明确与结构x/y/z方向刚体平动相关联的不可约表示，求出两组节点相应的投影算子矩阵，继而求出两组与节点坐标相关联的对称子空间的基向量，得到对应的力密度分块子矩阵，分析分块子矩阵的秩亏条件及零空间，直接求得结构的可行力密度及节点坐标。对于任何结构最多只需求出六个对称子空间和三个分块力密度矩阵，避免了求对称群中所有对称操作下的转换矩阵，减少了程序的计算量和复杂度，显著提高了索杆结构找形分析的计算效率，且本发明所述方法不仅适用于有约束的索杆结构，对于自平衡张拉整体结构同样适用。

附图说明

图1为基于群论的对称索杆结构找形方法流程示意。

图2为一个c2v对称张拉整体结构。

图2中，细实线表示拉索构件，粗实线表示压杆构件。第一节点a、第二节点b、第三节点c、第四节点d表示结构的四个自由节点，第五节点e、第六节点f表示结构的两个约束节点；第一拉索1连接第一节点a和第五节点e、第二拉索2连接第二节点b和第六节点f、第三拉索3连接第三节点c和第六节点f、第四拉索4连接第四节点d和第五节点e，第五拉索5连接第一节点a和第二节点b，第六拉索6连接第三节点c和第四节点d，第一压杆7连接第一节点a和第四节点d，第二压杆8连接第二节点b和第三节点c，第一拉索1、第二拉索2、第三拉索3、第四拉索4为斜向拉索，第五拉索5、第六拉索6为水平拉索。

具体实施方式

下面结合实施例和说明书附图对本发明作进一步的说明。

基于matlab编程平台，对图2结构采用本发明所述方法进行求解。图2所示结构有6个节点，8根杆件，其中节点1、2、3、4为自由节点，节点5、6为约束节点，属于c2v对称结构，根据对称性，将杆件分为三组，水平拉索、斜向拉索和压杆，分别用th、tv和td表示水平拉索、斜向拉索和压杆的力密度。

如图2所示结构有四个对称操作：恒等变换(e)，绕z轴旋转180°(c2)、相对于水平面的镜像对称(σv)和相对于竖向平面的镜像对称(σ′v)。

当引入对称群的相关概念后，计算基于的坐标系自动从笛卡尔坐标系转换为对称坐标系。从结构所有节点的角度出发，在对称操作c2(即绕z轴旋转180°)作用下，第一节点a转换到原先第三节点c的位置，第二节点b转换到原先第四节点d的位置，第五节点e转换到原先第六节点f的位置，第三节点c转换到原先第一节点a的位置，第四节点d转换到原先第二节点b的位置，第六节点f转换到原先第五节点e的位置，可写成如下形式：

同理，可写出结构所有节点在其他对称操作下的置换矩阵。

从结构自由节点的角度出发，在对称操作c2(即绕z轴旋转180°)作用下，第一节点a转换到原先第三节点c的位置，第二节点b转换到原先第四节点d的位置，第三节点c转换到原先第一节点a的位置，第四节点d转换到原先第二节点b的位置，可写成如下形式：

同理，可写出结构自由节点在其他对称操作下的置换矩阵。

c2v对称结构不可约表示有4个：a1、a2、b1、b2，其节点坐标z对应a1表示的一维子空间x对应于对称子空间y对应于对称子空间完整的力密度矩阵由四个小矩阵块对角形成，这里只需求解对称子空间和对应的三个转换小矩阵

则

同理可求出和

同理，可求出节点坐标x，y对应的力密度分块子矩阵：

结构维度为2，完整的力密度矩阵需满足秩亏至少为d+1＝3的要求，分析上述三个力密度分块子矩阵，使得每个矩阵分别有至少一个秩亏，和均是维度1′2的非方阵矩阵，自动满足秩亏要求，无需进行计算，计算方阵的行列式等于

即可得出各个单元力密度之间的解析关系：

td＝-0.5*th

根据拉索受拉的原则即可指定单元力密度。

确定单元力密度后，在matlab平台中，对单元力密度赋予数值，将符号矩阵转换为数值矩阵求零空间，其中的零空间是对称坐标系下节点的z方向坐标的零空间是对称坐标系下节点的x方向坐标的零空间是对称坐标系下节点的y方向坐标通过坐标转换矩阵v^(μ)i求出笛卡尔坐标系下结构某尺寸为单位长度的节点坐标，并且通过已知的约束节点坐标求出比例放大系数，最后通过平移求出真实结构在笛卡尔坐标系中的节点坐标。

因此，一种对称索杆结构找形方法从有支座边界约束索杆结构的所有节点和自由节点两个角度出发，求出与结构x/y/z方向刚体平动相关联的不可约表示的对称子空间，根据这些不可约表示对应的力密度分块子矩阵的秩亏条件及零空间，对维度较小的力密度分块子矩阵进行分析，直接求得结构的可行力密度及节点坐标。

上述实施例仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和等同替换，这些对本发明权利要求进行改进和等同替换后的技术方案，均落入本发明的保护范围。