|
% 目的：驗(yàn)證ANSYS中輸出整體剛度矩陣K和整體質(zhì)量矩陣M % % % 驗(yàn)證基礎(chǔ)： % (1) 基于模態(tài)分析中，驗(yàn)證ANSYS振型是否關(guān)于質(zhì)量矩陣M歸一化 % (2) 特征值分析，由K,M計(jì)算特征值，即頻率和振型，與ANSYS計(jì)算結(jié)果對(duì)比 % % 注意事項(xiàng)： % （1） ANSYS輸出矩陣以Harwell-Boeing格式存儲(chǔ)，并且由于K,M的對(duì)稱性 % % % （2） % % % % clc;clear; % 提取K,M K1=full(hb_to_msm('stiffness')); M1=full(hb_to_msm('mass')); % 還原對(duì)角線以上元素，使之為對(duì)稱陣 K =diag(diag(K1)/2)+K1-diag(diag(K1)); K = K+K'; M =diag(diag(M1)/2)+M1-diag(diag(M1)); M = M+M'; % 特征值分析，即計(jì)算頻率Freq和振型Phi [X1,Lamda]=eig(K,M); [nd,nd]=size(K); fdof=[1:18]'; pdof=[1,2,17]';% 約束自由度 fdof(pdof)=[]; [nfdof,nfdof]=size(X1); for j=1:nfdof; end [L,i]=sort(diag(Lamda)); X2=X1(:,i); Phi=zeros(18,nfdof); Phi(fdof,:)=X2; Freq = sqrt(L)/(2*pi); % 讀取ANSYS輸出的頻率與振型; Freq_ansys = load('Freq.txt'); for i =1:6 end % 計(jì)算兩者誤差 Err_F = norm(Freq(1:6)-Freq_ansys(1:6)); Err_Phi = norm(Phi(:,1:6)-Phi_ansys(:,1:6)); % 驗(yàn)證歸一化 Phi_ansys(pdof,:)=[]; test = Phi_ansys'*M*Phi_ansys; bar3(test) |
相關(guān)文章
- 2021-09-08BIM技術(shù)叢書Revit軟件應(yīng)用系列Autodesk Revit族詳解 [
- 2021-09-08全國(guó)專業(yè)技術(shù)人員計(jì)算機(jī)應(yīng)用能力考試用書 AutoCAD2004
- 2021-09-08EXCEL在工作中的應(yīng)用 制表、數(shù)據(jù)處理及宏應(yīng)用PDF下載
- 2021-08-30從零開始AutoCAD 2014中文版機(jī)械制圖基礎(chǔ)培訓(xùn)教程 [李
- 2021-08-30從零開始AutoCAD 2014中文版建筑制圖基礎(chǔ)培訓(xùn)教程 [朱
- 2021-08-30電氣CAD實(shí)例教程AutoCAD 2010中文版 [左昉 等編著] 20
- 2021-08-30電影風(fēng)暴2：Maya影像實(shí)拍與三維合成攻略PDF下載
- 2021-08-30高等院校藝術(shù)設(shè)計(jì)案例教程中文版AutoCAD 建筑設(shè)計(jì)案例
- 2021-08-29環(huán)境藝術(shù)制圖AutoCAD [徐幼光 編著] 2013年P(guān)DF下載
- 2021-08-29機(jī)械A(chǔ)utoCAD 項(xiàng)目教程 第3版 [繆希偉 主編] 2012年P(guān)DF