㈠ 動態裂紋擴展模擬
一、數值流形方法的動力學求解格式
前面所討論的裂紋擴展問題,都只限於靜態載入或緩慢載入下的斷裂問題和裂紋的穩態擴展問題,它們都屬於斷裂靜力學的范疇。但是在實際工程中,會經常遇到含裂紋的結構在高載入速率下的斷裂問題或裂紋在失穩後的快速擴展問題,而這兩類問題都屬於斷裂動力學問題。斷裂動力學是研究慣性效應不能忽略的斷裂力學問題(范天佑,1990),斷裂動力學問題的求解方法明顯不同於斷裂靜力學問題,在動態載入時,試件除產生彈塑性變形外,內部各質點的自由振動獲得一定的加速度,從而產生慣性力,即所謂的動態載入時的慣性效應(Chen E P et al.,1976)。
數值流形方法在處理動力學問題時,在當前步的計算中,各單元繼承了前一時間步的速度,而不是像處理靜力學問題一樣置當前時間步的速度為零。這種處理方法考慮了單元的慣性效應,但沒有考慮單元的阻尼,對一般動力學問題精度不夠,本節在數值流形方法的動力學求解格式分析的基礎上,進行改進,以更好地模擬一般動力學問題。
在數值流形方法的求解公式中,動力問題與靜力問題的主要區別在於質量矩陣,在動力學問題中,這是一個很重要的矩陣,當時間步比較小時,慣性力矩陣將控制著整個材料體所有各點的運動和穩定。
考慮當前時間步時,(u(x,y,t)v(x,y,t))T表示為單元e的任一點(x,y)與時間相關的位移,假設M表示單位面積的質量,e是q個覆蓋(Ue(1),Ue(2),…,Ue(q))的交集。單位面積上的慣性力為
岩石斷裂與損傷
式中:
[T(e)(x,y)]=(Te(1) Te(2) Te(3)…Te(q))
{D(e)(t)}={De(1)(t)De(2)(t)De(3)(t)…De(q)(t)}T
在單元e中,由慣性力引起的勢能為
岩石斷裂與損傷
設{D(e)}={0}是在時間步起始時的單元位移,{D(e)(Δ)}={D(e)}是在時間步終了時的位移,Δ是時間步長。則
岩石斷裂與損傷
所以
岩石斷裂與損傷
式中:,是時間步開始時的單元速度。該時間步終了時的速度{V(Δ)}是
岩石斷裂與損傷
則勢能變為
岩石斷裂與損傷
上式可分解為兩部分,第一部分對應於剛度矩陣,為
岩石斷裂與損傷
第二部分對應於荷載向量,為
岩石斷裂與損傷
其中:r,s=1,…,12。
由上可知,數值流形方法對動力學問題的求解是顯式演算法,即利用當前時間步的位移和速度求下一時步的位移和速度,在每一時間開始時,繼承了前一時步的速度和位移,且位移是用泰勒級數展開的,其截斷誤差是三階的,當時間步取得較大時,誤差也較大。另外該演算法也是有條件的,在計算過程中時間步長的選擇必須要小於某一個值,即(王勖成等,1995):,其中:Tn為求解系統中最小尺寸單元的最小固有振動周期,因此求解系統中最小單元的尺寸將決定其時間步長的選擇。
這種演算法比較適合於波動問題的求解,一方面是因為這種方法的求解特點正好和波的傳播特點相一致;另一方面,研究波的傳播過程需要採用小的時間步長,也符合該演算法的要求。
一般的結構動力學問題結構的動力相應中低頻成分佔主導,從計算精度上考慮,允許採用較大的時間步長,通常採用Newmark無條件穩定隱式演算法,它是借用動力有限元法的求解思想,採用動力有限元中的Newmark方法來求解動力學問題。
在有限元方法中,對動力學問題的模擬採用的方程為結構動力學方程(張國新等,2002;張子明等,2001):
岩石斷裂與損傷
其中:M為質量矩陣;C為阻尼矩陣;d為位移增量;,分別是位移速度和加速度。K=Ke+Kcn+Kcs+Kf,Ke是剛度矩陣,Kcn、Kcs分別為塊體及不連續面之間的接觸矩陣,Kf是約束矩陣;F為總載荷向量,F=Fp+Fb+Ff-F0+Fcn+Fcs+Ffr,Fp是外載荷向量,Fb是體積力向量,Ff是已知約束位移引起的等效載荷向量,F0 是初應力向量,Fcn、Fcs分別為法向和切向接觸引起的等效載荷向量,Ffr為接觸面之間的摩擦力引起的等效載荷向量。
對一般的結構動力學問題,動力方程式(12-26)通常採用Newmark解法。Newmark積分法實質上是線性加速度法的一種推廣,它採用如下假設(王勖成等,1995):
岩石斷裂與損傷
其中:α、β是按積分精度和穩定性要求而決定的兩個參數。當β≥0.5,α≥0.25(0.5+β)2時,Newmark法是一種無條件穩定演算法,而且求解時可以採用比一般差分法大得多的時間步長,它比較適合於求解結構動力學問題。在這種方法中t+Δt時刻的位移解答dt+Δt是通過滿足t+Δt時刻的運動方程:而得到的。
利用Newmark法求解運動方程的演算法步驟如下:
1.初始計算
(1)形成剛度矩陣K,質量矩陣M和阻尼矩陣C。
(2)給定d0、和。
(3)選擇時間步長Δt、參數α和β,並計算積分常數:
岩石斷裂與損傷
(4)形成有效剛度矩陣:
岩石斷裂與損傷
(5)三角分解:
岩石斷裂與損傷
2.對於每一時間步長
(1)計算時間t+Δt的有效載荷:
岩石斷裂與損傷
(2)求解時間t+Δt的位移:
岩石斷裂與損傷
(3)計算時間t+Δt的速度和加速度:
岩石斷裂與損傷
綜上可見,當選擇合適的參數後,Newmark法是一種無條件穩定的隱式演算法,即時間步長Δt的大小不影響解的穩定性,此時Δt的選擇主要根據解的精度要求而選擇。與數值流形方法中最初的動力學求解方法相比,這種演算法是以對求逆獲取較大的時間步長(彭自強,2003)。Newmark演算法考慮了動力求解問題中的阻尼效應。本章中動力學問題求解採用Newmark演算法。
二、動態裂紋計算模型
與靜態載入時一樣,當動態載荷作用於含裂紋的物體時,在裂紋尖端也將產生應力集中現象,這種裂紋尖端的應力集中現象用動態應力強度因子來描述,而動態應力強度因子不但是載荷、裂紋長度和帶裂紋體幾何尺寸的函數,而且是時間的函數,有文獻提出了動態應力強度因子K(t)與相應的靜態應力強度因子K(0)的關系式:
岩石斷裂與損傷
式中k(V)是裂紋擴展速度的函數,而與幾何形狀無關。其中k(V)可近似表達為
岩石斷裂與損傷
式中:CR為Rayleigh表面波速;h是一個彈性波速的函數,有
岩石斷裂與損傷
其中:c1、c2分別為材料的膨脹波速和剪切波速,有
岩石斷裂與損傷
其中:ρ、G、ν分別為材料的密度、剪切模量和泊松比。
根據以上的討論,可以由材料的性質求出裂紋的動態斷裂應力強度因子,並根據相應的准則來對動態裂紋的擴展進行判斷。動態裂紋擴展的判據很多,最簡單實用的是應力強度因子判據,對Ⅰ型裂紋:K(t)≥KID(V),KID(V)是純Ⅰ型裂紋的動態斷裂韌性,它是一個和裂紋擴展速度有關的函數,但其變化不是很大,因而在實際應用中,通常認為它是一個不變的材料常數(劉再華等,1996)。
對於平面Ⅰ、Ⅱ型混合動態裂紋問題來說,裂紋尖端的應力場和位移場的形式與式(12-3)完全一樣,所不同的就是把式(12-3)中的靜態應力強度因子換為動態應力強度因子即可(劉再華等,1996)。盡管形式一樣,但是由於動態應力強度因子是時間的函數,所以動態問題的應力場和位移場也都是時間的函數,因而與靜態問題有著本質的不同。下面就根據以上動態裂紋擴展的計算原理對動態裂紋擴展問題進行模擬。
三、算例分析
沖擊載荷是工程實踐中常見的一種動態載荷,也是引起大多數工程材料破壞的主要外載形式。當沖擊載荷作用於含裂紋等缺陷的結構物時,這些結構物的動態斷裂行為與靜態情況下有很大的差別。沖擊載荷作用於物體會產生應力波在物體中傳播,當應力波作用於裂紋時,裂紋尖端的應力強度因子將因材料結構和裂紋模式的變化而具有不同的動態響應,並具有不同的起裂和擴展行為。同時,裂紋的動態行為也對應力波的傳播具有不同的散射作用,因此存在著各種波與裂紋間的相互作用關系。所以對含裂紋缺陷的結構在沖擊載荷作用下動態行為的研究,不僅要研究材料的動態力學性能,而且要對工程結構的安全評估等進行研究。下面利用數值流形方法對沖擊載荷作用下裂紋的動態擴展行為進行模擬。
為了與相關試驗結果進行對比分析,採用姚學鋒等(1996)中的試驗模型進行模擬。計算模型為如圖12-7所示的含偏置裂紋的簡支梁,為了能夠更好地觀察到裂紋的曲裂現象,簡支梁的寬度取得較大,具體尺寸如圖所示,所受的載荷為一階躍的沖擊載荷,其大小為107N。梁的彈性常數為:彈性模量E=100GPa,泊松比ν=0.2。對於不同的偏離距離a進行模擬,模擬結果如圖12-8所示。
圖12-7 含偏置裂紋的簡支梁(單位:cm)
由圖12-8的模擬結果可以看出,在不同偏置距離的情況下,裂紋的擴展路徑表現出不同的曲裂現象,其現象和原因分析如下:
(1)裂紋開始起裂時擴展方向沿著初始裂紋方向,擴展達到一定長度後,開始偏離原來的方向發生曲裂,其方向是朝著沖擊載荷作用點的方向。
(2)初始裂紋的偏置距離越大,曲裂發生的時間越早,沿初始裂紋方向擴展的長度也越小,曲裂現象越嚴重。
(3)裂紋向上擴展到一定長度後,又突然偏離了原來的方向而發生第二次曲裂,最終與上邊界貫通。
圖12-8 不同偏置距離時的裂紋擴展路徑
產生曲裂現象的主要原因是:當沖擊載荷產生的彎曲波作用於裂紋時,裂紋尖端作用的力是Ⅰ型和Ⅱ型復合載荷,裂紋的擴展屬於Ⅰ-Ⅱ型復合裂紋。偏置距離越小,Ⅰ型應力分量所佔有的比重越大,所以剛開始起裂時,裂紋沿著原來的方向擴展。隨著偏置距離的增加,裂紋的曲裂程度也有所增加。當裂紋擴展到梁的上半部時,裂紋在慣性作用下繼續作曲裂運動,但受到壓應力的作用,使得裂紋的擴展速度有所減慢。由於裂紋的擴展導致梁內的應力重新分布,在裂紋尖端的下部出現了一個新的中心點,並在該點的上部和下部分別形成了相應的拉應力分量和壓應力分量,這些分量構成了新的彎矩使裂紋產生了第二次曲裂,並最終導致裂紋到達上邊界而產生結構的整體破壞。
圖12-9為不同偏置距離情況下的裂紋尖端兩類應力強度因子隨時間的變化關系。從圖中可以看到:
(1)KⅠ均大於零,而KⅡ有正有負。正KⅡ和負KⅡ的意義體現裂紋表面相對位移的方向相反。
(2)應力強度因子在開始的時候比較小,隨著裂紋擴展長度的增加,裂紋尖端的應力強度因子也迅速增加,最後導致裂紋貫穿,在即將貫穿的最後時刻,應力強度因子出現了波動。但在整個過程中,KⅡ始終是比較小的,這說明在曲裂過程中Ⅰ型裂紋擴展佔主導作用,而Ⅱ型裂紋僅佔次要地位。
圖12-9 不同偏置距離時的應力強度因子隨時間的變化曲線圖(時間單位:s;應力強度因子單位:104N/m3/2)
(3)從裂紋擴展時間上來看,當偏置距離較小時,裂紋擴展到梁頂部即貫通的時間較短,而隨著偏置距離的增加,裂紋擴展的路徑也隨著加長,裂紋的貫通時間也增大。
(4)從裂紋擴展過程中應力強度因子的大小來看,在同一時刻,隨著裂紋偏置距離的增加,裂紋尖端的應力強度因子是逐漸減小的,與理論分析結果一致,因為隨著偏置距離的增加,在同樣的外力作用下,裂紋受力減小,因而裂紋尖端的應力強度因子也相應地有所減小。
㈡ 哪些模擬演算法是穩定的
標 題: 電力系統電磁暫態分析與機電暫態分析的功能區別和模型特點
發信站: BBS 水木清華站 (Thu Jul 22 16:46:31 2004), 站內
電磁暫態過程數字模擬是用數值計算方法對電力系統中從數微秒至數秒之間的電磁暫
態過程進行模擬模擬。電磁暫態過程模擬必須考慮輸電線路分布參數特性和參數的頻
率特性、發電機的電磁和機電暫態過程以及一系列元件(避雷器、變壓器、電抗器等
)的非線性特性。因此,電磁暫態模擬的數學模型必須建立這些元件和系統的代數或
微分、偏微分方程。一般採用的數值積分方法為隱式積分法。
由於電磁暫態模擬不僅要求對電力系統的動態元件採用詳細的非線性模型,還要計及
網路的暫態過程,也需採用微分方程描述,使得電磁暫態模擬程序的模擬規模受到了
限制。一般進行電磁暫態模擬時,都要對電力系統進行等值化簡。
電磁暫態模擬程序目前普遍採用的是電磁暫態程序(electromagnetic transients program
,簡稱為EMTP),1987年以來,EMTP的版本更新工作在多國合作的基礎上繼續發展,
中國電力科學研究院(簡稱電科院)在EMTP的基礎上開發了EMTPE。具有與EMTP相似功
能的程序還有加拿大Manitoba直流研究中心的EMTDC/PSCAD、加拿大哥倫比亞大學的
MicroTran、德國西門子的NETOMAC等。
機電暫態過程的模擬,主要研究電力系統受到大擾動後的暫態穩定和受到小擾動後的
靜態穩定性能。其中暫態穩定分析是研究電力系統受到諸如短路故障,切除線路、發
電機、負荷,發電機失去勵磁或者沖擊性負荷等大擾動作用下,電力系統的動態行為
和保持同步穩定運行的能力。
電力系統機電暫態模擬的演算法是聯立求解電力系統微分方程組和代數方程組,以獲得
物理量的時域解。微分方程組的求解方法主要有隱式梯形積分法、改進尤拉法、龍格
-庫塔法等,其中隱式梯形積分法由於數值穩定性好而得到越來越多的應用。代數方
程組的求解方法主要採用適於求解非線性代數方程組的牛頓法。按照微分方程和代數
方程的求解順序可分為交替解法和聯立解法。
目前,國內常用的機電暫態模擬程序是電力系統綜合程序(PSASP)和中國版BPA電力系
統分析程序。國際上常用的有美國PTI公司的PSS/E,美國EPRI的ETMSP,以及國際電
氣產業公司開發的程序如:ABB的SIMPOW程序、德國西門子的NETOMAC也有機電暫態仿
真功能。
㈢ gproms軟體里有什麼優化演算法,可以求解穩態MINLP問題嗎
gproms的主要功能是流程工業工藝單元、過程設計、模擬和優化。本人判斷嵌入的優化演算法應該不能解MINLP,既是能解,效果也不能保證。建議用專業gams,或者lingo等。詳細的可以跟gproms的代理商聯系詢問。
㈣ 穩態與非穩態有什麼區別
穩態與非穩態求解的區別
1,流體在管道中流動時,若與流動有關的參數(例如速度)隨時間而變化,這時流體的流動稱為非穩態流動。
截面上的速度雖不相同,但所有各截面上的速度均不隨時間而變,這時,速度僅隨空間位置的改變而變化,而與時間無關,這種流動情況稱之為穩態流動。
2,兩種求解的結果對比而言:
一種情況是要求解的物理問題是steady(穩態)的,從理論上來說那麼兩種求解方式收斂之後的結果都是一樣或者近似的。只不過需要注意的是,在 unsteady(非穩態)的求解中dt的選取會影響計算結果,有可能會計算發散,而且還必須要足夠的計算步達到收斂才能和steady的結果進行比較
換句話 說,如果用unsteady(非穩態)的方法去求解steady(穩態)流,如果計算本身就沒收斂,就取結果進行比較,那麼肯定是不行的。一般而言,在進行 unsteady(非穩態)求解的時候,前面一段時間的計算結果基本上是不予採用的,因為有一個數值收斂的過程。
3,是你要求解的物理問題是unsteady(非穩態)的,那麼用steady(穩態)的求解方法得出的結果就是一堆垃圾了,沒有任何價值。
(4)符合求解穩態演算法擴展閱讀:
非穩態的特點:
1 熱導體內的溫度場分布隨時間的變化而變化
2在垂直於熱流量的方向上,每一截面的熱流量都不相等。
3非穩態導熱可以分為周期性和非周期性兩種類型。
穩態隨著控制論和其他生命科學的發展,穩態已不僅指內環境的穩定狀態,也擴展到有機體內極多的保持協調、穩定的生理過程,例如生命活動功能以及正常姿勢(直立以及行路姿勢)的維持等;也用於機體的不同層次或水平(細胞、組織器官、系統、整體、社會群體)的穩定狀態。
以及在特定時間內(由幾毫秒直至若干萬年)保持的特定狀態。穩態不僅是生理學,也是當今生命科學的一大基本概念。它對控制論、遺傳學(基因的穩態調節)、心理學(情緒穩態等)、病理學、臨床醫學等多種學科都有重要意義。
㈤ 萬百五的科學成就
中國古代自動裝置及其原理分析的研究
1963年起萬百五在西安交通大學學報上以「我國古代在自動調整系統方面的成就」為題發表兩篇論文,介紹初步收集到的不少珍貴材料及其成就探討。1965年他匯合成論文「我國古代自動裝置的原理分析及其成就的探討」,發表於自動化學報。論文認為我國古代勞動人民在自動裝置上的創造和成就是極其突出的。1991年他又補充新材料為《中國大網路全書:自動控制與系統工程卷》寫成新條目「我國古代自動裝置」。文中例舉:指南車是採用擾動補償原理的方向開環自動調整系統;銅壺滴漏計時裝置是採用非線性限制器的多級阻容濾波;浮子式閥門是用於銅壺滴漏計時裝置中保持水位恆定的閉環自動調節系統,又用於飲酒速度自動調節器;記里鼓車是備有路程自動測量裝置的車;漏水轉渾天儀是天文表現儀器,採用模擬原理的水運渾象;候風地動儀是觀測地震用的自動檢測儀器;水運儀象台採用模擬原理演示或觀測天象的水力天文裝置,內有樞輪轉速恆定系統採用內部負反饋並進行自振盪的系統。
大工業過程遞階穩態模型辯識
1991年萬百五提出利用工業過程在隨機雜訊下的動態信息辨識過程的穩態模型。
他和課題組以設定點的正常階躍變動作為輸入激勵信號,取得過程的輸出數據,進行對線性近似動態模型的開環辨識,採用最小二乘估計模型的參數。最後由線性近似動態模型建立穩態模型。重要的是,他們證明:階躍激勵信號對於穩態模型是充分激勵的。課題組將該方法推廣到多輸入多輸出工業過程以及一類非線性工業過程。對於雙線性、Hammerstein型和Wiener型非線性過程,採用穩態辨識方法和動態辨識方法結合的集成辨識技術,先求整體穩態模型結構,然後再用數次設定點變動所得數據進行模型參數估計,所得的整體穩態模型以及過程導數都是強一致性的。
非凸靜態大系統的凸化
靜態大系統在分解後採用目標協調法進行協調,這就要求它的目標函數具有嚴格凸性、存在鞍點,並Lagrange函數具有可分性。對於大工業過程常常採用一些經濟指標(如能耗、成本、利潤等)作為目標函數,它是非凸的。1986年萬百五和課題組採用二次等價性原理及點凸化技術法、梯度偏移,來解決線性系統具有線性約束及目標函數具有非嚴格凸性的凸化問題,以及輸出預估和平衡凸化法來解決非線性系統並具有不等式約束問題。但對一般可分問題具有等式和不等式約束,他們採用乘子及部分罰函數凸化方法。上述諸方法求得的非凸靜態大系統最優解都是局部最優解。1998年萬百五和課題組把此問題嵌入到一個多目標優化問題中,然後通過p次冪凸化非劣前沿,再從非劣解中挑出原問題的全局最優解。
基於模型未知的大工業過程穩態優化演算法
對於模型未知的(大)工業過程,英國P.Roberts教授提出基於系統優化與參數估計集成研究法(ISOPE)的修正兩步法。萬百五和課題組採用ISOPE演算法與分解—協調相結合的自適應法,可以解決工業過程含有與輸出相關不等式約束和非凸目標函數問題。課題組提出的雙模型法也能解決含有與輸出相關不等式約束。他和課題組提出的利用變數增廣技術的凸化目標函數法,可以應用於上述諸方法,從而增加它們的在線迭代效率和便於使用者選擇增益因子。後他們提出的三環迭代法,一般能進一步減少設定點變動次數,當工業過程呈線性或以線性為主時,效果更為明顯。他和課題組將序列模型逼法與分解—協調結合推廣應用到大工業過程並形成單迭代、雙迭代兩種演算法,為大工業過程建模與穩態優化開辟新途徑。
工業過程的廣義穩態優化
工業過程的穩態優化過去都指常值穩態的優化控制。1998 年萬百五和課題組提出穩態的概念除常值穩態外,應予以擴充即包括隨機穩態、周期穩態、擬周期穩態和混沌穩態。1989年他們對大工業系統建立一套隨機穩態優化理論與方法。利用動態信息辨識目標函數及其導數的估價值的基礎上,將隨機穩態優化問題等價為一個確定性優化問題,子系統模型則採用遞階穩態辨識方法。辨識和優化可以集成在一起形成單迭代及雙迭代演算法。他們將二次方差分析引入到上述問題中使原最優解獲得改善。對廣義穩態(特別是混沌穩態),他們利用遍歷理論的Borel不變測度積分定義廣義穩態集合的中心名義值,並據此對廣義穩態進行描述。並給出等效優化問題存在條件及其求解演算法步驟。對混沌穩態的優化控制,他和課題組研究和提出等效優化問題及其求解演算法步驟。
大工業過程的智能控制
1993年萬百五將它定義為應用人工智慧的概念和方法,來解決大系統的辨識、
鎮定、優化、協調、決策以及故障檢測、診斷和定位、處理等問題,其中特別是總結和模擬人在控制和決策進程中的經驗和規律。具體而言就是運用神經網路、迭代學習、模糊控制以及由知識庫、資料庫、學習機、推理機所組成的智能決策單元(專家系統)等來解決大系統的上述問題。他和合作者在大工業過程智能穩態優化方面作了開創性工作。諸如利用(大)工業過程動態信息基於三層或四層前向神經網路建立(大)工業過程穩態模型,可以離線訓練、在線更新,優化時可仍使用像SUMT那樣的優化演算法。基於規則的智能推理決策機構的利用,將大系統理論和知識工程相結合,可以解決問題智能協調、智能穩態優化等問題,簡便而實用。具有模糊參數的大工業過程智能穩態優化問題,也是採用基於規則的智能推理決策機構求解。迭代學習控制的引用,使得大工業過程的多次迭代而引起的設定點變動,減輕對工業過程的擾動,避免工業過程的振盪或失穩,縮短動態時間。基於Hopfield神經網路的大系統各種遞階優化智能智能演算法,具有非迭代性和快速性。特別是將大系統穩態優化問題通過Hopfield神經網路轉化為具有非迭代性和快速性的網路求解一組微分方程問題,對於存在模型與實際差異採用閉環控制的穩態優化問題,一次迭代就可求得最優解。
㈥ 求解穩態電壓調整問題
課本例題沒有考慮功率損耗,直接把末端功率當始端功率用,但就功率上不差啥,但是線路阻抗在那兒放著呢,算功率損耗和不算功率損耗能差個幾千伏,楊淑英那本書上給了兩種解法:一,設末端電壓(歸算至高壓側)為高壓側額定電壓,從後往前推,求出始端功率,用該始端功率和已知始端電壓求電壓降;二,就是您說的用一元二次方程解。這兩種演算法之間也有誤差,真的很糾結該用著三種方法里的哪一種。鑒於整道題的重點考察點不是這,為了考試里節省時間,用課本上那種不計損耗計算最方便,真這樣寫上去考試時候不會扣分吧?就怕圖方便最後扣分就哭了 查看原帖>>
㈦ 為什麼冬季可以採用穩態演算法計算採暖負荷而夏天卻一定要採用動態演算法計算空調負荷
穩態演算法即Q=KF(T1-T2),K和F就不去講,T1-T2即室內外的設計時的溫度,這種演算法簡單,你都可以用這個公式估算,但T1和T2是兩個不定的值,如果波動小則誤差倒小,在我國北方的冬季,室外溫度的波動幅度遠小於室內外的溫差,因此在做採暖負荷計算時,採用的是日平均溫差的穩態計演算法。
若南方的夏季也採用穩態演算法,則算出來的負荷就會小很多,所以只能用不穩定傳熱演算法,通過計算逐時溫度來確定室外的干球溫度.
㈧ 在數值計算中實現演算法的數值穩定性的若干方法
拿華騰的交通卡系統做例子吧
先說數據結構:
程序裡面有很多用戶信息,這些信息是用結構體存放,是用鏈表的形式,還是用哈希表等等的方法進行排列,這個就是數據結構
演算法:
要從這些數據結構里找出一個結點,是用遍歷搜索呢還是用折半搜索之類的方法,叫做演算法。
LZ可以體會下。
㈨ 電力系統穩態,式子演算法
事先聲明,我不是專家,就是知道點。很復雜,所以只能簡單一說,對不住啊。
穩態和暫態,是由於電路中有電磁振盪。首先一定是交流電路,暫態分析是指電路接通瞬間,由於電壓跳變,會使電路原本狀態發生改變,由0-時的狀態變為0+時的狀態,在這瞬間的改變是暫態分析的范圍。在0+之後,電路中的電壓電流會隨著時間變化,其變化方式由電路元件決定,當滿足一定條件時,時間足夠長後,電路中的電壓電流會趨近於穩定,這時的分析就是穩態分析了。
㈩ 穩態誤差的計算方法有哪兩種
1)靜態誤差系數法(s趨於0)
Kp=limG(s) R(t)=V1(t) ess=V/1+Kp
Kv=limsG(s) R(t)=Vt ess=V/Kv
Ka=lims^2G(s) R(t)=V(1/2)t^2 ess=V/Ka
2)終值定理
ess=limsE(s)(s趨於0)
=limsΦe(s)R(s)
注意:無論哪一種做法,你都要先用勞斯判據判斷一下穩定性再求解穩態誤差