龍格庫塔程序源代碼，用于解微分方程或微分方程組，給初學者作為參考。
上傳時間： 2013-12-25
上傳用戶：gyq
Delaunay三角剖分是1934年發明的將空間點連接為三角形，使得所有三角形中最小角最大的一個技術。
標簽： 剖分算法
上傳時間： 2015-07-13
上傳用戶：xxcl1980
每種程序設計語言都有描述程序語法結構的規則。例如，Pascal程序由程序塊（又叫分程序）構成，程序塊由語句組成，語句由表達式組成，表達式由記號組成等等。這些規則可以用上下文無關文法或BNF范式（Backus-Naur Form）描述。 編譯器常用的文法分析方法有自上而下和自下而上兩種。正如它們的名字所示，自上而下分析器建立分析樹是從根結點到葉結點，而自下而上分析器恰好反過來。它們的共同點是從左向右地掃描輸入，每次一個符號。 最有效的自上而下和自下而上的分析法都只能處理上下文無關文法的子類。這些子類足以描述程序設計語言的大多數語法結構，其中LL文法的分析器通常用手工實現，而LR文法的分析器通常利用自動工具構造。 本章致力于編譯器采用的典型語法分析方法。我們首先提出有關上下文無關文法的基本概念，然后介紹適合于手工實現的預測分析技術，最后給出自動工具用的LR分析算法。由于程序員準備的代碼經常會出現一些語法錯誤，因此我們還擴展所介紹的分析方法，使之能從常見的錯誤中恢復過來。 3.1 上下文無關文法 ..... 3.4 自下而上分析 3.6 二義文法的應用
上傳時間： 2013-12-25
上傳用戶：hanli8870
一階有限元通用程序，邊界條件可變,可以用于機械應力分析，電磁場數值計算等領域。而且可以自編自東剖份程序作為入口。
上傳時間： 2015-10-18
上傳用戶：chenxichenyue
這是我學習并行編程的第一個程序， 原來時fortran的， 我將其改寫為c的， 兩分程序都在包中， 程序比較簡單，但很有助于理解mpi
上傳時間： 2015-12-11
上傳用戶：xuanjie
經常看到很多朋友問定步長的龍格庫塔法設置問題，下面吧定步長四階龍格庫塔程序貼出來，有需要的可以看看
標簽：
上傳時間： 2016-07-26
上傳用戶：rocwangdp
這是FreeFem2D的姊妹軟件，對于學習有限元程序設計的人來說有很好的參考價值，我在Fedora Core9下編譯運行通過(安裝很容易)，編譯后會生成ff3d可執行程序，通過讀取描述問題的文件就可以自動剖分／求解等。
上傳時間： 2013-12-17
上傳用戶：thinode
function [alpha,N,U]=youxianchafen2(r1,r2,up,under,num,deta) %[alpha,N,U]=youxianchafen2(a,r1,r2,up,under,num,deta) %該函數用有限差分法求解有兩種介質的正方形區域的二維拉普拉斯方程的數值解 %函數返回迭代因子、迭代次數以及迭代完成后所求區域內網格節點處的值 %a為正方形求解區域的邊長 %r1,r2分別表示兩種介質的電導率 %up,under分別為上下邊界值 %num表示將區域每邊的網格剖分個數 %deta為迭代過程中所允許的相對誤差限 n=num+1; %每邊節點數 U(n,n)=0; %節點處數值矩陣 N=0; %迭代次數初值 alpha=2/(1+sin(pi/num));%超松弛迭代因子 k=r1/r2; %兩介質電導率之比 U(1,1:n)=up; %求解區域上邊界第一類邊界條件 U(n,1:n)=under; %求解區域下邊界第一類邊界條件 U(2:num,1)=0;U(2:num,n)=0; for i=2:num U(i,2:num)=up-(up-under)/num*(i-1);%采用線性賦值對上下邊界之間的節點賦迭代初值 end G=1; while G>0 %迭代條件：不滿足相對誤差限要求的節點數目G不為零 Un=U; %完成第n次迭代后所有節點處的值 G=0; %每完成一次迭代將不滿足相對誤差限要求的節點數目歸零 for j=1:n for i=2:num U1=U(i,j); %第n次迭代時網格節點處的值 if j==1 %第n+1次迭代左邊界第二類邊界條件 U(i,j)=1/4*(2*U(i,j+1)+U(i-1,j)+U(i+1,j)); end if (j>1)&&(j U2=1/4*(U(i,j+1)+ U(i-1,j)+ U(i,j-1)+ U(i+1,j)); U(i,j)=U1+alpha*(U2-U1); %引入超松弛迭代因子后的網格節點處的值 end if i==n+1-j %第n+1次迭代兩介質分界面（與網格對角線重合）第二類邊界條件 U(i,j)=1/4*(2/(1+k)*(U(i,j+1)+U(i+1,j))+2*k/(1+k)*(U(i-1,j)+U(i,j-1))); end if j==n %第n+1次迭代右邊界第二類邊界條件 U(i,n)=1/4*(2*U(i,j-1)+U(i-1,j)+U(i+1,j)); end end end N=N+1 %顯示迭代次數 Un1=U; %完成第n+1次迭代后所有節點處的值 err=abs((Un1-Un)./Un1);%第n+1次迭代與第n次迭代所有節點值的相對誤差 err(1,1:n)=0; %上邊界節點相對誤差置零 err(n,1:n)=0; %下邊界節點相對誤差置零 G=sum(sum(err>deta))%顯示每次迭代后不滿足相對誤差限要求的節點數目G end
標簽： 有限差分
上傳時間： 2018-07-13
上傳用戶：Kemin
目前,在電壓互感器設計中,雖有人進行過可靠性設計利優化設計方面的研究,但采用的方法仍為傳統方法.本文采用現代設計方法,它將有限元分析、可靠性設計技術利優化設計技術有機的結合起來,因此采用現代設計方法得到的方案比利用傳統設計方法設計出的方案更加經濟合理.首先,本文簡單介紹了電壓互感器的原理,描述了電壓互感器的分類、基本參數和誤差分析.第二,本文研究了電磁場有限元分析原理,介紹了麥克斯韋方程和電磁場微分方程.本文采用大型通用有限元分析軟件ANSYS對電壓互感器進行二維電磁場有限元分析,對電壓互感器建立了有限元數學模型和網格剖分,對有限元模型加載了邊界條件并進行了求解.研究了二維磁場分析單元PLANE53單元利電路模擬單元CIRCU124單元的特點及使用方法.第三,對電壓互感器的瓷套部分進行了可靠性設計.瓷套所受的彎曲負荷應力很多,主要包括:風力負荷產生的彎曲應力,地震負荷產生的彎曲應力,產品運輸中傾斜產生的彎曲應力.本文研究了瓷套的應力分布的確定方法,將多種應力疊加在一起,推出了應力分布參數的計算公式.瓷套的應力、強度利各設計變量均可認為服從正態分布,在設計時作為正態分布變量處理.本文應用應力-強度干涉理論,對電壓互感器瓷套的可靠性設計方法進行了研究.第四,研究了ANSYS軟件的優化設計模塊,研究了采用ANSYS軟件進行優化設計的步驟和優化工具及方法.利用ANSYS軟件的參數化設計語言與其OPT模塊,實現了有限元數值計算與優化設計的有機結合.并以額定一次電壓35KV,額定二次電壓100V,額定頻率50HZ的電壓互感器為例,進行了有限元分析計算利優化設計.根據電壓互感器產品設計的實際情況,確定設計變量為繞組導線規格和鐵心結構尺寸.優化循環結束以后,可以選擇列出所有參數的數值,也可以只列出優化變量,可以用圖顯示指定的參數隨序列號的變化情況,通過多方案的比較,得到最優方案.將現代設計方法應用于生產廠家,可節省研究開支,大大縮短開發周期,減少計算誤差,減少試驗費用,降低成本,提高產品的可靠性,因此本項目的研究具有良好的經濟效益和社會效益.
上傳時間： 2013-06-10
上傳用戶：tuilp1a
該文介紹了一種新型高壓發電機電力發生器,它無需升壓變壓器即可直接連接到電網,其定子采用多層同心式繞組,槽內導體為高壓電纜,高壓電纜的引入克服了傳統發電機輸出電壓不能高于36kV的限制;并簡要介紹了這種發電機的全新設計與應用前景;最后針對電力發生器不同于傳統發電機的結構,借助有限元分析軟件進行了端部的建模、端部磁場、端部漏抗與端部電磁力的求解.文中圍繞一模型樣機,首先介紹了三維渦流場計算與利用磁場儲能進行參數計算的理論基礎.之后進行了對定子端部區域的建模,由于電力發生器采用多層同心式繞組,其端部結構較為復雜,這對模型的建立、剖分都帶來了相當大的難度.為了達到簡化分析計算的目的,我們對所求解的實際模型進行了簡化處理,并闡述了簡化的理論根據.在此基礎上,詳細介紹了如何利用有限元分析軟件ANSYS進行具體分析計算,包括網格剖分、電流加載及邊界條件的處理.最后得出了端部磁場矢量分布圖,端部漏抗值及端部繞組的電磁力分布規律.該文采用了簡化模型的方法進行計算,為了驗證簡化的合理性,我們進行了實例計算驗證.結果表明,文中所采用的簡化方法是合理的.該文所進行端部磁場、端部漏抗及端部電磁力計算,為進一步分析其他工況下電力發生器端部電磁力及振動提供了參考.
上傳時間： 2013-06-26
上傳用戶：zhanditian
