matlab 流固耦合,关于Abaqus在流固耦合方面的应用研究.pdf

本文介绍了ABAQUS在流固耦合中的应用,特别是其处理理想气体状态方程的能力。通过对气体压缩过程的仿真与Matlab理论数值解的比较,验证了ABAQUS软件的可靠性。文中还展示了在特定产品分析中的实际应用。
摘要由CSDN通过智能技术生成

关于ABAQUS 在流固耦合方面的应用

左 辉 阚玉红 黄彦

中国兵器工业第二零八研究所 北京 102202

摘 要:针对 ABAQUS 有限元分析软件中的流固耦合功能,简述了其中理

想气体状态方程的各参数含义以及流固耦合的分析要点。文章通过ABAQUS 仿

真分析和理论数值解的对比,证明了ABAQUS 软件计算理想气体状态方程的可

信性,最后介绍其在某产品上的分析应用。

关键词: 理想气体方程 流固耦合 ABAQUS

前言

随着有限元技术的发展和用户要求的提高,各大有限元软件都含有流固耦合

模块,其主要用于液体、理想气体和JWL 的模拟,本文着重介绍ABAQUS 中理

想气体状态方程的功能和应用。为了验证ABAQUS 理想气体状态方程模拟气体

压缩的正确性,首先利用其模拟简单的气体压缩过程,并获得该过程中气体的状

态变化曲线 (仿真曲线);然后通过Matlab 求解该模型理论上的气体状态方程,

并获得数值解(理论曲线)。将仿真曲线和理论曲线进行对比,发现二者非常吻

合,证明了ABAQUS 模拟理想气体状态的可信性。在此基础上,将其用在某产

品上的流固耦合分析。

1 理想气体方程的参数含义

在ABAQUS 有限元分析软件中,气体压缩方程为:

p p R ( ) (1)

a z

其中: :气体压强的增量,初始增量为零,ODB 文件输出的压强

p

5

p a :

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
目 录 第1章 导论 1 1.1 引言 1 1.2 有限元分析与用户 1 1.3 本书的目的 1 1.4 本书的结构 2 第2章 杆单元 4 2.1 引言 4 2.2 一维杆单元 4 2.2.1 刚度矩阵:直接法 4 2.2.2 二维杆单元 6 2.3 整体刚度矩阵的组集 8 2.3.1 离散化 8 2.3.2 局部坐标系下的单元刚度矩阵 8 2.3.3 整体坐标系下的单元刚度矩阵 9 2.3.4 整体刚度矩阵的组集 10 2.3.5 整体力矢量 11 2.4 边界条件 12 2.4.1 一般情况 12 2.5 方程组求解 14 2.6 支反力 14 2.7 杆件的轴力 15 2.8 计算机程序:truss.m 16 2.8.1 数据准备 16 2.8.2 单元矩阵 18 2.8.3 整体刚度矩阵的组集过程 19 2.8.4 整体坐标系下力矢量的组集 19 2.8.5 整体方程组的解 19 2.8.6 节点位移 19 2.8.7 单元力 20 2.8.8 程序代码 20 2.9 问题求解 24 2.9.1 问题2.1 24 2.9.2 问题2.2 28 2.10 基于Abaqus的简单桁架分析 31 2.10.1 Abaqus概述 31 2.10.2 使用Abaqus的用户交互版本分析桁架 32 2.10.3 使用Abaqus的关键字版本分析桁架 47 第3章 梁单元 51 3.1 引言 51 3.2 刚度矩阵 52 3.3 均布载荷 55 3.4 中间铰链 58 3.5 计算机程序:beam.m 60 3.5.1 数据准备 60 3.5.2 整体刚度矩阵的组集和求解 63 3.5.3 节点位移 63 3.5.4 单元力 63 3.6 问题求解 67 3.6.1 问题3.1 67 3.6.2 问题3.2 70 3.6.3 问题3.3 73 3.7 基于Abaqus的简单梁模型分析 76 3.7.1 用户交互版本 76 3.7.2 关键字版本 85 第4章 刚架 89 4.1 引言 89 4.2 梁柱单元的刚度矩阵 89 4.3 两端铰接的梁柱单元的刚度矩阵 90 4.4 整体坐标系和局部坐标系 90 4.5 整体刚度矩阵的组集和未知位移求解 91 4.6 计算机程序:frame.m 91 4.6.1 数据准备 91 4.6.2 单元矩阵 93 4.6.3 整体刚度矩阵的组集 95 4.6.4 线性方程组的求解 95 4.6.5 节点位移 95 4.6.6 单元力 95 4.7 基于Abaqus的简单刚架分析 105 4.7.1 用户交互版本 105 4.7.2 关键字版本 111 第5章 应力和应变分析 114 5.1 引言 114 5.2 应力张量 114 5.2.1 定义 114 5.2.2 应力张量和应力矢量之间的关系 116 5.2.3 应力张量的变换 117 5.2.4 平衡方程 117 5.2.5 主应力 119 5.2.6 Mises应力 119 5.2.7 应力矢量的法向和切向分量 120 5.2.8 莫尔应力圆 121 5.2.9 应力的工程表示方法 121 5.3 变形和应变 122 5.3.1 定义 122 5.3.2 拉格朗日描述和欧拉描述 123 5.3.3 位移矢量 123 5.3.4 位移和变形梯度 124 5.3.5 格林应变矩阵 125 5.3.6 小变形理论 127 5.3.7 主应变 129 5.3.8 应变张量的变换 129 5.3.9 应变的工程表示方法 130 5.4 应力-应变的本构关系 130 5.4.1 广义胡克定律 130 5.4.2 材料的对称性 132 5.4.3 各向同性材料 134 5.4.4 平面应力与平面应变 137 5.5 问题求解 139 5.5.1 问题5.1 139 5.5.2 问题5.2 139 5.5.3 问题5.3 141 5.5.4 问题5.4 143 5.5.5 问题5.5 144 5.5.6 问题5.6 145 5.5.7 问题5.7 145 5.5.8 问题5.8 147 第6章 加权残值法 148 6.1 引言 148 6.2 基本方程 148 6.3 伽辽金法 148 6.4 伽辽金法的弱形式 150 6.5 二维或三维问题的分部积分法(格林公式) 151 6.6 瑞利-里兹法 154 6.6.1 定义 154 6.6.2 积分形式的函数表达式 155 6.6.3 瑞利-里兹法简介 155 6.6.4 自然泛
在固耦合问题中,可以使用Matlab来开发数值模拟软件,以解决航空发动机叶片颤振等具体问题。在开发这样的程序时,可以基于流固耦合的数学模型和有限元数值模型,加入非线性的影响,使模拟结果更贴近实际情况。然而,需要注意的是,目前国内在流固耦合技术方面仍处于初级阶段,还需要不断完善与发展。另外,还有一些流固耦合计算软件可供使用,如ADINA、CFX WB ANSYS、comsol、fluent、mpcci、abaqus等,大部分软件都支持弱耦合,部分软件还可以进行双向耦合。因此,如果您需要编写流固耦合Matlab程序,可以参考已有的流固耦合数学模型和有限元数值模型,结合Matlab的编程能力,开发出适用于具体问题的数值模拟软件。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [谈谈流固耦合计算](https://blog.csdn.net/weixin_39739234/article/details/116123658)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 50%"] - *2* [详讲流固耦合](https://blog.csdn.net/weixin_30651345/article/details/115851814)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值