Matlab有限元结构动力学分析与工程应用（源文件）

ansys1986

Matlab有限元结构动力学分析与工程应用 徐斌著 源文件。
希望对大家有帮助

电脑克星LK

嗯，这里的大多程序运行都会遇到问题，需要修改一些地方才能顺利计算出结果。但这都是小问题，有了程序源文件，不用再自己编那么多行代码，还是非常感谢作者的辛勤付出。看懂这本书是需要对有限元和结构动力学有一些基本掌握，新手急着看懂它是不现实的，得好好学习了结构动力学和有限元之后才能渐渐明白这本书的精彩之处。总之，这是一本好书，在杆系结构的有动力分析方面，是非常好的学习资料，并且具有相当的工程应用价值。

wuluo1213

徐斌老师的这本书很不错的，我打印了一本，后来有买个一本。呵呵，还发邮件给老师要回来了程序。很不错的老师！

zhwaing

ok非常感谢

rocketma

就是程序运行不了啊

rocketma

2# wuluo1213
里面的程序能运行吗

ck2006

恩，确实是一本好书！

zhenjunwei1987

为什么我程序不能运行呢

zhenjunwei1987

我照着书里面的程序自己不程序写入也不能运行啊 这是为什么呢

zhenjunwei1987

现在写书的人太不责任了

only8810205028

那个程序是有点问题，请问大家有没有用EX543求振型啊，我的怎么求出来是错误的啊，有两条命令，我不是很懂
index=zeros(2*3,1);                         % 初始化指数矢量
%  (2) 计算约束
%--------------------------------------------------------------------------
ConNode=[ 1,  1,  1,  1;...                            % 约束的节点
         17,  1,  1,  1];
ConVal =[ 1,  0,  0,  0;...                              %约束自由度的值
         17,  0,  0,  0];
bcdof=zeros(51,1);                               % 初始化自由度矢量
bcval=zeros(51,1);                                % 初始化自由度矢量值
[n1,n2]=size(ConNode);
                                                
for ni=1:n1
  ki=ConNode(ni,1);
  bcdof((ki-1)*3+1:ki*3)=ConNode(ni,2:3+1);                                             
  bcval((ki-1)*3+1:ki*3)=ConVal(ni,2:3+1);                                          
end
%--------------------------------------------------------------------------
%  (3) 施加节点力
%--------------------------------------------------------------------------
P=[-500,10,2];                        % load applied at node 21 in the negative y direction施加载荷
ff=zeros(51,1);                            % 初始化力矢量
ff(3*(P(2)-1)+P(3))=P(1);
%  (4) 坐标变换
%--------------------------------------------------------------------------
% x, y, z 整体坐标系
xx=zeros(17,1); yy=zeros(17,1); zz=zeros(17,1);
dx=0.011;
xx=0:dx:(17-1)*dx; xx=xx';
gcoord=[xx  yy  zz];
for iel=1:16                           
  nd(1)=iel;                    
  nd(2)=iel+1;                    
  x(1)=gcoord(nd(1),1); y(1)=gcoord(nd(1),2); z(1)=gcoord(nd(1),3);
                                                      
  x(2)=gcoord(nd(2),1); y(2)=gcoord(nd(2),2); z(2)=gcoord(nd(2),3);
                                                      
  leng=sqrt((x(2)-x(1))^2+(y(2)-y(1))^2+(z(2)-z(1))^2);
end
K=sysK;
M=sysM;
%  (5) 施加约束
%--------------------------------------------------------------------------
% 检查总体刚度矩阵的主元是否为0
[kk0,mm0,ff0,bcdof0,bcval0,sdof0]=kkCheck1(K,M,ff,bcdof,bcval);
% 对结构系统施加边界条件
%[kk1,mm1,ff1,sdof1]=femApplybc1(K,M,ff0,bcdof0,bcval0);  
ni=length(bcdof);
sdof=size(K);
for ii=1:ni
   if bcdof(ii)==1
     for j=1:sdof
      K(ii,j)=0;
       K(j,ii)=0;
       M(ii,j)=0;
       M(j,ii)=0;
       ff(j)=ff(j)-bcval(ii)*K(j,ii);
     end
     K(ii,ii)=1;
     ff(ii)=bcval(ii);
   end
end
% 检查边界条件，消去相应的自由度
[kk2,mm2,ff2,sdof2]=bcCheck1(kk0,mm0,ff0,bcdof0,bcval0);
                     
[V,D]=eig(kk2,mm2);                        
[lambda,ki]=sort(diag(D));                      % 排序
omega=sqrt(lambda);                           
omega1=sqrt(lambda)/(2*pi);                           
V=V(:,ki);                                                   
%--------------------------------------------------------------------------
jk=3;                               %选择显示第几阶模态
  Vi=[V(:,jk)];                        
  ik=0;                     

  for ii=1:sdof0                                
    if bcdof0(ii)==0
      mcoord(ii,1)=Vi(ii-ik);
    else
      mcoord(ii,1)=0;
      ik=ik+1;
    end
  end

  for ii=1:17                           
    for ij=1:3
      mcoordA(ii,ij)=mcoord((ii-1)*3+ij,1);
    end
  end

  nv=20;
  for i=1:nv+1
    t=(i-1)*(2*pi)/20;
    mcoordB=gcoord+[zeros(17,1),mcoordA(:,1),zeros(17,1)]*cos(t);这条命令是什么意思啊
    for iel=1:16                          
      nd(1)=iel;                              
      nd(2)=iel+1;                           
      x(1)=mcoordB(nd(1),1); y(1)=mcoordB(nd(1),2); z(1)=mcoordB(nd(1),3);
                                                        
      x(2)=mcoordB(nd(2),1); y(2)=mcoordB(nd(2),2); z(2)=mcoordB(nd(2),3);
                                                      
        plot(x,y,'b');
        xlabel('x');
        ylabel('y');
        title('2-D frame  element')
      end
  end
%--------------------------------------------------------------------------
%  打印结果
%--------------------------------------------------------------------------
  disp('The calculation is use of:')
  disp('2-D frame  element')
  disp('and consistent mass matrix')
disp(' ')
disp('     mode   numrical   ')

num=1:1:10;                                          
frequency=[num' omega1(1:10)  ]

ghostzer

不错，早买了徐斌老师这本书，一直看不懂呀，感谢楼主，

许伍sky

wuluo1213 发表于 2011-3-18 16:57
徐斌老师的这本书很不错的，我打印了一本，后来有买个一本。呵呵，还发邮件给老师要回来了程序。很不错的老 ...

你好，我的matlab有限元结构动力学分析与工程应用，例子2.4.2运行不了，我想问一下你的怎么样，谢谢

yujie021639

会有程序运行不了？

lg80224

买本书学习一下！

guchengansimwe

不错，早买了徐斌老师这本书，一直看不懂呀，感谢楼主，

yuanmeng2015

先看看学习下

qinhaoxing

HAO DONGXI

yuxiaojian01

多谢了楼主！

yuansen

刚找到电子书，这个源程序收藏啦，谢谢！:)