学习opensees时的几个问题

xiaomm09

学习opensees时的几个问题，欢迎讨论指点，谢谢.1.对自由度编号的RCM算法不是很明白，手册上说是使带宽最小化，但没看到相关的文献。
2.system中提供了很多的方法BandGeneral 、BandSPD、ProfileSPD、SparseGeneral 、
UmfPack、SparseSPD，相互之间有什么区别及作用条件？
3.定义截面后，构件的B和H是如何放置的？几何变换应该只是局部到整体坐标的一种坐标变换。

bob704

学习opensees时的几个问题，欢迎讨论指点，谢谢.1.对自由度编号的RCM算法不是很明白，手册上说是使带宽最小化，但没看到相关的文献。
2.system中提供了很多的方法BandGeneral 、BandSPD、ProfileSPD、SparseGeneral ...
xiaomm09 发表于 2010-10-10 14:37

实际上对第1、2个问题，用户不用彻底弄的很明白（从应用角度来说），因为关于自由度编号优化有很多种方法，甚至有些求解方程组方法不用优化节点编号（因为其不存储所有0元素）；对于system，其主要是用于求解方程组，解法也很多，例如硕士生课程里有数值方法这门课，里面也讲了很多方法。对于自己编过程序的网友，可能理解的更深一点，当然若有研究的网友给出相应的参考文献更好了。
第3个问题我目前也正在研究，有时候可能是局部坐标与整体坐标关系，在用ANSYS做杆系结构时，有时是试着来的，看看图形是什么，若能把OS里数据写成ANSYS的APDL，导入到ANSYS里直观看看就什么问题也没有了，也有助于理解。

xiaomm09

感谢bob704，确实不必要弄太清楚，属于不同的算法，只是觉得有些问题不弄明白总感觉像黑匣子，积累多了就糊涂了，呵呵。
期待有朋友上传相关文献。谢过。

chb1973

支持，帮顶

linguan118

对于2D梁单元，BH方向已经定死了，对于3D梁单元，只需给出与xz平面平行的方向
http://opensees.berkeley.edu/wiki/index.php/Linear_Transformation

snowbird

多多学习！

cobalt

对于杆系结构，当对刚度矩阵采用一维变带宽存储时，使用RCM（Reverse-Cuthill-Mckee）算法对节点编号进行优化，可以很大程度上减小刚度矩阵的最大半带宽和存储元素总个数。
我以前编制过RCM算法，对于813个节点、1090个单元的桥梁杆系有限元模型，最大半带宽只有未优化时的4.39%，一维存储元素总个数只有未优化时的11.06%。效果非常明显。
关于RCM方法的详细介绍，可以参见下面的论文：
[1] 黄志超,包忠祤,周天瑞. 有限元结点编号优化[J]. 南昌大学学报(理科版),Vol.28(3),2004.
[2] 贾建军,彭颖红. 三种基于图论的有限元结点编号优化算法 [J]. 机械科学与技术, Vol.17(5),1998.