復(fù)習(xí)-ees2015課件第二講

單位系統(tǒng)EES是一個(gè)對(duì)單位要求嚴(yán)格的軟件單位系統(tǒng)是進(jìn)行設(shè)計(jì)時(shí)需要設(shè)定的首要因素在

Options菜單中設(shè)定用指令明確設(shè)定單位更合理

（將覆蓋對(duì)話框設(shè)置）$UnitSystemSIMASSDEGKPACKJ新建一個(gè)EES文檔保存為

EESLecture2.1–UnitSystem.ees添加

$TabStops

和

$UnitSystem

指令單位系統(tǒng)單獨(dú)的常數(shù)可以設(shè)定單位如下:m=25[kg]a=2.5[m/s^2]F=m*a不能給

方程設(shè)定單位,但EES

能自動(dòng)生成F的單位(在結(jié)果窗口中以紫色顯示)單位系統(tǒng)

在Optionsmenu>Preferences>Options選中

“Checkunitsautomatically”

和

“Setunitsautomaticallywhenpossible”,點(diǎn)擊

Store

然后選擇

EES.PRF文件應(yīng)用，也可以按照需要選擇其他選項(xiàng)單位系統(tǒng)變量在Options>VariableInfo

菜單中列出F

的單位是[kg-m/s^2]將其改為[N]

可以觀察到EES

沒(méi)有給出警告然后改為[kN]

再檢查有沒(méi)有警告!將m的單位改為[lbm]

同時(shí)將F的單位改為[N]

再檢查有沒(méi)有警告!單位系統(tǒng)預(yù)測(cè)值上下限紫色顯示自動(dòng)生成單位單位系統(tǒng)顯示單位改為

[N]沒(méi)有單位問(wèn)題改為

[kN]出現(xiàn)問(wèn)題單位系統(tǒng)EES允許單位轉(zhuǎn)換給出方程

F=ma,要求

F

的單位是

kN.如果在單位設(shè)置中將其設(shè)為[kN],會(huì)出現(xiàn)警告

因此可以用單位轉(zhuǎn)換函數(shù):

F=m*a*convert(N,kN) “Alternativelyyoucandothis,butthenyouneedtoknowtheconversionconstant”（或者也可以如下，但必須已知轉(zhuǎn)換常數(shù)） F_1=(m*a)

/1000[N/kN] F_2=m*a*0.001[kN/N]單位系統(tǒng)英制單位與標(biāo)準(zhǔn)單位的相互轉(zhuǎn)換:

m_3=10[lbm] a_3=3.5[m/s^2] {F_3=m_3*a_3}“Thiswouldgiveanerror”（錯(cuò)誤）

F_3=(m_3*convert(lbm,kg))*a_3

單位系統(tǒng)為了使某個(gè)常數(shù)的含義更加清晰，可以直接在后面設(shè)置單位，例如:"Thisisclearerthanthenext..."time=3.5[h]*3600[s/h]"Thefactthatthisis3.5hoursisnotasapparent!"time=12600[s]后者顯示不出是3.5hEES相關(guān)例子:Examples/Unitsconversion/Checkingunitsandunitconversion(HeatEx.EES)Examples/GettingStartedwithEES/Convertingunitsfunction(Convert.EES).單位系統(tǒng)EES的內(nèi)部功能有以下幾種數(shù)學(xué)功能流體物性計(jì)算固體/液體物性計(jì)算EES程序庫(kù)外部程序Example編碼可以被粘貼FunctionInfo(Help)新建一個(gè)EES文檔保存為 EESLecture2.2–BuildInFunction.ees內(nèi)部功能內(nèi)部功能數(shù)學(xué)算例x=cos(Value)

這是粘貼過(guò)來(lái)的函數(shù)現(xiàn)在可以隨意改變這個(gè)函數(shù)，例如：theta=30[deg]x_coordinate=cos(theta)或者z=cos(33)數(shù)學(xué)算例LogValue=log10(100)

Thelog10(Value)是直接粘貼的加下劃線防止出錯(cuò)gamma=gamma_(Value)T=140[C]E=E_(Aluminum,T)注意為英式英語(yǔ)積分方程EES能進(jìn)行積分和微分計(jì)算.例如下面的方程：EES函數(shù)：y=Integral(x^3,x,0,3,0.06)內(nèi)部功能物性算例對(duì)于物性計(jì)算，通常都需要確定某些狀態(tài)，例如壓力、溫度等還需要確定某種材料，一般為固體或者流體例如求空氣的密度，直接從FluidPropertiesFunctionInfo對(duì)話框中粘貼此函數(shù):rho_1=Density(Air,T=T_1,P=P_1)內(nèi)部功能rho_1=Density(Air,T=T_1,P=P_1)注意以下幾點(diǎn):變量的名稱是體統(tǒng)給定的，但是可以修改對(duì)于流體的名稱，它是一個(gè)帶有單引號(hào)的字符串，其中不可含有空格接著壓力和溫度狀態(tài)，可以寫(xiě)程序或者賦值內(nèi)部功能流體物性算例T和P在語(yǔ)句中,它們的單位?rho_1=

Density(Air,T=25,P=100)這樣編寫(xiě)程序更合理

T1=25[C]P1=100[kPa]rho_2=

Density(Air,T=T1,P=P1)"T1,P1是變量!"h=Enthalpy(Water,T=T1,P=P1)v_specific=Volume(Water,T=T1,P=P1)內(nèi)部功能固體物性算例求鋁的彈性模量T=140[C]E=E_(Aluminum,T)注意英美拼寫(xiě)差異內(nèi)部功能編程習(xí)慣編寫(xiě)出來(lái)的語(yǔ)句可讀性非常重要，平時(shí)應(yīng)該養(yǎng)成好的編程習(xí)慣。直接粘貼過(guò)來(lái)的代碼是:rho_1=Density(Air,T=T_1,P=P_1)改為下面的語(yǔ)句可更清晰rho_1=Density(Air,T=T_1,P=P_1)等號(hào)兩邊使用空格,逗號(hào)后面加上空格使分隔更為明顯.尤其當(dāng)方程比較復(fù)雜的時(shí)候要更加注意內(nèi)部功能N=5 Duplicatej=1,5 A[j]=h[j]-T*s[j] Endareequivalentto: A[1]=h[1]-T*s[1] A[2]=h[2]-T*s[2] A[3]=h[3]-T*s[3] A[4]=h[4]-T*s[4] A[5]=h[5]-T*s[5]運(yùn)行該程序?qū)⒔⑷缦聢D所示表格，可參考Exzamples/Duplicatecommand給出的三個(gè)例子Thelowerandupperlooplimitscanbeanyofthefollowing:i)anumericalquantityorexpressionii)anEESvariablethathasbeensettoavalueaheadoftheDuplicatestatement.iii)anEESvariablethatissettoavalueintheParametrictableiv)theTableRun#keyword.EEScompilestheequationsbeforeexecutionoccurs.ThatisthereasontheindexlimitsmustbeknownbeforetheDuplicatestatementisprocessed.Thisisoneoffewsituationsinwhichtheorderoftheequationsisimportant.NotethatEESdoesnotallowmorethanoneDuplicateclausetoappearwithintherangeofanouterDuplicateclause.Thefollowingcode,forexample,wouldresultinacompilationerror.duplicatei=1,5 duplicatej=1,6 x[i,j]=i*j end duplicatek=1,5 "thiswillnotwork" y[i,k]=i+k endendHowever,theintentoftheseequationscanbeenteredbyplacingeachinnerDuplicateclausewithinitsownotherclause,asshownbelow.Theseequationswillcompilewithnoproblems.duplicatei=1,5 duplicatej=1,6 x[i,j]=i*j endendduplicatei=1,5 duplicatek=1,5 y[i,k]=i+k endend正常運(yùn)行TheDuplicatecommandshouldgenerallyNOTbeusedinFunctionsandProceduressinceaRepeat-Untilclauseismoreefficientandprovidesmorecapability.Forexample,theequations Duplicatej=1,N A[j]:=h[j]-T*s[j] Endshouldbewrittenas j:=0 repeat j:=j+1; A[j]:=h[j]-T*s[j] until(j>=N)REPEAT-UNTILstatementsprovidelogicaloperationsinthesamemannerasinIF-THEN-ELSEstatements,buttheyareoftensomewhatmoreconvenient,especiallyforloopingoperations.TheREPEAT-UNTILstatementsarewritteninthesamemannerasinPascal.HereisasmallfunctionwhichusestheREPEAT-UNTILstatements.FUNCTIONSUMINV(X)SUMX:=0;i:=0REPEATi:=i+1SUMX:=SUMX+1/X^iUNTIL(i>10)SUMINV:=SUMXEND下圖為X=5時(shí)，函數(shù)的運(yùn)行結(jié)果NotethatREPEAT-UNTILstatementscanonlybeusedininternalfunctionsorprocedures.NotethatDuplicateloopsshouldgenerallynotbeusedininternalfunctionsorprocedures.TheDuplicatestatementliterallyduplicatestheequation,substitutingtheduplicateindexasappropriate.ThisconstructionisdoneduringthecompilingstagesoitisthennotpossibletopassthelowerorupperboundoftheDuplicateloopasa