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