再论ANSYS中的总体矩阵提取(In Python)

菲尼克斯方 · 发表于 2013-3-31 11:11

你好，我看了你发的帖子，但我现在遇到一个问题不知道该怎么解决。我用ANSYS的hbmat方法提取出了结构的总刚度矩阵，但是是Harwell-boeing格式的矩阵，有没有什么方法或者程序能让它转换成普通的矩阵。我看了王新敏编的那本ansys书里的例题，但是运行有错误

JLBhaidao · 发表于 2013-4-20 00:01

技术牛贴，必须顶

yunze · 发表于 2013-5-5 23:00

特定用户是什么，我下不了呀，同志们

qinzh_86 · 发表于 2014-5-27 10:09

好东西谢谢楼主的分析

yhchen · 发表于 2014-12-24 16:27

mark 一下，感谢楼主分享
顺便请教一下，我导出来的质量阵有为0的列和行，这是为什么呢

yhchen · 发表于 2014-12-24 16:32

Rainyboy 发表于 2011-5-17 17:25
回复 6 # 心灯的帖子

哦，这样啊，看看这个帖子呗：

楼主，这个链接失效了，能在共享一个吗

雨之冥 · 发表于 2015-4-3 16:51

好贴，必须学习之

Jmantan · 发表于 2015-5-29 19:51

Rainyboy 发表于 2011-4-12 18:22
补充一下，如何得到矩阵的各行与模型各自由度之间的映射关系：
在HBMAT命令中打开mapping开关，即形如（最 ...

请问：怎么看mapping生成的文件啊？比如我想把HBMAT生成的矩阵怎么处理成安节点编号排序的矩阵？

Rainyboy · 发表于 2015-5-30 16:52

本帖最后由 Rainyboy 于 2015-5-30 10:15 编辑

Jmantan 发表于 2015-5-29 12:51
请问：怎么看mapping生成的文件啊？比如我想把HBMAT生成的矩阵怎么处理成安节点编号排序的矩阵？

例如，KK_RHS.mapping文件内容为：

Matrix Eqn       Node DOF
         1          3 UX
         2          3 UY
         3          3 ROTZ
         4          2 UX
         5          2 UY
         6          2 ROTZ

数据的意思是：
总体自由度向量第1行表示 3节点 UX位移
总体自由度向量第2行表示 3节点 UY位移
总体自由度向量第3行表示 3节点 ROTZ位移
总体自由度向量第4行表示 2节点 UX位移

以此类推。
关于你的第二个问题，以上面的数据为例，原始的动力方程（以静力学为例）为：

{\bf K}\left(\begin{array}{c}{\bf u}_3\\ {\bf u}_2\end{array}\right) =\left(\begin{array}{c}{\bf f}_3\\ {\bf f}_2\end{array}\right)

根据你想要的按照节点顺序排列，可以写出坐标变换：

\left(\begin{array}{c}{\bf u}_3\\ {\bf u}_2\end{array}\right)=\begin{bmatrix}{\bf 0} & {\bf I} \\{\bf I} & {\bf 0} \end{bmatrix}\left(\begin{array}{c}{\bf u}_2\\ {\bf u}_3\end{array}\right) ={\bf T}\left(\begin{array}{c}{\bf u}_2\\ {\bf u}_3\end{array}\right)

以及

\left(\begin{array}{c}{\bf f}_3\\ {\bf f}_2\end{array}\right)={\bf T}\left(\begin{array}{c}{\bf f}_2\\ {\bf f}_3\end{array}\right)

将这个变换代入初始的动力方程，可得：

{\bf T}^{-1}{\bf K}{\bf T}\left(\begin{array}{c}{\bf u}_2\\ {\bf u}_3\end{array}\right) =\left(\begin{array}{c}{\bf f}_2\\ {\bf f}_3\end{array}\right)

那么{\bf T}^{-1}{\bf K}{\bf T}就是你要的新刚度矩阵了。

关键是构造出{\bf T}矩阵，但这个也不难，只要知道新旧自由度向量中的对应位置就可以了。

Jmantan · 发表于 2015-5-30 23:52

Rainyboy 发表于 2015-5-30 16:52
例如，KK_RHS.mapping文件内容为：

Matrix Eqn Node DOF

先谢谢详细的回答，我等会再仔细看看，

Amonia · 发表于 2016-11-25 11:04

权限不够，下载不了package，好桑梓

think2015 · 发表于 2016-11-28 08:47

Amonia 发表于 2016-11-25 11:04
权限不够，下载不了package，好桑梓

多参与论坛讨论升级就好了

Amonia · 发表于 2016-11-28 11:05

think2015 发表于 2016-11-28 08:47
多参与论坛讨论升级就好了

嗯嗯

Edinburgh · 发表于 2016-11-29 08:47

楼主真用心谢谢

ESCI_HEU · 发表于 2017-7-1 17:29

谢谢楼主的好贴

账号		自动登录	找回密码
密码			我要加入

[前后处理] 再论ANSYS中的总体矩阵提取(In Python)

点评

点评