中国科学院流固耦合系统力学重点实验室

合集下载

输流波纹管流固耦合振动的初步研究

输流波纹管流固耦合振动的初步研究
水泵 , 一个 压 缩 机 , 路 , 纹 管 连 接 真 空 室外 部 管路 管 波
和真空 室 内部 管路 。通 过 冷却 水 与 晶体基 础 的对 流换 热, 达到冷却效果。其 中 , 水泵 由变频 器调节流量 , 工 作 时冷 却循 环 系统 处 于 连 续 运 转 状 态 , 路 中 的 阀 门 管
波纹管 元 件作 为输 流 管 道 的一 种 , 面 呈 波 纹状 , 壁 管 径 不断变 化 , 点 是 薄 壁 、 性 大 , 变 形 , 度 差 , 特 弹 易 刚
力、 水利 、 航空航 天、 核能等各领域 的管道系统都 采用 波纹 管作 为 管 路 连 接 件 。此 外 , 于 波 纹 管 的表 面 积 由 大 , 被 广泛用 于 换热 器 , 还 起到 加强 换 热 的效 果 。 输流 管道 的 振 动 是 一 个 典 型 的流 固耦 合 问 题 , 所 谓 流 固耦合 是 流 体 与 固体 之 间 的 交 互 作 用 不 能 分 开 : 变 形 固体在 流 体 载荷 作 用 下 会 产 生 变 形 或 运 动 , 变 而 形 或运 动又 反 过 来 影 响 流 场 , 而 改 变 流 体 载 荷 的分 从 布 和 大小 。 因而 , 法 显 式 地 消 去 描 述 流 体 运 动 的独 无 立 变量 或描 述 固体 运 动 的独 立 变 量 , 体 域 或 固体 域 流 均 不 可能单 独地 求解 。输 流 管道 的振 动 不是 单 纯 的流 体或其他 因素造成管道 的振动 , 而是管流之 间形成一 个整 体 发生振 动 , 为流 固耦 合振 动 。 故 对 于一般 的输 流 管 道振 动 , 引起 振 动 的原 因很 多 , 泵 的周期 性 运 动 , 门 的 开 闭 , 液 系 统 中吸 入 气 体 , 阀 输 高速 流 动 的失 稳 、 体 的可 压 缩 性 等 都 会 诱 发 管 路 振 流 动 。另外 , 境 中 的 振 动 , 泵 的机 械 振 动 , 动 机 转 环 如 电 子偏 心产 生 的 振 动 等 , 接 作 用 于 管路 和 整 个 系统 也 直 会导 致振 动 。通 常将输 流 管道 振 动 分 为 三 ’ : 引起 流 的管 道振 动及 稳 定性 , 究 具 有 定 常 流 速 的 不 可 压 缩 研 流体 与管 道弯 曲振 动 的耦 合 以及 在 此 流 速 下 的 管 道稳 定性 , 可称 为液 . 弹耦合 振 动 ; 非流 动 因素 引起 的管 道振 动 , 指忽 略 流速 的影 响 , 管道 振 动 归 因于 压 力 流量 是 将 的脉 动 、 锤及 机 械激 励 等 因素 , 类 耦 合 振 动 的 机理 水 这 有 摩擦 耦合 、 松耦 合 、 接 耦 合 ; 压 流 体 中 的声 波 泊 连 可 与 管道 振 动 的耦 合 , 种 耦 合 振 动 可 称 为 声一 耦 合 这 弹 振动。 对 输流 管道 的 振 动 研 究 的 常 用 计 算 方 法 有 : 限 有 元 法 、 征线 法 、 递 矩 阵 法 、 抗 分 析 法 和模 态 综 合 特 传 阻 法 。对 复杂 空 间 管 系 动 态 特 性 分 析 的 方 法 中 , 限元 有

流固耦合系统计算方法及应用软件概述

流固耦合系统计算方法及应用软件概述
种 各 样 的 液 固 耦 合 现 象 。 流 固 耦 合 问 题 可 由 其 耦 合 方 程 来 定 义。这组方程的定 义域 同时有 流体域 和 固体域 , 未知 变量含 而
显地分开 , 使描述物理现象 的方程 , 特别是本构 方程需要 针对具 体的物理 现象来建立 , 其耦 合效应 通过 描述 问题的微 分 方程而 体现。第 二大类 问题 的特征是耦合作 用仅仅发 生在两相 交界向 上, 在方程上 耦合 是 由两 相耦 合 面的 平衡 及协调 关 系 引人 的。 这类 问题 又 可 以具 体 的 分 为 三 种 情 况 。第 一 种 足 流 体 与固 体 结
问题 。张庆华 。 不考虑液体 的波动所 产生 的动压力 , 液体静 将
压力 作 为外 载 荷 通 过 A S S模 拟 加 载 , 用 有 限 元 模 拟 出 弹 性 NY 利 罐 体 的 变形 情 况 。王 华 , 刚 讨 论 了含 液 容 器 流 固耦 合 动 力 曹 特性 的有 限元 方 程 , A S S软 件 对 正 方 形 充 液 容 器 进 行 了模 用 NY 态分 析 , 出 其 理 论 模 态 特 性 , 果 表 明 该 模 型 有 效 可 靠 。徐 得 结
有描述流体现象 的变量及 描述 固体现象 的变 量。一般 而 言 , 具 有 以下 两 点 特 征 : () 1 流体域 或固体域均不 可能单独求解 ; () 2 无法显式地 消去描述 流体运 动的独立 变量 或描述 固体
运 动 的 独 立变 量 。

3・ 6
广州化 工
21 年 3 0 1 9卷第 1 期
1 2 化工机械充液系统流 固耦合计算方法的发展概述 .
充液系统 是 由固体和液体混合组成 的具有无 限多 自由度 的力学 系统 。在工程 实际 中, 化工 机械充液 系统 中液体 晃动 及 与 结 构 的耦 合 振 动 是 一 个 非 常 常 见 而 且 重 要 的 问 题 。这 类 耦 合 问 题 由 于 涉 及到 流 体 力 学 、 体 力 学 、 力 学 、 算 力 学 以及 高 动 计 等数学等学科 , 其机理 和本 质的探讨非 常复 杂和 困难。一般 根 据腔体旋转方式 不同 , 液体 晃动分为 绕对称轴 旋转运 动和 三 把 维 液 体 运 动 。如 航 天 器 上 的 液 体 燃 料 贮 箱 、 体 火 箭 等 属 于 前 液 者 , 对 于 像 油 罐 车 、 液 汽 车 、 机 油 箱 等 这 类 储 运 系 统 , 于 而 载 飞 属

基于 CFD 降阶的非线性气动弹性稳定性分析

基于 CFD 降阶的非线性气动弹性稳定性分析

基于 CFD 降阶的非线性气动弹性稳定性分析周强;陈刚;李跃明【摘要】为了快速寻找基于 CFD/CSD 的非线性气弹系统的颤振边界,根据Lyapunov 稳定性理论对非线性流固耦合系统进行稳定性分析。

首先通过微扰理论对非线性流固耦合系统处理建立近似线性化状态空间方程,再通过 POD (Proper Orthogonal Decomposition)方法将高维状态空间方程降阶为便于分析的降阶系统,根据系统所有特征值即可判定原始非线性系统稳定性。

Lyapunov 稳定性理论主要针对非线性系统,在实现过程中采用了 POD 降阶的方法,与大多数对降阶系统稳定性判定不同,其数学理论上是反映原始非线性流固耦合系统稳定性。

POD 降阶方法从系统内部流场出发,可以较好反映系统内部特性。

二维三维算例仿真结果验证了该方法预测颤振边界的正确性,分析发现,在亚音速阶段,稳定性主要由于结构模态主导;而在跨音速和超音速流动阶段,气弹稳定性主要由受流体特性影响。

%In order to quickly find the flutter boundary of a nonlinear aeroelstic system based on CFD /CSD,this paper used the Lyapunov stability theory to analyze the stability of a nonlinear fluid-structure coupling system.Firstly,a linearized state space model was obtained based on a nonlinear fluid-structure coupling system through the small perturbation theory;then a reduced order model was achieved by reducing the high dimensional linearized model through a POD(proper orthogonal decomposition)method.According to the eigenvalues of the ROM system,the stability of the original nonlinear system could be determined.Lyapunov stability theory was mainly aimed at the nonlinear system,and the POD method was adopted to realize the process.Differentfrom other stability methods,the mathematical theory reflected the stability of the original nonlinear fluid-structure coupling system.As the POD /ROM used in this paper was based on the CFD fluid field,so it could better reflect the internal characteristics of the original nonlinear coupling system.Numerical cases including two-dimensional aerofoil and three-dimensional wing's models were used to validate the accuracy of the stability method.The results show that in the subsonic domain,the stability is determined by the structural model;but in transonic and supersonic domain,the aeroelastic stability is mainly affected by fluid characteristics.【期刊名称】《振动与冲击》【年(卷),期】2016(035)016【总页数】7页(P17-23)【关键词】流固耦合;李雅普洛夫;正则正交分解;状态空间;颤振【作者】周强;陈刚;李跃明【作者单位】西安交通大学机械结构强度与振动国家重点实验室,西安 710049;西安交通大学机械结构强度与振动国家重点实验室,西安 710049;西安交通大学机械结构强度与振动国家重点实验室,西安 710049【正文语种】中文【中图分类】V211;V215.3飞行器气动弹性是研究飞行器结构与气动力相互耦合关系,是典型的流固耦合振动问题。

光滑粒子动力学方法在复杂流动中的研究进展

光滑粒子动力学方法在复杂流动中的研究进展

l f u i d s ;f l u i d — s t uc r ur t e i n t e r a c t i o n; me s h — l e s s me t h o d


精 度 ,而 网格 的重构 过程通 常又 异 常复杂 。由于无
网格方 法避 免 了复杂 的 网格 相关 操作 ,在 这方 面具
基于 网格 的传统 数值 方法 ,如有 限差 分 、有 限
体 积 、有 限元等 ,广泛 应用 于各 类流体 动 力学与 固
有 独特 优势 , 因此 自 2 0世 纪 9 O年代 以来 被充 分重 视 。较 常 见 的 无 网 格 方 法 包 括 光 滑 粒 子 动 力 学
第6 5卷 第 4期 f ¨ 2 0 1 4年 4月




Vl 0 1 . 6 5 No . 4 Ap r i l 2 01 4
CI ESC J o ur n a l

一 综一
光滑粒子动力学方法在 复杂流动 中的研究进 展
周光正, 葛蔚
( 中国科学院过程 工程研 究所 多相复杂系统 国家重点实验室,北京 1 0 0 1 9 0 )
B e i j i n g 1 0 0 1 9 0 , C h i n a )
A b s t r a c t :As a f u l l y me s h — l e s s , L a g r a n g i a n n u me r i c a l me t h o d , s mo o t h e d p a r t i c l e h y d r o d y n a mi c s( S P H)h a s
DOI :1 0 . 3 9 6 9  ̄ . i s s n . 0 4 3 8 — 1 1 5 7 . 2 0 1 4 . 0 4 . 0 0 1

隧道开挖过程中复杂裂隙围岩的固流耦合分析

隧道开挖过程中复杂裂隙围岩的固流耦合分析
第3 2卷 第 2期
21 O O年 O 月 4
J u n l fC vl o r a ii,Ar ht cu a o c i t r l& En io m e t l gn e ig e vr n n a En i e rn
土 木 建 筑 与 环 境 工 程
Vol 3 No. _2 2
隙岩体 水 力 学模 型 , 采用有 限元 法对 复 杂 裂 隙岩 体 中开 挖 隧道 的 固流 耦 合 过 程进 行 了数 值 模 拟 ,
模 拟 结果体 现 了主干 裂 隙在 渗 流 中的 强 导 水作 用 和 网络 状 裂 隙 的 贮水 功 能 与渗 流 滞后 效 应 , 开挖 过 程 中复 杂裂 隙岩体 渗流 场与 应 力场的耦 合作 用 显著 的增加 了隧道 围岩屈服 区。 关键 词 : 隙岩 体 ; 裂 固流耦 合 ; 隧道 中 图分 类号 : TU9 . 41 文献 标 志码 : A 文章编 号 :6 44 6 (O 0 0 —0 50 1 7 -7 4 2 1 ) 20 3 —5
Ab t a t s r c :W h n t ne s r hr ug wa e e un l a e t o h t 卜be rng e to f f a t e oc a s,a c a i s c i n o r c ur d r k m s ompl x i e a ton e nt r c i pr c s e we n r c a s de o m a i nd fui r n p r a i n i e o he m o ti p r a e s ns f r o e s b t e o k m s f r ton a l d t a s o t to s on f t s m o t nt r a o o

深海采矿扬矿管道工作特性的流固耦合分析与综合评价研究

深海采矿扬矿管道工作特性的流固耦合分析与综合评价研究

参考内容二
随着工业化的快速发展,大量的流体被用于各种工业生产过程中,而输流管 道作为传输流体的主要工具,其安全性与稳定性对工业生产具有重要意义。在输 流管道系统中,流固耦合振动现象是普遍存在的,它对管道的安全运行和稳定性 产生重要影响。因此,对输流管道流固耦合振动特性进行分析,对于预防和减少 管道事故,提高管道运行效率,具有重要的理论和实践价值。
管道系统中流固耦合机制尚不明确;(2)缺乏针对深海采矿扬矿管道工作 特性的综合评价方法。
研究方法
为了深入探讨深海采矿扬矿管道工作特性,本次演示采用流固耦合分析方法, 从理论和实验两个方面进行研究。首先,基于流体力学和固体力学的基本理论, 建立深海采矿扬矿管道的流固耦合模型;然后,利用数值计算方法对模型进行求 解,并设计实验进行验证。此外,本次演示还收集了大量实际工况数据,对模型 进行进一步优化和验证。
一、输流管道流固耦合振动的基 本概念
输流管道流固耦合振动是指管道内流体与管道壁之间的相互作用,导致管道 发生振动。这种振动不仅会对管道本身的结构产生影响,还会对管道内的流体产 生影响,反过来进一步影响管道的振动。这种相互影响、相互作用的物理现象, 即为输流管道的流固耦合振动。
二、输流管道流固耦合振动特性 的分析方法
感谢观看
深海采矿扬矿管道工作特性的 流固耦合分析与综合评价研究
01 引言
03 研究方法 05 结论
目录
02 文献综述 04 综合评价研究 06 参考内容
引言
随着海洋资源的日益枯竭,深海采矿作为一种新型的矿产资源开发方式,逐 渐引起了人们的。而在深海采矿过程中,扬矿管道的工作特性直接影响到整个采 矿系统的效率和安全性。因此,本次演示将针对深海采矿扬矿管道工作特性进行 流固耦合分析与综合评价研究,以期为提高深海采矿效率和安全性提供理论支持。

多物质流体动力学方法与结构动力学方法结合的流固耦合计算技术

多物质流体动力学方法与结构动力学方法结合的流固耦合计算技术
使 用 网格离散 , 或基 于点 离散 的方 法 . 系统之 间可 以使用 不 同 的相 互作 用 方 式 , 用 动 量 守恒 方 式 在 流 固 点 使 界 面处 理 非 滑 移 边 界 条 件 以 实 现 耦 合 . [ ] 展 了 一 种 S T S ( t izd S aeTme FudSrc- 文 5 发 S F I Sa l e p c・i li-t t bi u It atn 方法 , 于计算 血 液流 动过 程 中的 流 固耦 合 , n rc ) e o 用 其耦 合 原理也 是 基于 边界 条件 的处 理. 实 际上 L gag 过 程 和 E l 过 程 的 自由结合 可 以视 为一种 重 叠 网格 技术 . 网格 划 分 角度 看 , 种技 arn e ue r 从 这
提 高 几何 判 断 效 率 、避 免 祸 合 遗 漏 、 流 固计 算 时 间 步 协 调 等 处 理 技 术 . 应 用 建 立 的 流 固耦 合 计 算 程 序 模 拟 爆 炸 容器 中爆 炸 流 场 对 容 器 内部 结 构 、容 器 壳 体 的 作 用 过 程 . 计 算 结 果 表 明 其 可 以反 映 流 场 发 展 过 程 、流 场 结 构 相
互作用 和结构的运动变形和破坏过程 .
关键 词 : 固耦 合 ; 值 模 拟 ;多 物 质 流 体 动 力 学 程 序 ; 构 动 力 学 程 序 ; 炸 容 器 流 数 结 爆
中 图 分 类 号 :03 3:1 8 文 献 标 识 码 :A
0 引 言
在工 程实 际 中广泛 涉及 的流 体 和结 构相 互作 用 问题 , 可作 为流 固耦 合 问题进 行分 析 , : 立 一 种 基 于 重 叠 网格 空 间 覆 盖 判 断 的流 固耦 合 计 算 技 术 ,把 流 场 计 算 的 E l 多 物 质 流 体 动 力 学 计 算 建 ur e

基于相互作用积分方法的裂纹扩展分析

基于相互作用积分方法的裂纹扩展分析

基于相互作用积分方法的裂纹扩展分析陈旻炜;李敏;陈伟民【摘要】In the process of analyzing crack growth with traditional finite element method,crack tip singularity is an urgent problem that must be faced with,and then the direction and length of the propaga-tion are closely related to the mesh of model.In this paper,a model that is not thoroughly dependent on the mesh and can be used for analyzing problems of crack growth is set up.Because interaction integral method can relieve the effects of singularity to some extent,severe requirements for meshing caused by singularity can be avoided.Stress intensity factor at the crack tip is solved by interaction integral via finite element method,and direction of the propagation is determined by the maximum circumferential tensile stress theory.When cracks propagate,only the meshes of concerned areas are regenerated in order to ob-tain the final extension path in a relatively precise way.At last,several typical examples,including growth of mode-I crack,cracks with a hole plate and cracks subj ected to compression,are compared with refer-ences.Reliability of the method in this paper is verified and applicability in many damaged patterns is shown as well.%使用传统有限元方法对裂纹扩展问题进行分析时,裂尖奇异性是首先必须要面对的问题,其次裂纹扩展的方向、长度与模型的网格密切相关.本文通过建立一个不完全依赖于网格的模型来分析裂纹扩展问题,利用相互作用积分方法来求解裂尖参数,积分方法可以一定程度上消除奇异性影响,从而避免了裂纹尖端的奇异性对网格划分的苛刻要求.采用相互作用积分法则求解裂纹尖端的应力强度因子,由最大周向拉应力理论来判断裂纹的扩展方向.模拟裂纹扩展的过程中,只对涉及扩展区域的网格进行重新剖分来获得最终的扩展路径以保证计算的精度.最后通过I型裂纹扩展、带孔板裂纹扩展及受压裂纹扩展等几个典型算例与参考文献进行对比,证实了本文方法具有良好的可靠性,并展示了在多种破坏模式下方法的适用性.【期刊名称】《湖南大学学报(自然科学版)》【年(卷),期】2018(045)004【总页数】8页(P74-81)【关键词】裂尖奇异性;相互作用积分;应力强度因子;有限元方法;裂纹扩展【作者】陈旻炜;李敏;陈伟民【作者单位】北京航空航天大学航空科学与工程学院,北京 100191;北京航空航天大学航空科学与工程学院,北京 100191;中国科学院力学研究所流固耦合力学重点实验室,北京 100190【正文语种】中文【中图分类】O346.1关于裂纹扩展问题的研究,对于预防工程结构的断裂失效、评价结构的完整性与安全性具有非常重要的意义.在非常规油气资源的开采中,往往通过水力压裂技术增大岩层中裂纹的表面积来实现提高产量的目的,这也促进了众学者对裂纹扩展问题的研究.裂纹扩展问题属于不定边界问题,具有高度的非线性[1].它是一个连续面与非连续面不断转换的过程,而该问题的关键是如何实现对不连续面的准确刻画.通常用来分析裂纹扩展的方法主要有边界元方法和有限单元方法.边界元法(BEM)[2-3]在模拟复杂的模型时,严重依赖于问题的基本解,大大降低了处理非线性、非均质等问题时的优势.有限元方法中的单元消除(杀死)技术也可以简单、有效地实现有限单元法对于断裂破坏过程的模拟[4].为了解决有限元模型中重划分的问题[5],学者们在传统有限元基础上引入了结点释放技术[6]以及扩展有限元法[7-8].结点释放技术需要预先知道裂纹潜在的扩展方向,因此在模拟复杂情况下的裂纹扩展问题时会比较麻烦.扩展有限元法将裂缝独立于计算网格之外,利用扩充位移项来描述不连续的位移场从而表征裂缝特征.该方法模拟精度高,裂缝处理方便,但针对不同的裂纹问题需要选取不同的富集函数.Belytschko和Black[9]针对二维裂纹问题引入Westergaard函数作为富集函数(Enriched functions);Moes等[10]则引入阶跃函数作为富集函数;对于一些尚无解析解的问题,Waisman和Belytschko[11]还提出了一种带未知参数的富集函数方法来提高计算精度.另外,裂缝条数过多时会引入附加自由度,计算速度慢,收敛性较差.本文采用有限元方法来进行裂纹扩展的数值模拟,无论采用何种模式,对于裂纹扩展数值模拟问题,裂尖应力场的分析均是关键问题.由于裂纹尖端是奇异点,所以在利用有限元法求解裂尖应力场时,其结果的准确性与单元类型和网格尺寸都密切相关,最好在特定位置采用奇异单元与足够小的单元尺寸.在模拟裂纹扩展的过程中完全满足此类要求是非常困难的,因为计算时需要跟踪裂纹几何位置的变化不断重构网格并变更单元类型.本文不需要直接面对有限元算不准的裂尖奇异应力场,而是通过裂纹尖端的应力强度因子衡量裂纹的开裂特性,并且它还可以判断出裂纹的开裂方向,避免了需要预先定义裂纹扩展路径的问题.应力强度因子是采用相互作用积分方法计算得到,而积分方法可以减轻计算结果对裂尖处单元类型与单元尺寸的依赖.该方法只需要通过远离裂尖处某一积分路径上的应力场解与位移场解,便可计算出裂尖的应力强度因子.远场处的有限元结果往往具有很高的精度,所以该方法能够获得较为满意的结果,这样就给裂纹扩展的分析带来了极大的方便.本文给出了进行裂纹扩展分析所需使用的策略,力求在计算含奇异点的裂纹扩展问题时寻求一种较为可靠的方法.为了提高可用与易用性,在通用有限元软件上进行简单的改进,实现了ANSYS软件的二次开发,计算了几个比较有代表性的断裂力学问题,比较单元杀死方法和扩展有限元的结果以及网格划分的影响,验证了本文方法的正确性与有效性.1 裂纹扩展的理论基础1.1 应力强度因子的计算应力强度因子常见的计算方法有相互作用积分方法、结点位移外推法[12]以及复变函数方法[13].本文的应力强度因子是通过相互作用积分方法得到的.然而,相互作用积分方法是基于准静态条件下的计算方法,所以无法求解裂纹扩展速度等动态参数.具体的基本理论如下[14].如图1所示,相互作用积分是包含裂纹尖端回路G的能量积分.它的定义如下[7]: I=(1)式中考虑了两个独立的平衡状态,状态为真实应力-变形场的变量,状态为附加应力-变形场的变量;为相互作用应变能密度;nj为积分回路上的外法线单位矢量.附加场可以是任何满足变形方程和本构方程的应力-变形场,一般情况下取为裂尖渐近场的解析解[15].图1 相互作用积分示意图Fig.1 Schematic diagram of interaction integration 当回路G 接近裂纹尖端时,相互作用积分与真实变形场和附加变形场的应力强度因子之间有如下关系:(2)式中E*为平面应变或平面应力下的弹性模量.在实际应用中,为方便有限元的计算,通常把相互作用积分转换成一个面积分,在回路G 的外面再设置另外一个回路C,两个回路之间所围成的区域为S,通过散度定理将式(1)的线积分转换为面积分,另外再考虑裂纹面上的面力作用,则相互作用积分就可以写为[16]I=(3)式中:q是平滑的权函数,在区域S内取值为1,其余地方为0;my是裂纹面上的外法线.通过以上分析可以知道,选取特定的辅助场就可以得到真实状态下的应力强度因子.为了求解参数可以假设辅助场为纯I型的状态,并且满足如下条件(4)由于纯I型状态下的应力场和位移场都有理论解,因此很容易通过式(3)得到相互作用积分的值,最后结合式(2)和式(4)就可以得到参数的表达式(5)同理假设辅助场为纯II型状态就可以得到参数的结果.1.2 裂纹扩展准则典型的断裂准则包括能量释放率、J积分和应力强度因子等相关准则.本文采用应力强度因子准则,在相互作用积分方法所求得裂纹尖端i应力强度因子和的基础上,定义等效应力强度因子[17]:(6)当裂尖的等效应力强度因子超过临界应力强度因子时,裂尖就开始扩展.对于复合型的裂纹扩展,应用最大周向拉应力强度因子理论来确定出裂纹扩展的方向角[18]: θi=(7)1.3 模拟裂纹扩展的过程在裂纹开裂过程的计算中,主要是通过计算确定出裂纹的起裂条件与裂纹的开裂角.另外一个问题就是裂纹扩展时模型网格的重新划分,重新划分的区域只需要选择所关心的扩展部分即可,因此只是增加了部分的结点自由度,详细的划分尺寸可以由下文的算例得到.由于利用相互作用积分方法求解裂纹尖端的应力强度因子,所以在网格划分时不需要采用特殊单元,这样就降低了网格重新划分的难度.具体的计算流程如下:1) 通过相互作用积分方法计算出裂纹尖端i处的应力强度因子和然后利用式(5)得到裂纹开裂方向角θi,最后根据复合裂纹断裂准则确定开裂点.2) 设置裂纹尖端的开裂步长为Δl,沿裂纹扩展方向角劈裂模型.关于开裂步长的设置,参考文献[19]中给出的建议增长步长为其中h为单元尺寸.3) 重新划分模型网格,需要特别注意裂纹尖端附近的网格尺寸,并防止重新划分的网格出现畸变.4) 重复过程1) ~3),进行下一个开裂点的操作.如果开裂过程中不满足开裂准则那么裂纹扩展过程中止.2 典型问题的分析与比较在ANSYS有限元程序中,模型均采用八结点的PLANE183平面应变单元.分析接触问题时,接触单元采用三结点CONTA172单元,目标单元选择三结点TARGE169单元.在选择的这些典型问题中,将单元杀死方法、相互作用积分方法和扩展有限元方法进行了比较.首先比较了单元杀死方法和文章中的相互作用积分方法,目的是说明本文方法要比单元杀死方法计算更正确.其次比较了相互作用积分方法与扩展有限元方法,进一步证明了本文方法是有效可靠的.下面详细讨论这些算例.2.1 应力强度因子的计算应力强度因子的计算是本文方法中判断裂纹开裂的主要参数,其对模型网格的敏感性是首先必须要说明的问题.只有部分裂纹的应力强度因子有封闭解,所以本文选择有限宽中心裂纹板的解析解为例来说明应力强度因子的敏感性计算问题.该计算模型如图2所示,单向受均匀拉应力σ的有限尺寸板,宽为2W,高为2H,中心存在一条长为2a的裂纹.假设板为各向均匀的线弹性材料,裂纹尖端的应力强度因子表达式为[20]:(8)图2 含中心裂纹的有限大板Fig.2 The finite plate with a central crack其中F为有限宽度的修正系数,代表性系数包括以下3种.Isida近似系数:对a/W<0.7时误差为0.3%.修正的Koiter公式:误差为0.3 %.修正的Feddersen公式:误差为0.1 %.本算例中a/W的取值为0.2.对计算出来的应力强度因子进行无量纲化,记为Kn,Kn=KI/(σ表1中列出了不同的应力强度因子修正结果与本文计算得到的结果.从表1中可以看出,计算得到的结果与手册中的解析结果非常接近.表明了本文求解应力强度因子的方法是可靠的.表1 裂纹尖端的无量纲应力强度因子Tab.1 Normalized stress intensity factors at the crack tip求解方法Kn相对误差/%Isida方法1.02520.18Koiter方法1.0224-0.10Feddersen方法1.02440.10本文方法1.0234—通常在采用奇异单元求解应力强度因子的方案时,其计算结果对网格的尺寸有较高的要求[21],这样在某些问题中就会因追求精度而降低了计算效率,并且不利于对裂纹扩展的计算.因此表2针对不同网格尺寸与类型所得到的结果进行了对比.结果表明网格的改变对最终结果的影响不大,并且对计算值的影响也很小,建议网格尺寸的取值为0.1a.由此可见网格划分对该方法的影响相比加密网格的传统有限元方法来说是不敏感的,并且不依赖于裂纹尖端近场应力的数据,所以本文方法具有一定的优越性.2.2 I型单裂纹的扩展模拟例证2为最简单的I型裂纹扩展问题.图3为该问题的模型简图,板的边缘有一条长度为I的初始裂纹,板的左右两侧承受σ的拉应力作用.最终的计算结果还与利用单元杀死方法的结果进行了对比并绘制出图4,试图通过这两种方法的对比说明本文方法预测裂纹扩展路径的可靠性.其中的实线为使用本文方法所计算出来的结果,而带三角符号的线(映射划分网格)、带四边形符号的线(自由划分网格1)和带十字星符号的线(自由划分网格2)均为单元杀死方法的结果.自由划分网格1与网格2的区别是网格尺寸的不同,网格1的尺寸是网格2尺寸的二倍.其余结果的网格尺寸与自由划分网格1的尺寸相同.表2 不同网格尺寸下裂纹尖端的无量纲应力强度因子Tab.2 Normalized stressintensity factors at the crack tip with different mesh size网格类型求解方法应力强度因子参考值0.02a0.05a0.10a0.20a计算值相对误差/%计算值相对误差/%计算值相对误差/%计算值相对误差/%三角形网格Isida方法1.02521.02450.071.02310.201.02340.181.0220.31Koiter方法1.02241.0245-0.211.0231-0.071.0234-0.101.0220.04Feddersen方法1.02441.0245-0.011.02310.131.02340.101.0220.23四边形网格Isida方法1.02521.02440.081.02240.271.02310.201.02150.36Koiter方法1.02241.0244-0.201.022401.0231-0.071.02150.09Feddersen方法1.02441.024401.02240.201.02310.131.02150.28注:a为裂纹半长.图3 I型裂纹拉伸的模型示意图Fig.3 Model of mode-I crack subjected to tension图4 I型裂纹拉伸的扩展结果图Fig.4 Propagation results of mode-I subjected to tension从图4中可以看出虽然各条曲线的趋势都是向前扩展,但是具体细节上还是存在差异.裂纹扩展过程中带四边形的线与带十字星的线都出现了偏折,这是因为单元杀死方法只能按照预先设定好的单元分布情况进行扩展,在自由划分的网格中往往不能保证网格线平直.然而带四边形的线与带十字星的线的扩展路径并不重合,这说明单元杀死方法在网格尺寸不同的情况下计算结果出现了差异.图中带三角形的线与实线基本重合,在映射网格下单元杀死方法可以捕捉到简单的裂纹开裂规律.虽然单元杀死方法可以预测出裂纹扩展的总体趋势,但是它严重依赖于网格划分的方式,并且由于网格划分的问题会使裂纹扩展路径不同,不利于后续问题的研究,因此不是一种稳定可靠的计算方法.利用本文方法计算出来的虚线不受网格划分的影响,才最接近真实的裂纹扩展结果.2.3 带孔板边缘多裂纹扩展模拟在验证完简单的裂纹扩展模型后,通过较复杂的例证对本文方法进行深度的校验.用来验证的模型示意图如图 5所示,旨在研究受圆孔缺陷的影响下裂纹的扩展规律.具体描述为含两个中心对称圆孔的矩形平板,底部固定,顶端承受预先给定的位移加载,使它们只能在竖直方向上有移动.两条长度均为1 mm的预制裂纹在孔洞离边界的较远一侧,并且呈中心对称分布.模型的材料参数如下:弹性模量E=200 GPa,泊松比ν=0.3,材料的临界应力强度因子KIC =1 500 N/mm3/2. 图5 带孔板边缘多裂纹扩展的模型示意图(单位:mm)Fig.5 Model of multiple edge cracks propagation with hole plate为了验证本文方法的可靠性,本文将通过改变单元类型的相互作用积分方法计算得到的裂纹扩展路径与文献中Khoei等人[22]的路径进行对比.裂纹的扩展均需要预先假定一个阈值,但是由于裂纹尖端具有奇异性,因此一般方法在处理裂纹尖端问题时对裂尖网格的划分要求比较高,即使单元类型相同,而尺寸的差异也会导致裂纹扩展路径的变化.依据之前介绍的计算步骤,本文方法所计算出来的裂纹扩展径图如图 6中带四边形实线与带三角形实线所示,而Khoei等人的计算结果则用实线标识.从图 6中可以看出本文的计算结果与Khoei的结果[22]都显示出相似的裂纹扩展路径.裂纹均因圆孔的影响而在开裂初期朝圆孔的方向扩展,在后期由于两条裂纹之间的相互作用而互相靠拢.由此可知本文的计算结果与文献的计算结果吻合较好,并且不同的网格划分方式对于路径的影响小.裂纹扩展不同阶段的网格划分图如图 7所示,裂纹每次向前扩展一定距离,在裂纹附近区域的网格都会进行重构,并且裂纹附近的网格尺寸比其他区域的尺寸更小,原因除了保证一定的精度外还因为能较好地保证网格的划分质量.图6 裂纹扩展路径对比图Fig.6 Graph of crack growth with different methods综上所述,本文利用的相互作用积分方法,对网格的划分方式并没有特别严格的要求,既适用于三角形网格又适用于四边形网格,而且扩展路径也与文献中的结果相差不大,因此是一种更可靠的方法.2.4 受压裂纹扩展的模拟在地层中有很多岩层都处于受压状态,然而其中的裂纹依然会产生错动而发生开裂.例如在水力压裂工程中,由于其具有围压的特殊性,天然裂纹会受水压裂纹以及地应力的影响而扩展,若在进行数值模拟时不考虑天然裂纹的压裂纹情况则会造成偏差.本文从这一角度出发,探讨利用同样的方法是否可以解决裂纹在受压条件下的扩展问题.图7 裂纹扩展不同阶段的单元划分图Fig.7 Element partition plot of crack growth path at different stage裂纹面在受压时会相互接触并且产生摩擦力作用,因此必须考虑接触力作用,而且为了防止裂纹面之间发生相互渗透,模拟时需要在裂纹面之间设置接触分析,一般来说就是用接触单元来模拟接触面的作用.其中较为常见的分析方法有罚函数法和增广型拉格朗日乘子法等.但是罚函数法存在收敛性差的问题[23],因此本文采用的是增广型拉格朗日乘子法,通过判断双面约束来进行计算.计算过程中实时监测了接触面上的接触压力,保证接触面上不会产生破坏.模型采用与岩层类似的材料参数,具体数据为:弹性模量E =0.33 GPa,泊松比ν= 0.19,Keq=1.4N/mm3/2.模型的几何尺寸如图 8所示,模型长60 mm,高120 mm,裂纹长为a=12 mm,裂纹之间的竖向距离为b=20 mm,裂纹的倾斜角均为45°,在模型的底部和顶端均有压应力的作用.为了验证本文方法的正确性,将计算结果与石路杨等人[17]利用扩展有限元方法所得到的结果进行对比.图9(a)中的黑折线为本文方法计算出来的裂纹扩展路径,而图9(b)中的白折线为文献中的计算结果,可以看出两种方法所得到的结果基本一致.从图9中的计算路径可以得出,模型在受压状态下,裂纹的扩展过程始终是沿着模型的最大主应力进行的,这与初始裂纹考虑了摩擦力有关[18].除了初始裂纹面处于受压力作用外,扩展生成的新裂纹面是处于受拉作用的,因此这种破坏形式属于劈裂破坏.图 8 受压裂纹扩展的模型示意图Fig. 8 Model of crack subjected to compression(a)本文计算结果 (b)文献计算结果图9 裂纹扩展路径图Fig.9 Crack growth path如果受压的裂纹面之间没有摩擦,则裂纹的开裂路径则与带摩擦的情形有显著区别,角度大小与摩擦系数有关.裂纹的起裂方向不会沿着模型的最大主应力方向,在本例中与竖直方向有一个大约35°的夹角.在扩展一段时间后,裂纹才会逐渐转向并沿着最大主应力方向继续扩展.裂纹扩展路径图可详见图10,其中虚线代表无摩擦条件下的扩展路径,而实线则表示有摩擦的情况.该方法可以嵌入通用的有限元软件中,例如ANSYS,ABAQUS等.由于ANSYS擅长进行多物理场的分析,而水力压裂工程的模拟就是流固耦合的计算,因此基于本文所述的方法结合通用有限元平台进行流固耦合分析就较容易实现水力压裂中的裂纹扩展.这也是笔者今后要进一步深入研究的方向,因为依托大型通用有限元软件,所以计算方法也很容易推广.图10 无摩擦与有摩擦下的裂纹扩展路径图Fig.10 Crack growth paths under friction and frictionless condition3 结论本文提供了一种模拟裂纹扩展的计算方法.按照本文提出的计算流程与策略,便可以较正确地得到裂纹的扩展路径.利用相互作用积分方法来求解裂尖参数可以一定程度上消除奇异性影响,从而避免了裂纹尖端的奇异性对网格划分的苛刻要求.在传统有限元方法的基础上,通过相互作用积分法则求解裂纹尖端的应力强度因子,由最大周向拉应力理论来判断裂纹的扩展方向.模拟裂纹扩展的过程中,对涉及扩展区域的网格进行重新剖分并获得了较正确的扩展路径与较高的计算精度.应用该方法在模拟裂纹扩展的过程中,兼顾计算精度与效率的情况下建议网格尺寸取0.1倍的裂纹半长.该方法既能模拟裂纹拉伸时的情况也可以模拟受压时的扩展.相比裂纹受拉伸情况下的扩展,裂纹受压时的扩展需要设置接触,因此在计算时间上相对较长,但在网格尺寸方面并没有差别.在分析岩石破坏的过程中,由于材料存在较大的分散性,因此一般只关心岩层中裂纹的开裂方向,而对精度往往没有过高的要求.今后希望可以应用这种可靠性的方法来做进一步的研究.参考文献[1] 陶伟明, 班勇婷, 王慰军,等. 一种裂纹扩展的有限元仿真分析方法及其实现[J].力学季刊,2009, 30(4): 612-617.TAO Weiming,BAN Yongting,WANG Weijun,et al. A finite element simulation procedure and its realization for crack growth[J].Chinese Quarterly of Mechanics,2009,30(4):612-617. (In Chinese)[2] PORTELA A, ALIABADI M H, ROOKE D P. The dual boundary element method: effective implementation for crack problems[J]. International Journal for Numerical Methods in Engineering, 1992,33(6):1269-1287. [3] 雷钧, 黄祎丰, 杨庆生. 裂纹动态扩展的边界元模拟[J] . 北京工业大学学报, 2013, 39(6): 806-810.LEI Jun,HUANG Yifeng,YANG Qingsheng. Dynamic crack propagation by boundary element method[J].Journal of Beijing University of Technology,2013,39(6):806-810. (In Chinese)[4] 王强, 宋雪迪, 郝中华,等. 基于ABAQUS纤维梁单元的钢筋混凝土柱受力破坏全过程数值模拟[J]. 土木工程学报, 2014, 47(12):16-26.WANG Qiang, SONG Xuedi, HAO Zhonghua,et al.Numerical simulation of failure process of RC columns based on fiber beam element in ABAQUS[J]. China Civil Engineering Journal, 2014,47(12):16-26. (In Chinese)[5] BAO J Q, FATHI E, AMERI S. A coupled finite element method for the numerical simulation of hydraulic fracturing with a condensation technique[J]. Engineering Fracture Mechanics, 2014,131:269-281.[6] FU P, JOHNSON S M, CARRIGAN C R. An explicitly coupled hydro-geomechanical model for simulating hydraulic fracturing in arbitrary discrete fracture networks[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2012,37:278-300.[7] DOLBOW J, BELYTSCHKO T. A finite element method for crack growth without remeshing[J].International Journal for Numerical Methods in Engineering,1999,46:131-150.[8] 曾青冬, 姚军. 水平井多裂缝同步扩展数值模拟[J]. 石油学报,2015,36(12):1571-1579.ZENG Qingdong,YAO Jun.Numerical simulation of multiple fractures simultaneous propagation in horizontal well[J].Acta PetroleiSinica,2015,36(12):1571-1579. (In Chinese)[9] BELYTSCHKO T, BLACK T. Elastic crack growth in finite elements with minimal remeshing[J]. International Journal for Numerical Methods in Engineering, 1999, 45: 601-620.[10] MOES N, DOLBOW J, BELYTSCHKO T. A finite element method for crack growth without remeshing[J]. International Journal for NumericalMethods in Engineering, 1999, 46: 131-150.[11] WAISMAN H, BELYTSCHKO T. Parametric enrichment adaptivity by the extended finite element method[J]. International Journal for Numerical Methods in Engineering, 2008, 73: 1671-1692.[12] 肖汉宁, 彭苏华, 高朋召. 支柱瓷绝缘子表面裂纹应力强度因子的有限元分析[J]. 湖南大学学报(自然科学版), 2014, 41(12):53-58.XIAO Hanning, PENG Suhua, GAO Pengzhao. Stress intensity factor analysis of surface crack on rod porcelain insulator through finite element method[J]. Journal of Hunan University(Natural Sciences), 2014, 41(12):53-58. (In Chinese)[13] 肖万伸, 欧阳成, 张春雨. 一维六方准晶中界面裂纹与位错的交互作用[J]. 湖南大学学报(自然科学版), 2016, 43(12):71-75.XIAO Wanshen, OUYANG Cheng, ZHANG Chunyu. Interaction between circular interfacial crack and screw dislocation in 1d hexagonal quasicrystals[J]. Journal of Hunan University(Natural Sciences), 2016,43(12):71-75. (In Chinese)[14] 庄茁,柳占立,成斌斌,等. 扩展有限单元法[M]. 北京:清华大学出版社, 2012:20-21.ZHUANG Zhuo,LIU Zhanli,CHENG Binbin,et al.Extended finite element method[M].Beijing:Tsinghua University Press,2012:20-21.(In Chinese) [15] YAU J F, WANG S S, CORTEN H T. A mixed-mode crack analysis of isotropic solids using conservation laws of elasticity[J]. Journal of Applied Mechanics, 1980, 47(2):335.[16] DOLBOW J, N, BELYTSCHKO T. An extended finite element method。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

中国科学院流固耦合系统力学重点实验室Key Laboratory for Mechanics in Fluid Solid Coupling Systems Institute of Mechanics, Chinese Academy of Sciences季报2019年第1期(总第17期)目录中科院流固耦合系统力学重点实验室现场评估工作顺利完成 (2)中科院流固耦合系统力学重点实验室召开2019年室务会 (3)中国航空学会空气动力学分会飞行载荷专业工作会在扬州召开 (6)圆柱阵列波浪力幅值的波动现象和预报公式 (8)轻质金属点阵圆柱壳结构制备与力学性能研究进展 (9)力学所提出一种大幅提升3D打印点阵结构力学性能的新方法 (11)雾化稠油掺稀降粘技术研究进展 (12)南海天然气水合物试采安全评价研究进展 (14)油气水多相流量计研究进展 (15)空化致板间液滴界面稳定性研究获得多个奖项 (16)空泡与柔性膜的流固耦合研究获得2019度中国力学大会优秀墙报奖. 18中科院流固耦合系统力学重点实验室现场评估工作顺利完成7月15日,中科院前沿科学与教育局、中科院重点实验室现场评估专家组一行14人莅临中科院力学所,对依托力学所建设的流固耦合系统力学重点实验室进行现场评估。

专家组组长顾逸东院士主持了评估会议并宣布了现场评估的议程安排。

力学所所长秦伟,党委书记、副所长刘桂菊,副所长魏宇杰,副所长尹明及流固耦合系统力学重点实验室学术委员会主任、实验室主任参加会议。

实验室主任黄晨光做实验室主任工作报告,围绕发展定位与研究方向、科研任务与代表性成果、队伍建设与人才培养、开放交流与运行管理等方面,向专家组汇报了评估期内的发展成果和工作成效。

杨国伟研究员、王展研究员分别做“高速列车气动设计与流固耦合动力学特性研究”和“极端海洋环境及其与工程结构的流固耦合理论”代表性成果报告。

专家组肯定了实验室取得的成绩以及工作亮点,并就汇报和自评估报告中的存疑事项进行了交流。

现场评估专家组还查看了高速列车动模型试验平台、海洋流固土耦合实验室、多相流体力学实验室、冲击与耦合效应实验室的科研仪器建设、大型科研仪器设备使用共享等情况,同时,参观了实验室的展板窗口。

在此基础上,专家组召开会议,根据现场考核情况对实验室进行打分,并初步形成了评估意见。

经过努力,实验室顺利完成了此次中科院重点实验室现场评估工作,并在评估中充分展现了自身的优势和特色,最终取得良好的评估成绩。

在国家科技创新基地优化整合的背景下,实验室将积极适应新形势和新要求,进一步加强实验室建设和运行管理工作,全面提升科研平台建设水平和运行效率,为加快科技创新提供良好的条件支撑。

(流固耦合系统力学重点实验室供稿)中科院流固耦合系统力学重点实验室召开2019年室务会9月11日至9月12日,中科院流固耦合系统力学重点实验室(LMFS)在力学所怀柔园区召开2019年度室战略研讨会暨室基金评选及中期进展评审会议。

实验室学术委员会主任李家春院士,实验室副主任杨国伟、周济福、高福平及实验室科研人员等60余人出席了本次会议。

9月11日下午举行了实验室发展战略研讨会,李家春院士首先致开幕词。

他指出,召开战略研讨会是实验室的一项重要工作,它有助于更好地协调不同课题组之间的选题方向,梳理、聚焦实验室的学科发展,以期更好地服务于国家的重大战略需求,在海洋工程、环境工程和交通工程等应用领域做出实验室应有的贡献。

在下一步的工作中,还是要紧紧围绕研究流体与固体之间的相互作用规律这根主线不动摇,重点研究气、水、土与固体结构之间的相互作用。

在研究人员配置方面,要形成一支老中青相结合、团结协作的高水平研究队伍,更要重点关注年轻人的成长,为实验室未来的发展储备力量。

随后,各个课题组的组长汇报了各组近期取得的成果及下一步的研究规划,与会领导及科研人员对此进行了讨论。

大家一致表示,下一步的工作中,在解决重大工程问题的同时更要注重提炼其中所蕴含的科学问题,从学科发展的角度促进基础理论的进步,与此同时引领相关工程技术的进步。

在保持住实验室在高速列车空气动力学与结构、极端海洋环境与海洋工程、环境工程、非常规能源力学理论和开采技术以及高速水动力学和高超热流固耦合等方面优势地位的基础上,不同课题组之间要学会取长补短、协同作战,并加强与其他单位的沟通和交流,以此来寻求新的学科增长点。

在促进力学学科发展和服务国家安全与国民经济发展等方面走在前列。

次日举行了实验室科研基金申请及中期进展评审会议。

室科研基金从成立伊始至今,一直为众多刚刚独立开展科研工作的年轻人保驾护航,激发了他们的创新活力,培育和孵化了一大批具有较大发展潜力的科研方向。

基金的获得者在国家级科研项目申请、原创成果发表及专利申请方面都取得了不俗的成绩,有力地支撑了实验室的整体发展。

经过评选,今年有12项课题获得了室科研基金的资助。

同时,去年获得资助的10项课题进行了中期进展汇报,各个课题都如期甚至超额完成了申请时所定下的研究目标,并有多个项目于今年获得了国家自然科学基金的资助。

流固耦合系统力学重点实验室2019年室务会于9月12日下午结束,圆满完成了预定的各项会议议程。

今后,实验室将继续定期组织会议共同探讨实验室未来发展,凝聚共识,团结一致,为祖国科技事业的发展贡献力量!图1.实验室发展战略研讨图2.实验室基金评选(流固耦合系统力学重点实验室供稿)2019 APCATS暨AJSAE会议总结2019年08月28日至08月31日,第十届亚太航天技术和科学会议(The 10th Asian-Pacific Conference on Aerospace Technology and Science,简称APCATS)暨第四届亚洲航天工程联合研讨会(The 4th Asian Joint Symposium on Aerospace Engineering简称AJSAE)在中国台湾新竹市隆重召开。

本次学术会议由台湾交通大学承办,来自中国大陆、中国台湾、韩国、日本、印度、新加坡等地区的100多名航天领域著名学者出席了会议。

APCATS与AJSAE重点关注航空航天技术和科学领域的最新进展与创新,致力于为科学家和工程师提供一个讨论最新理念、方法和成果的平台,促进以亚太地区为主的航天科技工作者在不同领域展开交流与合作。

中国科学院力学研究所流固耦合与系统力学实验室的杨国伟研究员受邀为大会委员会成员出席了本次会议,同行的还有赵伟研究员、申义庆研究员、郭迪龙高级工程师、郑冠男高级工程师、孙振旭副研究员、银波副研究员、汪球副研究员,分别作了题为“Moving Model Rig and Its Application for High-Speed Trains”,“Study on a New Gas Wave Refrigerator with Rotating Gas Distrib utor”,“A New Fifth-order WENO-Z Scheme and Its Applications”,“Effects of bogies on the wake flow of a high-speed train”,“The Investigation into the flutter of missile rudder’s surface under disturbance flow”,“Aero-Acoustic optimization of the wing mock-up configuration”,“Investigation of Flexibility Effect on the Aerodynamics ofInsect Wing in Complicated Environment”和“Influence of model material on the accuracy of transient heat transfer in impulse facilities”的学术报告,并和与会代表深入交流,进行了热烈的讨论,建立了良好的合作关系,会议取得圆满成功。

(流固耦合与数值计算课题组供稿)中国航空学会空气动力学分会飞行载荷专业工作会在扬州召开由中国航空学会空气动力学分会、沈阳飞机设计研究所扬州协同创新研究院、中国科学院力学研究所共同主办的中国航空学会空气动力学分会飞行载荷专业工作会议于2019年10月9-11日在沈阳飞机设计研究所扬州协同创新研究院召开。

来自中国科学院、清华大学、北京航空航天大学以及航空工业空气动力研究院、沈阳所、直升机所、上飞院、航天十一院等32所国内外知名院校、科研单位的飞行载荷领域专家教授参会,共同探讨飞行载荷专业的发展方向,并针对飞行载荷领域的设计难题开展了技术交流。

随着我国航空航天领域科学技术的迅猛发展,飞行载荷设计在有/无人驾驶的军用、民用和通用飞行器设计领域的重要性与日俱增。

此次会议针对模飞试验、动载荷设计、高校飞行载荷专业设置以及非固定翼飞行器飞行载荷设计等议题开展了深入探讨,致力于提升飞行载荷设计水平,不断加大科技创新投入力度,为航空航天领域飞行载荷设计技术的发展提供了有益参考。

图1.会议代表合影图2.会场照片图3.会场照片(流固耦合与数值计算课题组供稿)圆柱阵列波浪力幅值的波动现象和预报公式圆柱阵列与波浪的相互作用是海洋工程中的重要问题。

在考虑入射波和绕射波的情况下,对一长排穿透水面的大直径圆柱阵列中每个圆柱所受波浪力变化规律进行研究,有助于深刻理解波浪与大数量圆柱的相互作用机理,而且所得结果对于波浪力变化曲线的高效计算、提高疲劳载荷谱预报精度也有实际意义。

结构力学课题组在上述研究中取得重要进展,相关结果以“Fluctuation of magnitude of wave loads for a long array of bottom-mounted cylinders”为题发表于流体力学杂志:Journal of Fluid Mechanics, 2019, 868: 244-285。

结构力学课题组通过理论分析和数值计算,发现在波浪力随无量纲波数变化曲线的非捕获区中,波浪力幅值的波动间距不随无量纲波数变化,该间距只与圆柱阵列总数、圆柱在阵列中的位置、波浪入射角有关;基于水波绕射的相互作用理论和波的干涉理论这两种方法,通过理论研究和数学推导获得了波动间距、非捕获区峰谷值横坐标的预报公式,所获的多个理论公式给出的预报结果与实验和数值计算结果吻合良好;揭示了该波动现象和理论预报公式所蕴含的内在物理机制;最后,讨论了上述工作的工程实际意义。

图1. 波浪力幅值曲线非捕获区的极值点和波动间距(结构力学课题组供稿)轻质金属点阵圆柱壳结构制备与力学性能研究进展柱壳结构是广泛应用于火箭、导弹及飞行器的主承力构件,在服役过程中往往承受较大的轴向过载,如何既满足轻量化指标、又提高轴压屈曲性能是一项有挑战性的工作。

相关文档
最新文档