08.02.2013 Views

JULIO.SEPTIEMBRE 2007 NUM. 27 - CIVILTEC

JULIO.SEPTIEMBRE 2007 NUM. 27 - CIVILTEC

JULIO.SEPTIEMBRE 2007 NUM. 27 - CIVILTEC

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

formaS de reSolver<br />

laS eCuaCioneS<br />

1. Método de Newton: a. Supone unos valores iniciales para las cargas<br />

desconocidas, H A , H B , H C . b. Plantea un sistema de ecuaciones lineales<br />

en donde los coeficientes son las derivadas parciales de cada una de<br />

las ecuaciones con respecto a las cargas desconocidas y las incógnitas<br />

de este sistema lineal, son las modificaciones a las cargas supuestas.<br />

Los términos independientes son las mismas ecuaciones (su diferencia<br />

es cero) con signo negativo, c. Sistema lineal de Newton:<br />

(∂F(A)/∂HA)XA + (∂F(A)/∂HB)XB + (∂F(A)/∂HC)XC = -F(A)<br />

(∂F(B)/∂HA)XA + (∂F(B)/∂HB)XB + (∂F(B)/∂HC)XC = -F(B)<br />

(∂F(C)/∂HA)XA + (∂F(C)/∂HB)XB + (∂F(C)/∂HC)XC = -F(C)<br />

o en otras palabras<br />

IJI IXI = -IFI Resolver por iteraciones<br />

2. Método orientado a la computadora.<br />

El sistema de ecuaciones:<br />

Ya(HD-HA) 0.54 + Yb(HB-HA) 0.54 + Yc(HC-HA) 0.54 + QA = 0<br />

Yb(HA-HB) 0.54 + Yd(HD-HB) 0.54 + Ye(HC-HB) 0.54 + Yf(HE-HB) 0.54 + QB = 0<br />

Yc(HA-HC) 0.54 + Ye(HB-HC) 0.54 + Yg(HE-HC) 0.54 + QC = 0<br />

Al igual que el método de Newton y Cross, se supondrán valores de<br />

H A , H B , H C , pero con los valores supuestos se calcularán constantes<br />

temporales para el sistema de ecuaciones no lineales. Estas serán de<br />

la forma:<br />

Ka = Ya / (abs[HD-HA] 0.46 )<br />

Kb = Yb / (abs[HB-HA] 0.46 )<br />

Kc = Yc / (abs[HC-HA] 0.46 )<br />

Con estas constantes temporales<br />

se tendrá un sistema lineal:<br />

Ka(HD-HA) + Kb(HB-HA) + Kc(HC-HA) + QA = 0<br />

Kb(HA-HB) + Kb(HA-HB) + Ke(HC-HB) + Kf(HE-HB) + QB = 0<br />

Kc(HA-HC) + Ke(HB-HC) + Kg(HE-HC) + QC = 0<br />

Estas ecuaciones se pueden simplificar:<br />

-(Ka+Kb+Kc)HA + KbHb + KcHC = -(QA + KaHD)<br />

KbHA - (Kd+Kb+Ke+Kf) HB + KeHC = -(QB+KfHE+KdHD)<br />

KCHA + KeHB - (Kc+Ke+Kg)HC = -(QC+KgHE)<br />

Esto es un sistema lineal simétrico con diagonal principal dominante.<br />

Los posibles métodos de solución: Gauss sustitución hacia atrás,<br />

Gauss Jordan, Gauss Seydel, Cholesky.<br />

SeCuela:<br />

1. Suponer cargas en los nudos de carga desconocida.<br />

2. Cálculos de constantes Ki temporales.<br />

3. Resolver el sistema de ecuaciones<br />

4. Conocidos valores de (H) despejadas y más aproximadas se regresa<br />

al paso 2, hasta obtener la precisión deseada.<br />

nota: Al manejar el denominador de los K temporales en valor absoluto,<br />

se arreglan los signos del sistema lineal.<br />

referenCiaS<br />

Treviño Garza Rodolfo, “Programa de cálculo electrónico para redes de<br />

abastecimiento de agua” Apuntes de la clase, ITESM, Departamento<br />

de Ingeniería Civil.<br />

a rt Í C u l o t É C n i C o<br />

<strong>JULIO</strong>.SEP TIEMBRE <strong>27</strong>.20 07 <strong>CIVILTEC</strong> 9

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!