CN107958124A - 一种基于应力分布模式优化的混合尺度耦合方法 - Google Patents
一种基于应力分布模式优化的混合尺度耦合方法 Download PDFInfo
- Publication number
- CN107958124A CN107958124A CN201711281595.2A CN201711281595A CN107958124A CN 107958124 A CN107958124 A CN 107958124A CN 201711281595 A CN201711281595 A CN 201711281595A CN 107958124 A CN107958124 A CN 107958124A
- Authority
- CN
- China
- Prior art keywords
- stress distribution
- stress
- model
- node
- distribution mode
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000010168 coupling process Methods 0.000 title claims abstract description 40
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 70
- 230000008878 coupling Effects 0.000 claims abstract description 16
- 238000005859 coupling reaction Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000004044 response Effects 0.000 claims description 16
- 230000009471 action Effects 0.000 claims description 10
- 238000004134 energy conservation Methods 0.000 claims description 10
- 238000005452 bending Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 8
- 230000003938 response to stress Effects 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000007670 refining Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 8
- 230000036541 health Effects 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 abstract description 2
- 239000013598 vector Substances 0.000 description 13
- 239000007787 solid Substances 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及基于应力分布模式优化的混合尺度耦合方法，包括S1应力分布计算模型的建模；S2单位荷载作用下的应力分布计算模型中部截面各节点应力和坐标值提取；S3初步确定连接界面处的应力分布模式；S4基于应力分布计算模型进行应力分布模式优化的迭代计算，并以可决系数R2作为评判指标来确定连接界面处的最优应力分布模式；S5多点约束方程系数矩阵所需单元形函数的选取；S6建立多点约束方程，以实现多尺度模型中不同尺度单元在连接界面处的耦合。采用本发明提高了界面应力分布模式的准确性，解决了界面剪应力分布粗糙假定导致的多尺度模型在连接界面处应力集中，模型可靠性降低等问题。从而建立起高精度的多尺度模型，以辅助结构健康监测系统。
Description
技术领域
本发明涉及土木工程结构健康监测领域，具体涉及一种基于应力分布模 式优化的混合尺度耦合方法。
背景技术
建筑结构朝着新颖、复杂的方向发展，尤其是新型建筑结构体系的研究 与运用，具有超高层、超跨度、大悬挑等超出规范限制要求的特点，这对结 构设计参数选取，结构分析模拟计算方法、结构施工方法和结构服役行为估 计等方面带来了挑战。针对此挑战，结构健康监测通过在结构上布置传感器 实时获得结构响应，运用结构响应监测数据体现结构受力状态和工作性能。 但由于布置在结构上的传感器数量有限，结构受力复杂部位小空间的传感器 安装位置受限，导致结构响应监测数据不完备。而相较于监测数据对结构响应获取不完备，多尺度有限元模型则通过对关键结构构件和局部复杂部位建 立细观尺度模型实现在细观尺度方面的结构响应获取，同时通过细观尺度和 宏观尺度模型的综合运用，有效兼顾计算精度和效率。
多尺度有限元模型根据具体建模方法的不同分为信息传递多尺度模型和 一致多尺度模型。其中，信息传递多尺度模型建模时需进行尺度分离，这使 得模型的信息仅由低层尺度模型向高层尺度模型单向传递。显然，信息传递 多尺度模型的计算精度主要取决于模型边界条件的选取是否准确，然而大型 土木工程结构受力复杂，准确选取施加在局部模型上的边界条件较为困难， 而不精确的边界条件将直接导致在此基础上进行的局部模型分析结果的不准 确，甚至导致分析结果不可用。一致多尺度模型把不同尺度下的模型通过恰 当的连接方式耦合起来，使得不同尺度模型协同计算，且模拟结果内涵丰富。 因此，一致多尺度模型在土木工程领域的应用受到了越来越多的关注。
相关技术中，混合尺度耦合对于结构一致多尺度模型的准确性能分析具 有重要意义，而应力分布模式的准确性又是影响混合尺度耦合方法性能的关 键之所在。现有的基于能量守恒的混合尺度耦合方法对截面切应力常进行简 化假设，如假定切应力沿截面壁厚为定值，但实际上切应力沿壁厚方向并非 恒为定值；对于材料力学并未给出的如箱型截面等较复杂截面类型的切应力 公式，采用材料力学已有公式来近似替代原截面上的切应力分布也是较粗糙 的。这种不准确的应力分布模式的引入将直接导致一致多尺度模型在连接界 面处局部应力集中，结构局部细节性能分析受到影响，进而降低了结构多尺 度模型的可靠性。因此，有必要对现有的基于能量守恒的混合尺度耦合方法 进行改进，提高假定的界面应力分布模式的精度，从而建立起更为真实，可 信的结构一致多尺度模型。
发明内容
为至少在一定程度上克服相关技术中存在的问题，本发明提供一种能够 有效提高假定的界面应力分布模式的精度，并建立起高精度的结构一致多尺 度模型的基于应力分布模式优化的混合尺度耦合方法，以解决现有基于能量 守恒的多点约束方程法对界面剪应力分布的假定较粗糙，进而导致结构一致 多尺度模型计算精度降低的问题。
根据本申请的实施例提供的技术方案，为实现以上目的，本发明采用如 下技术方案：一种基于应力分布模式优化的混合尺度耦合方法，其特征在于： 该方法包括以下步骤：
S1、选取模型参数，完成对应力分布计算模型的建模；
S2、单位荷载作用下的应力分布计算模型中部截面各节点应力和坐标值 提取；
S3、基于提取的节点响应信息，对应力沿横截面分布特征进行判断，初 步确定连接界面处的应力分布模式；
S4、基于应力分布计算模型进行应力分布模式优化的迭代计算，并以可 决系数R2作为评判指标来确定连接界面处的最优应力分布模式；
S5、多点约束方程系数矩阵所需单元形函数的选取；
S6、基于能量守恒原理、最优应力分布模式和所选单元形函数建立多点 约束方程，以实现多尺度模型中不同尺度单元在连接界面处的耦合。
进一步的，所述步骤S1中应力分布计算模型包括基于高维单元类型的主 体部分、位于主体部分两端截面形心的约束点和加载点，以及约束点和加载 点与主体部分间的刚性连接。
进一步的，所述应力分布计算模型的主体部分相当于多尺度模型中细观 尺度部分的缩影，其所用模型参数应与多尺度模型中的细观尺度部分所用相 一致，具体的模型参数为单元类型、形函数类型和横截面尺寸；应力分布计 算模型主体部分的长度根据圣维南原理确定；所述位于主体部分两端截面形 心的约束点和加载点由通用有限元软件中的质量单元模拟；所述约束点和加 载点与主体部分间的刚性连接由通用有限元软件中的建立刚性面法来实现。
进一步的，所述步骤S2中的单位荷载包括轴力F1、两个方向的剪力F2和F3、扭矩M1和两个方向的弯矩M2和M3；相应地，在不同荷载工况下需提 取的应力分布计算模型中部截面各节点应力响应分别为由轴力引起的正应力 由两个方向剪力引起的切应力和由扭矩引起的两个方向切应力 和以及由两个方向弯矩引起的正应力和提取的节点坐标值 为应力分布计算模型主体部分横截面所在坐标平面内的具体坐标值。
进一步的，所述步骤S3中，利用提取的节点应力值和坐标值绘制出不同 荷载工况下应力沿横截面的分布图，从而初步掌握各类型应力的分布特征； 并基于此，对连接界面处备选应力分布模式的变量阶次进行限定，同时明确 应力分布在横截面纵横两个方向上的内在关系。
进一步的，在各单位荷载工况下，基于限定的变量阶次范围和明确的应 力在横截面纵横向上的分布特点，对需提取应力响应的分布模式的具体函数 形式进行初选，并将应力分布模式具体为若干种函数形式，而该选定的函数 形式即构成了当前单位荷载工况下的一系列备选应力分布模式。
进一步的，所述步骤S4中，可决系数R2作为评判各备选应力分布模式 优劣的指标，从中初选出较优的应力分布模式；通过不断细化应力分布计算 模型的网格尺寸来更新各节点应力和坐标值，进而通过迭代计算来不断优化 初选应力分布模式；最终，以可决系数R2>0.9作为终止条件来确定连接界面 处的最优应力分布模式；
所述可决系数R2的计算公式为：
进一步的，应力分布模式由各单位荷载工况下提取的相应节点应力和坐 标数据，利用基于非线性最小二乘法的曲面拟合技术拟合得到；其中，节点 应力值为应力分布模式的因变量，即输出项；节点坐标值为应力分布模式的 自变量，即输入项。
进一步的，所述步骤S5中，计算多点约束方程系数矩阵所需单元形函数 由连接界面维度和多尺度模型中细观尺度部分所用单元类型来确定；所述连 接界面维度包括两维多尺度模型中的一维连接界线和三维多尺度模型中的两 维连接界面；所述细观尺度部分所用单元类型，具体考察的是单元是否含有 中间节点以及单元形状。
进一步的，所述步骤S6中，对单位荷载作用下，优化得到的应力分布模 式和选取的单元形函数，利用多尺度模型在连接界面处做功相等来求解约束 方程系数矩阵和建立多点约束方程，以实现多尺度模型中不同尺度单元在连 接界面处的耦合；
进一步的，多点约束方程以多尺度模型连接界面一侧节点为主节点，另 一侧节点为从节点，通过指定从节点某一方向自由度与主节点相应自由度间 的内在关系，来实现不同尺度单元在连接界面处的自由度耦合，其表达式为
式(2)中，us_i为从节点i的某一方向自由度；um_j为主节点j的相应自 由度，且j≠i；n为主节点自由度数；Cj为系数项；C0为常数项。
本申请的实施例提供的技术方案可以包括以下有益效果：考虑混合尺度 耦合对结构一致多尺度模型准确性能分析的重要性，提出基于应力分布模式 优化的混合尺度耦合方法，以提高界面应力分布模式的准确性，解决界面剪 应力分布粗糙假定导致的一致多尺度模型在连接界面处局部应力集中，结构 局部细节性能分析受到影响，结构多尺度模型可靠性降低等问题。从而建立 起可模拟结构真实行为的高精度一致多尺度模型，以辅助结构健康监测系统。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性 的，并不能限制本申请。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申 请的实施例，并与说明书一起用于解释本申请的原理。
图1是本发明基于应力分布模式优化的混合尺度耦合方法流程图；
图2是本发明多尺度模型在连接界面处梁单元-实体单元力传递示意图；
图3是本发明应力分布计算模型示意图之侧视图；
图4是本发明应力分布计算模型示意图之横截面示意图；
图5是本发明两维多尺度模型建模时单元形函数选取示意图；
图6是本发明三维多尺度模型建模时单元形函数选取示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的 描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的 要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所 有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一 些方面相一致的方法的例子。
如图1所示，本发明提供一种基于应力分布模式优化的混合尺度耦合方 法，该方法包括以下步骤：
S1、合理选取模型参数，完成对应力分布计算模型的建模；
S2、单位荷载作用下的应力分布计算模型中部截面各节点应力和坐标值 提取；
S3、基于提取的节点响应信息，对应力沿横截面分布特征进行判断，初 步确定连接界面处的应力分布模式；
S4、基于应力分布计算模型进行应力分布模式优化的迭代计算，并以可 决系数R2作为评判指标来确定连接界面处的最优应力分布模式；
S5、多点约束方程系数矩阵所需单元形函数的选取；
S6、基于能量守恒原理、最优应力分布模式和所选单元形函数建立多点 约束方程，以实现多尺度模型中不同尺度单元在连接界面处的耦合。
由于本实施例中所提出的基于应力分布模式优化的混合尺度耦合方法最 终落脚于多尺度模型在连接界面处做功相等，来实现不同尺度单元间多点约 束方程的建立，故首先以基于梁单元和实体单元的一致多尺度模型为例，说 明多尺度模型在连接界面处的能量守恒(即做功相等)，具体在连接界面处 的力传递示意图如图2所示。其中，梁单元节点有6自由度{u，v，w，θ1，θ2， θ3}，实体单元节点有3自由度{u'，v'，w'}；F1为梁端轴力，F2、F3为梁端沿 2、3轴方向的剪力，M1为梁端扭矩，M2、M3为梁端沿2、3轴方向的弯矩； σ1为连接界面处实体单元正应力，τ2、τ3为连接界面处实体单元沿2、3轴方 向的切应力；A为实体单元在耦合界面处的面积。J为第J个实体单元。
根据能量守恒原理，图2中不同尺度单元在连接界面处应做功相等，即 等式(1)成立
式(1)左侧为梁单元在连接界面处做功，右侧为实体单元在连接界面处 做功。其中，u′J、v′J、w′J为实体单元J沿1、2和3轴的节点位移向量；K为 连接界面处实体单元总数；AJ为第J个实体单元的面积；NJ为第J个实体单 元的形函数向量。对式(1)进行简化、合并得如下的矩阵形式
式(2)中，FB为连接界面处梁端节点力向量；UB为连接界面处梁端节 点位移向量；US为连接界面处实体单元节点位移向量；为多点约束方程系 数矩阵。
结构一致多尺度模型同时包含宏观尺度单元和细观尺度单元，但不同尺 度单元节点自由度数往往不一致，故需对不同尺度单元在连接界面处的节点 建立自由度耦合关系，以确保力和位移在不同尺度单元间正确传递。
多点约束方程以连接界面一侧节点为主节点，另一侧节点为从节点，通 过指定从节点某一方向自由度与主节点相应自由度间的内在关系，来实现不 同尺度单元在连接界面处的自由度耦合，其表达式为
式(3)中，us_i为从节点i的某一方向自由度；um_j为主节点j的相应自 由度(j≠i)；n为主节点自由度数；Cj为系数项。
式(3)的矩阵简化形式为
us+Cum＝0(4)
式中，us为从节点自由度向量；um为主节点自由度向量；C为多点约束 方程系数矩阵。
对比式(2)与式(4)可知，若取式(2)中梁端节点力向量为单位向量， 并将梁端节点作为从节点，实体模型在连接界面处各节点作为主节点来考虑， 则式(2)与式(4)等价，此时式(2)中的矩阵即为建立多点约束方程所 需的系数矩阵。故可考虑梁端分别作用有单位集中力或力矩的情况，通过计 算各单位荷载作用下的系数向量，即可对梁端节点的6个自由度分别建立起 各自的多点约束方程。如当梁端仅作用有单位轴力时，梁端节点力向量 FB＝[1,0,0,0,0,0]T，则式(2)可简化为
式(5)即为梁端节点与连接界面处实体单元节点关于1轴方向平动自由 度间的多点约束方程。其中，即为建立多点约束方程所需 的系数向量。同理，当梁端分别仅作用有单位剪力、扭矩和弯矩时，关于梁 端节点其余5个自由度方向的多点约束方程也可通过能量守恒等式(2)获得
整合梁单元六个自由度方向的多点约束方程系数向量得完整的系数矩阵 为
从式(11)可知，多点约束方程系数矩阵与梁端单位荷载在实体单元上 产生的应力分布模式和单元形函数相关，故求取梁端节点在单位荷载作用下 传递至连接界面处实体单元上的真实应力分布模式，以及恰当选取单元形函 数是建立合理、有效的多点约束方程的关键。
下面对步骤S1的实施做进一步的详述：
为获取梁端节点在单位荷载作用下传递至实体单元的真实应力分布模 式，建立如图3所示的应力分布计算模型，其横截面示意图如图4所示。该 应力分布计算模型由基于高维单元类型的主体部分、位于主体部分两端截面 形心的约束点和加载点，以及约束点和加载点与主体部分间的刚性连接三部 分组成。所述应力分布计算模型的主体部分其相当于多尺度模型中细观尺度 部分的缩影，故该应力分布计算模型的单元类型、形函数类型和横截面尺寸 应与多尺度模型中的细观尺度部分所用相一致。
此外，应力分布计算模型主体部分的长度根据圣维南原理确定；所述位 于主体部分两端截面形心的约束点和加载点由通用有限元软件中的质量单元 模拟；所述约束点和加载点与主体部分间的刚性连接由通用有限元软件中的 建立刚性面法来实现。
此外，在建立应力分布计算模型时还需考虑以下两点内容：
(1)由圣维南原理可知，在距离边界一定长度范围内，结构内部的应力 分布会受到边界效应的影响，进而导致该区域内的结构应力分布复杂。因此， 在建立应力分布计算模型时应充分考虑这一点，选取的计算模型长度不宜过 小，以确保模型中部一定长度范围内的应力分布不受边界效应影响，同时， 选取的计算模型长度也不宜过长，以免造成较大的计算量。为此，可采用含 连接界面构件长度的1/3作为应力分布计算模型长度。
(2)考虑到细观尺度模型其横截面上的应力分布在一定程度上受划分网 格尺寸的影响，为确保运用该计算模型所得的应力分布模式真实有效，应对 其采用比多尺度模型中细观尺度部分所用单元尺寸更细化的尺寸进行网格划 分，如对应力分布计算模型网格尺寸采用多尺度模型中细观尺度部分网格尺 寸的1/2。
下面对步骤S2和S3的实施做进一步的详述：
基于图3所示应力分布计算模型，在其左侧约束点处施加固定端约束， 右侧加载点处分别施加单位轴力F1、剪力F2、剪力F3、扭矩M1、弯矩M2以 及弯矩M3，并进行求解计算以获取应力分布计算模型的全部节点响应值，如 节点位移、应力和应变等。具体地：首先提取图3虚线框所示应力分布计算 模型中间截面的各节点编号；然后通过节点编号提取这些节点的响应值用于 应力分布模式优化，需注意的是在不同的荷载工况下所需的节点响应会各有 不同，如在梁端单位轴力F1作用下需提取响应为节点正应力在梁端单位 剪力F2和F3作用下需提取响应为节点切应力和在梁端单位扭矩M1作用下需提取响应为节点切应力和在梁端单位弯矩M2和M3作用下 需提取响应为节点正应力和
如图3-图4所示，本实施例中所提方法采用基于非线性最小二乘法的曲 面拟合技术进行连接界面处高维单元应力分布模式优化，为此需对应力分布 计算模型中间截面各节点坐标值进行提取，同时提取的节点坐标应属于计算 模型横截面所处坐标平面。考虑到本发明所提方法建立的应力分布计算模型， 其横截面始终处于y-z坐标平面内(如图4所示)，故针对本发明所提方法需 提取的节点坐标为y轴和z轴坐标值，并记作Node(y,z)。
作为一种优选的实施方式，本实施例所述步骤S3中，
在进行曲面拟合之前，首先利用提取的节点应力值和坐标值绘制出不同 荷载工况下应力沿横截面的分布图，然后通过观察各应力分布图的趋势，对 应力沿横截面纵横向的分布特点有初步地了解，以便在拟合应力分布的回归 曲面模型时对可选的变量阶次加以限制，同时明确y、z两个方向对应力分布 的相互影响关系。
需要补充说明的是，在各单位荷载工况下，基于限定的变量阶次范围和 明确的应力在横截面纵横向上的分布特点，对需提取应力响应的分布模式的 具体函数形式进行初选，并将应力分布模式具体为若干种函数形式，而该选 定的函数形式即构成了当前单位荷载工况下的一系列备选应力分布模式。
下面对步骤S4的实施做进一步的详述：
根据各类节点应力值和坐标值，通过非线性最小二乘法来拟合得到各备 选函数的具体参数值，从而得到了多组具体的应力分布回归曲面函数，并通 过定义评价拟合质量的拟合优度指标来选定最优的应力分布模式。在此，选 取可决系数R2作为评价指标，其具体的数学表达式为
式(12)中，SSR为回归平方和；SST为总平方和。为节点i的应力拟 合值；fi为节点i的应力真值；为应力分布计算模型中间截面各节点应力平 均值。
需要进一步说明的是，可决系数R2的值在0～1之间变动，其值越接近于 1则表明拟合函数对应力真值的描述越准确，获得的应力分布模式质量越高。 确定最优应力分布模式的步骤具体为：
(1)对于每一种应力响应，计算其所有备选函数的可决系数R2，并按其 数值大小进行排序。然后选取可决系数R2最大的函数作为初始应力分布模式；
(2)给定可决系数的阈值为0.9，则当初始应力分布模式的可决系数大 于该阈值时，初始应力分布模式即为最优应力分布模式。若初始应力分布模 式的可决系数小于阈值0.9，则需进行以下各步骤；
(3)对应力分布计算模型的网格尺寸进行细化，再次将单位集中力和力 矩分别施加于计算模型的右侧加载点上，并通过求解计算，提取模型中间截 面新的节点应力和坐标值；
(4)基于步骤(3)提取的中间截面节点应力和坐标值，对初始应力分 布模式的函数形式运用非线性最小二乘法来确定其新的参数值，从而实现对 回归曲面模型的更新。
(5)计算更新的回归曲面模型的可决系数R2，并与阈值0.9进行比较。 若更新模型的可决系数R2值大于阈值，则该更新模型即为最优应力分布模式。 若更新模型的可决系数R2小于阈值，则重复步骤(3)-(5)直至更新模型的 可决系数值大于给定的阈值0.9。
将通过上述步骤获取的最优应力分布模式分别记作 和需要说明的是，本发明所提方法中最优应力分 布模式的具体数学表达式是基于应力分布计算模型的坐标系统推导得到。如 前所述，本发明所提方法中建立的应力分布计算模型的坐标系统始终保持x 轴平行于横截面形心线，y-z坐标平面为计算模型横截面所处平面。但对于大 型土木工程结构的多尺度模型而言，其结构构件具有各自的空间位形，即本 发明所提方法中定义的应力分布计算模型的坐标系统并不适用于所有需进行 界面耦合的构件。此时，需对不满足要求的结构构件根据应力分布计算模型 的坐标系统，重新指定构件的局部坐标系统，并将连接界面处高维单元的各 节点坐标系旋转至相同的局部坐标系下。
下面对步骤S5的实施做进一步的详述，
所述步骤S5中，计算多点约束方程系数矩阵所需单元形函数由连接界面 维度和多尺度模型中细观尺度部分所用单元类型来确定；所述连接界面维度 包括两维多尺度模型中的一维连接界线和三维多尺度模型中的两维连接界 面；所述细观尺度部分所用单元类型，具体考察的是单元是否含有中间节点 以及单元形状。
计算多点约束方程系数矩阵所需的单元形函数与连接界面维度和多尺度 模型中细观尺度部分所用单元类型相关，具体地：
如图5所示，具体参阅如图5a)和图5b)中的梁-板单元耦合情况，对该 两维多尺度模型，其连接界面为一维的线，故用于计算多点约束方程系数矩 阵的单元形函数应与一维单元相对应。而具体选用的形函数形式则由多尺度 模型中所用高维单元是否含有中间节点来最终确定。部分形函数选取情况如 图5所示，其中i、j、k分别表示一维单元的起始节点、中间节点和终止节点；
如图6所示，具体参阅如图6a)和图6b)中的梁-实体单元耦合情况，对 该三维多尺度模型，其连接界面为两维的面，故用于计算多点约束方程系数 矩阵的单元形函数应与两维平面单元相对应。而具体选用的形函数形式则由 多尺度模型中所用高维单元的单元形状(如六面体、五面体和四面体单元) 以及单元是否含有中间节点来最终确定。如选取四节点四边形单元的形函数 (具体如图6a)所示)或八节点四边形单元的形函数(具体如图6b)所示)。
对步骤S6的实施补充说明如下：
由式(11)可知，将拟合得到的高维单元的最优应力分布模式与相应的 单元形函数结合，通过积分运算即可求得多点约束方程的系数矩阵其中， 高维单元的最优应力分布模式通过步骤S1～S4获得；相应的单元形函数通过 步骤S5获得。当该系数矩阵已知时，直接运用通用有限元软件的相关命令即 可实现不同尺度单元在连接界面处的节点自由度耦合，如中的CE命 令和中的EQUATION命令。
本发明从辅助结构健康监测角度出发，针对现有基于能量守恒的多点约 束方程法对界面剪应力分布的假定较粗糙，导致结构一致多尺度模型计算精 度降低这一问题，提出了基于应力分布模式优化的混合尺度耦合方法，以提 高假定的界面应力分布模式的精度，并建立起高精度的结构一致多尺度模型。
可以理解的是，上述各实施例中相同或相似部分可以相互参考，在一些实 施例中未详细说明的内容可以参见其他实施例中相同或相似的内容。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表 示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码 的模块、片段或部分，并且本申请的优选实施方式的范围包括另外的实现，其 中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或 按相反的顺序，来执行功能，这应被本申请的实施例所属技术领域的技术人员 所理解。
尽管上面已经示出和描述了本申请的实施例，可以理解的是，上述实施例 是示例性的，不能理解为对本申请的限制，本领域的普通技术人员在本申请的 范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种基于应力分布模式优化的混合尺度耦合方法，其特征在于：该方法包括以下步骤：
S1、选取模型参数，完成对应力分布计算模型的建模；
S2、单位荷载作用下的应力分布计算模型中部截面各节点应力和坐标值提取；
S3、基于提取的节点响应信息，对应力沿横截面分布特征进行判断，初步确定连接界面处的应力分布模式；
S4、基于应力分布计算模型进行应力分布模式优化的迭代计算，并以可决系数R2作为评判指标来确定连接界面处的最优应力分布模式；
S5、多点约束方程系数矩阵所需单元形函数的选取；
S6、基于能量守恒原理、最优应力分布模式和所选单元形函数建立多点约束方程，以实现多尺度模型中不同尺度单元在连接界面处的耦合。
2.根据权利要求1所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述步骤S1中应力分布计算模型包括基于高维单元类型的主体部分、位于主体部分两端截面形心的约束点和加载点，以及约束点和加载点与主体部分间的刚性连接。
3.根据权利要求2所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述应力分布计算模型的主体部分相当于多尺度模型中细观尺度部分的缩影，其所用模型参数应与多尺度模型中的细观尺度部分所用相一致，具体的模型参数为单元类型、形函数类型和横截面尺寸；
应力分布计算模型主体部分的长度根据圣维南原理确定；所述位于主体部分两端截面形心的约束点和加载点由通用有限元软件中的质量单元模拟；所述约束点和加载点与主体部分间的刚性连接由通用有限元软件中的建立刚性面法来实现。
4.根据权利要求1所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述步骤S2中的单位荷载包括轴力F1、两个方向的剪力F2和F3、扭矩M1和两个方向的弯矩M2和M3；相应地，在不同荷载工况下需提取的应力分布计算模型中部截面各节点应力响应分别为由轴力引起的正应力由两个方向剪力引起的切应力和由扭矩引起的两个方向切应力和以及由两个方向弯矩引起的正应力和提取的节点坐标值为应力分布计算模型主体部分横截面所在坐标平面内的具体坐标值。
5.根据权利要求1所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述步骤S3中，利用提取的节点应力值和坐标值绘制出不同荷载工况下应力沿横截面的分布图，从而初步掌握各类型应力的分布特征；并基于此，对连接界面处备选应力分布模式的变量阶次进行限定，同时明确应力分布在横截面纵横两个方向上的内在关系。
6.根据权利要求5所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：在各单位荷载工况下，基于限定的变量阶次范围和明确的应力在横截面纵横向上的分布特点，对需提取应力响应的分布模式的具体函数形式进行初选，并将应力分布模式具体为若干种函数形式，而该选定的函数形式即构成了当前单位荷载工况下的一系列备选应力分布模式。
7.根据权利要求1至6任一项所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述步骤S4中，可决系数R2作为评判各备选应力分布模式优劣的指标，从中初选出较优的应力分布模式；通过不断细化应力分布计算模型的网格尺寸来更新各节点应力和坐标值，进而通过迭代计算来不断优化初选应力分布模式；最终，以可决系数R2>0.9作为终止条件来确定连接界面处的最优应力分布模式；
所述可决系数R2的计算公式为：
<mrow>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
应力分布模式由各单位荷载工况下提取的相应节点应力和坐标数据，利用基于非线性最小二乘法的曲面拟合技术拟合得到；其中，节点应力值为应力分布模式的因变量，即输出项；节点坐标值为应力分布模式的自变量，即输入项。
8.根据权利要求7所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：可决系数R2的值在0～1之间变动，其值越接近于1则表明拟合函数对应力真值的描述越准确，获得的应力分布模式质量越高；确定最优应力分布模式的步骤具体为：
(1)对于每一种应力响应，计算其所有备选函数即选应力分布模式的可决系数R2，并按其数值大小进行排序；然后选取可决系数R2最大的函数作为初始应力分布模式；
(2)给定可决系数的阈值为0.9，则当初始应力分布模式的可决系数大于该阈值时，初始应力分布模式即为最优应力分布模式；若初始应力分布模式的可决系数小于阈值0.9，则需进行以下各步骤；
(3)对应力分布计算模型的网格尺寸进行细化，再次将单位集中力和力矩分别施加于计算模型的右侧加载点上，并通过求解计算，提取模型中间截面新的节点应力和坐标值；
(4)基于步骤(3)提取的中间截面节点应力和坐标值，对初始应力分布模式的函数形式运用非线性最小二乘法来确定其新的参数值，从而实现对回归曲面模型的更新；
(5)计算更新的回归曲面模型的可决系数R2，并与阈值0.9进行比较。若更新模型的可决系数R2值大于阈值，则该更新模型即为最优应力分布模式。若更新模型的可决系数R2小于阈值，则重复步骤(3)-(5)直至更新模型的可决系数值大于给定的阈值0.9。
9.根据权利要求1所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述步骤S5中，计算多点约束方程系数矩阵所需单元形函数由连接界面维度和多尺度模型中细观尺度部分所用单元类型来确定；所述连接界面维度包括两维多尺度模型中的一维连接界线和三维多尺度模型中的两维连接界面；所述细观尺度部分所用单元类型，具体考察的是单元是否含有中间节点以及单元形状。
10.根据权利要求1所述的基于应力分布模式优化的混合尺度耦合方法，其特征在于：所述步骤S6中，对单位荷载作用下，优化得到的应力分布模式和选取的单元形函数，利用多尺度模型在连接界面处做功相等来求解约束方程系数矩阵和建立多点约束方程，以实现多尺度模型中不同尺度单元在连接界面处的耦合；
多点约束方程以多尺度模型连接界面一侧节点为主节点，另一侧节点为从节点，通过指定从节点某一方向自由度与主节点相应自由度间的内在关系，来实现不同尺度单元在连接界面处的自由度耦合，其表达式为：
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>C</mi>
<mi>j</mi>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>m</mi>
<mo>_</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，us_i为从节点i的某一方向自由度；um_j为主节点j的相应自由度，且j≠i；n为主节点自由度数；Cj为系数项；C0为常数项。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711281595.2A CN107958124B (zh) | 2017-12-06 | 2017-12-06 | 一种基于应力分布模式优化的混合尺度耦合方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711281595.2A CN107958124B (zh) | 2017-12-06 | 2017-12-06 | 一种基于应力分布模式优化的混合尺度耦合方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107958124A true CN107958124A (zh) | 2018-04-24 |
| CN107958124B CN107958124B (zh) | 2021-04-02 |
Family
ID=61958130
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711281595.2A Active CN107958124B (zh) | 2017-12-06 | 2017-12-06 | 一种基于应力分布模式优化的混合尺度耦合方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107958124B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110090887A (zh) * | 2019-04-09 | 2019-08-06 | 哈尔滨理工大学 | 一种基于数值模拟的充液拉深成形周向加压方法 |
| CN110209057A (zh) * | 2019-06-26 | 2019-09-06 | 哈尔滨工业大学（深圳） | 一种最优状态反馈保性能的amd控制器的设计方法 |
| CN111090906A (zh) * | 2019-10-10 | 2020-05-01 | 惠州市德赛西威汽车电子股份有限公司 | 一种以非线性弹簧代替胶水实体建模的方法 |
| CN111861176A (zh) * | 2020-07-10 | 2020-10-30 | 哈尔滨工业大学（深圳） | 空间杆系结构重要监测杆件确定方法、装置及设备 |
| CN111859745A (zh) * | 2020-07-10 | 2020-10-30 | 哈尔滨工业大学（深圳） | 型钢混凝土结构响应分布的获取方法、装置和设备 |
| CN113255196A (zh) * | 2021-07-05 | 2021-08-13 | 广州中望龙腾软件股份有限公司 | 一种网格优化方法、网格生成器、存储介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102819632A (zh) * | 2012-07-26 | 2012-12-12 | 苏州工业园区设计研究院股份有限公司 | 一种钢框架结构节点多尺度有限元模型建模方法 |
| CN104730596A (zh) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
| US20150205006A1 (en) * | 2010-03-25 | 2015-07-23 | Schlumberger Technology Corporation | Downhole modeling using inverted pressure and regional stress |
| EP2975437A1 (en) * | 2014-07-15 | 2016-01-20 | Services Petroliers Schlumberger | Method and system to invert for fault activity and tectonic stress |
| CN106934185A (zh) * | 2017-04-27 | 2017-07-07 | 中国石油大学(华东) | 一种弹性介质的流固耦合多尺度流动模拟方法 |
| CN107368642A (zh) * | 2017-07-13 | 2017-11-21 | 武汉大学 | 金属增材制造多尺度多物理场耦合仿真方法 |
-
2017
- 2017-12-06 CN CN201711281595.2A patent/CN107958124B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150205006A1 (en) * | 2010-03-25 | 2015-07-23 | Schlumberger Technology Corporation | Downhole modeling using inverted pressure and regional stress |
| CN102819632A (zh) * | 2012-07-26 | 2012-12-12 | 苏州工业园区设计研究院股份有限公司 | 一种钢框架结构节点多尺度有限元模型建模方法 |
| EP2975437A1 (en) * | 2014-07-15 | 2016-01-20 | Services Petroliers Schlumberger | Method and system to invert for fault activity and tectonic stress |
| CN104730596A (zh) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
| CN106934185A (zh) * | 2017-04-27 | 2017-07-07 | 中国石油大学(华东) | 一种弹性介质的流固耦合多尺度流动模拟方法 |
| CN107368642A (zh) * | 2017-07-13 | 2017-11-21 | 武汉大学 | 金属增材制造多尺度多物理场耦合仿真方法 |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110090887A (zh) * | 2019-04-09 | 2019-08-06 | 哈尔滨理工大学 | 一种基于数值模拟的充液拉深成形周向加压方法 |
| CN110209057A (zh) * | 2019-06-26 | 2019-09-06 | 哈尔滨工业大学（深圳） | 一种最优状态反馈保性能的amd控制器的设计方法 |
| CN111090906A (zh) * | 2019-10-10 | 2020-05-01 | 惠州市德赛西威汽车电子股份有限公司 | 一种以非线性弹簧代替胶水实体建模的方法 |
| CN111090906B (zh) * | 2019-10-10 | 2023-07-07 | 惠州市德赛西威汽车电子股份有限公司 | 一种以非线性弹簧代替胶水实体建模的方法 |
| CN111861176A (zh) * | 2020-07-10 | 2020-10-30 | 哈尔滨工业大学（深圳） | 空间杆系结构重要监测杆件确定方法、装置及设备 |
| CN111859745A (zh) * | 2020-07-10 | 2020-10-30 | 哈尔滨工业大学（深圳） | 型钢混凝土结构响应分布的获取方法、装置和设备 |
| CN111861176B (zh) * | 2020-07-10 | 2024-03-22 | 哈尔滨工业大学（深圳） | 空间杆系结构重要监测杆件确定方法、装置及设备 |
| CN111859745B (zh) * | 2020-07-10 | 2024-09-20 | 哈尔滨工业大学（深圳） | 型钢混凝土结构响应分布的获取方法、装置和设备 |
| CN113255196A (zh) * | 2021-07-05 | 2021-08-13 | 广州中望龙腾软件股份有限公司 | 一种网格优化方法、网格生成器、存储介质 |
| CN113255196B (zh) * | 2021-07-05 | 2021-11-19 | 广州中望龙腾软件股份有限公司 | 一种网格优化方法、网格生成器、存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107958124B (zh) | 2021-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107958124B (zh) | 一种基于应力分布模式优化的混合尺度耦合方法 | |
| Lai et al. | Designing a shape–performance integrated digital twin based on multiple models and dynamic data: a boom crane example | |
| JP6832475B1 (ja) | ランダムアイソジオメトリック解析に基づくブレードの高剛性の設計方法 | |
| CN108959686A (zh) | 一种基于灵敏度分析的有限元模型修正方法 | |
| JP2022521907A (ja) | ハイブリッド繊維複合材料の板巻きシェル構造に対する高速協調最適化方法 | |
| CN112016167B (zh) | 基于仿真和优化耦合的飞行器气动外形设计方法及系统 | |
| CN103902785B (zh) | 一种基于多元不确定性的结构有限元模型修正方法 | |
| JP6784780B2 (ja) | 大規模再生可能エネルギーのデータについて確率モデルを構築する方法 | |
| CN106021644B (zh) | 确定混合维模型界面约束方程系数的方法 | |
| CN105631136B (zh) | 一种考虑封闭框架剪切变形和刚性节点的内力计算方法 | |
| CN113221475B (zh) | 一种用于高精度流场分析的网格自适应方法 | |
| CN103106305A (zh) | 一种基于实测模态的空间网格结构模型分步修正方法 | |
| CN112214843B (zh) | 一种用于风洞试验机翼模型有限元刚度修正方法及装置 | |
| CN103838913A (zh) | 曲线箱梁弯桥的有限单元法 | |
| CN113343427A (zh) | 一种基于卷积神经网络的结构拓扑构型预测方法 | |
| CN112414668A (zh) | 一种风洞试验数据静气弹修正方法、装置、设备及介质 | |
| CN109271655A (zh) | 一种基于非对称有限元算法的材料尺度效应分析方法 | |
| CN117725859A (zh) | 飞机疲劳损伤关键部位的疲劳裂纹扩展预测方法及系统 | |
| CN109977467A (zh) | 一种机翼结构可靠性灵敏度分析方法 | |
| CN109101731A (zh) | 基于B-distance的时变可靠性序列优化方法 | |
| Długosz et al. | The optimal design of UAV wing structure | |
| CN115221622A (zh) | 一种大尺寸复合材料机身壁板装配定位夹持布局优化方法 | |
| Kaminsky et al. | Adaptive sampling techniques for surrogate modeling to create high-dimension aerodynamic loading response surfaces | |
| CN110083956B (zh) | 机床运动机构可靠性灵敏度分析方法 | |
| CN106709167A (zh) | 一种基于相似度理论的装配有效性评估方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |