一种基于强耦合整体技术的索膜结构抗风设计方法技术

本发明专利技术涉及一种基于强耦合整体技术的索膜结构抗风设计方法，建立索膜结构的初始形态；建立流固耦合体系控制方程和流固交界面处的耦合条件，对其进行时空离散，得到流固耦合体系的强耦合整体方程；采用SSTK-ω湍流模型模拟湍流，得到考虑湍流模型的流固耦合体系的强耦合整体方程，并采用Newton方法将其进行线性化；构造线性化后的强耦合整体方程矩阵的预处理器，得到带有预处理器的流固耦合体系的强耦合整体方程，采用Krylov子空间投影法对其进行求解，得到强耦合整体方程流体压力、流体速度、索膜结构位移和线弹性有限元模型的位移；从而进行索膜结构的抗风设计，得到抗风设计参数。

【技术实现步骤摘要】
一种基于强耦合整体技术的索膜结构抗风设计方法
本专利技术属于建筑行业索膜结构抗风分析与设计领域，具体涉及一种基于强耦合整体技术的索膜结构抗风设计方法。
技术介绍
目前国内外索膜结构抗风设计理论尚不成熟。风与索膜结构的流固耦合作用是索膜结构抗风设计中的难点问题，流固耦合作用研究的滞后已成为制约索膜结构风振理论完善和发展的主要因素。沈世钊院士也曾指出，对于索膜结构的抗风设计，重点是解决风与索膜结构的流固耦合问题。因此采用合理方法研究风与索膜结构的流固耦合问题对于完善索膜结构的风振设计理论、指导工程实践具有重要的理论意义和实用价值。目前工程界已经可以通过使用相对成熟的有限元软件对索膜结构在风荷载下的应力、变形效应进行精确计算。同时，也能够通过使用计算流体动力学软件对索膜结构所受的风荷载进行准确估计。然而，对于流固耦合计算的数值方法和应用却少之又少。流固耦合作用分析的数值模拟技术包括弱耦合分区技术、强耦合分区技术和强耦合整体技术。弱耦合分区技术是在每一时间步内先对流体控制方程和结构控制方程分别独立求解，然后将作用在流体域模型上的气动力荷载通过流体和结构的交界面传递给结构域模型，用来预测结构的位移，然后再将结构的位移作为新的荷载传递给流体域，此过程如此反复，直到结果收敛于指定值。强耦合分区技术就是在流体域和结构域各自求解器的基础上再增加一次迭代循环，在每一个时间步求解非线性方程组，计算出全场变量的值。强耦合整体技术是指在每一时间步内对流体控制方程和结构控制方程同时联立求解，即在每一时间步内流体域和结构域互相传递气动力和结构位移的过程是同时进行的；狭义的说，就是对整个流固耦合问题用一个方程组表示，并进行求解。目前对于索膜结构抗风设计中的流固耦合作用的分析主要采用弱耦合分区和强耦合分区技术，原因是它们可以利用现有的流体或结构计算软件，但这两种方法由于流体网格和结构网格的划分不一致，结构的反应要滞后于流体域，从而导致交界面处的计算能量耗散，引起数值计算不稳定，最终影响结果的准确性。并且不论流体和结构单个求解器的准确性有多高，这种方法也只能达到一阶准确性。而且对于索膜结构这类流体和结构的耦合作用很强的情况，弱耦合分区和强耦合分区技术非常容易发散，计算准确性会受到较大影响。尽管采用强耦合整体技术计算流固耦合问题的研究还非常有限，但强耦合整体技术在稳定性和准确性方面却显示出其优势。目前国内外对于索膜结构抗风设计中流固耦合问题研究的主要发展趋势是力争采用更准确高效的计算方法和手段解决复杂的流固耦合计算，为索膜结构的抗风设计提供更可靠的理论或方法，因此强耦合整体技术可以成为研究风与索膜结构流固耦合作用的有利工具。采用强耦合整体技术对索膜结构进行抗风设计，首先需要重点解决的问题之一就是处理空气流体域和结构域交界面处的数据信息传递，即流体域计算域和结构域计算域的计算结果的相互交换，完成下一时间步计算的网格更新过程。需要注意的是，对交界面处信息的处理要适应强耦合整体技术的实施。并在处理流固交界面处的数据信息传递之后，采用强耦合整体技术需要采用一个非线性方程组表示空气流体和索膜结构形成的流固耦合体系，即需要推导出耦合体系的非线性强耦合整体方程，强耦合整体技术所采用的非线性方程组通常需要采用Newton方法线性化后求解。在求解该方程过程中，会耗费在对雅可比矩阵的反复集成和相应线性体系解的Newton修正中，这往往会导致该问题的计算量是不可接受的。因此在使用强耦合整体技术时对方程的高效求解也是采用强耦合整体技术需要解决的重点问题。构造预处理器是求解线性化后的非线性方程的一种有效方法，因此在强耦合整体方程的基础上构造合适的预处理器，使得方程的求解过程简化，并节约大量计算机时也是非常关键的步骤。目前的流固耦合计算的预处理器都是分别单独的针对流体域和结构域的，在使用过程中分别对流体域和结构域使用这些预处理器，这样在很大程度上就会降低计算效率。因此需要针对所要采用的强耦合整体技术的特点，构造相应的预处理器，以实现整体隐式方程的求解的高效性和稳定性。索膜结构属于低矮的大跨度结构，位于大气边界层中湍流强度较高的近地区域，索膜结构周围的流场会使其发生变形的缘故，索膜结构常常会遇到包括分离、再附、漩涡脱落大尺度流瞬时流动过程。因此需要采用湍流模型以便更准确的进行索膜结构的抗风设计。这里选用合适的湍流模型并将其融入提出的整体隐式方程中，也是非常重要的环节，也是索膜结构抗风计算的主要技术难点之一。
技术实现思路
针对现有技术的不足，目前索膜结构抗风设计中的流固耦合分析的分区数值模拟技术存在易发散、准确性差的缺点，为索膜结构抗风设计带来安全隐患，本专利技术提出一种基于强耦合整体技术的索膜结构抗风设计方法。将线弹性有限元模型应用于空气流体与索膜结构形成的耦合系统中，采用线弹性有限元模型处理流体和结构交界面处的数据信息交换传递，并基于线弹性有限元模型方程提出一种适应于强耦合整体技术的网格更新方法；采用流体方程、结构方程和线弹性有限元模型方程的变分形式推导出强耦合整体方程；采用Newton方法线性化强耦合整体方程，构造强耦合整体方程的预处理器，采用两方程的SSTk-ω湍流模型模拟湍流，并将其融入整体隐式方程中，给出了计算膜结构风振中流固耦合计算的实现过程和步骤，为索膜结构抗风设计提供更准确可靠的理论方法和造作简便的计算平台。一种基于强耦合整体技术的索膜结构抗风设计方法，包括以下步骤：步骤1：根据实际膜结构材料的物理性质和实际膜结构的初始预应力，通过找形方法建立索膜结构的初始形态；所述的通过找形方法建立索膜结构有限元模型具体方法为：采用小杨氏模量逐步刚度逼近的方法建立索膜结构的初始形态；步骤2：建立流固耦合体系控制方程和流固交界面处的耦合条件，流固耦合体系控制方程包括：流体域控制方程、结构域控制方程和线弹性有限元模型控制方程；流固交界面处的耦合条件包括流固交界面线弹性有限元模型和结构的位移耦合条件、流固交界面处流体和结构的应力耦合条件、流固交界面处流体和结构的速度耦合条件；所述的建立流固耦合体系控制方程的具体方法为：采用纳维叶-斯托克斯方程建立流体域控制方程，采用拉格朗日方程建立结构域控制方程，采用动力平衡方程的半离散有限元方程的准静态形式建立线弹性有限元模型控制方程；所述的采用线弹性有限元模型处理流体域变形的具体方法为：将流体网格视为线弹性固体，利用流体单元刚化矩阵中的刚化函数得到线弹性有限元模型单元的单元刚化矩阵，得到流固耦合体系的网格更新公式，对流固耦合体系网格进行实时更新，采用控制网格单元的纵横比控制更新后网格质量；步骤3：采用线弹性有限元模型处理流体域变形，利用线弹性有限元模型控制方程对流固耦合体系网格进行实时更新；步骤4：对流固耦合体系的控制方程和流固交界面处的耦合条件进行时空离散，得到流固耦合体系的各个控制方程的变分弱形式和流固交界面处的各耦合条件的变分弱形式，合并各个控制方程和耦合条件的变分弱形式，得到流固耦合体系的强耦合整体方程；所述的对流固耦合体系的控制方程和流固交界面处的耦合条件进行时空离散的具体方法为：采用伽辽金有限元法对流固耦合体系控制方程、流固交界面处线弹性有限元模型和结构的位移耦合条件、流固交界面处流体和结构的应力耦合条件进行空间离散，本文档来自技高网...

【技术保护点】
一种基于强耦合整体技术的索膜结构抗风设计方法，其特征在于，包括以下步骤：步骤1：根据实际膜结构材料的物理性质和实际膜结构的初始预应力，通过找形方法建立索膜结构的初始形态；步骤2：建立流固耦合体系控制方程和流固交界面处的耦合条件，流固耦合体系控制方程包括：流体域控制方程、结构域控制方程和线弹性有限元模型控制方程；流固交界面处的耦合条件包括流固交界面线弹性有限元模型和结构的位移耦合条件、流固交界面处流体和结构的应力耦合条件、流固交界面处流体和结构的速度耦合条件；步骤3：采用线弹性有限元模型处理流体域变形，利用线弹性有限元模型控制方程对流固耦合体系网格进行实时更新；步骤4：对流固耦合体系的控制方程和流固交界面处的耦合条件进行时空离散，得到流固耦合体系的各个控制方程的变分弱形式和流固交界面处的各耦合条件的变分弱形式，合并各个控制方程和耦合条件的变分弱形式，得到流固耦合体系的强耦合整体方程；步骤5：采用SST K‑ω湍流模型模拟湍流，建立SST K‑ω湍流模型控制方程，对SST K‑ω湍流模型控制方程进行空间离散，得到SST K‑ω湍流模型控制方程的变分弱形式，用SST K‑ω湍流模型的变分弱形式替换流固耦合体系的强耦合整体方程中的流体域控制方程的变分弱形式，得到考虑湍流模型的流固耦合体系的强耦合整体方程；步骤6：采用Newton方法将考虑湍流模型的流固耦合体系的强耦合整体方程进行线性化，得到线性化的流固耦合体系的强耦合整体方程；步骤7：构造线性化后的强耦合整体方程矩阵的预处理器，得到带有预处理器的流固耦合体系的强耦合整体方程；步骤8：采用Krylov子空间投影法求解带有预处理器的流固耦合体系的强耦合整体方程，得到强耦合整体方程流体压力、流体速度、索膜结构位移和线弹性有限元模型的位移；步骤9：利用流体压力、流体速度和索膜结构位移进行索膜结构的抗风设计，得到抗风设计参数：索膜结构风压系数、索膜结构风振响应和索膜结构气动特性。...

【技术特征摘要】
1.一种基于强耦合整体技术的索膜结构抗风设计方法，其特征在于，包括以下步骤：步骤1：根据实际膜结构材料的物理性质和实际膜结构的初始预应力，通过找形方法建立索膜结构的初始形态；步骤2：建立流固耦合体系控制方程和流固交界面处的耦合条件，流固耦合体系控制方程包括：流体域控制方程、结构域控制方程和线弹性有限元模型控制方程；流固交界面处的耦合条件包括流固交界面线弹性有限元模型和结构的位移耦合条件、流固交界面处流体和结构的应力耦合条件、流固交界面处流体和结构的速度耦合条件；步骤3：采用线弹性有限元模型处理流体域变形，利用线弹性有限元模型控制方程对流固耦合体系网格进行实时更新；步骤4：对流固耦合体系的控制方程和流固交界面处的耦合条件进行时空离散，得到流固耦合体系的各个控制方程的变分弱形式和流固交界面处的各耦合条件的变分弱形式，合并各个控制方程和耦合条件的变分弱形式，得到流固耦合体系的强耦合整体方程；步骤5：采用SSTK-ω湍流模型模拟湍流，建立SSTK-ω湍流模型控制方程，对SSTK-ω湍流模型控制方程进行空间离散，得到SSTK-ω湍流模型控制方程的变分弱形式，用SSTK-ω湍流模型的变分弱形式替换流固耦合体系的强耦合整体方程中的流体域控制方程的变分弱形式，得到考虑湍流模型的流固耦合体系的强耦合整体方程；所述的对SSTK-ω湍流模型控制方程进行空间离散的具体方法为：采用伽辽金有限元法对SSTK-ω湍流模型的控制方程进行空间离散，即将滑移边界条件作为湍流的边界，将结构壁面的法线方向速度作为垂直于壁面的流体速度，得到离散后的SSTK-ω湍流模型的控制方程；所述SSTK-ω湍流模型模拟湍流的公式如下所示：其中，vf为流体流动速度，ρf为流体密度，σf为流体的全应力张量，k为湍动能，ε为湍动能的耗散率，ω＝ε/k，vfi为流体在水平方向xi上的速度，vfj为流体在竖直方向xj上的速度，xi为水平方向，xj为竖直方向，μp为分子扩散所造成的动力粘性，μt为流体脉动速度对应的动力粘性，φ＝{α，β}为SST模型中的常数，β*＝0.09；φ＝F1φ1+(1-F1)φ2，φ1＝{α1，β1，σk1，σω1}为原k-ω模型的常数，φ2＝{α2，β2，σk2，σω2}为变形后k-ω模型的常数，F1为第一个混合函数；1h是膜结构中心距离最远壁面的距离；采用伽辽金有限元法对SSTK-ω湍流模型的控制方程进行空间离散，将滑移边界条件作为湍流的边界，垂直于壁面的流体速度由固体壁面的法线方向速度给定，结合的壁面摩擦边界条件；得到SSTK-ω湍流模型控制方程的变分弱形式如下所示：其中，λ和γ为检验函数，ftur1和ftur2分别表示SSTK-ω湍流模型控制方程的变分弱形式；用SSTK-ω湍流模型的变分弱形式替换流固耦合体系的强耦合整体方程中的流体域控制方程的变分弱形式，即得到考虑湍流模型的流固耦合体系的强耦合整体方程，如下所示：

【专利技术属性】
技术研发人员：孙芳锦，梁爽，张大明，冯旭，张辉，詹宇楠，王岩露，张尚祥，余磊，
申请(专利权)人：辽宁工程技术大学，
类型：发明
国别省市：辽宁;21