搜索
整体刚度阵提取方法方法一1.利用ansys建立模型，施加约束与载荷，并进行一次模态分析2.利用如下命令生成单元刚度阵、质量阵、阻尼阵/solantype,7seopt,matname,1nsel,allm,all,allsolveselist,matname,3seopt,matname,1：1-生成刚度阵；2-生成刚度阵质量阵；3-生成刚度阵、质量阵、阻尼阵selist,matname,3：0-列出summary；1-列出排序；2-列出排序+载荷右端项；3-列出排序+载荷右端项+各矩阵；3.此方法适用于规模较小的模型。4.此方法总刚度阵是按节点从小到大顺序排序的。方法二1.利用ansys建立模型，施加约束与载荷，并进行一次模态分析。FINISH/SOLUEMATWRITE,YESANTYPE,2MODOPT,LANB,2FINISH2.利用如下命令生成整体刚度阵、质量阵、阻尼阵的HB格式文件，其中HB文件格式说明见附注。FINISH/AUX2FILE,'file','full',''HBMAT,'k','txt','',ASCII,STIFF,YESFINISHHBMAT,'m','txt','',ASCII,MASS,YESHBMAT,'c','txt','',ASCII,DAMP,NOYES-列出载荷右端项；NO-不列出载荷右端项；3.利用matlab程序“hb_to_msm.m”即可得到matlab稀疏矩阵形式的MCK阵，程序见附录2function[a,rhsval]=hb_to_msm(input_file)M、C、K：稀疏矩阵；rhsval：载荷右端项input_file：步骤2导出的数据文件4.该方法适用于大型结构，复杂模型。5.该方法生成的整体刚度阵，是按优化后节点序列排序的，需要找到各节点在刚度阵中的对应的位置，此处通过EMAT文件获取，其中EMAT文件的格式说明见附录。FINISH/AUX2FORM,LONGFILEAUX2,'file','emat',''DUMP,5,FINISH6.利用matlab程序“readAnsysSort2.m”即可读取EMAT文件，程序见附录3functionNodeDofs=readAnsysSort2(file1,file2)file1：通过5导出的EMAT文件，即“DUMP.lis”file2：模型约束文件节点号约束自由度111111121212121212UXUYUZROTXROTYROTZUXUYUZROTXROTYROTZCONST.lis文件格式