喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預覽哦。。。下載后都有，，請放心下載，，文件全都包含在內，，【有疑問咨詢QQ:414951605 或 1304139763】
=============================================
喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預覽哦。。。下載后都有，，請放心下載，，文件全都包含在內，，【有疑問咨詢QQ:414951605 或 1304139763】
=============================================
無錫太湖學院
信 機 系 機械工程及自動化 專業(yè)
畢 業(yè) 設 計論 文 任 務 書
一、題目及專題：
1、題目 灌裝生產(chǎn)線上灌裝閥的設計
2、專題
二、課題來源及選題依據(jù)
在現(xiàn)代灌裝廠中，灌裝機和灌裝閥已經(jīng)成為了一個廠的命脈。一個高效率高精度的灌裝閥可以為一個灌裝長帶來強大的活力。灌裝閥的設計需要了解灌裝生產(chǎn)工藝的各項需求，且設計到不少專業(yè)知識。此設計難度中等，設計量合適，可以很大程度提升學生的專業(yè)水準。隨著技術的發(fā)展，生活水準的提高，灌裝技術在未來必將發(fā)揮出巨大的作用。
三、本設計（論文或其他）應達到的要求：
① 能夠詳細的了解灌裝生產(chǎn)線的灌裝工序，熟悉灌裝的每一步的過程，以及各部件的功效；
② 詳細的了解灌裝機的全局機構和布置；
③ 詳細的了解灌裝機的傳動機構，以及撥瓶機構的傳動機構；
④詳細的了解灌裝閥工作的原理，灌裝閥的組成，以及灌裝閥工作時的各個步驟；
四、接受任務學生：
機械97 班 姓名 杭建國
五、開始及完成日期：
自2012年11月7日 至2013年5月25日
六、設計（論文）指導（或顧問）：
指導教師 簽名
簽名
簽名
教研室主任
〔學科組組長研究所所長〕 簽名
系主任 簽名
2012年11月12日
無錫太湖學院
畢業(yè)設計（論文）
開題報告
題目： 液體灌裝生產(chǎn)線灌裝閥的設計
信機 系 機械工程及自動化專業(yè)
學 號： 0923821
學生姓名： 杭建國
指導教師： 何雪明（職稱：副教授 ）
（職稱： ）
2012年11月25日
課題來源
在現(xiàn)代灌裝廠中，灌裝機和灌裝閥已經(jīng)成為了一個廠的命脈。一個高效率高精度的灌裝閥可以為一個灌裝長帶來強大的活力。是對機械專業(yè)學生所學專業(yè)知識的綜合應用。該課題難度適中，設計量合適，能提升學生對知識的應用能力。
科學依據(jù)（包括課題的科學意義；國內外研究概況、水平和發(fā)展趨勢；應用前景等）
灌裝閥是灌裝系統(tǒng)的主要零部件之一，其工作的效率以及灌裝精度，不但影響到灌裝生產(chǎn)線的整體效率，而且影響到整個工廠的效益。在現(xiàn)代隨著國內外的科技的發(fā)展，灌裝生產(chǎn)在人們的生活中發(fā)揮著越來越重大的作用。為了能夠更好的進行灌裝的生產(chǎn)，我們對灌裝閥做了一定的設計改動以便于其能夠更好的工作。以及提高其的生產(chǎn)效率經(jīng)濟效益。
研究內容
1. 根據(jù)灌裝生產(chǎn)工藝，分析零件的灌裝生產(chǎn)的詳細步驟，了解灌裝生產(chǎn)線的總體結構與布置；
2. 分析灌裝生產(chǎn)線的傳動系統(tǒng)，以及各系統(tǒng)如何進行傳送；
3. 設計灌裝閥，并進行分析其是否合理，以及如何運作，是否符合灌裝工藝的各種要求；
擬采取的研究方法、技術路線、實驗方案及可行性分析
先下廠調研，分析現(xiàn)有的一些現(xiàn)有的灌裝閥，然后設計路線。確定灌裝生產(chǎn)線的工作原理，以及各部件的各種功能，確定灌裝閥在灌裝生產(chǎn)線上的作用。了解灌裝閥的工作原理，對其進行分析，并對其進行合理的改動設計，保證其的灌裝性能，以及灌裝效率，提高其的經(jīng)濟效益。
研究計劃及預期成果
研究計劃：
2012年10月12日-2012年12月25日：按照任務書要求查閱論文相關參考資料，填寫畢業(yè)設計開題報告書。
2013年1月11日-2013年3月5日：填寫畢業(yè)實習報告。
2013年3月8日-2013年3月14日：按照要求修改畢業(yè)設計開題報告。
2013年3月15日-2013年3月21日：學習并翻譯一篇與畢業(yè)設計相關的英文材料。
2013年3月22日-2013年4月11日：分析灌裝閥。
2013年4月12日-2013年4月25日：灌裝閥的整體設計。
2013年4月26日-2013年5月20日：畢業(yè)論文撰寫和修改工作。
預期成果：
形成一套完整的加工工藝，可以進行高效且精準的灌裝，擁有合理的結構整體與傳動系統(tǒng)，最后有一套完整的說明書。
特色或創(chuàng)新之處
簡易、造價較低、能夠滿足灌裝效率以及精度要求，提高經(jīng)濟效益。
已具備的條件和尚需解決的問題
已具備：設計思路明確，有著清晰的條理。具備一定的UG圖形處理。
尚需解決：設計細節(jié)還需提升。UG能力還需加強。
指導教師意見
指導教師簽名：
年 月 日
教研室（學科組、研究所）意見
教研室主任簽名：
年 月 日
系意見
主管領導簽名：
年 月 日
英文原文
Screw Compressors
The direction normal to the helicoids, can be used to calculate the coordinates of the rotor helicoids and from x and y to which the clearance is added as:
， ， （2.19）
where the denominator D is given as：
（2.20）
and serve to calculate new rotor end plane coordinates, x0n and y0n,with the clearances obtained for angles θ = /p and τ respectively. These and now serve to calculate the transverse clearance δ0 as the difference between them, as well as the original rotor coordinates and .
If by any means, the rotors change their relative position, the clearance distribution at one end of the rotors may be reduced to zero on the flat side of the rotor lobes. In such a case, rotor contact will be prohibitively long on the flat side of the profile, where the dominant relative rotor motion is sliding, as shown in Fig. 2.29. This indicates that rotor seizure will almost certainly occur in that region if the rotors come into contact with each other.
Fig. 2.29. Clearance distribution between the rotors: at suction, mid rotors, and discharge with possible rotor contact at the discharge
Fig. 2.30. Variable clearance distribution applied to the rotors
It follows that the clearance distribution should be non-uniform to avoid hard rotor contact in rotor areas where sliding motion between the rotors is dominant.
In Fig. 2.30, a reduced clearance of 65 μm is presented, which is now applied in rotor regions close to the rotor pitch circles, while in other regions it is kept at 85 μm, as was done by Edstroem, 1992. As can be seen in Fig. 2.31, the situation regarding rotor contact is now quite different. This is maintained along the rotor contact belt close to the rotor pitch circles and fully avoided at other locations. It follows that if contact occurred, it would be of a rolling character rather than a combination of rolling and sliding or even pure sliding. Such contact will not generate excessive heat and could therefore be maintained for a longer period without damaging the rotors until contact ceases or the compressor is stopped.
2.6 Tools for Rotor Manufacture
This section describes the generation of formed tools for screw compressor hobbing, milling and grinding based on the envelope gearing procedure.
2.6.1 Hobbing Tools
A screw compressor rotor and its formed hobbing tool are equivalent to a pair of meshing crossed helical gears with nonparallel and nonintersecting axes. Their general meshing condition is given in Appendix A. Apart from the gashes forming the cutter faces, the hob is simply a helical gear in which.
Fig. 2.31. Clearance distribution between the rotors: at suction, mid of rotor and discharge with a possible rotor contact at the discharge
Each referred to as a thread, Colburne, 1987. Owing to their axes not being parallel, there is only point contact between them whereas there is line contact between the screw machine rotors. The need to satisfy the meshing equation given in Appendix A, leads to the rotor – hob meshing requirement for the given rotor transverse coordinate points and and their first derivative.The hob transverse coordinate points and can then be calculated. These are sufficient to obtain the coordinate The axial coordinate , calculated directly, and are hob axial plane coordinates which define the hob geometry.
The transverse coordinates of the screw machine rotors, described in the previous section, are used as an example here to produce hob coordinates. he rotor unit leadsare 48.754mm for the main and ?58.504mm for the ate rotor. Single lobe hobs are generated for unit leads :6.291mm for the main rotor and ?6.291mm for the gate rotor. The corresponding hob helix angles ψ are 85? and 95?. The same rotor-to-hob centre distance C = 110mm and the shaft angle Σ = 50? are given for both rotors. Figure 2.32 contains a view to the hob.
Reverse calculation of the hob – screw rotor transformation, also given in Appendix A permits the determination of the transverse rotor profile coordinates which will be obtained as a result of the manufacturing process. These ay be compared with those originally specified to determine the effect of
Fig. 2.32. Rotor manufacturing: hobbing tool left, right milling tool
manufacturing errors such as imperfect tool setting or tool and rotor deformation upon the final rotor profile.
For the purpose of reverse transformation, the hob longitudinal plane coordinatesand andshould be given. The axial coordinate is used to calculate , which is then used to calculate the hob transverse coordinates:
， （2.21）
These are then used as the given coordinates to produce a meshing criterionand the transverse plane coordinates of the “manufactured” rotors.
A comparison between the original rotors and the manufactured rotors is given in Fig. 2.33 with the difference between them scaled 100 times. Two types of error are considered. The left gate rotor, is produced with 30um offset in the centre distance between the rotor and the tool, and the main rotor with
Fig. 2.33. Manufacturing imperfections
0.2? offset in the tool shaft angle Σ. Details of this particular meshing method are given by Stosic 1998.
2.6.2 Milling and Grinding Tools
Formed milling and grinding tools may also be generated by placing in the general meshing equation, given in Appendix A, and then following the procedure of this section. The resulting meshing condition now reads as:
(2.22)
However in this case, when one expects to obtain screw rotor coordinates from the tool coordinates, the singularity imposed does not permit the calculation of the tool transverse plane coordinates. The main meshing condition cannot therefore be applied. For this purpose another condition is derived for the reverse milling tool to rotor transformation from which the meshing angle τ is calculated:
(2.23)
Once obtained, τ will serve to calculate the rotor coordinates after the “manufacturing” process. The obtained rotor coordinates will contain all manufacturing imperfections, like mismatch of the rotor – tool centre distance, error in the rotor – tool shaft angle, axial shift of the tool or tool deformation during the process as they are input to the calculation process. A full account of this useful procedure is given by Stosic 1998.
2.6.3 Quantification of Manufacturing Imperfections
The rotor – tool transformation is used here for milling tool profile generation. The reverse procedure is used to calculate the “manufactured” rotors. The rack generated 5-6 128mm rotors described by Stosic, 1997a are used as given profiles: x(t) and y(t). Then a tool – rotor transformation is used to quantify the influence of manufacturing imperfections upon the quality of the produced rotor profile. Both, linear and angular offset were considered.
Figure 2.33 presents the rotors, the main manufactured with the shaft angle offset 0.5? and the gate with the centre distance offset 40 μm from that of the original rotors given by the dashed line on the left. On the right, the rotors are manufactured with imperfections, the main with a tool axial offset of 40 μm and the gate with a certain tool body deformation which resulted in 0.5? offset of the relative motion angle θ. The original rotors are given by the dashed line.
3Calculation of Screw Compressor Performance
Screw compressor performance is governed by the interactive effects of thermodynamic and fluid flow processes and the machine geometry and thus can be calculated reliably only by their simultaneous consideration. This may be chieved by mathematical modelling in one or more dimensions. For most applications, a one dimensional model is sufficient and this is described in full. 3-D modelling is more complex and is presented here only in outline. A more detailed presentation of this will be made in a separate publication.
3.1 One Dimensional Mathematical Model
The algorithm used to describe the thermodynamic and fluid flow processes in a screw compressor is based on a mathematical model. This defines the instantaneous volume of the working chamber and its change with rotational angle or time, to which the conservation equations of energy and mass continuity are applied, together with a set of algebraic relationships used to define various phenomena related to the suction, compression and discharge of the working fluid. These form a set of simultaneous non-linear differential equations which cannot be solved in closed form.
The solution of the equation set is performed numerically by means of the Runge-Kutta 4th order method, with appropriate initial and boundary conditions.
The model accounts for a number of “real-life” effects, which may significantly influence the performance of a real compressor. These make it suitable for a wide range of applications and include the following:
– The working fluid compressed can be any gas or liquid-gas mixture for which an equation of state and internal energy-enthalpy relation is known, i.e. any ideal or real gas or liquid-gas mixture of known properties.
– The model accounts for heat transfer between the gas and the compressor rotors or its casing in a form, which though approximate, reproduces the overall effect to a good first order level of accuracy.
– The model accounts for leakage of the working medium through the clearances between the two rotors and between the rotors and the stationary parts of the compressor.
– The process equations and the subroutines for their solution are independent of those which define the compressor geometry. Hence, the model can be readily adapted to estimate the performance of any geometry or type of positive displacement machine.
– The effects of liquid injection, including that of oil, water, or refrigerant can be accounted for during the suction, compression and discharge stages.
– A set of subroutines to estimate the thermodynamic properties and changes of state of the working fluid during the entire compressor cycle of operations completes the equation set and thereby enables it to be solved.
Certain assumptions had to be introduced to ensure efficient computation.These do not impose any limitations on the model nor cause significant departures from the real processes and are as follows:
– The fluid flow in the model is assumed to be quasi one-dimensional.
– Kinetic energy changes of the working fluid within the working chamber are negligible compared to internal energy changes.
– Gas or gas-liquid inflow to and outflow from the compressor ports is assumed to be isentropic.
– Leakage flow of the fluid through the clearances is assumed to be adiabatic.
3.1.1 Conservation Equations
For Control Volume and Auxiliary Relationships
The working chamber of a screw machine is the space within it that contains the working fluid. This is a typical example of an open thermodynamic system in which the mass flow varies with time. This, as well as the suction and discharge plenums, can be defined by a control volume for which the differential equations of the conservation laws for energy and mass are written. These are derived in Appendix B, using Reynolds Transport Theorem.
A feature of the model is the use of the non-steady flow energy equation to compute the thermodynamic and flow processes in a screw machine in terms of rotational angle or time and how these are affected by rotor profile modifications. Internal energy, rather than enthalpy, is then the derived variable. This is computationally more convenient than using enthalpy as the derived
Variable since, even in the case of real fluids, it may be derived, without reference to pressure. Computation is then carried out through a series of iterative cycles until the solution converges. Pressure, which is the desired output variable, can then be derived directly from it, together with the remaining required thermodynamic properties.
The following forms of the conservation equations have been employed in the model:
中文翻譯
螺桿式壓縮機
幾何的法線方向的螺旋，可以用來計算的坐標轉子螺旋和的從x和y的間隙加入如：
， ， (2.19)
其中分母D被給定為：
(2.20)
，服務來計算新的轉子端的平面的坐標， 和，得到的間隙角θ =鋅/ p和τ 。這些，現(xiàn)在的差額計算的橫向間隙δ0在它們之間，以及原來的轉子坐標和 。
如果以任何方式，轉子的改變它們的相對位置，該間隙的平側面上分布在轉子的一端，也可以減少到零的轉子葉片。在這樣的情況下，轉子的接觸將是令人望而卻步長側扁的檔案中，其中占主導地位的相對滑動轉子運動，如示于圖。 2.29。這表明，轉子扣押幾乎肯定會
如果轉子進入彼此接觸，發(fā)生在該區(qū)域。
圖.29。：吸力，中間轉子和轉子之間的間隙分布可能轉子接觸放電在放電
圖 2.30?？勺冮g隙分布應用到轉子
如下的間隙分布應該非均勻，以避免在轉子轉子之間的滑動運動的地方是硬轉子的接觸占主導地位。
另外，在圖2.30 ，清除率降低65微米，這是現(xiàn)在應用在轉子靠近轉子節(jié)圓的區(qū)域，而在其他區(qū)域是保持在85微米所做的那樣， 1992年由Edstroem 。正如在圖中可以看出的。 2.31，現(xiàn)在的情況，轉子的接觸是完全不同的。這是保持靠近轉子的節(jié)圓沿轉子的接觸帶，并完全避免在其他位置。因此，如果發(fā)生接觸，這將是一個滾動字符，而不是相結合的滾動和滑動，甚至是純滑動。這樣的接觸不會產(chǎn)生過多的熱量，因此可以保持一段較長時間，而不會損壞轉子直到接觸終止或使壓縮機停止。
2.6 為轉子制造的工具
本節(jié)描述了一代形成的工具螺桿壓縮機滾齒，銑床和磨床的基礎上信封資產(chǎn)負債程序。
2.6.1 滾齒機工具
螺桿壓縮機轉子和其形成的滾齒機工具相當于一個對相互嚙合的交錯軸斜齒輪與非平行不相交軸。他們的嚙合條件一般除了見附錄A。張裂縫形成刀具的面孔時，僅僅是一個螺旋齒輪的滾刀。
圖2.31。轉子之間的間隙分布：在抽吸，中期的轉子和可能轉子接觸放電在放電
每個齒被稱為作為一個線程， Colburne ，1987 。由于其自身的軸線不平行，它們之間的唯一的點接觸，而有行螺桿機轉子之間的接觸。需要滿足的嚙合在附錄A中給出的公式，導致轉子A “滾刀嚙合要求對于給定的轉子橫向坐標點X01和Y01和他們的第一次衍生。滾刀橫向坐標點X02和Y02計算出來的。這是足夠的，以獲得坐標軸向坐標z2的，直接計算，和R2是滾刀軸向平面坐標它定義滾刀的幾何形狀。
軸向坐標z2的，直接計算，和R2是滾刀軸向平面坐標它定義滾刀的幾何形狀。螺桿機轉子的橫向坐標，描述在前面的部分，被用來作為一個例子在這里產(chǎn)生滾刀坐標。轉子單元導致p1的有48.754毫米的主和- 58.504毫米的門轉子。單葉爐產(chǎn)生單位領導P2 ： 6.291毫米的主轉子和閘轉子- 6.291毫米為。相應的滾刀螺旋角ψ分別為85 ?和95 ? 。相同的轉子滾刀中心距離C = 110毫米和軸角Σ = 50 ?給出了兩個轉子。圖2.32中包含一個查看的爐灶。
反向計算的爐灶 - 螺桿轉子的改造，也給出了附錄A，允許確定轉子型線的橫向坐標這將在制造過程中得到的結果。這些可能與最初指定的進行比較，以確定影響
圖。 2.32。轉子制造業(yè)：滾齒刀具左，右銑刀
制造的錯誤，如完美的工具設置或工具，轉子變形后，最終的轉子型線。
如果在反向變換的目的，滾刀的縱向平面坐標R2和Z2和應給予。軸向坐標z2的使用計算τ = z2/p2 ，然后將其用于計算滾刀橫向的坐標：
， (2.21)
然后用這些作為給定的坐標，以產(chǎn)生一個嚙合判據(jù)的橫向平面上的坐標的“人造”轉子。
原來的轉子之間的比較和所制造的轉子給出圖。 2.33與它們之間的區(qū)別縮放100倍。二被認為是類型的錯誤。左邊的門轉子，30微米抵消在該轉子與該工具，和主旋翼之間的中心距離0.2 ?偏移刀具軸角Σ 。這個特殊的網(wǎng)格劃分方法的詳情給出由Stosic 1998。
圖。 2.33。制造缺陷
2.6.2 銑削和磨削工具
形成銑削和磨削工具也可以通過放置p2= 0產(chǎn)生嚙合方程，一般在附錄A中，然后按照本節(jié)的程序的?，F(xiàn)在的嚙合條件內容：
(2.22)
然而，在這種情況下，當一個人希望獲得螺桿轉子坐標從工具坐標，所施加的奇異性，不允許計算該工具的橫向平面坐標。主要的嚙合條件不能因此其應用。為此目的，另一個條件推導了扭轉銑刀到轉子的變換從該嚙合角τ計算方法是:
(2.23)
一旦獲得，τ后，將用來計算轉子坐標“制造”的過程。得到的轉子的坐標將包含所有制造不完善的地方，如不匹配的轉子 - 刀具中心的距離，在轉子中的誤差 - 工具軸角度，軸向移位的工具或工具變形在這個過程中，因為它們是輸入到計算處理。一個完整的帳戶這個有用的程序是Stosic于 1998年定義的。
2.6.3 量化的制造缺陷
轉子 - 這里使用的工具轉換生成銑削刀具輪廓。相反的步驟是用來計算“制成品”轉子。機架產(chǎn)生的Stosic的5-6 128毫米轉子的描述， 1997年a使用給定的概況：X （t）和Y（T） 。然后一個工具 - 轉子的變換是用來量化后，所生成的質量的影響的制造缺陷轉子型線。這兩種，線性和角度偏移量進行了考慮。
圖2.33轉子，主要制造與軸角度偏移0.5 ?門與中心的距離偏移量40微米原來的轉子由在左邊的虛線給出。在屏幕的右側，轉子制造的缺陷，主要與刀具軸向偏移為40μm，具有一定的工具主體變形而導致的柵極0.5 ?偏移量的相對運動的角度θ 。原來的轉子由下式給出的虛線。
3螺桿壓縮機性能的計算
螺桿式壓縮機的性能是受熱力學的相互影響和流體流動過程和機器的幾何形狀，從而可以只能由他們同時考慮可靠地計量。這可能是實現(xiàn)在一個或多個維度的數(shù)學建模。對于大多數(shù)應用，一個三維模型是足夠的，這充分說明。3-D模型比較復雜，而且只在這里介紹的輪廓。一個更詳細介紹了這個將在一個單獨的出版物。
3.1 一維數(shù)學模型
來描述的熱力學和流體的流動過程中所使用的算法的螺桿壓縮機的基礎上的數(shù)學模型。這個定義的瞬時工作腔的體積，其變化與旋轉角或時間，能量和質量的連續(xù)性的守恒方程施加，連同一組用于定義各種代數(shù)關系的吸入，壓縮和排出的工作有關的現(xiàn)象流體。這些形成了一套同時非線性微分方程不能得到解決封閉形式。
方程組的溶液通過數(shù)值進行龍格庫塔四階的方法，用適當?shù)某跏紬l件和邊界條件。
模型考慮了一些“現(xiàn)實生活”的影響，這可能顯著一個真正的壓縮機的性能產(chǎn)生影響。這使得它適合一個廣泛的應用范圍，包括以下內容：
-工作流體壓縮可以是任何氣體或液 - 氣混合物的方程的狀態(tài)和內能，焓關系是已知的，即任何理想的或真正的氣體或液 - 氣混合物的已知屬性。
- 的氣體之間的熱傳遞和壓縮機的模型占轉子或外殼的一種形式，雖然近似，再現(xiàn)一個良好的第一順序的精度水平的整體。
-工作介質通過間隙泄漏的模型占兩個轉子之間的和之間的轉子和固定的壓縮機部分。
- 他們的解決方案的過程方程和子程序是獨立的的那些定義壓縮機的幾何形狀。因此，該模型可以估計很容易適應任何幾何形狀的性能或類型容積式機器。
- 液體注入的影響，包括油，水，制冷劑或可以占的吸入，壓縮及排出階段期間。
- A組的子程序估計的熱力學性質和變化的工作流體的狀態(tài)在整個壓縮機的操作循環(huán)期間完成方程組，從而使其能夠解決。
若干假設被引入，以確保有效的計算。這些不施加任何限制的模式，也沒有造成顯著偏離真實進程如下：
-在模型中的流體流被假定為準一維。
-內的工作腔的工作流體的動能變化是相比微不足道內部能量的變化。
-假定氣體或氣 - 液的流入和流出，從壓縮機端口等熵。
-流體通過間隙泄漏流被假定為絕熱。
3.1.1 守恒方程
控制音量和輔助的關系
工作腔的螺桿機內的空間，它包含工作流體。這是一個開放的熱力學系統(tǒng)的一個典型例子在其中的質量流量隨時間變化。這一點，以及吸入管和排出增壓室，可以定義由一個控制體積的差動方程，能量和質量守恒定律被寫入。這些都是在附錄B中派生，使用雷諾運輸定理。
一個模型的特征是使用的非定常流動能量方程計算的螺桿機方面的熱力學和流動過程旋轉角度或時間，以及如何，這些受影響的的轉子配置文件修改。內部的能量，而不是焓，然后導出變量。這是計算更方便，比焓派生變量，因為，即使在實際流體的情況下，它可以導出，而不參考的壓力。計算，然后通過一系列的迭代進行周期，直到該溶液收斂。的壓力，這是所需的輸出變量，然后，可以直接從它派生，連同其余必需的熱力學性質。
已使用了以下形式的守恒方程模型：
15