Problemas EDoS SRW y PR CO2

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 8

PROBLEMA.

Calcular la compresibilidades, fugacidad del líquido y vapor del CO2 (Dióxido de Carbono) a 300
K y 6.7131 MPa usando la ecuación de Soave Redlich-Kwong y Peng-Robinson.
a) Calcular la fugacidad del CO2 (Dióxido de Carbono) a 300 K y 6.7131 MPa usando la ecuación de
Soave Redlich-Kwong.
La ecuación de Soave - Redlich - Kwong se define como:
𝑅𝑇 𝑎(𝑇)
𝑃= − (1)
𝑣 − 𝑏 𝑣(𝑣 + 𝑏)
Dicha ecuación en función del factor de compresibilidad se escribe como:
𝑓(𝑍) = 𝑍 3 − 𝑍 2 + ( 𝐴 − 𝐵 − 𝐵 2 )𝑍 − 𝐴𝐵 (2)
Derivando la ecuación 2 con respecto a Z,
𝑓´(𝑍) = 3𝑍 2 − 2𝑍 + ( 𝐴 − 𝐵 − 𝐵 2 ) (3)
Donde:
2
𝑎𝑃 [1 + (0.480 + 1.574𝜔 − 0.176𝜔2 )(1 − 𝑇𝑟 1⁄2 )] 𝑇𝑐 2 𝑃
𝐴 = 2 2 = 0.42748 (4)
𝑅 𝑇 𝑃𝑐𝑇 2
𝑏𝑃 𝑇𝑐𝑃
𝐵= = 0.08664 (5)
𝑅𝑇 𝑃𝑐𝑇
Partiendo de la definición de fugacidad, la expresión matemática para determinar la fugacidad de la ecuación de
Soave-Redlich-Kwong es:
𝑓 𝐴 𝐵
ln ( ) = (𝑍 − 1) − ln(𝑍 − 𝐵) − ln (1 + ) (6)
𝑃 𝐵 𝑍
Para poder encontrar el valor correcto del factor de compresibilidad con mayor rapidez podemos emplear un
método numérico, en este caso utilizamos el método de Newton-Raphson el cual está definido como:
𝑛+1 𝑛
𝑓( 𝑥 𝑛 )
𝑋 =𝑋 − ′ 𝑛 (7)
𝑓 (𝑥 )
Escribiendo la ec.7 en función de Z tenemos como resultado:
𝑛+1 𝑛
𝑓( 𝑍 𝑛 )
𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍 =𝑍 − ′ 𝑛 (8)
𝑓 (𝑍 )
Reescribiendo la ec.8 utilizando la ec.3 y ec.3,
𝑍 3 − 𝑍 2 + ( 𝐴 − 𝐵 − 𝐵 2 )𝑍 − 𝐴𝐵
𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍𝑠𝑢𝑝 − (9)
3𝑍 2 − 2𝑍 + ( 𝐴 − 𝐵 − 𝐵 2 )
Para comenzar la solution del problema debemos definir los datos proporcionados. Estos datos se encuantran en
la tabla 1.
Tabla 1. Datos proporcionados en el problema.
Datos Proporcionados
𝑻 [𝑲] 300
𝑷 [𝑴𝒑𝒂] 6.7131
Sustancia CO2

Al analizar los parametros A y B notamos que necesitamos conocer los valores de la temperatura crítica, la presión
crítica y el factor acéntrico para el CO2, en la table 2 se reportan los valores. Para encontrar los valores de estas
variables consultamos la base de datos del NIST.
Tabla 2. Datos requeridos para la solución del problema.
Datos Requeridos
𝑻𝑪 [𝑲] 304.1282
𝑷𝑪 [𝑴𝒑𝒂] 7.3773
𝝎 0.22394
Debido a que se desconoce el valor de 𝑍 podemos suponer un valor (𝑍𝑠𝑢𝑝 ) para iniciar el cálculo. De acuerdo a
(Smith Van Ness & Abbott 2007); en cálculos del factor de comprensibilidad de la fase vapor y fase líquida, el
valor inicial sugerido de 𝑍𝑠𝑢𝑝 , puede ser:
𝑍𝑠𝑢𝑝 = 1 (𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟)
𝑍𝑠𝑢𝑝 = 0 (𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎)
Procedimiento iterativo (fase vapor):
En la tabla 3, se muestran todos los resultados de cada iteración para determinar 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 .
Tabla 3. Resultados de las iteraciones para determinar 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 .
Para la fase vapor
No. Iteración 𝒁𝒔𝒖𝒑 𝒇(𝒁) 𝒇´(𝒁) 𝒁𝒏𝒖𝒆𝒗𝒂 𝑬𝒓𝒓𝒐𝒓
0 1.000000 0.285646 1.317957 0.783266 2.44E-01
1 0.783266 0.083766 0.591942 0.641755 1.99E-01
2 0.641755 0.024196 0.269995 0.552137 1.50E-01
3 0.552137 0.006711 0.128249 0.499806 9.96E-02
4 0.499806 0.001654 0.067763 0.475392 5.01E-02
5 0.475392 0.000283 0.045166 0.469124 1.33E-02
6 0.469124 0.000016 0.039941 0.468711 8.81E-04
7 0.468711 0.000000 0.039605 0.468709 3.74E-06
8 0.468709 0.000000 0.039604 0.468709 6.72E-11

𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 = 0.468709


El método termina en la iteración número ocho con un valor de 0.468709 y un error de 6.72E-11. Con este valor
del factor de compresibilidad podemos calcular:
𝑃𝑣 𝑍𝑅𝑇
𝑍= ⟹𝑣=
𝑅𝑇 𝑃
−6 𝑀𝑃𝑎𝑚3
(0.488709) (8.314𝑥10 ) (300𝐾) 𝑚3
𝑚𝑜𝑙 𝐾
𝑣= = 0.000174145
6.7131 𝑀𝑃𝑎 𝑚𝑜𝑙

𝑚3
0.000174145 3
v=( 𝑚𝑜𝑙 ) (1000 g) = 0.003956993 m
0095 g 1 Kg Kg
44.
mol
El valor verdadero (encontrado en la base de datos del NIST) es 0.003723 m3 /Kg. sto genera un error del 6.2865
% con respecto a la realidad. Por lo tanto, la ecuación de Soave - Redlich - Kwong promete resultados muy
próximos a la realidad. Utilizando la ec.6., se calcula fugacidad de la fase vapor:
𝑓 𝐴 𝐵
ln ( ) = (𝑍 − 1) − ln(𝑍 − 𝐵) − ln (1 + )
𝑃 𝐵 𝑍
Donde : 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 = 0.468709
𝐴 = 0.404269355
𝐵 = 0.079924431
Sustituyendo los datos, se obtienen:
𝑓 0.404269355 0.079924431
ln ( ) = (0.468709 − 1) − ln(0.468709 − 0.079924431) − ln (1 + )
𝑃 0.079924431 0.468709
𝑓
ln ( ) = −0.382957973
𝑃
Resolviendo para 𝑓,
𝑓 𝑣 = 𝑃𝑒 −0.382957973 = (6.7131 𝑀𝑃𝑎)𝑒 −0.382957973
∴ 𝑓 𝑣 = 4.577270537 𝑀𝑃𝑎
Procedimiento iterativo (fase líquida):
1. Determinar los valores de A y B, utilizando la ec.4 y ec.5, respectivamente.
2
𝑎𝑃 [1 + (0.480 + 1.574𝜔 − 0.176𝜔2 )(1 − 𝑇𝑟 1⁄2 )] 𝑇𝑐 2 𝑃
𝐴 = 2 2 = 0.42748
𝑅 𝑇 𝑃𝑐𝑇 2
𝑇 300 𝐾
𝑇𝑟 = = = 0.9864
𝑇𝐶 304.1282 𝐾
2
[1 + (0.480 + 1.574(0.22394) − 0.176(0.22394)2 )(1 − (0.9864)1⁄2 )] (304.1282 𝐾)2 (6.7131 𝑀𝑃𝑎)
𝐴 = 0.42748
(7.3773 𝑀𝑃𝑎)(300 𝐾)2

𝐴 = 0.404269355

𝑇𝑐𝑃
𝐵 = 0.08664
𝑃𝑐𝑇
(304.1282 𝐾)(6.7131 𝑀𝑃𝑎)
𝐵 = 0.08664 = 0.079924431
(7.3773 𝑀𝑃𝑎)(300 𝐾)

En la tabla 4, se muestran los resultados de cada iteración para determinar 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 .

Tabla 4. Resultados de las iteraciones para determinar 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 .

Para la fase liquida


No. Iteración 𝒁𝒔𝒖𝒑 𝒇(𝒁) 𝒇´(𝒁) 𝒁𝒏𝒖𝒆𝒗𝒂 𝑬𝒓𝒓𝒐𝒓
0 0.000000 -0.032311 0.317957 0.101621 -
1 0.101621 -0.009277 0.145696 0.165297 4.86E-01
2 0.165297 -0.002560 0.069333 0.202225 2.02E-01
3 0.202225 -0.000637 0.036192 0.219829 8.35E-02
4 0.219829 -0.000116 0.023274 0.224832 2.25E-02
5 0.224832 -0.000008 0.019942 0.225253 1.87E-03
6 0.225253 0.000000 0.019668 0.225256 1.30E-05
7 0.225256 0.000000 0.019666 0.225256 6.29E-10

𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 = 0.225256

El método termina en la iteración número siete con un valor de 0.225256 y un error de 6.29E-10. Con este valor
del factor de compresibilidad podemos calcular:
𝑃𝑣 𝑍𝑅𝑇
𝑍= ⟹𝑣=
𝑅𝑇 𝑃
𝑀𝑃𝑎𝑚3
(0.225256) (8.314𝑥10−6 ) (300𝐾) 𝑚3
𝑚𝑜𝑙 𝐾
𝑣= = 8.36919𝐸 − 05
6.7131 𝑀𝑃𝑎 𝑚𝑜𝑙
𝑚3
8.36919𝐸 − 05 3
v=( 𝑚𝑜𝑙 ) (1000 g) = 0.001901679 m
0095 g 1 Kg Kg
44.
mol
El valor verdadero (encontrado en la base de datos del NIST) es 0.0014722 m3 /Kg. sto genera un error del
22.5841 % con respecto a la realidad. Por lo tanto, la ecuación de Soave - Redlich - Kwong no promete resultados
muy próximos a la realidad en la fase líquida.
Utilizando la ec.6., se calcula fugacidad de la fase líquida:
𝑓 𝐴 𝐵
ln ( ) = (𝑍 − 1) − ln(𝑍 − 𝐵) − ln (1 + )
𝑃 𝐵 𝑍
Donde : 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 = 0.225256

𝐴 = 0.404269355
𝐵 = 0.079924431
Sustituyendo los datos, se obtienen:
𝑓 0.404269355 0.079924431
ln ( ) = (0.225256 − 1) − ln(0.225256 − 0.079924431) − ln (1 + )
𝑃 0.079924431 0.225256
𝑓
ln ( ) = −0.381991381
𝑃
Resolviendo para 𝑓,

𝑓 𝑙 = 𝑃𝑒 −0.381991381 = (6.7131 𝑀𝑃𝑎)𝑒 −0.381991381

∴ 𝑓 𝑙 = 4.581697032 𝑀𝑃𝑎
b) Calcular la fugacidad del CO2 (Dióxido de Carbono) a 300 K y 6.7131 MPa usando la ecuación de
Peng - Robinson.
La ecuación de Peng - Robinson se define como:
𝑅𝑇 𝑎
𝑃= − (1)
(𝑣 − 𝑏) 𝑣(𝑣 + 𝑏) + 𝑏(𝑣 − 𝑏)
Dicha ecuación en función del factor de compresibilidad se escribe como:
𝑓(𝑍) = 𝑍 3 − (1 − 𝐵)𝑍 2 + (𝐴 − 2𝐵 − 3𝐵 2 )𝑍 − (𝐴𝐵 − 𝐵 2 − 𝐵 3 ) (2)
Derivando la ecuación 2 con respecto a Z,
𝑓´(𝑍) = 3𝑍 2 − 2(1 − 𝐵)𝑍 + (𝐴 − 2𝐵 − 3𝐵 2 ) (3)
Donde:
2
𝑎𝑃 [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 )(1 − 𝑇𝑟 1⁄2 )] 𝑇𝑐 2 𝑃
𝐴 = 2 2 = 0.45725 (4)
𝑅 𝑇 𝑃𝑐𝑇 2

𝑏𝑃 𝑇𝑐𝑃
𝐵= = 0.07780 (5)
𝑅𝑇 𝑃𝑐𝑇
Partiendo de la definición de fugacidad, la expresión matemática para determinar la fugacidad de la ecuación de
Peng - Robinson es:
𝑓 𝐴 𝑍 + 2.414𝐵
ln ( ) = (𝑍 − 1) − ln(𝑍 − 𝐵) − ln ( ) (6)
𝑃 2√2𝐵 𝑍 − 0.414𝐵
Para poder encontrar el valor correcto del factor de compresibilidad con mayor rapidez podemos emplear un
método numérico, en este caso utilizamos el método de Newton-Raphson el cual está definido como:
𝑛+1 𝑛
𝑓( 𝑥 𝑛 )
𝑋 =𝑋 − ′ 𝑛 (7)
𝑓 (𝑥 )
Escribiendo la ec.7 en función de 𝑍 tenemos como resultado:
𝑛+1 𝑛
𝑓( 𝑍 𝑛 )
𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍 =𝑍 − ′ 𝑛 (8)
𝑓 (𝑍 )
Reescribiendo la ec.8 utilizando la ec.3 y ec.3,
𝑍 3 − (1 − 𝐵)𝑍 2 + (𝐴 − 2𝐵 − 3𝐵 2 )𝑍 − (𝐴𝐵 − 𝐵 2 − 𝐵 3 )
𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍𝑠𝑢𝑝 − (9)
3𝑍 2 − 2(1 − 𝐵)𝑍 + (𝐴 − 2𝐵 − 3𝐵 2 )
Para comenzar la solution del problema debemos definir los datos proporcionados. Estos datos se encuantran en
la tabla 1.
Tabla 1. Datos proporcionados en el problema.
Datos Proporcionados
𝑻 [𝑲] 300
𝑷 [𝑴𝒑𝒂] 6.7131
Sustancia CO2

Tabla 2. Datos requeridos para la solución del problema.


Datos Requeridos
𝑻𝑪 [𝑲] 304.1282
𝑷𝑪 [𝑴𝒑𝒂] 7.3773
𝝎 0.22394

De acuerdo a (Smith Van Ness & Abbott 2007); en cálculos del factor de comprensibilidad de la fase vapor y fase
líquida, el valor inicial sugerido de 𝑍𝑠𝑢𝑝 , puede ser:
𝑍𝑠𝑢𝑝 = 1 (𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟)
𝑍𝑠𝑢𝑝 = 0 (𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎)
Procedimiento iterativo (fase vapor):
1. Determinar los valores de A y B, utilizando la ec.4 y ec.5, respectivamente.
2
𝑎𝑃 [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 )(1 − 𝑇𝑟 1⁄2 )] 𝑇𝑐 2 𝑃
𝐴= = 0.45725
𝑅2𝑇 2 𝑃𝑐𝑇 2
𝑇 300 𝐾
𝑇𝑟 = = = 0.9864
𝑇𝐶 304.1282 𝐾
2
[1 + (0.37464 + 1.54226(0.22394) − 0.26992(0.22394)2 )(1 − (0.9864)1⁄2 )] (304.1282 𝐾)2 (6.7131 𝑀𝑃𝑎)
𝐴 = 0.45725
(7.3773 𝑀𝑃𝑎)(300 𝐾)2

𝐴 = 0.431736928

𝑇𝑐𝑃
𝐵 = 0.07780
𝑃𝑐𝑇
(304.1282 𝐾)(6.7131 𝑀𝑃𝑎)
𝐵 = 0.07780
(7.3773 𝑀𝑃𝑎)(300 𝐾)
𝐵 = 0.071769630
Tabla 5. Resultados de las iteraciones para determinar 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 .

Para la fase vapor


No. Iteración 𝒁𝒔𝒖𝒑 𝒇(𝒁) 𝒇´(𝒁) 𝒁𝒏𝒖𝒆𝒗𝒂 𝑬𝒓𝒓𝒐𝒓
0 1.000000 0.319050 1.421435 0.775544 2.54E-01
1 0.775544 0.094225 0.642534 0.628899 2.10E-01
2 0.628899 0.027674 0.296911 0.535692 1.60E-01
3 0.535692 0.007997 0.144302 0.480272 1.09E-01
4 0.480272 0.002200 0.078274 0.452164 6.03E-02
5 0.452164 0.000528 0.051828 0.441984 2.28E-02
6 0.441984 0.000096 0.043420 0.439779 5.00E-03
7 0.439779 0.000013 0.041680 0.439460 7.25E-04
8 0.439460 0.000002 0.041431 0.439420 9.23E-05
9 0.439420 0.000000 0.041400 0.439415 1.15E-05
10 0.439415 0.000000 0.041396 0.439414 1.43E-06
11 0.439414 0.000000 0.041395 0.439414 1.79E-07

𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 = 0.439414

El método termina en la iteración número once con un valor de 0.439414 y un error de 1.79E-07. Con este valor
del factor de compresibilidad podemos calcular:
𝑃𝑣 𝑍𝑅𝑇
𝑍= ⟹𝑣=
𝑅𝑇 𝑃
𝑀𝑃𝑎𝑚3
(0.439414 ) (8.314𝑥10−6 ) (300𝐾) 𝑚3
𝑚𝑜𝑙 𝐾
𝑣= = 0.000163261
6.7131 𝑀𝑃𝑎 𝑚𝑜𝑙
𝑚3
0.000163261 3
v=( 𝑚𝑜𝑙 ) (1000 g) = 0.003709674 m
0095 g 1 Kg Kg
44.
mol
El valor verdadero (encontrado en la base de datos del NIST) es 0.003723 m3 /Kg. sto genera un error del 0.3592
% con respecto a la realidad. Por lo tanto, la ecuación de Peng – Robinson promete resultados muy próximos a la
realidad.
Utilizando la ec.6., se calcula fugacidad de la fase vapor:
𝑓 𝐴 𝑍 + 2.414𝐵
ln ( ) = (𝑍 − 1) − ln(𝑍 − 𝐵) − ln ( )
𝑃 2√2𝐵 𝑍 − 0.414𝐵
Donde : 𝑍𝑓𝑎𝑠𝑒 𝑣𝑎𝑝𝑜𝑟 = 0.439414

𝐴 = 0.431736928
𝐵 = 0.071769630
Sustituyendo los datos, se obtienen:
𝑓 0.431736928 0.439414 + 2.414(0.071769630)
ln ( ) = (0.439414 − 1) − ln(0.439414 − 0.071769630) − ln [ ]
𝑃 2√2(0.071769630) 0.439414 − 0.414(0.071769630)

𝑓
ln ( ) = −0.41576591
𝑃
Resolviendo para 𝑓,

𝑓 𝑣 = 𝑃𝑒 −0.41576591 = (6.7131 𝑀𝑃𝑎)𝑒 −0.41576591

∴ 𝑓 𝑣 = 4.429536414 𝑀𝑃𝑎

Procedimiento iterativo (fase líquida):


1. Determinar los valores de A y B, utilizando la ec.4 y ec.5, respectivamente.
2
𝑎𝑃 [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 )(1 − 𝑇𝑟 1⁄2 )] 𝑇𝑐 2 𝑃
𝐴 = 2 2 = 0.45725
𝑅 𝑇 𝑃𝑐𝑇 2
𝑇 300 𝐾
𝑇𝑟 = = = 0.9864
𝑇𝐶 304.1282 𝐾
2
[1 + (0.37464 + 1.54226(0.22394) − 0.26992(0.22394)2 )(1 − (0.9864)1⁄2 )] (304.1282 𝐾)2 (6.7131 𝑀𝑃𝑎)
𝐴 = 0.45725
(7.3773 𝑀𝑃𝑎)(300 𝐾)2

𝐴 = 0.431736928

𝑇𝑐𝑃
𝐵 = 0.07780
𝑃𝑐𝑇
(304.1282 𝐾)(6.7131 𝑀𝑃𝑎)
𝐵 = 0.07780
(7.3773 𝑀𝑃𝑎)(300 𝐾)
𝐵 = 0.071769630
Tabla 6. Resultados de las iteraciones para determinar 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 .

Para la fase liquida


No. Iteración 𝒁𝒔𝒖𝒑 𝒇(𝒁) 𝒇´(𝒁) 𝒁𝒏𝒖𝒆𝒗𝒂 𝑬𝒓𝒓𝒐𝒓
0 0.000000 -0.025465 0.272745 0.093366 -
1 0.093366 -0.007278 0.125567 0.151324 3.83E-01
2 0.151324 -0.001983 0.060515 0.184085 1.78E-01
3 0.184085 -0.000474 0.032660 0.198593 7.31E-02
4 0.198593 -0.000076 0.022382 0.201992 1.68E-02
5 0.201992 -0.000004 0.020157 0.202181 9.33E-04
6 0.202181 0.000000 0.020035 0.202182 2.83E-06
7 0.202182 0.000000 0.020035 0.202182 2.60E-11

𝑍𝑛𝑢𝑒𝑣𝑎 = 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 = 0.202182


El método termina en la iteración número siete con un valor de 0.202182 y un error de 2.60E-11. Con este valor
del factor de compresibilidad podemos calcular:
𝑃𝑣 𝑍𝑅𝑇
𝑍= ⟹𝑣=
𝑅𝑇 𝑃
−6 𝑀𝑃𝑎𝑚3
(0.202182) (8.314𝑥10 ) (300𝐾) 𝑚3
𝑚𝑜𝑙 𝐾
𝑣= = 7.5119𝐸 − 05
6.7131 𝑀𝑃𝑎 𝑚𝑜𝑙
𝑚3
7.5119𝐸 − 05 1000 g m3
v=( 𝑚𝑜𝑙 )( ) = 0.001706882
0095 g 1 Kg Kg
44.
mol

El valor verdadero (encontrado en la base de datos del NIST) es 0.0014722 m3 /Kg. sto genera un error del 13.7491
% con respecto a la realidad. Por lo tanto, la ecuación de Peng - Robinson no promete resultados muy próximos
a la realidad en la fase líquida.
Utilizando la ec.6., se calcula fugacidad de la fase vapor:
𝑓 𝐴 𝑍 + 2.414𝐵
ln ( ) = (𝑍 − 1) − ln(𝑍 − 𝐵) − ln ( )
𝑃 2√2𝐵 𝑍 − 0.414𝐵
Donde : 𝑍𝑓𝑎𝑠𝑒 𝑙í𝑞𝑢𝑖𝑑𝑎 = 0.202182
𝐴 = 0.431736928
𝐵 = 0.071769630
Sustituyendo los datos, se obtienen:
𝑓 0.431736928 0.202182 + 2.414(0.071769630)
ln ( ) = (0.202182 − 1) − ln(0.202182 − 0.071769630) − ln [ ]
𝑃 2√2(0.071769630) 0.202182 − 0.414(0.071769630)
𝑓
ln ( ) = −0.415148725
𝑃

Resolviendo para 𝑓,
𝑓 𝑙 = 𝑃𝑒 −0.415148725 = (6.7131 𝑀𝑃𝑎)𝑒 −0.415148725
∴ 𝑓 𝑙 = 4.432271101 𝑀𝑃𝑎

También podría gustarte