计算物理 第9章 有限单元法（附录）.docVIP

PAGE PAGE 3 第9章 附录 9.2.2 输入数据 %demo_fem_2.m function demo_fem_2 clc;clear; nel=15; nnel=2; ndof=1; nnode=nel+1; sdof=nnode*ndof; gcoord(1:nnode) = 0.0:1.0/nel:1.0; nodes(1:nel,1)=1:nel; nodes(1:nel,2)=2:nel+1; acoef=1; bcoef=-3; ccoef=2; % ?￠·?·?3ì?μêy bcdof(1)=1;bcval(1)=0; % ±??? ff=zeros(sdof,1); kk=zeros(sdof,sdof); index=zeros(nnel*ndof,1); for iel=1:nel; nl=nodes(iel,1);nr=nodes(iel,2); xl=gcoord(nl);xr=gcoord(nr); eleng=xr-xl; index=feeldof1(iel,nnel,ndof); k=feode2l(acoef,bcoef,ccoef,eleng); f=fef1l(xl,xr); [kk,ff]=feasmbl2(kk,ff,k,f,index); end