偏微分方程的有限元法求解_第1頁
偏微分方程的有限元法求解_第2頁
偏微分方程的有限元法求解_第3頁
偏微分方程的有限元法求解_第4頁
偏微分方程的有限元法求解_第5頁
已閱讀5頁，還剩17頁未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、 %對于d2u/dx2=f的FEM解算器，其中f=x*(1-x)%邊界條件u(0)=0, u(1)=0.%精確解用以比對xx=linspace(0,1,101);%產(chǎn)生0-1之間的均分指令，101為元素個數(shù)uex=(1/6).*xx.3-(1/12).*xx.4-(1/12).*xx;%對力項設(shè)置高斯點的數(shù)目NGf=2;if (NGf=2)xiGf=-1/sqrt(3);1/sqrt(3);%1、2的值aGf=1 1;else,NGf=1;xiGf=0.0;aGf=2.0;end%單元數(shù)目Ne=5;%建立網(wǎng)格節(jié)點x=linspace(0,1,Ne+1);%零剛性矩陣K=zeros(Ne+1,N

2、e+1);b=zeros(Ne+1,1);%對所有單元循環(huán)計算剛性和殘差for ii=1:Ne,kn1=ii;kn2=ii+1;x1=x(kn1);x2=x(kn2);dx=x2-x1;%每一個單元的長度dxidx=2/dx;%d/dxdxdxi=1/dxidx;%dx/ddN1dxi=-1/2;%d1/ddN2dxi=1/2;%d2/ddN1dx=dN1dxi*dxidx;%-1/(xj-xj-1)dN2dx=dN2dxi*dxidx;%1/(xj-xj-1)K(kn1,kn1)=K(kn1,kn1)-2*dN1dx*dN1dx*dxdxi;%Rj的第二項K(kn1,kn2)=K(kn1,k

3、n2)-2*dN1dx*dN2dx*dxdxi;K(kn2,kn1)=K(kn2,kn1)-2*dN2dx*dN1dx*dxdxi;K(kn2,kn2)=K(kn2,kn2)-2*dN2dx*dN2dx*dxdxi;%用高斯積分估計力項的積分for nn=1:NGf%NGf=2xiG=xiGf(nn);%得到高斯點的N1=0.5*(1-xiG);%求N1和N2(即在xiG的權(quán)重/插值) 形狀函數(shù)在的值N2=0.5*(1+xiG);%值fG=xiG*(1-xiG);%對點求fgG1=N1*fG*dxdxi;%在節(jié)點處估計權(quán)函數(shù)在高斯點的被積函數(shù)gG2=N2*fG*dxdxi;%估計是個積分值b(kn1)=b(kn1)+aGf(nn)*gG1;% aGf為1b(kn2)=b(kn1)+aGf(nn)*gG2;endend%在x=0處設(shè)置Dirichlet條件kn1=1;K(kn1,:)=zeros(size(1,Ne+1);K(kn1,kn1)=1;b(kn1)=0;%在x=1處設(shè)置Dirichlet條件kn1=1;K(kn1,:)=zeros(size(1,Ne+1);K(kn1,kn1)=1;b(kn1)=0;%求解方程v=Kb;%v為Kx=b的解plot(x,v,'*-');%畫圖并比較hold on;p

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論