JULIO.SEPTIEMBRE 2007 NUM. 27 - CIVILTEC
JULIO.SEPTIEMBRE 2007 NUM. 27 - CIVILTEC
JULIO.SEPTIEMBRE 2007 NUM. 27 - CIVILTEC
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