|
% 目的：驗證ANSYS中輸出整體剛度矩陣K和整體質量矩陣M % % % 驗證基礎： % (1) 基于模態分析中，驗證ANSYS振型是否關于質量矩陣M歸一化 % (2) 特征值分析，由K,M計算特征值，即頻率和振型，與ANSYS計算結果對比 % % 注意事項： % （1） ANSYS輸出矩陣以Harwell-Boeing格式存儲，并且由于K,M的對稱性 % % % （2） % % % % clc;clear; % 提取K,M K1=full(hb_to_msm('stiffness')); M1=full(hb_to_msm('mass')); % 還原對角線以上元素，使之為對稱陣 K =diag(diag(K1)/2)+K1-diag(diag(K1)); K = K+K'; M =diag(diag(M1)/2)+M1-diag(diag(M1)); M = M+M'; % 特征值分析，即計算頻率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 % 計算兩者誤差 Err_F = norm(Freq(1:6)-Freq_ansys(1:6)); Err_Phi = norm(Phi(:,1:6)-Phi_ansys(:,1:6)); % 驗證歸一化 Phi_ansys(pdof,:)=[]; test = Phi_ansys'*M*Phi_ansys; bar3(test) |
相關文章
- 2021-09-08BIM技術叢書Revit軟件應用系列Autodesk Revit族詳解 [
- 2021-09-08全國專業技術人員計算機應用能力考試用書 AutoCAD2004
- 2021-09-08EXCEL在工作中的應用 制表、數據處理及宏應用PDF下載
- 2021-08-30從零開始AutoCAD 2014中文版機械制圖基礎培訓教程 [李
- 2021-08-30從零開始AutoCAD 2014中文版建筑制圖基礎培訓教程 [朱
- 2021-08-30電氣CAD實例教程AutoCAD 2010中文版 [左昉 等編著] 20
- 2021-08-30電影風暴2：Maya影像實拍與三維合成攻略PDF下載
- 2021-08-30高等院校藝術設計案例教程中文版AutoCAD 建筑設計案例
- 2021-08-29環境藝術制圖AutoCAD [徐幼光 編著] 2013年PDF下載
- 2021-08-29機械AutoCAD 項目教程 第3版 [繆希偉 主編] 2012年PDF