馬后炮化工論壇-dynamic

AD、ACM初級(jí)培訓(xùn)游存芳2005.6.202023/2/61IntroductiontoAspenDynamicsAD:AspenDynamics

ACM:AspenCustomModeler2023/2/62IntroductiontoAspenDynamics課程安排-Day11. AspenDynamics概述

例子-AspenDynamics的應(yīng)用

練習(xí)-脫乙烷塔2. 創(chuàng)建一個(gè)動(dòng)態(tài)模擬過(guò)程

練習(xí)-加動(dòng)態(tài)數(shù)據(jù)3. 運(yùn)行一個(gè)動(dòng)態(tài)模擬過(guò)程

練習(xí)-動(dòng)態(tài)模擬2023/2/63IntroductiontoAspenDynamics4. 建模要點(diǎn)

練習(xí)-RPlugThermalInertia

練習(xí)-OverfilledVessel2023/2/64IntroductiontoAspenDynamics課程安排-Day25.腳本

練習(xí)-Scripts6. 任務(wù)語(yǔ)言

練習(xí)-Tasks7. 壓力驅(qū)動(dòng)的過(guò)程模擬

練習(xí)-PressureDrivenSimulation

2023/2/65IntroductiontoAspenDynamics8.過(guò)程控制

練習(xí)-PIDControllerTuning

練習(xí)-串級(jí)控制(可選)2023/2/66IntroductiontoAspenDynamicsAspenDynamics簡(jiǎn)介學(xué)習(xí)AspenDynamics的目的Aspendynamics的應(yīng)用實(shí)例AspenDynamics用戶(hù)界面2023/2/67IntroductiontoAspenDynamicsAspenDynamics和AspenCustomModelerFeatureACMADACM+ADRunGUIand

calculationsxxxCreatenewmodelsxxUsetheDynamics

libraryxxCallPropertiesPlusyesyesyes2023/2/68IntroductiontoAspenDynamics工作流AspenPropertiesPlusAspenPlusGUICreateandEditCustomModelsBrowseModelsandModifyFlowsheetAspenDynamicsAspenCustomModelerExportTransferModels2023/2/69IntroductiontoAspenDynamics什么是AspenDynamics?用戶(hù)界面自動(dòng)生成動(dòng)態(tài)模擬輸入定義壓力驅(qū)動(dòng)的流程模擬流量驅(qū)動(dòng)的流程模擬

可以瀏覽動(dòng)態(tài)模型數(shù)學(xué)方程式必要時(shí)可以利用圖形界面修改流程引入已有的流程模塊動(dòng)態(tài)模擬的自動(dòng)初始化利用AspenPlus的計(jì)算結(jié)果2023/2/610IntroductiontoAspenDynamics自動(dòng)插入內(nèi)置控制器可用GUI組態(tài)用戶(hù)控制器可以引入已有的控制模塊三相(vapor-liquid-liquid平衡)計(jì)算能力電解質(zhì)處理能力可與PolymersPlus集成支持石油餾分的分析/調(diào)和功能,包括石油物性的計(jì)算什么是AspenDynamics?2023/2/611IntroductiontoAspenDynamics什么是AspenDynamics?支持以下模擬運(yùn)行模式SteadyStateInitializationDynamicOptimizationSteadyStateDynamic(OptimalControl)EstimationHomotopy（同倫拓展法）Note:本課程不涵蓋optimization,estimationorhomotopy等運(yùn)行模式2023/2/612IntroductiontoAspenDynamics什么是AspenDynamics利用熱量計(jì)的結(jié)果可以預(yù)估動(dòng)力學(xué)參數(shù)SimulationAccesseXtensions(SAX)可以利用外部程序接口訪問(wèn)和控制模擬過(guò)程WindowsOLEAutomationNote:Estimation,andSAX不包括在本課程內(nèi)

2023/2/613IntroductiontoAspenDynamics什么是AspenDynamics利用OLEforProcessControl(OPC)實(shí)現(xiàn)在線連接標(biāo)準(zhǔn)Windows系統(tǒng)可以與過(guò)程控制系統(tǒng)通訊.能夠與控制系統(tǒng)或擁有OPCserver的其他應(yīng)用系統(tǒng)交換數(shù)據(jù)，不需要專(zhuān)門(mén)的用戶(hù)接口程序.動(dòng)態(tài)模擬可用于:培訓(xùn)軟測(cè)量和數(shù)據(jù)預(yù)估以當(dāng)前裝置條件自動(dòng)初始化的動(dòng)態(tài)預(yù)測(cè)模型2023/2/614IntroductiontoAspenDynamics什么是AspenDynamics自動(dòng)生成DMCPlus線性模型與DMCPlus多變量控制器的接口可在AspenDynamics內(nèi)部使用DMCPlus控制器可在AspenDynamics內(nèi)進(jìn)行階躍實(shí)驗(yàn)可將結(jié)果以DMCPlusCollect文件格式輸出來(lái)進(jìn)行模型識(shí)別Note:DMCPlus接口內(nèi)容不在本課程2023/2/615IntroductiontoAspenDynamics什么是AspenDynamicsControlDesignInterface(CDI)GeneratelinearizedmodelinformationforuseincontrolpackagessuchasMATLAB?State-space(A,B,C,D)matricesSteady-stategainmatricesRelativeGainArray(RGA)CalculatetheNiederlinskiIndex(NI)CalculatetheMorariResiliencyIndex(MRI)Note:CDIisnotcoveredinthiscourse2023/2/616IntroductiontoAspenDynamics目標(biāo):ShowhowtoapplyAspenDynamics,includingthekeyfeaturesoftheuserinterface.例子:AspenDynamics的應(yīng)用2023/2/617IntroductiontoAspenDynamics例子:AspenDynamics的應(yīng)用水-甲醇的閃蒸過(guò)程Flowrate100Kmol/hrTemperature50CPressure2BarMole-Fractions:Water0.5Methanol0.5Vaporfraction0.5Pressuredrop0.0atmVerticalVesselLength3.0mDiameter2.0mConstantdutyheattransferInitialliquidfillagefraction0.5FEEDVAPLIQFLASH2023/2/618IntroductiontoAspenDynamics初始狀態(tài)為收斂的流程模擬:StartExample.bkp利用AspenPlus增加動(dòng)態(tài)數(shù)據(jù)利用動(dòng)態(tài)按鈕訪問(wèn)動(dòng)態(tài)數(shù)據(jù)文件夾動(dòng)態(tài)數(shù)據(jù)表單需要如下數(shù)據(jù)：設(shè)備幾何尺寸

過(guò)程熱傳遞方法初始持液量設(shè)備熱傳遞包括對(duì)環(huán)境的熱損失自動(dòng)生成AspenDynamics的問(wèn)題文件并帶入AspenDynamics中例子:AspenDynamics的應(yīng)用2023/2/619IntroductiontoAspenDynamics用戶(hù)界面介紹

ProcessFlowsheetwindowSimulationExplorerAllItemspaneContentspaneSimulationMessageswindowMenubarToolsbuttonsStatusBarSpecificationstatuswindowSpecificationanalysiswindow例子:AspenDynamics的應(yīng)用2023/2/620IntroductiontoAspenDynamics動(dòng)態(tài)模擬的初始化使用最新的“照相”數(shù)據(jù)運(yùn)行模擬計(jì)算并從預(yù)制好的表格和曲線分析結(jié)果生成定制的表格和曲線可以通過(guò)增加新的PID控制器模塊來(lái)修改流程控制方案例子:AspenDynamics的應(yīng)用2023/2/621IntroductiontoAspenDynamics用戶(hù)界面TitlebarMenubarToolbarsSimulationExplorerStatusbarFlowsheetWindowSimulationMessagesWindow2023/2/622IntroductiontoAspenDynamics用戶(hù)界面帶有通常下拉菜單的菜單欄Context-sensitivemenus可以快速執(zhí)行某項(xiàng)動(dòng)作的工具欄Run,Pause(tostop),Restart,Rewind,NewPlot,NewTable,etc.包括在線的和隨機(jī)的幫助關(guān)于一些功能和特點(diǎn)的幫助關(guān)于建模的假定和模型框架的幫助關(guān)于AspenDynamics提供的例子的幫助2023/2/623IntroductiontoAspenDynamics幫助主題窗口2023/2/624IntroductiontoAspenDynamics用戶(hù)界面工藝流程窗口:設(shè)備模塊和物流的圖形化連接設(shè)備模塊通過(guò)物流相連接物流連接到接口接口附著在設(shè)備模塊上2023/2/625IntroductiontoAspenDynamics工藝流程(PFS)窗口流程圖窗口中可以點(diǎn)擊鼠標(biāo)右鍵彈出菜單2023/2/626IntroductiontoAspenDynamics工藝流程窗口流程通過(guò)模塊的“拖和放”來(lái)構(gòu)建使用鼠標(biāo)左鍵將圖例從模擬管理器中拖出放到工藝流程窗口點(diǎn)擊鼠標(biāo)右鍵彈出流程圖操作菜單FlowsheetRMBactionsinclude:Zoom,Pan,Redraw,Print,FindObject,etc…塊和物流同樣擁有自己的彈出菜單2023/2/627IntroductiontoAspenDynamics工藝流程窗口塊和物流的右鍵彈出菜單：訪問(wèn)表格模塊和物流的屬性

模塊和物流的重新命名和刪除對(duì)圖例進(jìn)行移動(dòng),旋轉(zhuǎn),放大縮小或更換圖例重新連接物流的源頭和目的地調(diào)用計(jì)算結(jié)果表格,曲線或其他窗體調(diào)用變量尋找窗口等等...2023/2/628IntroductiontoAspenDynamics工藝流程窗口接口物流通過(guò)連接到設(shè)備模塊上的接口將模塊連接起來(lái)通過(guò)拖動(dòng)接口可將接口放在圖例周邊的任何位置紅色接口時(shí)必需的

蘭色接口是可選的單個(gè)接口允許多股物流的連接2023/2/629IntroductiontoAspenDynamics模擬管理器窗口項(xiàng)目總匯窗格包含了在內(nèi)容窗格中顯示的目標(biāo)所在的文件夾內(nèi)容窗格中顯示的是在項(xiàng)目總匯窗格中選中的文件夾中的內(nèi)容2023/2/630IntroductiontoAspenDynamicsSimulationExplorerWindowSimulationExplorerwithRMBMenuAllItemsPaneContentsPaneExpand/CollapseContents2023/2/631IntroductiontoAspenDynamics模擬管理窗口模擬文件夾包括:庫(kù)動(dòng)態(tài)庫(kù),包括模型物流類(lèi)型其他庫(kù)(如果有許可證的話)用戶(hù)模型庫(kù)

由ACM用戶(hù)創(chuàng)建的用戶(hù)庫(kù)TRAFLOW(two-phasepipelinedynamicsimulationmodule)Note:TRAFLOWisnotcoveredinthiscourse2023/2/632IntroductiontoAspenDynamics模擬管理窗口模擬文件夾包括:流程文件夾

樹(shù)狀顯示用于當(dāng)前模擬的設(shè)備模塊和物流以及相關(guān)的預(yù)先編制好的腳本、表格和曲線表格(用戶(hù)的和預(yù)先設(shè)置的)包括歷史表格在內(nèi)的結(jié)果表格連續(xù)和分布的曲線局部變量表(局部變量的應(yīng)用限定在本流程的編輯器內(nèi))任務(wù)腳本2023/2/633IntroductiontoAspenDynamics模擬管理窗口模擬文件夾包括:解法器選擇窗口可以選擇公式的解法,積分器,優(yōu)化器,診斷輸出,誤差以及其他的模擬選擇2023/2/634IntroductiontoAspenDynamicsSolverPropertiesWindow2023/2/635IntroductiontoAspenDynamics模擬管理窗口模擬文件夾包括:全局變量表應(yīng)用在模擬中的全局變量的詳細(xì)內(nèi)容2023/2/636IntroductiontoAspenDynamicsGlobalsTableForm全局變量表2023/2/637IntroductiontoAspenDynamics模擬管理窗口模擬文件夾包括:組分表流程中所定義的組分名稱(chēng)和物性方法表2023/2/638IntroductiontoAspenDynamics組分表窗口訪問(wèn)物性選擇窗口的按鈕組分表窗口物性選擇窗口2023/2/639IntroductiontoAspenDynamics模擬管理窗口模擬文件夾包括:診斷

報(bào)告計(jì)算公式和變量在模擬計(jì)算中的進(jìn)展信息(只在運(yùn)行計(jì)算中或計(jì)算后得到)2023/2/640IntroductiontoAspenDynamics模擬信息窗口顯示來(lái)自解法器的診斷信息

顯示來(lái)自腳本,任務(wù)的輸出,編輯等的信息…右鍵彈出菜單動(dòng)作:改變解法器的報(bào)告等級(jí)

選擇信息輸出地點(diǎn)屏幕文件清除模擬信息窗口的內(nèi)容2023/2/641IntroductiontoAspenDynamics模擬信息窗口SimulationMessagesWindowwithRMBMenu2023/2/642IntroductiontoAspenDynamics定義狀態(tài)條和窗口定義狀態(tài)條和窗口

顏色會(huì)告訴使用者定義的狀態(tài)GreenforGo當(dāng)流程被修改時(shí)會(huì)自動(dòng)更新定義雙擊狀態(tài)按鈕將訪問(wèn)定義狀態(tài)窗口激活定義分析窗口記錄定義的變化激活變量尋找窗口(將在后面討論)創(chuàng)建流程腳本(腳本將在后面討論)2023/2/643IntroductiontoAspenDynamics定義狀態(tài)窗口SpecificationStatusWindow-General2023/2/644IntroductiontoAspenDynamics定義狀態(tài)窗口SpecificationStatusWindow-FixedChanges2023/2/645IntroductiontoAspenDynamics作業(yè)(15分鐘):脫乙烷塔

目的:用AspenPlus

對(duì)脫乙烷塔進(jìn)行穩(wěn)態(tài)模擬

完成塔的設(shè)計(jì)計(jì)算，用嚴(yán)格精餾法確定塔徑2023/2/646IntroductiontoAspenDynamics作業(yè):脫乙烷塔

(1)

工藝描述這是對(duì)來(lái)自乙烯廠的一個(gè)脫乙烷塔的模擬.進(jìn)料來(lái)自脫甲烷塔且溫度較脫乙烷塔低.因此進(jìn)料首先用C2預(yù)熱以適應(yīng)脫乙烷塔的工況.閃蒸后的氣相和液相分別送到塔的最合適的進(jìn)料位置.再沸器用急冷水加熱，冷凝器用C3冷卻.控制目標(biāo)是要分開(kāi)C2s和C3s

.盡管塔頂?shù)腃3s超標(biāo)會(huì)導(dǎo)致下游C2分離出現(xiàn)問(wèn)題，但仍將塔底的C2s濃度作為硬性指標(biāo).2023/2/647IntroductiontoAspenDynamicsWorkshop:DeethanizerTower

(2)E101C2INHCFEEDC2OUTFLASHFDT101VAPFEEDLIQFEEDC2PRODC3PRODD101C2OUT2023/2/648IntroductiontoAspenDynamics作業(yè):脫乙烷塔(3)

已經(jīng)準(zhǔn)備了一個(gè)backup文件Start-Deethanizer.bkp.工藝流程已經(jīng)構(gòu)建好,組分表,物性方法,進(jìn)料物流(C2INandHCFEED)組分都已規(guī)定完畢.剩下的工作是要補(bǔ)充其他所需數(shù)據(jù),并且進(jìn)行塔的設(shè)計(jì)計(jì)算(求出塔徑).1.打開(kāi)ASPENPLUS并打開(kāi)Start-Deethanizer.bkp文件2.用下頁(yè)提供的數(shù)據(jù)完成模擬所需的條件輸入3.運(yùn)行模擬計(jì)算2023/2/649IntroductiontoAspenDynamics作業(yè):脫乙烷塔(4)PropertyMethod:PSRKFeedStreams

HCFEED C2INMassFlowrateKg/hr: 92,000 150,000TemperatureC: -17 0PressureBar: 21.0 25Molefractions: Methane: 0.003 Ethane: 0.145 1.000 Propane: 0.032 Ethylene: 0.615 Propylene: 0.2052023/2/650IntroductiontoAspenDynamics作業(yè):脫乙烷塔(5)

FlashdrumD101:FLASH2blockwithadiabaticflashHeatDuty=0.0Pressuredrop=-0.1Bar(value=<0.0ispressuredrop)2023/2/651IntroductiontoAspenDynamics作業(yè):脫乙烷塔(6)

Pre-heaterE101:ShortcutHEATXblockwithcold-sideoutlettemperaturespecified.Hotsideinletstream:C2INColdsideinletstream:HCFEEDSpecification:Coldstreamoutlettemperature=-12.8CHotsidepressuredrop=-0.2Bar(value=<0.0ispressuredrop)Coldsidepressuredrop=-0.2Bar(value=<0.0ispressuredrop)2023/2/652IntroductiontoAspenDynamics作業(yè):脫乙烷塔(7)

De-ethanizerColumnT101:Numberofstages=21TotalcondenserDistillaterate(mass)=62541Kg/hrRefluxratio(molar)=0.6Feeds:Liquidtostage13,ON-STAGE;Vaportostage17,ON-STAGEPressures:Stage1/Condenser=20bar;Stage2pressure=20.1barStagepressuredrop=0.01bar2023/2/653IntroductiontoAspenDynamics作業(yè):脫乙烷塔(8)

ColumnResultsSummary-Condenser/Topstage2023/2/654IntroductiontoAspenDynamics作業(yè):脫乙烷塔(9)

ColumnResultsSummary-Reboiler/Bottomstage2023/2/655IntroductiontoAspenDynamics作業(yè):脫乙烷塔(10)

De-ethanizerColumnTraySizing:4.FromtheDatabrowserexplorerview,expandtheblocksobject.UnderBlocks,selectblockT101objectandexpand.SelectTraySizingandcreateanewTraySizingobjectwiththefollowingdata;Startingstage:2,Endingstage:20Traytype:SieveNumberofpasses=25.重新運(yùn)行計(jì)算–注意計(jì)算出的塔板直徑是多少?6.保存你的模擬為backup文件，下次作業(yè)要用.2023/2/656IntroductiontoAspenDynamics作業(yè):脫乙烷塔(11)

TraySizingResultsSummary2023/2/657IntroductiontoAspenDynamics作業(yè):脫乙烷塔Review2023/2/658IntroductiontoAspenDynamics創(chuàng)建一個(gè)動(dòng)態(tài)模擬Objective:DescribeandunderstandthedatarequiredtocreateadynamicsimulationinputfromaconvergedAspenPlusflowsheet2023/2/659IntroductiontoAspenDynamics完成動(dòng)態(tài)模擬輸入這部分將要介紹以下問(wèn)題：訪問(wèn)動(dòng)態(tài)數(shù)據(jù)表格增加動(dòng)態(tài)模擬所需的數(shù)據(jù)

創(chuàng)建動(dòng)態(tài)問(wèn)題文件2023/2/660IntroductiontoAspenDynamicsToviewDynamicToolbarmakesureDynamiccheckboxisselectedfromToolbarsdialogueboxunderViewmenu訪問(wèn)動(dòng)態(tài)數(shù)據(jù)表格先從動(dòng)態(tài)工具條中調(diào)出動(dòng)態(tài)按鈕然后按下2023/2/661IntroductiontoAspenDynamics增加動(dòng)態(tài)數(shù)據(jù)動(dòng)態(tài)數(shù)據(jù)用于以下計(jì)算:容器的幾何尺寸(計(jì)算容器體積時(shí)所需)容器的初始持液狀態(tài)(用于計(jì)算持液量)過(guò)程的熱傳遞方法設(shè)備熱傳遞選擇設(shè)備的熱容

環(huán)境的熱損失嚴(yán)格精餾塔的水利學(xué)和壓降計(jì)算2023/2/662IntroductiontoAspenDynamics容器的幾何動(dòng)態(tài)數(shù)據(jù)容器類(lèi)別瞬間平衡

對(duì)大多數(shù)容器的缺省類(lèi)型–不需要容器的幾何尺寸立式臥式容器幾何尺寸封頭類(lèi)別橢圓半球平面長(zhǎng)度直徑2023/2/663IntroductiontoAspenDynamics容器的幾何數(shù)據(jù)表VesselGeometryDataForm-FlashDrum2023/2/664IntroductiontoAspenDynamics容器的幾何結(jié)構(gòu)2023/2/665IntroductiontoAspenDynamics容器的初始條件表2023/2/666IntroductiontoAspenDynamics過(guò)程熱傳遞選擇恒定熱負(fù)荷

缺省定義(不需要增加進(jìn)一步數(shù)據(jù))恒定的熱/冷媒溫度熱負(fù)荷取決于工藝介質(zhì)溫度和熱/冷媒溫度的溫差當(dāng)穩(wěn)態(tài)熱負(fù)荷Q等于零時(shí)，缺省為恒定熱負(fù)荷計(jì)算LMTD(logmeantemperaturedifference)熱負(fù)荷取決于工藝介質(zhì)溫度和熱/冷媒溫度的對(duì)數(shù)平均溫差

2023/2/667IntroductiontoAspenDynamics過(guò)程熱傳遞選擇表HeatTransferOptionForm2023/2/668IntroductiontoAspenDynamics恒定熱負(fù)荷選擇在動(dòng)態(tài)模擬中負(fù)荷是一個(gè)“fixed”類(lèi)型的變量初始值為在穩(wěn)態(tài)AspenPlus計(jì)算的結(jié)果在動(dòng)態(tài)模擬中熱負(fù)荷能夠被控制

可以通過(guò)人為改變數(shù)值直接控制也可以通過(guò)PID調(diào)節(jié)器控制2023/2/669IntroductiontoAspenDynamics恒定溫度選擇在動(dòng)態(tài)輸入表格中定義熱/冷媒的溫度在動(dòng)態(tài)模擬中可以直接或用調(diào)節(jié)器改變熱/冷媒的溫度計(jì)算方程:Q=UA.(T_process–T_medium)這里:Q =熱負(fù)荷UA =總傳熱系數(shù)與傳熱面積乘積T_process =工藝介質(zhì)溫度T_medium =熱/冷媒溫度2023/2/670IntroductiontoAspenDynamicsLMTD選擇[1]Q=UA.LMTD [2]Q=Fmmed.Cpmed.(Tmed_Out–Tmed_In)變量

說(shuō)明LMTD 對(duì)數(shù)平均溫差Fmmed 熱/冷媒質(zhì)量流量Tmed_In 熱/冷媒進(jìn)口溫度Tmed_Out 熱/冷媒出口溫度Cpmed 熱/冷媒定義熱容UA 總傳熱系數(shù)與傳熱面積乘積2023/2/671IntroductiontoAspenDynamicsLMTD選擇Thisvariable HasavaluethatTmed_in “fixed”型變量,但是能在動(dòng)態(tài)模擬中控制T_approach 在動(dòng)態(tài)模擬中是一個(gè)變化的值Cpmed 在動(dòng)態(tài)模擬中為“fixed”型變量2023/2/672IntroductiontoAspenDynamics設(shè)備傳熱當(dāng)設(shè)備溫度發(fā)生較大變化時(shí)設(shè)備的熱容就顯得十分重要，這種情況一般發(fā)生在開(kāi)車(chē)、停車(chē)或泄壓時(shí)。所需數(shù)據(jù)設(shè)備質(zhì)量

設(shè)備熱容2023/2/673IntroductiontoAspenDynamics設(shè)備傳熱環(huán)境傳熱

當(dāng)工藝過(guò)程對(duì)環(huán)境溫度的變化敏感時(shí)，環(huán)境的傳熱就十分重要了所需數(shù)據(jù)設(shè)備重量

設(shè)備熱容

總傳熱系數(shù)環(huán)境溫度在全局定義表單上輸入環(huán)境溫度(或壓力)2023/2/674IntroductiontoAspenDynamics設(shè)備熱傳遞表EquipmentHeatTransferForm-FlashDrum2023/2/675IntroductiontoAspenDynamicsSetupSpecificationsGlobalSheet2023/2/676IntroductiontoAspenDynamics嚴(yán)格精餾法的水利學(xué)計(jì)算簡(jiǎn)單計(jì)算(缺省)TrayPacking帶有壓力修正的嚴(yán)格計(jì)算TraysPacking2023/2/677IntroductiontoAspenDynamics嚴(yán)格水利學(xué)計(jì)算表格2023/2/678IntroductiontoAspenDynamics嚴(yán)格精餾核算定義表單RadFracRatingSpecificationForm2023/2/679IntroductiontoAspenDynamics嚴(yán)格精餾核算的壓力修正表單壓力修正確認(rèn)框2023/2/680IntroductiontoAspenDynamics生成AspenDynamics文件輸出生成并存為Aspen動(dòng)態(tài)問(wèn)題文件runid.dynf生成并存為計(jì)算物性所需的AspenPropertiesPlusDefinition文件(APPDF)runiddyn.appdf發(fā)送具有輸出功能且自動(dòng)開(kāi)始Aspen動(dòng)態(tài)并調(diào)入輸出問(wèn)題文件runid#.dynfrunid#.appdf2023/2/681IntroductiontoAspenDynamics輸出/發(fā)送動(dòng)態(tài)模擬文件輸出

菜單欄發(fā)送到AspenDynamics2023/2/682IntroductiontoAspenDynamics輸出動(dòng)態(tài)模擬文件注意:你只能輸出/發(fā)送以下?tīng)顟B(tài)的模擬:有結(jié)果的模擬有警告結(jié)果的模擬2023/2/683IntroductiontoAspenDynamics作業(yè)(15分鐘):增加動(dòng)態(tài)數(shù)據(jù)目的:用AspenPlus把動(dòng)態(tài)數(shù)據(jù)增加到穩(wěn)態(tài)模擬中去.完成

Radfrac

塊的塔板核算，進(jìn)行塔的水利學(xué)計(jì)算.2023/2/684IntroductiontoAspenDynamics作業(yè):增加動(dòng)態(tài)數(shù)據(jù)

(1)這是前一個(gè)作業(yè)的延續(xù).你可以使用提供的文件也可以使用上一個(gè)作業(yè)的結(jié)果.1.調(diào)用文件Start-DynData.bkp–這是上一個(gè)作業(yè)的結(jié)果2.在上一個(gè)作業(yè)中增加所需的動(dòng)態(tài)數(shù)據(jù)(數(shù)據(jù)如下)到穩(wěn)態(tài)計(jì)算中–在View/Toolbar菜單中調(diào)入動(dòng)態(tài)按鈕E101Pre-heaterAssumeinstantaneous(steady-state)operationD101Flashdrum:

Verticalvessel,withconstantheatdutyheatingoption

Length=3m,Diameter=2m

Initialliquidvolumefraction=0.52023/2/685IntroductiontoAspenDynamics作業(yè):增加動(dòng)態(tài)數(shù)據(jù)

(2)T101De-ethanizerColumn:Condenser: Constanttemperature(-45C)coolingmediumReboiler: LMTDheattransferoption, Quenchwatermediumtemperature:80

C Temperatureapproach:5

CRefluxDrum: Horizontalvessel Length:5m,Diameter:2mReboilerSump:Height=5m,Diameter=2.8m(fromtraysizingresult)Hydraulics: Rigorousoption 2023/2/686IntroductiontoAspenDynamics作業(yè):增加動(dòng)態(tài)數(shù)據(jù)

(3)3.Includeatrayrating(TrayRatingdataform)sectiontocalculatethecolumnhydraulicsTrayRatingforT101Column

Trayedsection:Stages2to20

SieveTray,2-pass

Diameter=2.8m(fromtraysizingresults)

Updatepressureprofile(checkononDesign/Pdroptab)4.運(yùn)行穩(wěn)態(tài)計(jì)算

5.保存backup文件.6.用File/Export生成.dynf和.appdf文件–選“FlowDrivenDynSimulation”文件類(lèi)型.2023/2/687IntroductiontoAspenDynamics作業(yè):增加動(dòng)態(tài)數(shù)據(jù)

(5)TrayRatingResultsSummary2023/2/688IntroductiontoAspenDynamics增加動(dòng)態(tài)數(shù)據(jù)Review2023/2/689IntroductiontoAspenDynamics運(yùn)行動(dòng)態(tài)模擬Objective:Describeandbecomefamiliarwiththeuserinterfacefeaturesrequiredforrunningdynamicsimulationsandviewingtheresults2023/2/690IntroductiontoAspenDynamics運(yùn)行動(dòng)態(tài)模擬這部分將要介紹以下問(wèn)題;運(yùn)行控制運(yùn)行動(dòng)態(tài)模擬的工具定義狀態(tài)和狀態(tài)欄水利壓頭計(jì)算模擬和系統(tǒng)工作文件2023/2/691IntroductiontoAspenDynamics運(yùn)行控制運(yùn)行模式初始化穩(wěn)態(tài)動(dòng)態(tài)運(yùn)行控制工具RunStepPauseRestartRewind運(yùn)行選擇通訊間隔2023/2/692IntroductiontoAspenDynamics運(yùn)行控制工具按鈕RunStepPauseRestartRewindReset運(yùn)行模式2023/2/693IntroductiontoAspenDynamics運(yùn)行控制菜單2023/2/694IntroductiontoAspenDynamics運(yùn)行選擇窗口2023/2/695IntroductiontoAspenDynamics所有變量的記錄時(shí)間歷史從運(yùn)行選擇對(duì)話框中可選擇保存所有變量的記錄時(shí)間歷史也可對(duì)某一表中所列變量保存記錄時(shí)間歷史從表格編輯屬性中選出記錄功能加到顯示欄中也可對(duì)變量尋找工具找到的變量選擇記錄時(shí)間歷史注意:當(dāng)顯示歷史表格或曲線時(shí)缺省選擇是變量的時(shí)間歷史不做記錄，在曲線顯示之前的數(shù)據(jù)不顯示2023/2/696IntroductiontoAspenDynamics變量的記錄時(shí)間歷史Plotcreatedaftertime0.1Hours-Nohistorydatarecordedpreviously2023/2/697IntroductiontoAspenDynamics解法選擇窗口2023/2/698IntroductiontoAspenDynamics工具設(shè)置–用戶(hù)缺省選擇照相和使用照相表單–查看結(jié)果測(cè)量單位-UOM變量尋找2023/2/699IntroductiontoAspenDynamics工具按鈕建立歷史表建立表建立分布曲線建立曲線打開(kāi)管理器變量尋找照相初始化模擬2023/2/6100IntroductiontoAspenDynamics工具菜單ToolsMenu2023/2/6101IntroductiontoAspenDynamics設(shè)置窗口SettingsWindow-UserDefaultPreferences2023/2/6102IntroductiontoAspenDynamics照相照相是保存在模擬過(guò)程中所有變量的結(jié)果（定義的或計(jì)算值）當(dāng)前照相適應(yīng)當(dāng)前流程檔案照相當(dāng)流程或模擬結(jié)構(gòu)發(fā)生任何變化時(shí)，照好的“底片”被打包并轉(zhuǎn)變成單精度的二進(jìn)制*.snp文件保存在問(wèn)題文件夾中。保存的結(jié)果

照相被保存為ASCII語(yǔ)言(dynfextension)文件.2023/2/6103IntroductiontoAspenDynamics照相自動(dòng)拍照初始定義(在運(yùn)行前)穩(wěn)態(tài)運(yùn)行初始化運(yùn)行預(yù)估運(yùn)行優(yōu)化Homotopy在零時(shí)間的動(dòng)態(tài)初始化(收斂的)預(yù)先定義的照相

在定義好的時(shí)間間隔在重新初始化安排的照相手工創(chuàng)建的照相2023/2/6104IntroductiontoAspenDynamics照相管理的一般窗口2023/2/6105IntroductiontoAspenDynamics照相管理的創(chuàng)建窗口2023/2/6106IntroductiontoAspenDynamics照相的使用模擬計(jì)算在零時(shí)間的重新開(kāi)始回到以前某一時(shí)刻的狀態(tài)選擇先前保存的拍照“底片”流程初始化可以從當(dāng)前結(jié)果選擇照相或從先前保存的照相中選擇狀態(tài)在各自的區(qū)域內(nèi)輸入所要求的變量名字形式(在變量尋找工具部分將討論變量名字形式)注意:對(duì)rewind或restart,如果模擬時(shí)間超出當(dāng)前時(shí)間范圍，“定時(shí)拍照將取消.2023/2/6107IntroductiontoAspenDynamics照相返回窗口2023/2/6108IntroductiontoAspenDynamics表格預(yù)定的表格類(lèi)型

表歷史表單/時(shí)間連續(xù)曲線輸出DMCPlus數(shù)據(jù)文件(*.clc)用于DMCPlus建模中的模型辯識(shí)分布曲線(stage-wiseprocessesonly)Microsoft?ActiveX?Controlform(例如PID控制器面板和組態(tài)表)流程用戶(hù)定制表格Note:CanAccessResultswithWindowsAutomation2023/2/6109IntroductiontoAspenDynamics預(yù)先設(shè)置的表格PredefinedForms-RadFrac2023/2/6110IntroductiontoAspenDynamics預(yù)先設(shè)置的表格–示例AllVariablesConfigureConfigureSensor(streamsonly)ManipulateResultsProfileTable(stage-wiseprocessesonly)TimeSeriesPlot/HistoryTable(streamsonly)ProfilePlot2023/2/6111IntroductiontoAspenDynamics預(yù)先設(shè)置表格-ConfigureSensorConfigureSensor表格（僅對(duì)物流）

可以訪問(wèn)用于物流附加物性計(jì)算的開(kāi)關(guān)相分率相組成相密度體積流率pH石油性質(zhì)從結(jié)果表格查看所有結(jié)果

2023/2/6112IntroductiontoAspenDynamicsConfigurSensorFormCheckSensorOn2023/2/6113IntroductiontoAspenDynamicsHistoryTable/TimeSeriesPlotFormAccessHistoryTableviaaPlot(andviceversa)RMBonPlottoaccesshistorytable2023/2/6114IntroductiontoAspenDynamics流程用戶(hù)表格NewPlotNewProfilePlotNewTableNewHistoryTableForms-CustomFormsToolButtons2023/2/6115IntroductiontoAspenDynamics從工具菜單選擇用戶(hù)表格Forms-CustomFormsTools2023/2/6116IntroductiontoAspenDynamics從變量尋找窗口定制表單制表按鈕2023/2/6117IntroductiontoAspenDynamics顯示變量屬性SelectRecordattributetoaddtotablecolumns2023/2/6118IntroductiontoAspenDynamics變量屬性2023/2/6119IntroductiontoAspenDynamics定義(SPEC)變量屬性SPEC屬性是識(shí)別變量的類(lèi)型Fixed變量能夠由用戶(hù)控制(階躍或斜坡)Free變量由解法器計(jì)算出來(lái)的

2023/2/6120IntroductiontoAspenDynamics單位變量屬性–測(cè)量單位單個(gè)變量的測(cè)量單位可在表格中改變?nèi)笔〉臏y(cè)量單位與AspenPlus模擬單位一致MetricUSENGMETMETCBARSISI-CBAR從Tools/UnitsofMeasuremenu選擇測(cè)量單位集2023/2/6121IntroductiontoAspenDynamics測(cè)量單位SelectUnitsforvariablefromdropdownlist2023/2/6122IntroductiontoAspenDynamics測(cè)量單位用戶(hù)可用腳本定義自己的測(cè)量單位集用戶(hù)可用腳本增加單個(gè)變量新的顯示測(cè)量單位集Note:MetricUOMsetisthebaseUOMsetforAspenDynamics.2023/2/6123IntroductiontoAspenDynamics水利學(xué)靜壓頭效果液相物流流出下列模塊時(shí)包括了液相壓頭效果Flash2Flash3DecanterMixerRadFracReboilerandCondenserRCSTR設(shè)置LIQUIDHEAD參數(shù)為“YES“-缺省為“NO”全局變量表(所有塊)單個(gè)塊的組態(tài)表單2023/2/6124IntroductiontoAspenDynamics水利學(xué)靜壓頭效果BlockConfigureFormwithLiquidHeadCalculationSwitchLiquidHeadParameter2023/2/6125IntroductiontoAspenDynamics變量尋找使用變量尋找工具包括;尋找和列出流程中塊變量名字

查看變量屬性的當(dāng)前值修改變量屬性生成用戶(hù)表格生成腳本(在后面討論)變量名格式的基礎(chǔ)為完全的“路徑名”路徑可以包括幾個(gè)“fields”2023/2/6126IntroductiontoAspenDynamics變量尋找窗口2023/2/6127IntroductiontoAspenDynamics變量名格式可以用包含統(tǒng)配符格式訪問(wèn)變量名變量名的路徑可以分成幾個(gè)段（fields）一個(gè)field是下列分界字符之間的文字:.(句號(hào))()(圓括號(hào))Note:ThepatternBlocks(“T101”).Stage(*).Tcontainsfourfields2023/2/6128IntroductiontoAspenDynamics變量名格式-統(tǒng)配符符號(hào)

意義* 查找當(dāng)前名字段的零個(gè)或多個(gè)字符~ 查找零個(gè)或多個(gè)名字段? 查找的指定的名字段中的任何單個(gè)字符可以混用統(tǒng)配符2023/2/6129IntroductiontoAspenDynamics變量名格式–示例*與當(dāng)前名字段的零個(gè)或多個(gè)字符相匹配X* 代表 X123X* 不代表 X(“WATER”)X* 不代表 X.Y2023/2/6130IntroductiontoAspenDynamics變量名格式-Examples*Tray.Temp匹配下面的名字;MyTray.TempHerTray.TempHisTray.Temp*Tray.Temp不能代表下面名字;MyColumn.MyTray.TempThesearchpatternhas*TrayinthefirstfieldbutthenameMyColumnisinadifferentfield2023/2/6131IntroductiontoAspenDynamics變量名格式-Examples~matcheszeroormorecharactersinanumberofnamefieldsinthepathX~matchesnamesoftheform X123X~matchesnamesoftheform X(“WATER”)X~matchesnamesoftheform X.YThesearchwillmatchallvariablesthathavethesamenameinthehierarchyforbothblocksandstreams2023/2/6132IntroductiontoAspenDynamics變量名格式-Examples~.*Tray.TempmatchesnamesofthefollowingformsMyTray.TempHerTray.TempHisTray.TempMyColumn.MyTray.TempMyProcess.MyColumn.MyTray.TempThesearchwillmatchallvariablesthathavethesamenameinthehierarchyforbothblocksandstreams2023/2/6133IntroductiontoAspenDynamics變量名格式-ExamplesAll4examplesbelowproducethesameresultBlocks(“T101”).Stage(*).X(*)Blocks(“T101”).S*(*).X(*)Blocks(“T101”).S*~(*).X(*)Blocks(“T101”).S*~(*).X~2023/2/6134IntroductiontoAspenDynamics變量名格式-Examples?matcheszeroorasinglecharacterintheCURRENTnamefieldTray?matchesnamesofformTray1Tray?matchesthenameformTray2Tray?doesnotmatchnamesoftheformTrayTray?doesnotmatchnamesoftheformTray12Note:Seetheon-linehelpformoreexamples2023/2/6135IntroductiontoAspenDynamics模擬和系統(tǒng)工作文件RUNID.BKP AspenPlusbackupRUNID.APW AspenPlus(binary)documentRUNIDDYN.APPDF AspenDynamicsproblem physicalpropertiesdefinitionfileRUNIDDYN.OBJ AspenDynamicsuser FORTRANobjectlibraryRUNIDDYN.HIS AspenDynamicsrunhistory2023/2/6136IntroductiontoAspenDynamics模擬和系統(tǒng)工作文件RUNID.DYNF AspenDynamicsproblemfile forflowdrivensimulationsRUNID.DYND AspenDynamicsproblem (binary)document..\RUNID AspenDynamicsworkingproblem systemfilessubdirectory,including snapshots(binary*.snpfiles)2023/2/6137IntroductiontoAspenDynamics作業(yè)(60分鐘):動(dòng)態(tài)模擬

目的:練習(xí)和熟悉AspenDynamics的特性運(yùn)行動(dòng)態(tài)模擬并實(shí)現(xiàn)一個(gè)階躍變化建立自己的結(jié)果表格、隨時(shí)間變化的曲線和分布曲線修改曲線屬性2023/2/6138IntroductiontoAspenDynamics在線幫助如果需要幫助

則有一個(gè)詳細(xì)的主題

從幫助主題對(duì)話框中點(diǎn)擊索引表.一個(gè)表或區(qū)域

在ASPENPLUS工具欄中, 點(diǎn)擊What‘sThis按鈕然后點(diǎn)擊需要幫助的區(qū)域或表格.一個(gè)對(duì)話框

在對(duì)話框中點(diǎn)擊幫助按鈕.光標(biāo)或鼠標(biāo)放置的位置按F12023/2/6139IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(1)將提供你一個(gè)Start-DynSim.dynfandStartDynSim.appdf文件.這個(gè)文件和上一次輸出的作業(yè)是等同的.1.打開(kāi)AspenDynamics并調(diào)入所給文件–或調(diào)入上一次作業(yè)所輸出的文件.2.從工具菜單下打開(kāi)變量尋找窗口3.確認(rèn)“IncludeParameters”框4.在Variable區(qū)域,輸入下列尋找格式以搜尋所有在這一層次的包含LiquidHead字符的變量~.LiquidHead*

5.從搜尋結(jié)果中,選中D101和T101的冷凝器和再沸器的液相壓頭變量.2023/2/6140IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(2)

2023/2/6141IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(3)6.點(diǎn)擊Table按鈕生成三個(gè)液相壓頭參數(shù)的用戶(hù)表格7.在表格中再增加相關(guān)的物流壓力變量,如下;8.將運(yùn)行方式設(shè)為初始化并運(yùn)行9.觀察物流壓力的變化2023/2/6142IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(4)10.通過(guò)改變液相壓頭參數(shù)將D101和T101冷凝器和再沸器的液相壓頭效果計(jì)算的開(kāi)關(guān)從“No”切到“Yes”11.重新進(jìn)行初始化計(jì)算12.注意物流壓力值是否有很大改變?13.建立一個(gè)新表來(lái)觀察下列變量;a)改變碳?xì)浠衔镞M(jìn)料(streamHCFEED)的質(zhì)量流率的定義(從Manipulate表中選擇FmR變量)b)D101出口物流(VAPFEED和LIQFEED)的質(zhì)量流率的變化(即各自結(jié)果表中的Fm變量)2023/2/6143IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(5)c)D101液位的變化(結(jié)果表中的液位變量)在物流HCFEED的結(jié)果表中哪一個(gè)變量和其控制表中的FmR變量是等同的?14.建立一個(gè)包含新建表格中所有變量的曲線圖15.用RMB的屬性/AxisMap來(lái)修改曲線圖使其所有的流量在同一個(gè)坐標(biāo)軸上而液位在第二個(gè)坐標(biāo)軸上16.打開(kāi)控制器LC1(D101液位控制器)的面板，用plot按鈕打開(kāi)相關(guān)的結(jié)果曲線17.建立一個(gè)新的歷史表格來(lái)觀察塔底乙烷摩爾分率和塔頂丙烷摩爾分率的變化2023/2/6144IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(6)18.用下列步驟建立一個(gè)新的分布曲線來(lái)顯示塔的乙烷和丙烷的液相組成.a)Tools菜單/NewProfilePlotb)輸入曲線的名字C2_C3_Split生成一個(gè)空的分布曲線c)在空?qǐng)D的任何位置,RMB點(diǎn)出Plot的菜單d)選出ProfileVariables彈出ProfileEditore)在Profile1的ProfileName區(qū)間輸入乙烷的名字f)點(diǎn)擊在Y-AxisVariables部分下的Add按鈕

g)輸入第一個(gè)分布的變量名字: Blocks(“T101”).stage(*).X(“Ethane”)Profileeditor如下頁(yè)顯示.2023/2/6145IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(7)

TheProfileEditor2023/2/6146IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(8)

h)將Profile的值改為2I)對(duì)Profile2,輸入丙烷名字j)點(diǎn)擊在Y-AxisVariables部分下的Add按鈕

k)輸入第二個(gè)分布的變量名字: Blocks(“T101”).stage(*).X(“Propane”)l)Ok.19.從Run菜單中選擇RunOptions,將“Communication”間隔設(shè)為0.1hours并將“Pauseat”時(shí)間設(shè)為0.1hours–必須先將運(yùn)行狀態(tài)設(shè)為動(dòng)態(tài)!20.在動(dòng)態(tài)模式下運(yùn)行模擬到0.1hours你能重現(xiàn)T101的下列組分分布嗎?2023/2/6147IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(9)T101CompositionProfilesatsimulationTime0.1Hours2023/2/6148IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(10)21.重新打開(kāi)ProfileEditor(右鍵在profileplot上打開(kāi)profilevariables).在時(shí)間設(shè)定部分,點(diǎn)擊時(shí)間設(shè)定單選按鈕將時(shí)間增至0.522.打開(kāi)碳?xì)浠衔镞M(jìn)料控制表將進(jìn)料的總質(zhì)量流率提高20%（從92,000到110,000kg/hr）.23.將pauseat時(shí)間改為0.5hours24.運(yùn)行模擬至完成你能得到下面的圖形嗎?2023/2/6149IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(11)FeedflowDisturbanceEffects2023/2/6150IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(12)T101CompositionHistory2023/2/6151IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(13)T101CompositionProfilesatsimulationTime0.5Hours2023/2/6152IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(14)ControllerLC1FacePlateatsimulationTime0.5Hours2023/2/6153IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(15)ControllerLC1ResultsPlotatsimulationTime0.5Hours2023/2/6154IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

(16/16)在裝置中這股進(jìn)料擾動(dòng)帶來(lái)的整個(gè)影響是什么?哪一個(gè)變量是最容易受到影響的?.閃蒸罐的液位控制器的調(diào)節(jié)正常嗎?為什么液位控制器不能回到它的設(shè)定值?25.將你的問(wèn)題文件保存為AspenDynamicsLanguage(*.dynf)文件類(lèi)型提示:為了更好的計(jì)算,請(qǐng)關(guān)閉你不需要的任何應(yīng)用程序(e.g.AspenPlus,InternetBrowsers,AspenDynamicsFormssuchastables,plots,etc….).2023/2/6155IntroductiontoAspenDynamics作業(yè):動(dòng)態(tài)模擬

Review2023/2/6156IntroductiontoAspenDynamics功能和主要建模特點(diǎn)目的:介紹AspenDynamics模型的功能及局限性.理解基本的建模機(jī)理及一般假設(shè)2023/2/6157IntroductiontoAspenDynamics功能及主要建模特點(diǎn)這部分將要討論以下問(wèn)題;支持功能物性組分類(lèi)型物流模型和物流類(lèi)型單元操作模型能處理滿(mǎn)或空的容器2023/2/6158IntroductiontoAspenDynamics物性完全三相計(jì)算能力(vapor-liquid-liquid平衡)

表觀方法的電解質(zhì)計(jì)算

物性計(jì)算選擇嚴(yán)格的物性計(jì)算“局部”物性方法(缺省)計(jì)算物性例如焓、熵、密度、K值2023/2/6159IntroductiontoAspenDynamics物性計(jì)算局部物性方法物性可以簡(jiǎn)單地?cái)M合成溫度和壓力的函數(shù)比完全的物性計(jì)算要快和完全的物性計(jì)算一樣準(zhǔn)確為求精確局部物性參數(shù)在每一個(gè)積分步長(zhǎng)都進(jìn)行更新2023/2/6160IntroductiontoAspenDynamics物性計(jì)算局部物性計(jì)算例子--液相摩爾焓計(jì)算方程

這里:hL

= 液相摩爾焓xi

=組分I的摩爾分率

AiBi = 組分I的局部物性參數(shù)T = 溫度nc = 組分?jǐn)?shù)2023/2/6161IntroductiontoAspenDynamics物性計(jì)算注意:電解質(zhì)的計(jì)算可能需要“嚴(yán)格”的物性計(jì)算選擇2023/2/6162IntroductiontoAspenDynamics聚合物模擬PolymersPlus是一個(gè)層次化的產(chǎn)品能處理聚合物的物性提供組分的特性描述包括聚合的動(dòng)力學(xué)模型LimitedtomaximumofonepolymerandoneZieglerNattacatalystOnlythefollowingkineticmodelsaresupportedFreeRadicalStepGrowthMulti-siteZieglerNattaNote:Examplesimulationsofpolymer