CN103177162A - 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法 - Google Patents

一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法 Download PDF

Info

Publication number
CN103177162A
CN103177162A CN2013101172068A CN201310117206A CN103177162A CN 103177162 A CN103177162 A CN 103177162A CN 2013101172068 A CN2013101172068 A CN 2013101172068A CN 201310117206 A CN201310117206 A CN 201310117206A CN 103177162 A CN103177162 A CN 103177162A
Authority
CN
China
Prior art keywords
thin
wall construction
fluid
node
model
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
Application number
CN2013101172068A
Other languages
English (en)
Other versions
CN103177162B (zh
Inventor
王文全
闫妍
张立翔
曾云
李泽
杨锡阶
杨洋
郭涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201310117206.8A priority Critical patent/CN103177162B/zh
Publication of CN103177162A publication Critical patent/CN103177162A/zh
Application granted granted Critical
Publication of CN103177162B publication Critical patent/CN103177162B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明属于水利水电建筑工程以及航天航空领域薄壁结构动力学特性分析和承载设计领域，特别涉及一种基于交错迭代耦合技术的薄壁结构动力学特性预测方法。根据薄壁结构的物理力学性质、约束条件和第一步计算得到的流体载荷，建立基于有限单元的薄壁结构计算模型；利用设计工况流量（流速）作为流体计算的进口边界条件；利用交错迭代耦合技术进行薄壁结构动力特性分析，从而为结构安全设计提供依据。该计算方法准确性较以往的材料力学方法和单向耦合方法有了显著的提高，而且计算速度快，使流固耦合从理论研究走上实际工程应用。

Description

一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法
技术领域
本发明属于水利水电建筑工程行业薄壁结构动力学特性分析和承载设计领域，而且特别涉及一种基于交错迭代耦合技术的薄壁结构动力学特性预测方法。 
背景技术
流体和薄壁结构的相互作用问题有着广泛的工程背景，普遍存在于人体和自然界中，如血管和血液的相互作用、血液流中的微型泵以及人工心脏、水轮发电机组叶片流激振动、飞机机翼和航空发动机叶片颤振问题以及飞行中的昆虫翅膀和空气等,且随着计算方法、计算手段和日益迫切的工程需要而越发受到学术界和工程界的广泛关注，也取得一定的成果。 
此外，工程界已经可以通过使用相对成熟的有限元软件对结构在流体载荷下的应力、变形等效应进行精确计算。同时，也能够通过使用计算流体动力学软件对结构所受的流体载荷进行准确估计。然而对于进行流固耦合计算的方法和工程应用却较少。其主要原因有：（1）简单的耦合方法由于流体的附加质量会导致求解不稳定（P. Causin, J. F. Gerbeau, F. Nobile, “Added-mass effect in the design of partitioned algorithms for fluid-structure problems”, Comput. Methods Appl. Mech.Engrg, 194 (2005), 4506-4527）。（2）过于复杂的界面插值方法计算量大，缺乏工程可操作性。过于简单界面插值方法导致计算误差大。（3）不同软件间缺乏有效的信息传递与交换接口。 
本方法针对以上问题，提出交错迭代耦合技术，界面信息交换采用移动最小二乘曲面拟合技术，使考虑流固耦合的薄壁结构动力学特性预测方法从理论变成现实。在此基础上，对一薄壁悬臂梁结构的流激振动进行实例分析，验证了上述方法的有效性和准确性。 
发明内容
由于薄壁结构具有较强的抗弯曲变形能力，又能参加总体承力，薄壁结构已取代部分杆系结构，成为水利水电建筑领域、流体机械以及航空结构的主要形式，广泛应用于流体机械的叶片、机翼、机身等。但近年来，薄壁结构的振动和开裂问题日益突出，越发受到学术界和工程界的关注。比如，水轮机导水机构翼型叶片、转轮叶片作为透平机械的核心部件，其水力振动也是典型的流固耦合振动。以往单一的流场计算或单一的结构计算对机组效率、空化预测或对结构设计虽起到了指导性作用，但对研究水力振动，叶片裂纹成因等现象无能为力，因此采用流固耦合计算的方法，来预测这一类复杂湍流与大型薄壁结构之间的耦合振动势在必行，其预测结果更符合真实的物理本质，能更好地指导薄壁结构设计。但目前薄壁结构动力学特性预测无法采用先进的流固耦合技术，就其原因主要有：流体或固体计算模块的稳定性，界面插值技术的准确性，数据交换与通信、计算效率等。因此，本发明提出一种基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其目的就是要通过解决以上问题，为薄壁结构动力学特性预测和后继的承载设计提供一种准确、高效的实用方法。 
为了达到上述目的，本发明提出一种基于交错迭代耦合技术的薄壁结构动力学特性预测方法，包括下列步骤： 
（1）根据薄壁结构的初始位置构型及其周围环境流场，建立适体的基于有限体积的流体计算模型；
（2）根据薄壁结构的物理力学性质、约束条件和步骤（1）计算得到的流体计算模型，建立基于有限单元的薄壁结构计算模型；
（3）将设计工况流量或流速作为流体计算的进口边界条件；
（4）首先用最小二乘曲面拟合技术将流体计算模型与薄壁结构计算模型信息进行界面信息交换（确保计算准确性），然后通过数据通讯管理和实施监控程序在每一个时间步内进行流体计算模型与薄壁结构计算模型的反复交错迭代（解决耦合计算中的准确性），获得薄壁结构的流激振动响应和结构主应力分布，实现薄壁结构动力学特性预测。从而为结构安全设计提供依据。
所述步骤（1）中流体计算模型的建立采用基于混合动态亚格子的大涡模拟模型、有限体积空间离散和全隐式时间离散格式。 
所述步骤（2）中的薄壁结构是指水利水电工程以及航天航空领域广泛采用的由薄型板件和加强构件组成的结构。板件有蒙皮、腹板、隔板等，加强构件有桁条和梁、肋、框的缘条等。不加强的薄壁壳体，如球形、柱形容器等，也属于薄壁结构的范围。 
所述步骤（2）中薄壁结构模型计算采用基于有限单元的薄壁计算模型和直接积分的纽马克法（Newmark）。 
所述步骤（4）中用最小二乘曲面拟合技术将流体计算模型与薄壁结构计算模型信息进行界面信息交换的方法为： 
A、求任意流体节点 
Figure 213786DEST_PATH_IMAGE001
处的位移值
Figure 406870DEST_PATH_IMAGE002
：在耦合界面
Figure 545728DEST_PATH_IMAGE003
上，有固体节点
Figure 191473DEST_PATH_IMAGE004
，见图2，节点
Figure 503505DEST_PATH_IMAGE005
相应的位置坐标为
Figure 488779DEST_PATH_IMAGE006
，其节点上的位移为
Figure 860854DEST_PATH_IMAGE007
，
Figure 931578DEST_PATH_IMAGE008
为任意流体节点
Figure 47302DEST_PATH_IMAGE001
的位置坐标，在拟合区域的一个局部子域上，拟合函数
Figure 887082DEST_PATH_IMAGE009
表示为
Figure 432989DEST_PATH_IMAGE010
                            （1）
式中
Figure 725430DEST_PATH_IMAGE011
为待求系数，它是坐标
Figure 644844DEST_PATH_IMAGE012
的函数，
Figure 401447DEST_PATH_IMAGE013
称为插值基函数，它是一个
Figure 53009DEST_PATH_IMAGE014
阶完备的多项式，
Figure 160642DEST_PATH_IMAGE015
是插值基函数的项数；例如对平面板与流体耦合，线性基
Figure 555851DEST_PATH_IMAGE016
，二次基
Figure 166961DEST_PATH_IMAGE017
，依次类推，还可取三次基等。加权离散范式，令
Figure 584353DEST_PATH_IMAGE019
           （2）
    式中
Figure 579991DEST_PATH_IMAGE020
是影响区域内节点的数目，
Figure 248870DEST_PATH_IMAGE021
是节点
Figure 570130DEST_PATH_IMAGE005
的权函数，常用的权函数是样
条函数，记
Figure 58880DEST_PATH_IMAGE022
，
Figure 126717DEST_PATH_IMAGE023
，则三次样条权函数
Figure 650103DEST_PATH_IMAGE024
                    (3.a)
Figure 876685DEST_PATH_IMAGE025
            (3.b)
Figure 180627DEST_PATH_IMAGE026
                                 (3.c)
为确定系数
Figure 721330DEST_PATH_IMAGE027
，式（2）应取最小值，将式（2）对
Figure 161538DEST_PATH_IMAGE028
求导，并令
Figure 762284DEST_PATH_IMAGE029
，得
Figure 287943DEST_PATH_IMAGE030
                            （4）
其中
Figure 632337DEST_PATH_IMAGE031
                          （5）
Figure 192631DEST_PATH_IMAGE032
               （6）
Figure 761016DEST_PATH_IMAGE033
                              （7）
将式（4）带入式（1），就可以得到MLS拟合函数
Figure 977233DEST_PATH_IMAGE034
                        （8）
式中
Figure 922056DEST_PATH_IMAGE035
称为形函数，
Figure 274540DEST_PATH_IMAGE014
表示基函数的阶数,
                   （9）
求得MLS拟合函数后,将耦合界面任意流体节点
Figure 203203DEST_PATH_IMAGE001
处位置坐标带入式(8),就可求得任意流体节点
Figure 472510DEST_PATH_IMAGE001
处的位移值
Figure 648277DEST_PATH_IMAGE002
；
B、求任意固体节点Si处的载荷值：根据步骤（A）中所述的方法，已知流体节点上的载荷值，同样可以获取耦合界面
Figure 573507DEST_PATH_IMAGE003
上的固体节点的载荷值，即可完成流体计算模型与薄壁结构计算模型进行的界面信息交换。
所述流体计算模型与薄壁结构计算模型在耦合界面满足应力和运动一致条件。 
如图3所示，所述步骤（4）中数据通讯管理和实施监控程序将流体计算模型、薄壁结构计算模型和界面信息交换分别设置为三个模块：流体计算模块，薄壁结构计算模块和信息交换模块，其中流体计算模块（CFD模块）实现用计算机编程语言对流体计算模型的数值计算，结构计算模块（CSD模块）实现用计算机编程语言对结构计算模型的数值计算，（界面）信息交换模块（插值程序）实现用计算机编程语言对最小二乘曲面拟合技术的数值计算。在每一个时间步内通过进行数据通讯管理和实施监控程序流体计算模型与薄壁结构计算模型的反复交错迭代的具体步骤包括：(Ⅰ)启动流体计算模块，执行信息交换模块（步骤A所述的方法），获取作用在结构上的流体载荷；(Ⅱ)启动薄壁结构计算模块，获得薄壁结构运动特征，执行信息交换模块，获取流固交界面上流体侧的位置，作为下一个交错迭代步流体计算的边界条件，继续执行步骤(Ⅰ)和(Ⅱ)，满足收敛标准后，跳到下一时间步，直到满足物理上所需的计算总时间，停止计算，即可完成流体计算模型与薄壁结构计算模型的反复交错迭代。 
所述收敛标准是指下一迭代步固体节点的位移值与上一迭代步固体节点位移值之比小于10-4。 
所述一个时间步是指把物理上连续时间分割成有限时间段，每一时间段在计算力学中就叫一个时间步；每一时间段的长短就叫时间步长；计算总时间就等于时间步长乘以总的时间步。 
所述薄壁结构动力学特性包括薄壁结构不同位置处的位移、应力、应变等物理量随时间的变化规律。 
本发明的优点和积极效果为：传统的流固耦合计算仅停留在小规模的科学研究，本方法提供的自动化技术手段和高效的信息交换技术，使这项技术在实际工程中的应用成为可能，并能改变过去只能采用简化手段进行薄壁结构动力学特性预测的现象。该计算方法准确性较以往的材料力学方法和单向耦合方法有了显著的提高，而且计算速度快，使流固耦合从理论研究走上实际工程应用。 
附图说明
图1本发明的流固耦合技术方案流程； 
图2本发明使用的界面插值技术；
图3本发明使用的交错迭代耦合技术方案流程；
图4本发明实施例采用的计算对象；
图5本发明实施例中监测点示意图；
图6本发明实施例中流场监测点湍涡频率演化图；
图7本发明实施例中结构监测点振动位移演化图。
图中：1-流体进口流向，2-流体出口流向，3-刚性方柱，4-耦合振动板。 
具体实施方式
下面结合附图对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。 
实施方式一：请参考图1，图1所示为本发明较佳实施例的基于交错迭代耦合技术的薄壁结构动力学特性预测方法流程图。本实施例包括下列步骤： 
（1）根据薄壁结构的初始位置构型及其周围环境流场，建立适体的基于有限体积的流体计算模型；流体计算模型的建立采用基于混合动态亚格子的大涡模拟模型、有限体积空间离散和全隐式时间离散格式。
（2）根据薄壁结构的物理力学性质、约束条件和步骤（1）计算得到的流体计算模型，建立基于有限单元的薄壁结构计算模型；薄壁结构模型计算采用基于有限单元的薄壁计算模型和直接积分的纽马克法（Newmark）。 
（3）将设计工况流量或流速作为流体计算的进口边界条件； 
（4）首先用最小二乘曲面拟合技术将流体计算模型与薄壁结构计算模型信息进行界面信息交换（确保计算准确性），然后通过数据通讯管理和实施监控程序在每一个时间步内进行流体计算模型与薄壁结构计算模型的反复交错迭代（解决耦合计算中的准确性），获得薄壁结构的流激振动响应和结构主应力分布，实现薄壁结构动力学特性预测。从而为结构安全设计提供依据。
步骤（4）中用最小二乘曲面拟合技术将流体计算模型与薄壁结构计算模型信息进行界面信息交换的方法为： 
A、求任意流体节点处的位移值：在耦合界面
Figure 533876DEST_PATH_IMAGE003
上，有固体节点
Figure 743141DEST_PATH_IMAGE004
，见图2，节点相应的位置坐标为
Figure 15039DEST_PATH_IMAGE006
，其节点上的位移为
Figure 470291DEST_PATH_IMAGE007
，
Figure 423642DEST_PATH_IMAGE008
为任意流体节点
Figure 52069DEST_PATH_IMAGE001
的位置坐标，在拟合区域的一个局部子域上，拟合函数表示为
Figure 776629DEST_PATH_IMAGE010
                            （1）
式中
Figure 960485DEST_PATH_IMAGE011
为待求系数，它是坐标
Figure 330287DEST_PATH_IMAGE012
的函数，
Figure 548778DEST_PATH_IMAGE013
称为插值基函数，它是一个
Figure 345833DEST_PATH_IMAGE014
阶完备的多项式，
Figure 16986DEST_PATH_IMAGE015
是插值基函数的项数；例如对平面板与流体耦合，线性基
Figure 987216DEST_PATH_IMAGE016
，二次基
Figure 997897DEST_PATH_IMAGE017
，依次类推，还可取三次基等。加权离散
Figure 28170DEST_PATH_IMAGE018
范式，令
Figure 124302DEST_PATH_IMAGE019
           （2）
    式中
Figure 901153DEST_PATH_IMAGE020
是影响区域内节点的数目，
Figure 94237DEST_PATH_IMAGE021
是节点
Figure 233094DEST_PATH_IMAGE005
的权函数，常用的权函数是样
条函数，记
Figure 613260DEST_PATH_IMAGE022
，
Figure 128555DEST_PATH_IMAGE023
，则三次样条权函数
Figure 176145DEST_PATH_IMAGE024
                    (3.a)
Figure 485904DEST_PATH_IMAGE025
            (3.b)
Figure 353366DEST_PATH_IMAGE026
                                 (3.c)
为确定系数
Figure 469089DEST_PATH_IMAGE027
，式（2）应取最小值，将式（2）对
Figure 308869DEST_PATH_IMAGE028
求导，并令
Figure 851846DEST_PATH_IMAGE029
，得
Figure 409866DEST_PATH_IMAGE030
                            （4）
其中
                          （5）
Figure 23567DEST_PATH_IMAGE032
               （6）
Figure 740375DEST_PATH_IMAGE033
                              （7）
将式（4）带入式（1），就可以得到MLS拟合函数
Figure 520112DEST_PATH_IMAGE034
                        （8）
式中
Figure 977638DEST_PATH_IMAGE035
称为形函数，
Figure 854328DEST_PATH_IMAGE014
表示基函数的阶数,
Figure 676790DEST_PATH_IMAGE036
                   （9）
求得MLS拟合函数后,将耦合界面
Figure 6140DEST_PATH_IMAGE003
任意流体节点
Figure 205040DEST_PATH_IMAGE001
处位置坐标带入式(8),就可求得任意流体节点
Figure 936236DEST_PATH_IMAGE001
处的位移值；
B、求任意固体节点Si处的载荷值：根据步骤（A）中所述的方法，已知流体节点上的载荷值，同样可以获取耦合界面
Figure 746246DEST_PATH_IMAGE003
上的固体节点的载荷值，即可完成流体计算模型与薄壁结构计算模型进行的界面信息交换。
流体计算模型与薄壁结构计算模型在耦合界面满足应力和运动一致条件。如图3所示步骤（4）中数据通讯管理和实施监控程序将流体计算模型、薄壁结构计算模型和界面信息交换分别设置为三个模块：流体计算模块，薄壁结构计算模块和信息交换模块，其中流体计算模块（CFD模块）实现用计算机编程语言对流体计算模型的数值计算，结构计算模块（CSD模块）实现用计算机编程语言对结构计算模型的数值计算，（界面）信息交换模块（插值程序）实现用计算机编程语言对最小二乘曲面拟合技术的数值计算。在每一个时间步内通过进行数据通讯管理和实施监控程序流体计算模型与薄壁结构计算模型的反复交错迭代的具体步骤包括：(Ⅰ)启动流体计算模块，执行信息交换模块（步骤A所述的方法），获取作用在结构上的流体载荷；(Ⅱ)启动薄壁结构计算模块，获得薄壁结构运动特征，执行信息交换模块，获取流固交界面上流体侧的位置，作为下一个交错迭代步流体计算的边界条件，继续执行步骤(Ⅰ)和(Ⅱ)，满足收敛标准后，跳到下一时间步，直到满足物理上所需的计算总时间，停止计算，即可完成流体计算模型与薄壁结构计算模型的反复交错迭代。收敛标准是指下一迭代步固体节点的位移值与上一迭代步固体节点位移值之比小于10-4。一个时间步是指把物理上连续时间分割成有限时间段，每一时间段在计算力学中就叫一个时间步；每一时间段的长短就叫时间步长；计算总时间就等于时间步长乘以总的时间步。 
作为流场中一类薄壁结构的典型代表，通过交错迭代耦合技术获取悬臂弹性薄板(长为300mm,宽为100mm,厚度为4mm)在方柱绕流形成的波动流场中的振动特性。计算模型如图4所示，流体及固体材料特性见表1。流动雷诺数定义为
Figure 748837DEST_PATH_IMAGE037
，其中
Figure 334539DEST_PATH_IMAGE038
为刚性方柱边长，为10cm。计算时。 
为了便于分析流固耦合中弹性薄板振动特性以及周围流场的流动特性，设置了4个监测点。薄板上有两个监测点，编号1-2，监测板的振动特性，薄板自由端附近设置了6个流场监测点，编号3-8，如图5所示。 
表1 流体及结构的材料属性
Figure 2013101172068100002DEST_PATH_IMAGE001
图7为结构上两个监测点的振动位移演化图。从图可见：（1）板自由端振动曲线呈现非周期、非对称的振动形态；（2）从测试或计算时间段内得到的振动信息看，都捕捉到了高位移出现的瞬间，此位移峰值比平均值要高出3倍左右，表明在流固耦合作用下，当某一瞬时流场扰动频率接近结构基频时，会产生共拍效应。如图6中点6的湍涡频率看, 在10.2s和21.7s时，湍涡频率接近60Hz，与结构的第二阶基频接近，发生“拍”现象，结构分别在12.9s和22.7s左右出现峰值振动位移；在34.3s时，湍涡频率接近10Hz，与结构的第一阶基频接近，再次发生“拍”现象，结构在34.6s左右出现最大峰值振动位移。同时也可以看出，这种瞬间的高位移振动并没有呈现明显的周期性，如自由端位移在12.9s出现峰值振动位移，而在22.7s时并没有出现峰值位移，而是出现突然升高－降低－再升高的振动现象，一方面，动态结构对流体激振荷载的响应产生“累积效应”的结果，同时结构振动反过来又会显著改变流体的运动特性，使流动频率背离结构的某一振动基频，从而使结构的振动位移发生突然改变；（3）实验测试得到的振动位移曲线和计算得到的振动位移曲线基本吻合，说明数值方法是可行的。（4）由上述较佳实施例实施结果看，对于薄壁结构流激振动特性，考虑流固耦合效应，会使结构的动力学特性更接近真实，能较好地预测到共拍效应，从而更好地指导结构设计。
虽然本发明已以较佳实施例揭露如上，然而并非限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。 

Claims (10)

1.一种基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于具体步骤包括如下：
（1）根据薄壁结构的初始位置构型及其周围环境流场，建立基于有限体积的流体计算模型；
（2）根据薄壁结构的物理力学性质、约束条件和步骤（1）的流体计算模型，建立基于有限单元的薄壁结构计算模型；
（3）将设计工况流量或流速作为流体计算的进口边界条件；
（4）首先用最小二乘曲面拟合技术将流体计算模型与薄壁结构计算模型信息进行界面信息交换，然后通过数据通讯管理和实施监控程序在每一个时间步内进行流体计算模型与薄壁结构计算模型的反复交错迭代，获得薄壁结构的流激振动响应和结构主应力分布，实现薄壁结构动力学特性预测。
2.根据权利要求1所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述步骤（1）中流体计算模型的建立采用基于混合动态亚格子的大涡模拟模型、有限体积空间离散和全隐式时间离散格式。
3.根据权利要求1所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述步骤（2）中的薄壁结构是水利水电工程以及航天航空领域广泛采用的由薄型板件和加强构件组成的结构。
4.根据权利要求1所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述步骤（2）中薄壁结构模型计算采用基于有限单元的薄壁计算模型和直接积分的纽马克法。
5.根据权利要求1所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于所述步骤（4）中用最小二乘曲面拟合技术将流体计算模型与薄壁结构计算模型信息进行界面信息交换的方法为：
A、求任意流体节点                                                
Figure 2013101172068100001DEST_PATH_IMAGE001
处的位移值
Figure 59845DEST_PATH_IMAGE002
：在耦合界面上，有固体节点
Figure 443553DEST_PATH_IMAGE004
，节点
Figure DEST_PATH_IMAGE005
相应的位置坐标为
Figure 698954DEST_PATH_IMAGE006
，其节点上的位移为
Figure DEST_PATH_IMAGE007
，
Figure 523297DEST_PATH_IMAGE008
为任意流体节点
Figure 106725DEST_PATH_IMAGE001
的位置坐标，在拟合区域的一个局部子域上，拟合函数
Figure DEST_PATH_IMAGE009
表示为
Figure 294124DEST_PATH_IMAGE010
                            （1）
式中
Figure DEST_PATH_IMAGE011
为待求系数，它是坐标
Figure 466348DEST_PATH_IMAGE012
的函数，
Figure DEST_PATH_IMAGE013
为插值基函数，
Figure 402206DEST_PATH_IMAGE014
是插值基函数的项数；
加权离散范式，令
           （2）
    式中是影响区域内节点的数目，
Figure 650970DEST_PATH_IMAGE018
是节点
Figure 490750DEST_PATH_IMAGE005
的权函数，记
Figure DEST_PATH_IMAGE019
，，则三次样条权函数
Figure DEST_PATH_IMAGE021
                    (3.a)
Figure 824703DEST_PATH_IMAGE022
            (3.b)
Figure DEST_PATH_IMAGE023
                                 (3.c)
    为确定系数
Figure 868752DEST_PATH_IMAGE024
，式（2）取最小值，将式（2）对求导，并令
Figure 500721DEST_PATH_IMAGE026
，得
Figure DEST_PATH_IMAGE027
                            （4）
    其中
Figure 840698DEST_PATH_IMAGE028
                          （5）
Figure DEST_PATH_IMAGE029
               （6）
Figure 10648DEST_PATH_IMAGE030
                              （7）
    将式（4）带入式（1），得到MLS拟合函数
Figure DEST_PATH_IMAGE031
                        （8）
    式中
Figure 29026DEST_PATH_IMAGE032
称为形函数，表示基函数的阶数,
Figure 515502DEST_PATH_IMAGE034
                   （9）
    求得MLS拟合函数后,将耦合界面任意流体节点
Figure 621310DEST_PATH_IMAGE001
处位置坐标带入式(8),求得任意流体节点
Figure 554631DEST_PATH_IMAGE001
处的位移值
Figure 410460DEST_PATH_IMAGE002
；
B、求任意固体节点处的载荷值：根据步骤（A）中所述的方法，已知流体节点上的载荷值，获取耦合界面上的固体节点的载荷值，完成流体计算模型与薄壁结构计算模型进行的界面信息交换。
6.根据权利要求1或5所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述流体计算模型与薄壁结构计算模型在耦合界面满足应力和运动一致条件。
7.根据权利要求1所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述步骤（4）中数据通讯管理和实施监控程序将流体计算模型、薄壁结构计算模型和界面信息交换分别设置为三个模块：流体计算模块，薄壁结构计算模块和信息交换模块，在每一个时间步内通过进行数据通讯管理和实施监控程序流体计算模型与薄壁结构计算模型的反复交错迭代的具体步骤包括：(Ⅰ)启动流体计算模块，执行信息交换模块，获取作用在结构上的流体载荷；(Ⅱ)启动薄壁结构计算模块，获得薄壁结构运动特征，执行信息交换模块，获取流固交界面上流体侧的位置，作为下一个交错迭代步流体计算的边界条件，继续执行步骤(Ⅰ)和(Ⅱ)，满足收敛标准后，跳到下一时间步，直到满足物理上所需的计算总时间，停止计算，即可完成流体计算模型与薄壁结构计算模型的反复交错迭代。
8.根据权利要求7所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述一个时间步是指把物理上连续的时间分割成有限时间段，每一时间段即为一个时间步；每一时间段的长短为时间步长；计算总时间为时间步长乘以总的时间步。
9.根据权利要求7所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述收敛标准是指下一迭代步固体节点的位移值与上一迭代步固体节点位移值之比小于10-4。
10.根据权利要求1所述的基于交错迭代耦合技术的薄壁结构动力学特性预测方法，其特征在于：所述薄壁结构动力学特性包括薄壁结构不同位置处的位移、应力、应变随时间的变化规律。
CN201310117206.8A 2013-04-07 2013-04-07 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法 Expired - Fee Related CN103177162B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310117206.8A CN103177162B (zh) 2013-04-07 2013-04-07 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310117206.8A CN103177162B (zh) 2013-04-07 2013-04-07 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法

Publications (2)

Publication Number Publication Date
CN103177162A true CN103177162A (zh) 2013-06-26
CN103177162B CN103177162B (zh) 2015-11-11

Family

ID=48637018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310117206.8A Expired - Fee Related CN103177162B (zh) 2013-04-07 2013-04-07 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法

Country Status (1)

Country Link
CN (1) CN103177162B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103853920A (zh) * 2014-02-24 2014-06-11 昆明理工大学 一种输流多壁碳纳米管动力学特性预测方法
CN105160092A (zh) * 2015-08-27 2015-12-16 中国运载火箭技术研究院 一种适用于热防护系统瞬态温度场计算的热环境插值方法
CN109241579A (zh) * 2018-08-15 2019-01-18 华北水利水电大学 一种考虑流体附加质量多维空间瞬变特性的快速建模方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2003524A2 (en) * 2007-06-15 2008-12-17 United Technologies Corporation Aeroelastic model
CN102663241A (zh) * 2012-03-29 2012-09-12 昆明理工大学 一种弹性水击下水轮机暂态功率的模拟计算方法
CN102913464A (zh) * 2012-09-05 2013-02-06 江苏大学 离心泵转子瞬态流固耦合特性预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2003524A2 (en) * 2007-06-15 2008-12-17 United Technologies Corporation Aeroelastic model
CN102663241A (zh) * 2012-03-29 2012-09-12 昆明理工大学 一种弹性水击下水轮机暂态功率的模拟计算方法
CN102913464A (zh) * 2012-09-05 2013-02-06 江苏大学 离心泵转子瞬态流固耦合特性预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
W.Q. WANG,ET AL: "Strongly coupling of partitioned fluid–solid interaction solvers using reduced-order models", 《APPLIED MATHEMATICAL MODELLING》 *
王文全: "薄壁结构流固耦合数值模拟及计算方法研究", 《中国博士学位论文全文数据库基础科学辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103853920A (zh) * 2014-02-24 2014-06-11 昆明理工大学 一种输流多壁碳纳米管动力学特性预测方法
CN105160092A (zh) * 2015-08-27 2015-12-16 中国运载火箭技术研究院 一种适用于热防护系统瞬态温度场计算的热环境插值方法
CN105160092B (zh) * 2015-08-27 2017-12-22 中国运载火箭技术研究院 一种适用于热防护系统瞬态温度场计算的热环境插值方法
CN109241579A (zh) * 2018-08-15 2019-01-18 华北水利水电大学 一种考虑流体附加质量多维空间瞬变特性的快速建模方法

Also Published As

Publication number Publication date
CN103177162B (zh) 2015-11-11

Similar Documents

Publication Publication Date Title
Wang et al. Turbulence modeling of deep dynamic stall at relatively low Reynolds number
Kang et al. Effects of flexibility on the aerodynamic performance of flapping wings
Zhang et al. Simulation of hydraulic transients in hydropower systems using the 1-D-3-D coupling approach
Cinnella et al. A numerical method for turbomachinery aeroelasticity
CN105631125A (zh) 基于降阶模型的气动-热-结构耦合分析方法
Guo et al. A CFD/CSD model for aeroelastic calculations of large-scale wind turbines
Xu et al. Correlation-based transition transport modeling for simulating crossflow instabilities
Shi et al. A coupling VWM/CFD/CSD method for rotor airload prediction
Min et al. Large-eddy simulation of corner separation in a compressor cascade
Rezvaya et al. Solving the hydrodynamical tasks using CFD programs
Prasad et al. Efficient prediction of classical flutter stability of turbomachinery blade using the boundary element type numerical method
CN103177162A (zh) 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法
Huang et al. Identification of flutter derivatives of bridge deck under multi-frequency vibration
Andersen et al. Simulation and prediction of wakes and wake interaction in wind farms
Ghommem Modeling and analysis for optimization of unsteady aeroelastic systems
Zhu et al. Identification of flutter derivatives of bridge decks using CFD-based discrete-time aerodynamic models
CN105138766B (zh) 基于模糊聚类的高超声速气动热降阶模型的加点方法
Rao Advances in aero structures
Douvi et al. Aerodynamic characteristics of s809 vs. NACA 0012 airfoil for wind turbine applications
Fabiano et al. Adjoint-based aerodynamic design on unstructured meshes
Fang et al. Design optimization of unsteady airfoils with continuous adjoint method
Fard et al. Energy harvesting performance of thick oscillating airfoils using a discrete vortex model
Nandi et al. Prediction and analysis of the nonsteady transitional boundary layer dynamics for flow over an oscillating wind turbine airfoil using the γ-reθ transition model
CN102663210B (zh) 利用建筑内表面设置空气薄层实现室内热环境数值模拟的方法
Yu et al. High-fidelity optimization of flapping airfoils for maximum propulsive efficiency

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151111

Termination date: 20210407

CF01 Termination of patent right due to non-payment of annual fee