abaqus dat文件 matlab_Abaqus/用Abaqus进行参数分析(python+maltab助攻)

(原文发在微信公众号“建筑工业产品经理”)


有时候需要对模型做大量的参数分析以尝试获得一些规律认识,这时若仅仅通过abaqus的cae界面手动修改参数以及前后处理,几乎不可能完成大量的分析(例如需要修改模型参数进行数百甚至上千次分析),这就需要通过写程序来提高效率。这个技能很实用,下面简述最近做此分析的过程,作为备忘。

step1:利用abaqus建立一个参数的模型并得到其inp文件。

step2:通过matlab批量修改inp文件中的参数,获得用于参数分析的批量的n个inp文件,例如要进行125个算例就可以获得125个inp文件。

参数分析一般需要选择好哪几个变量,这几个变量分别用那几个数值,例如三个参数各取值5个数,就有125个参数组合,这些参数组合提前准备好,并通过代码批量替换inp文件中的特定数据,得到125个inp文件。

step3:对批量的inp文件进行计算分析(在.bat文件中写入对n个inp的计算job,将bat文件放在abaqus工作目录中运行,abaqus可依次计算这n个job),得到n个odb文件;

step4:在abaqus的cae界面中对其中一个odb文件录制宏,得abaqus自动生成的.py文件(python语言代码写的后处理操作);

step5:修改与宏所对应的.py文件,加入循环语句,进而自动对n个odb文件进行后处理,分别提取所需要的结果。将所关心的数据写入到txt文件中,至此得到了含有n个算例下所关心的数据;

step6:在matlab中进一步分析n个算例的数据结果,得到参数分析所需要的图表,讨论与分析。

小结:整个过程似乎有点繁琐,是因为我对于python并不熟悉,所以有些工作通过matlab来完成。也许有更简单的办法,有朋友分享一下么?

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 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 自然泛
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值