驗(yàn)證ANSYS輸出的整體K和M

2017-01-16  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

%
% 目的:驗(yàn)證ANSYS中輸出整體剛度矩陣K和整體質(zhì)量矩陣M
% 以一個(gè)6單元簡支梁模型計(jì)算結(jié)果為例,相應(yīng)ANSYS代碼見例1所示。
%
% 驗(yàn)證基礎(chǔ):
% (1) 基于模態(tài)分析中,驗(yàn)證ANSYS振型是否關(guān)于質(zhì)量矩陣M歸一化
% (2) 特征值分析,由K,M計(jì)算特征值,即頻率和振型,與ANSYS計(jì)算結(jié)果對比
%
% 注意事項(xiàng):
% (1) ANSYS輸出矩陣以Harwell-Boeing格式存儲(chǔ),并且由于K,M的對稱性
% 僅存儲(chǔ)了對角線以下的元素,因此,在還原矩陣時(shí),還應(yīng)將對稱部分考慮進(jìn)去。
%
% (2) 關(guān)于如何輸出K,M的ANSYS命令流以及如何提取Harwell-Boeing格式矩陣,
% 可參考前期整理的資料。
%
% Author: J.H.Zhu
% 2009-06-09
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';
% 特征值分析,即計(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;
mnorm=sqrt(X1(:,j)'*M*X1(:,j));
X1(:,j)=X1(:,j)/mnorm;
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
name=['mod0',num2str(i),'.txt'];
Phi_ansys(:,i)=load(name);
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)





開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:驗(yàn)證ANSYS輸出的整體K和M Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn) 

編輯
在線報(bào)名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com