网站建设建站流程方案,北京最大的网站建设有限公司,网站虚拟主机哪个好,大气手机网站模板免费下载Abaqus 导出单元刚度矩阵和全局刚度矩阵 首次创建#xff1a;2023.7.29 最后更新#xff1a;2023.7.29 如有什么改进的地方#xff0c;欢迎大家讨论#xff01; 详细情况请查阅#xff1a;Abaqus Analysis User’s Guide 一、Abaqus 导出单元刚度矩阵
1.生成单元刚度矩阵…
Abaqus 导出单元刚度矩阵和全局刚度矩阵 首次创建2023.7.29 最后更新2023.7.29 如有什么改进的地方欢迎大家讨论 详细情况请查阅Abaqus Analysis User’s Guide 一、Abaqus 导出单元刚度矩阵
1.生成单元刚度矩阵方法
导出单元刚度矩阵需要修改inp文件先在abaqus中创建一个单元集合这个集合是你想导出的单元集合这里我创建的是Element_ALL。然后生成inp文件在inp文件中的step关键字对中添加如下一行指令 *Element Matrix Output, ElsetElement_ALL, File NameElementStiffness, Output FileUser Defined, stiffnessyes
.....
.....
*Restart, write, frequency0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variablePRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variablePRESELECT
** 添加下面这一行
*Element Matrix Output, ElsetElement_ALL, File NameElementStiffness, Output FileUser Defined, stiffnessyes
*End Step2.单元刚度矩阵文件格式分析
这样在abaqus提交任务计算后就会在abaqus 的当前项目工作目录下生成一个ElementStiffness.mtx文件其格式如下
**
** ELEMENT NUMBER 1 STEP NUMBER 1 INCREMENT NUMBER 1
** ELEMENT TYPE C3D4
*USER ELEMENT, NODES 4, LINEAR
** ELEMENT NODES
** 2, 1, 9, 81, 2, 3
*MATRIX,TYPESTIFFNESS127.06125836015 ,0.0000000000000 , 183.76386705487 0.0000000000000 , -121.64117806010 , 388.01179556583 -39.618685020706 , -22.681043477885 , 48.656471224038 , 60.711115702304 -15.120695651923 , -96.321293715418 , 80.880129060096 , 37.801739129808 92.212564977143 ,32.437647482692 , 60.499604560096 , -81.962788877760 , -13.159037040064 -19.738555560096 , 42.381078873928 -173.27005949011 , 0.0000000000000 , 0.0000000000000 , 28.116306917523 0.0000000000000 , -53.846153846154 , 287.62647738911 0.0000000000000 , -173.27005949011 , 80.769230769231 , 0.0000000000000 28.116306917523 , -80.769230769231 , 0.0000000000000 , 287.62647738911 0.0000000000000 , 121.15384615385 , -606.44520821538 , -80.769230769231 -121.15384615385 , 98.407074211329 , 0.0000000000000 , 0.0000000000000 1006.6926708619 ,85.827486150661 , 22.681043477885 , -48.656471224038 , -49.208737599120 -22.681043477885 , 34.567543403526 , -142.47272481653 , 0.0000000000000 80.769230769231 , 105.85397626499 15.120695651923 , 85.827486150661 , -40.008181769231 , -15.120695651923 -24.007578179248 , 40.008181769231 , 0.0000000000000 , -142.47272481653 0.0000000000000 , 0.0000000000000 , 80.652816845114 -32.437647482692 , -60.012272653846 , 300.39620152731 , 45.271796585256 60.012272653846 , -58.825364207497 , 53.846153846154 , 0.0000000000000 -498.65453685784 , -66.680302948718 , 0.0000000000000 , 257.08369953803 ........从内容可以看出*MATRIX,TYPESTIFFNESS后面的就是一个单元刚度矩阵其形式为下三角形式。上面是一个线性四面体的单元刚度矩阵由于线性四面体有4个节点每个节点有3个自由度所以单元刚度矩阵的维度是12x12的 k e 12 x 12 ke_{12x12} ke12x12 整理后如下所示
127.061258360150.0000000000000183.763867054870.0000000000000-121.64117806010388.01179556583-39.618685020706-22.68104347788548.65647122403860.711115702304……………
二、Abaqus 导出全局刚度矩阵
1.生成全局刚度矩阵方法
在inp文件最后添加一个输出全局刚度矩阵的step如下所示
....
....
*Output, field, variablePRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variablePRESELECT
** 导出单元刚度矩阵
*Element Matrix Output, ElsetElement_ALL, File NameElementStiffness, Output FileUser Defined, stiffnessyes
**
*End Step**
** 导出全局刚度矩阵
*Step, nameGlobalStiffnessMatrix
*MATRIX GENERATE, STIFFNESS
*MATRIX OUTPUT, STIFFNESS, FORMATCOORDINATE
*End Step
**2.全局刚度矩阵文件格式分析
同样提交任务后将在工程的工作目录下生成一个以你任务名称命名的mtx文件如我的任务名称为Job-exportStiffness-globalstiffness那么将生成Job-exportStiffness-globalstiffness_STIF2.mtx文件这个文件就存储着全局刚度矩阵。部分内容如下
1 1 3.405254068249726e02
1 2 6.072817378766025e01
2 1 6.072817378766025e01
1 3 1.269469880224359e02
3 1 1.269469880224359e02
1 4 -5.756613438782665e01
4 1 -5.756613438782665e01
1 5 -3.763929516105769e01
5 1 -3.763929516105769e01
1 6 8.076923076923076e01
6 1 8.076923076923076e01
1 7 6.181074378769931e01
7 1 6.181074378769931e01
....
....可以看出这是采用三元数组的方式来存储的即一行数据前面个为矩阵的x,y位置索引最后一个为该索引位置的值。
