FUNDACIONES ELASTICAS

EN MEDIO CONTINUO

Gustavo Arias Albán

FUNDACIONES ELASTICAS
EN MEDIO CONTINUO

PRIMERA EDICION

Gustavo Arias Albán

Depósito Legal

Diagramación: Ana Luisa Ces

Abaco Arte C.A.

Impreso en Venezuela
por

Caracas 2009

II
 PRESENTACION

El desarrollo de la teoría y solución del problema de fundaciones elásticas según el

planteamiento de M. Hetenyi [4], mostrado en esta publicación, es autoría del Profesor Alejandro
Segovia Gallegos [1] de la Universidad Central de Quito Ecuador y de la Escuela Politécnica del
Ejercito Ecuatoriano ESPE. Los apuntes relacionados con el problema de la flexión de vigas,
resuelto con su método de Rigideces Sucesivas [2] , [9] y editados por la ESPE como
poligrafiados, fueron generosamente cedidos al suscrito en la década del 70. Posteriormente,
recibimos manuscritos relacionados con la teoría de la torsión [1] que complementaban el
problema de retículos de fundación y facilitaban la solución integral.

Lamentablemente, la documentación citada no se encuentra disponible en las fuentes

originales, y la posibilidad de consulta, ampliación, aplicación o propuesta de soluciones alternas,
es prácticamente nula.

El suscrito, ex-alumno del Profesor Segovia Gallegos, siente la inquietud de rescatar tan
valioso trabajo de un tema muy poco tratado en la literatura especializada, mediante un resumen
convenientemente detallado de la solución de las ecuaciones diferenciales de la flexión y de la
torsión, conservando la nomenclatura, presentación de fórmulas, unidades y ejemplos del
original, añadiendo la aplicación del análisis matricial en la fase donde es aplicable.

La solución de las ecuaciones diferenciales se plantea con funciones trigonométricas e

hiperbólicas y la convención de signos de resistencia de materiales. La base del procedimiento
está en el cálculo de las constantes de integración de cada tramo que dependen de las rotaciones y
desplazamientos de las juntas que, a su vez, se determinan con el método de los desplazamientos
empleando las rigideces de barra en suelo elástico deducidas por el autor [1] y la convención de
signos de estructuras continuas. Las características en varias secciones de un tramo se calculan
con las constantes de integración mencionadas y la convención de signos de resistencia de
materiales.

Se anexan dos programas de computación: FUNDAVIGA para análisis de vigas

continuas y FUNDARETIC para análisis de retículos, escritos en Rocky Mountain HP Basic de
Hewlett Packard aplicable al HT Basic para Windows de TRANSERA [8], con los cuales se se
incrementa la utilidad de lo expuesto.

No podemos dejar de mencionar los trabajos sobre el tema, del Profesor Andrés Reti [11]
de 1951 limitado a vigas continuas y el del Profesor Celso Fortoul P. [12] de 1973 de amplia
aplicación al computador.

Gustavo Arias Albán

III
 RECONOCIMIENTOS

A la memoria del Profesor Alejandro Segovia Gallegos.

A la Empresa OTIP C.A. por el auspicio de esta versión.

A los colegas y amigos que estimularon la escritura de este resumen.

IV
 CONTENIDO

1. VIGAS DE FUNDACION ELASTICA

1.1. Ecuaciones Diferenciales de la Flexión 1

1.2. Clasificación de Vigas 2

1.3. Vigas de Longitud Intermedia 3

1.3.1. Solución Homogénea 3
1.3.2. Constantes de Integración 4

1.4. Vigas de Longitud Indefinida 4

1.4.1. Caso de Carga Concentrada 5
1.4.2. Caso de Momento Aislado 6
1.4.3. Caso de Carga Uniforme Parcial 7

1.5. Solución Particular en Vigas de Longitud Intermedia 9

1.5.1. Cargas Concentradas y Uniformes Parciales. Ejemplo Nº 1 9
1.5.2. Carga Uniforme Total 11

1.6 Estructuras Continuas 11

1.6.1 Ecuación Matricial de una Viga 11
1.6.2 Vigas Continuas. Ejemplo Nº 2 12
1.6.3. Vigas sobre Apoyos Fijos. Ejemplo Nº 3 15

1.7 Programa FUNDAVIGA 17

1.7.1 Codificación de Datos 18
Programa. 19
Ejemplo Nº 4 25

2. RETICULOS DE FUNDACION ELASTICA

2.1. Ecuaciones Diferenciales de la Torsión 27

2.2. Vigas de Longitud Intermedia 28

2.2.1. Solución Homogénea 28
2.2.2. Constantes de Integración 28

2.3. Vigas de Longitud Indefinida 29

2.3.1. Momento Torsor Aislado 29

V
 2.4 Solución Particular en Vigas de Longitud Intermedia 30
2.4.1. Momentos de Torsión Aislados. Ejemplo Nº 5 30
2.4.2. Momento Torsor Uniforme en la Viga 31

2.5. Ecuación Matricial de una Viga 32

2.6. Vigas Esviadas 33

2.7 Ecuación Matricial de un Retículo 34

2.8 Momento Máximo de Tramo 35

2.9 Programa FUNDARETIC 36

2.9.1 Preparación de Datos 36
2.9.2 Codificación de Datos 37
Programa. 38
Ejemplo Nº 6 48
Ejemplo Nº 7 51
2.9.3 Salida Resumida. Ejemplo Nº 8 54

3. ANEXO A 55

Propiedades Mecánicas del Concreto

Consideraciones sobre el Módulo de Balasto
Propiedades Geométricas de las Secciones

4. ANEXO B 56

Convención de signos de la Estructura y de la Fundació

Casos de Combinación de Solicitaciones

5. ANEXO C 57

Observaciones al Estudio de Suelos

6. BIBLIOGRAFIA. 58

VI
 1. VIGAS DE FUNDACION ELASTICA

1.1 Ecuaciones Diferenciales de la Flexión

En la Fig 1.1 se esquematiza el comportamiento de una viga prismática de pequeña altura

solicitada por una carga, sin deformación por corte. En las Figs 1.1(a,b,c) constan los vectores
positivos según la convención usada en resistencia de materiales: desplazamiento (w) y carga (pt)
hacia abajo, giro de la elástica (φ) en sentido horario, momento (m) positivo con tracción en la
fibra inferior y corte positivo (V) si el sector a la izquierda de la sección considerada se desplaza
hacia arriba.

(a) Elástica (b) Momentos (c) Cortes (d) Viga sobre Suelo Elástico

Fig 1.1 Convención de Signos según Resistencia de Materiales

Las ecuaciones diferenciales para este tipo de viga son:

dw

dx

d 2w
m   EI 2
dx
(1.1a, d)
3
d w
V   EI
dx3

d 4w
pt  EI
dx 4

Si la viga descansa sobre un medio elástico, Fig 1.1(d), la barra se deforma produciendo
una reacción (q) del suelo que actúa verticalmente en dirección opuesta a la deformación y que,
suponemos, es proporcional al desplazamiento de la viga en ese punto.

Fundaciones Elásticas 1
 Esta hipótesis implica la condición de que el soporte es elástico y que su elasticidad puede
caracterizarse por la carga unitaria que produce una deformación unitaria. Esta constante del
medio soportante (kgf/cm3) se denomina módulo de la fundación o coeficiente de balasto (β).
La reacción del suelo para la viga de ancho constante (r), será:

q   .r.w (1.2)

la carga efectiva: pt  p   .r.w

d 4w
y la (1.1d) queda: EI   .r.w  p (1.3)
dx 4

Ecuación diferencial lineal con coeficientes constantes cuya solución consta de dos partes:
la particular y la homogénea. La primera debe satisfacer la (1.3) sin necesidad de cumplir
condiciones de vínculo de los extremos de la pieza, la segunda satisfará la ecuación homogénea
que se obtiene anulando el miembro derecho de la (1.3) y usando el factor 4, por conveniencia de
la solución, o sea:
4 EI d 4 w
 4w  0 (1.4)
 .r dx 4

introduciendo el concepto de longitud elástica (λ) definida como:

1/ 4
 4 EI 
   (1.5)
  .r 
x
y adoptando la variable adimensional (u), tal que: u ; dx  .du

d 4w
la (1.4) queda:  4w  0 (1.6)
du 4

1.2 Clasificación de las Vigas

Matemáticamente se demuestra que la solución de la ecuación (1.6) representa 4 ondas

amortiguadas, 2 desde cada extremo de la viga, que decrecen conforme se acerca al otro extremo.
En la práctica se establece una clasificación orientada a considerar o no la influencia recíproca
de los dos extremos de la pieza de acuerdo con el valor de la relación (L/λ):

viga de longitud corta: (L/λ) ≤ 1.50, la viga se considera rígida y se analiza como tal.

viga de longitud intermedia: 1.50 < (L/λ) ≤ 5.0, se considera flexible con importante influencia
recíproca de las perturbaciones entre uno y otro extremo de la barra por lo que es necesario tomar
en cuenta las 4 ondas amortiguadas.

viga de longitud indefinida: (L/λ) > 5.0, se considera muy flexible con nula o escasa influencia
recíproca entre los extremos por lo se toman en cuenta 2 ondas amortiguadas.

2 Fundaciones Elásticas
 Esta clasificación no es invariable, existen criterios que señalan límites diferentes que
quedan a juicio del proyectista. El principal objetivo de esta clasificación es el de llamar la
atención sobre el comportamiento diferente de cada tipo de viga.

1.3 Vigas de Longitud Intermedia

1.3.1 Solución Homogénea

Brevemente se describe la solución de la (1.6) con base en una forma exponencial, con las
constantes F y α, cuya expresión y su respectiva derivada son:
d 4w
w  F.eu ;   4 F .e u
du 4
que reemplazadas en (1.6) dan origen a la ecuación característica: α 4 + 4 = 0
cuya solución, usando funciones hiperbólicas, resulta:

w = A1 C.c + A2 S.c + A3 C.s + A4 S.s (1.7)

con la siguiente nomenclatura:

c = cos u ; s = sen u ; u=x/λ

(1.8)
u –u u –u
C = cosh u = (e + e )/2 ; S = senh u = (e – e )/2

A1 , A2 , A3 , A4 = constantes de integración

Por derivación de la (1.7), según las (1.1a,b,c), con los términos correspondientes a la
solución particular se obtiene el sistema de ecuaciones del CUADRO Nº 1

CUADRO Nº 1 VIGAS DE LONGITUD INTERMEDIA (1.5 < L/λ < 5.0)

Valores en cualquier punto interior para u = x / λ)

Solución Factor C.c S.c C.s S.s

Particular
w= wo 1 A1 A2 A3 A4
φ = φo 1/λ A2 + A3 A1 + A 4 -A1 + A4 -A2 + A3
m= mo 2EI/λ2 -A4 -A3 A2 A1
V= Vo 2EI/λ3 A 2 – A3 A1 – A4 A 1 + A4 A 2+ A3

A modo de ejemplo se muestra la interpretación de la 3ª ecuación del CUADRO Nº 1:

m = mo + 2EI/λ2 (-A4 C.c – A3 S.c + A2 C.s + A1 S.s)

Según (1.2) la presión unitaria en el suelo es: q   .w (1.9)

Fundaciones Elásticas 3
 1.3.2 Constantes de Integración:

Las constantes de integración dependen de las condiciones de borde: rotaciones y

desplazamientos de los extremos de la viga, calculadas con las soluciones homogénea y
particular.

Denominando (w1,φ1) a las deformaciones según la solución homogénea y (w10,φ10) a las

de la solución particular del extremo izquierdo, y (w2,φ2 , w20 , φ20) las correspondientes al
extremo derecho -(ver Fig 1.5)- , se pueden deducir expresiones de las constantes de integración
a partir de las ecuaciones del CUADRO Nº 1. Así, para el extremo izquierdo : x = u = 0 , las
funciones (S , s ) se anulan y entonces, de las dos primeras ecuaciones del CUADRO Nº 1 se
obtienen:
w1 = w10 + A1

1  10  ( A2  A3 )

y para el extremo derecho: x = L, u = L/ λ, de las mismas ecuaciones del CUADRO Nº 1:

w2 = w20 + A1 C.c + A2 S.c + A3 C.s + A4 S.s

λ φ2 = λ φ20 + A1(S.c – C.s) + A2(C.c-S.s) + A3(C.c+S.s) + A4(S.c+C.s)

resolviendo el sistema que se obtiene al combinar las cuatro expresiones anteriores se obtiene el
CUADRO Nº 2

CUADRO Nº 2 VIGAS DE LONGITUD INTERMEDIA (1.5 < L/λ < 5.0)

Constantes de Integración con (u = L / λ)

Factor w 1 – w 10 λ(φ 1 – φ 10) w 2 – w 20 λ (φ 2 – φ 20)

A1 = 1 1 0 0 0
A2 = 1/(S 2 – s 2) -(C.S + c.s) -s 2 (S.c + C.s) -S.s
A3 = 1/(S 2 – s 2) (C.S + c.s) S2 -(S.c + C.s) S.s
A4 = 1/(S 2 – s 2) -(S 2 + s 2) -(C.S – c.s) 2 S.s (S.c – C.s)

1.4 Vigas de Longitud Indefinida

Se estudia este tipo de vigas porque aparte de su valor intrínseco, constituye un medio
muy útil para resolver el problema de la solución particular de vigas de longitud intermedia con
cargas entre los extremos.

Según la Sección 1.2, una viga es de sección indefinida si (L/λ) > 5.0 ; si la viga reposa
en un lecho elástico con una carga aplicada en una sección distante de los extremos y se elige el
punto de aplicación de la carga como origen de coordenadas, medidas hacia la derecha, los
desplazamientos, giros, momentos y cortes tienden a disminuir con el aumento de la abscisa
adimensional (u = x/λ).

4 Fundaciones Elásticas
 En tales casos, los términos con ( eu ) en las funciones hiperbólicas de (1.7) deben
anularse porque sus valores crecen con el incremento de la abscisa, lo cual es un absurdo físico, y
entonces la solución (1.7) con la nomenclatura (1.8) para la viga de longitud indefinida queda:

w  D1.eu .c  D2 .eu .s (1.10)

donde c = cos u ; s = sen u ; u = x/λ

D1 , D2 = constantes de integración que dependen del tipo de carga.

Con derivaciones sucesivas de (1.10), adicionando los términos correspondientes a la

solución particular se obtiene el sistema de ecuaciones del CUADRO Nº 3

CUADRO Nº 3 VIGAS DE LONGITUD INDEFINIDA ( L/λ > 5.0 )

Valores en cualquier punto para (u = x / λ)

Solución Factor e–u.c e–u.s

Particular
w wo 1 D1 D2
φ φo 1/λ -(D1-D2) -(D1+D2)
m mo 2EI/λ2 D2 -D1
3
V Vo 2EI/λ -(D1+D2) (D1-D2)

1.4.1 Caso de Carga Concentrada

En la Fig 1.4(a) se muestra una viga de longitud indefinida con una carga concentrada
aplicada en un punto alejado de los extremos. Tomando como origen de las abscisas el punto de
aplicación de la carga, se esquematizan los diagramas de desplazamientos (w), rotaciones (φ),
momentos (m) y cortes (V) con sus respectivas propiedades.

Fig 1.4(a)

Fundaciones Elásticas 5
 Constantes de Integración: Con carga uniforme nula (p = 0), resulta wo = 0 según (1.2),
y según (1.1a,d) también se anula el resto de la solución particular. Tomando en cuenta que en el
origen u = 0, de modo que: e-u.c = 1, e-u.s = 0, y reemplazando los valores de (φ , V) en la 2ª y
4ª ecuaciones del CUADRO Nº 3 se obtienen las constantes de integración y el CUADRO Nº 4.
P. 3
D1  D2 
8EI

CUADRO Nº 4. VIGAS DE LONGITUD INDEFINIDA (L/λ>5.0)

Con carga concentrada. Valores en cualquier punto con (u=x/λ)

Factor e-u.c e-u.s

w Pλ3/8EI 1 1
2
φ Pλ /4EI 0 -1
m Pλ/4 1 -1
V P/2 -1 0

1.4.2 Caso de Momento Aislado

En la Fig 1.4(b) se muestra una viga de longitud indefinida con un momento aislado
aplicado en un punto alejado de los extremos. Tomando como origen de las abscisas el punto de
aplicación del momento, se esquematizan los diagramas de desplazamientos (w), rotaciones (φ),
momentos (m) y cortes (V) con sus respectivas propiedades.

Fig 1.4(b)

con las mismas consideraciones hechas en el caso de carga concentrada y reemplazando los
valores de (w , m) en la 1ª y 3ª ecuaciones del CUADRO Nº 3 se obtienen las constantes de
integración D1 , D2 , y el CUADRO Nº 5:

6 Fundaciones Elásticas
 M 2
D1 = 0 ; D2 
4 EI

CUADRO Nº 5. VIGAS DE LONGITUD INDEFINIDA (L/λ>5.0)

Con momento aislado. Valores en cualquier punto con (u = x/λ)

Factor e-u.c e-u.s

2
w Mλ /4EI 0 1
Φ Mλ/4EI 1 -1
m M/2 1 0
V M/2λ -1 -1

1.4.3 Caso de Carga Uniforme Parcial

Las características en un punto dado de una viga de longitud indefinida con una carga
uniforme parcial que actúa en un sector determinado dependen de la ubicación del punto respecto
a la carga.

Fig 1.4(c)

Punto fuera de la carga: en la Fig 1.4(c) se muestra el punto (1) ubicado a la derecha de
la carga uniforme parcial aplicada en el sector AB = aλ de una viga de longitud indefinida. El
desplazamiento en dicho punto originado por una carga concentrada elemental: P = p.λ.du se
determina con el CUADRO Nº 3, y el efecto de toda la carga se obtiene por integración entre los
puntos B y A. Las demás características, de deducen por derivaciones.

p 4  2
w1 
8EI
 e .cos  2  e 1.cos 1 

p 3  1
1  e  cos 1  sin 1   e 2  cos  2  sin  2 
8EI 
(1.11a,d)
p 2
m1 
4
e  1
sin 1  e 2 sin  2 

p  1
V1  e  cos 1  sin 1   e 2  cos  2  sin  2 
4 

Fundaciones Elásticas 7
 Punto dentro de la carga: en la Fig 1.4(d) se muestra el punto (2) dentro de la carga. Para
determinar el desplazamiento en este punto, se superponen las influencias del sector izquierdo de
longitud (a1.λ) haciendo α1 = a1 y α2 = 0, y del sector derecho de longitud (a 2.λ) haciendo α1 = a2
y α2 = 0, utilizando la (1.11a)

Fig 1.4(d)

Las expresiones resultantes son las siguientes:

p 4
w2 
8EI
 2  e a1 cos a1  ea 2 cos a 2 

p 3  a1
2  e  cos a1  sin a1   ea 2  cos a 2  sin a 2 
8EI
(1.12a,d)
p 2
m2 
4
e  a1
sin a1  ea 2 sin a 2 

p  a1
V2  e  cos a1  sin a1   e a 2  cos a 2  sin a 2  
4

Punto en el centro de la carga: para el punto central (C) de la carga, que se muestra en la
Fig 1.4(d) se cumple a1 = a2 = a/2 , que reemplazadas en (1.12a,d) resultan:

p 4
wc 
4EI
1  e a / 2 .cos  a / 2  

c  0
(1.13a,d)
p  a / 2
2
mc  e sin  a / 2 
2

Vc  0

éstas propiedades son similares a las del caso de carga concentrada en cuanto a simetría y
antimetría respecto al centro de la carga repartida, excepto el valor del corte que en este caso es
nulo.

8 Fundaciones Elásticas
 1.5 Solución Particular en Vigas de Longitud Intermedia

1.5.1 Cargas Concentradas y Uniformes Parciales

En la Fig.(1.5) se muestra una viga de longitud intermedia y sección constante sobre un

suelo elástico de balasto también constante, con cargas arbitrarias entre apoyos que, para el
análisis de estructuras continuas, se suponen empotrados. El problema se resuelve en dos etapas.

Fig.(1.5) Deformaciones: a) Solución Particular , b) Solución Homogénea

c) Reacciones de Empotramiento

En primer lugar, se considera una viga sin restricciones en los apoyos, calculando para
las secciones S1, S2. …….Sn, convenientemente escogidas, las características wo, φo, mo, Vo
producidas por todas y cada una de las cargas presentes en la viga y luego superpuestas en una
solución total. Frecuentemente se acepta como solución particular la que corresponde a una viga
de longitud indefinida utilizando el CUADRO Nº 4 para cargas concentradas ó el Nº 5 para
momentos aislados, tomando como origen el punto de aplicación de la carga o momento con
abscisas crecientes (x) hacia la derecha, haciendo uso de las propiedades de simetría o antimetría
para las secciones hacia la izquierda del origen. En el caso de carga uniforme parcial se usará una
de las expresiones (1.11), (1.12), (1.13) según la ubicación de la sección: dentro o fuera de la
carga. Obsérvese la nomenclatura utilizada para las características en los extremos: w10, φ10, m10,
V10 en el extremo izquierdo (1) y w20, φ20, m20, V20 en el extremo derecho (2)

En segundo lugar, la viga real exige desplazamientos y rotaciones nulos en los apoyos,
por lo que con los valores de w10 , φ10 , w20 , φ20, haciendo w1= φ1 = w2 = φ2 = 0 porque la
solución homogénea es nula, se calculan las constantes de integración con el CUADRO Nº 2 y
las características adicionales con el CUADRO Nº 1 que permiten determinar las reacciones de
empotramiento en los apoyos.

Ejemplo N.1 [1]

En la figura anexa se muestra una viga con 3 tipos de cargas: uniforme parcial (p),
concentrada (P) y momento (M), con 5 secciones escogidas convenientemente, incluidas las de
los apoyos. En el cuadro anexo (X) representa la abscisa acumulada de cada sección a partir del
apoyo izquierdo mientras (x), no mostrada en el cuadro, representa la abscisa de una sección
respecto al origen de la carga en estudio. En la primera columna del cuadro se muestran las
cargas P, M con el número de la sección donde están aplicadas.

Fundaciones Elásticas 9
Características de la viga:

módulo de elasticidad E = 210 tf/cm2 ; Luz L = 6.00 m ; Balasto β = 2.0 kgf/cm3

sección: b = 0.60 m ; h = 0.40 m. ; base r = 0.60 m ; I = 32.0 dcm4

E.I = 6.72 x 103 tf.m2 ; β.r = 2 x 0.60 x 103 tf/m2 ; λ = (4 x 6.72 / 1.2 ) ¼ = 2.1755 m

X 0.00 1.50 3.00 4.50 6.00

p 0.9292 2.209 3.260 3.260 2.2090
M3 -1.3058 -1.686 0.000 1.686 1.3058
P4 0.1166 0.679 1.624 2.298 1.6236
wo (w10) -0.2600 1.202 4.884 7.244 (w20) 5.1384
w 0.0000 0.0000
p 0.7101 0.909 0.394 -0.394 -0.9091
M3 -0.4837 0.165 2.428 0.165 -0.4837
P4 0.2346 0.522 0.675 0.000 -0.6745
φo (φ10) 0.4610 1.596 3.497 -0.229 (φ20) -2.0673
φ 0.0000 0.0000
p -1.7482 0.687 3.351 3.351 0.6570
M3 -0.7201 -5.808 -15.00 / 15.00 5.808 0.7201
P4 -1.1185 -1.300 0.4435 6.527 0.4435
mo (m10) -3.5868 -6.451 -11.205 / 18.795 15.686 (m20) 1.8206
Mf -2.2177 -25.4097
p 0.4093 3.185 0.737 -0.727 -3.1854
M3 -2.0356 -4.871 -6.895 -4.871 -2.0356
P4 -0.3620 0.288 2.323 6.000 / -6.000 -2.3232
Vo (V10) -1.9883 -1.398 -3..835 -0.402/ -11.598 (V20) -7.5442
Vf -0.9744 -26.7466

Los números en negrilla representan los valores de cada item en el origen de las abscisas
parciales o sea la sección donde está aplicada la carga respectiva

El momento y el cortante presentan dos valores en las secciones donde están aplicados el
momento y la carga que corresponden a las magnitudes a la izquierda y derecha de la sección.

Con los valores w10 , φ10 , w20 , φ20 y el Cuadro Nº2 se obtienen las constantes de integración:

A1 = 0.2600 , A2 = -0.1131 , A3 = -0.8898 , A4 = -0.4821

y del CUADRO Nº 1 los valores finales de w , φ , Mf , Vf en los apoyos.

10 Fundaciones Elásticas
 1.5.2 Carga Uniforme Total

Generalmente, las cargas que actúan en el tramo de una viga de fundación se deben a su
peso propio y a una posible reacción de la losa de piso que constituyen el caso de carga uniforme
en el tramo. El problema puede resolverse aplicando la sección 1.4.3, sin embargo se presenta la
solución directa de inmediata aplicación.

Por acción de la carga uniforme (p) el desplazamiento (wo) es constante porque (p) es
constante, por lo tanto de (1.2) con (1.5) se tiene:

p p. 4
wo   (1.14)
 .r 4 EI

ahora w10 = w20 = wo , mientras w1 = φ1 = w2 = φ2 = φ10 = φ20 = 0, y del CUADRO Nº 2 las

constantes de integración resultan:

C c S s
A1   wo , A2  wo , A3   A2 , A4  wo (1.15)
Ss Ss

las reacciones de empotramiento, según el CUADRO Nº 1, haciendo m o = Vo = 0 tomando

en cuenta la (1.14) , quedan:
(S  s)
M f 1   M f 2  p 2
2S  s
(1.16a,b)
(C  c)
V f 1  V f 2  p
S  s
momentos con tracción en las fibras superiores de los apoyos y cortes con igual signo.

1.6 Estructuras continuas

Para calcular las deformaciones se aplica el método de los desplazamientos con la

nomenclatura y convención de signos mostrada en la Fig 1.6(1), que coincide con la convención
de signos de la Fig. 1.1 para las rotaciones y desplazamientos. El cambio de nomenclatura
adoptado permite diferenciar el cálculo matricial del análisis por resistencia de materiales.

1.6.1 Ecuación Matricial de una Viga:

La Fig.1.6(1) muestra una viga (m) entre las juntas i-j con coordenadas locales positivas,
numeradas convenientemente para el manejo de este tipo de problemas, con detalles de las
relaciones entre deformaciones y solicitaciones.

Fundaciones Elásticas 11
 Fig. 1.6 (1) Coordenadas locales de una viga

la ecuación matricial resulta:

(2i-1) (2i) (2j-1) (2j) * numeración de las casillas de la matriz

qm = km * dm (1.17)

con las rigideces de una barra en suelo elástico [1]

2EI CS  cs 2EI S 2  s 2 4EI CS  cs

k b 2 2 2 t
 S 2  s2  S s  3 S 2  s2
(1.18)
2EI Cs  Sc 4 EI Ss 4 EI Cs  Sc
a bo  to 
 S 2  s2  2
S  s2
2
 3 S 2  s2

y las funciones :

s = sen (L/λ) c = cos (L/λ) S = senh(L/λ) C = cosh(L/λ) (1.19)

1.6.2 Vigas Continuas

Una viga continua la forman tramos con longitudes fijadas por las distancias entre
columnas, cambios de sección o variación de balasto, es decir que el procedimiento propuesto se
limita a tramos con valores constantes de inercia, ancho de la base y balasto aunque pueden ser
diferentes de un tramo a otro.

En cada nodo de unión de dos tramos contiguos pueden actuar momentos y cargas
concentradas, y en cada tramo pueden actuar cargas que producen reacciones de empotramiento
según lo descrito en la Sección 1.5. Tales cargas producirán una rotación de la elástica y un
desplazamiento vertical que deben ser idénticos para los dos extremos que concurren al nodo y
que modifican los resultados de la etapa inicial.

12 Fundaciones Elásticas
 La ecuación matricial de la viga continua se construye directamente con las rigideces de
cada tramo siguiendo la numeración de las casillas indicadas en el sistema (1.17). En la Fig.
1.6(2) se muestra la ecuación matricial de una viga continua de 3 tramos:

Q = K * D

Fig 1.6(2) Ecuación Matricial de una Viga de 3 Tramos

Resuelto el sistema D = K-1.Q , se analiza cada tramo independientemente uno de otro

con la convención original de resistencia de materiales y los valores (-di , di+1) de cada extremo.

Ejemplo Nº 2: Sea la viga continua formada por un volado y dos tramos con las
características y cargas nodales que se muestran en la figura anexa, con un módulo de elasticidad
E = 230 tf/cm2 y balasto diferente para cada tramo.

Cargas y Propiedades de la Viga:

P

B,h .60 x .30 .60 x .30 .80 x .30 m

L 1.60 6.00 8.00 m
EI 3.105 3.105 4.140 (tf.m2).10 -3
β 2.00 1.80 1.50 kgf/cm3
λ 1.794 1.842 1.927 m
L/λ 0.892 3.26 4.15

Rigideces x 10-3 (según 1.18)

k 7.809 3.380 4.299
a 3.846 0.228 -0.043
b 7.438 1.832 2.232
bo 7.183 -0.033 -0.119
t 9.808 1.996. 2.317
to 8.851 -0.170 -0.101

Fundaciones Elásticas 13
 Matriz de Rigidez K
9.808 7.438 -8.851 7.183
7.438 7.809 -7.183 3.846
-8.851 -7.183 11.804 -5.606 0.170 -0.033
7.183 3.846 -5.606 11.189 0.033 0.228
0.170 0.033 4.313 0.401 0.101 -0.119
-0.033 0.228 0.401 7.679 0.119 -0.043
0.101 0.119 2.317 -2.232
-0.119 -0.043 -2.232 4.299

Vector de Cargas Q Vector de Desplazamientos D

q1 0.0 d1 0.34468
q2 0.0 d2 4.18758
q3 30.0 d3 6.86719
q4 = 20.0 d4 = 3.65566
q5 60.0 d5 13.83905
q6 - 40.0 d6 -6.33533
q7 40.0 d7 23.40265
q8 - 25.0 d8 6.65549

Con el par de deformaciones de cada extremo de un tramo (di , di+1) , (dj , dj+1) se
encuentran las constantes de integración del tramo con el CUADRO Nº 2 y luego las
propiedades en cualquier punto del tramo, con el CUADRO Nº 1 .

Por ejemplo, para el tramo 2, entre las juntas 2 y 3:

junta 2: w1 = d3 = 6.86719 , φ1 = d4 = 3.65566

junta 3: w2 = d5 = 13.83905 , φ2 = d6 = -6.33533

las constantes de integración se obtienen del CUADRO Nº 2, con u = L /λ = 3.26

A1 = 6.86719 ; A2 = - 8.18115 ; A3 = 14.91308 ; A4 = -13.07637

y, las características en los cuartos de la luz, del CUADRO Nº 1, con u = x /λ

x u = x /λ w φ m V Presión
m mm Rad*103 tf.m tf tf/m2
0.00 0.00 6.87 3.66 23.95 -22.97 12.36
1.50 0.814 7.28 -1.07 -1.23 -10.58 13.11
3.00 1.629 7052 1.98 -8.62 0.77 13.54
4.50 2.443 12.87 4.19 3.46 16.81 23.17
6.00 3.257 13.84 -6.34 46.15 40.52 24.91
De (1.9), la presión unitaria = β.w

14 Fundaciones Elásticas
 1.6.3 Vigas sobre Apoyos Fijos

Se consideran apoyos fijos aquellos que tienen una gran rigidez axial. En la fig 1.6(3) se
muestra una viga continua de tres tramos con apoyos fijos y tramos sobre suelo elástico. En tal
caso los desplazamientos wi son nulos y solo existe la solución particular con wo > 0, por lo que
la ecuación matricial se plantea en función de momentos y rotaciones en la forma tradicional.

La ecuación matricial (1.17) del tramo i-j queda:

(1.20)
en la que los valores de (mi , mj) corresponden a los de la solución particular. La ecuación
matricial para la viga de tres tramos, resulta:

Q = K * D (1.21)

Fig 1.6(3) Ecuación Matricial de una Viga de 3 Tramos con Apoyos Fijos

Si el apoyo es flexible, su rigidez a flexión se incluye en la diagonal de la matriz.

Conocido el valor de wo = p /β.r de cada tramo y las rotaciones φ de cada junta, entonces
w10 = w20 = wo del tramo respectivo mientras w1 = w2 = φ10 = φ20 = 0 y las constantes de
integración del tramo i-j , según el CUADRO Nº 2 , se presentan en el CUADRO Nº 6

CUADRO Nº 6 VIGAS DE LONGITUD INTERMEDIA CON APOYOS FIJOS

Constantes de Integración con (u = L / λ)

Factor – w 10 λφi – w 20 λ.φ j

A1 = 1 1 0 0 0
A2 = 1/(S 2 – s 2) -(CS + cs) -s 2 (Sc + Cs) -Ss
A3 = 1/(S 2 – s 2) (CS + cs) S2 -(Sc + Cs) Ss
A4 = 1/(S 2 – s 2) -(S 2 + s 2) -(CS – cs) 2 Ss (Sc – Cs)

Las características de cualquier punto del tramo se calculan con el CUADRO Nº 1.

Fundaciones Elásticas 15
 Ejemplo Nº 3 [1]: En el cuadro anexo se muestra una viga de tres tramos sobre apoyos
fijos, los dos primeros tramos descansan en el suelo mientras el tercero está libre con el extremo
derecho empotrado en el apoyo. Las cargas son del tipo uniforme, diferente de un tramo a otro. El
módulo de elasticidad es E = 210 tf/cm2.

Cargas y Propiedades de la Viga

b,h 60 x 50 60 x 50 60 x 50 m
L 4.00 4.50 4.00 m
EI 13.125 13.125 13.125 tf.m.10 -3
β 3.00 3.00 Kgf/cm3
λ 2.3239 2.3239 m
L/λ 1.7212 1.9364

Rigideces x 10-3 (según 1.18) (viga común)

k 14.1606 13.0936 13.125
a 5.7989 4.7917 6.563

Reac. Empotramiento (s/1.16a,b) (viga común)

M 3.766 -3.766 7.369 -7.369 2.667 -2.667 m
V 5.726 5.726 10.041 10.041 4.000 4.000 tf

Matriz de Rigidez (según 1.21)

14.1606 5.7989
5.7989 27.2542 4.7917
4.7917 26.2186

Vector de Cargas Vector de Rotaciones

φ1 0.216131
ΣM1 3.7660 φ2 = 0.121653
ΣM2 = 3.6030 φ3 -0.201572
ΣM3 -4.7020

Tramo 1: L = 4.00 λ = 2.3239 ; p = 3.00 ; wo = 1.6666

Constantes de integración: A1 = -1.6667 ; A2 = 1.1724 ; A3 = -0.6701 ; A4 = 0.0000
(según Cuadro Nº 6)

Características en los Cuartos de la Luz

(según Cuadro Nº 1)
x u = x/λ w φ m V Presión
0.00 0.00 0.0000 0.2161 0.0004 3.8537 0.0000
1.00 0.4303 0.1765 0.1062 2.4151 1.0296 0.5295
2.00 0.8606 0.1860 -0.0834 2.1209 -1.6153 0.5580
3.00 1.2909 0.0510 -0.1489 -0.8617 -4.3916 0.1530
4.00 1.7212 0.0000 0.1217 -6.7416 -7.3857 0.0000
Presión Unitaria (kgf/cm2)= β.w

16 Fundaciones Elásticas
Tramo 2: L = 4.50 λ = 2.3239 ; p = 4.80 ; wo = 2.6667
Constantes de integración: A1 = -2.6666 ; A2 = 2.5166 ; A3 = -2.2339 ; A4 = 1.3869

Características en los Cuartos de la Luz

x w φ m V Presión
0.000 0.0000 0.1217 -6.7415 9.9361 0.0000
1.125 0.3062 0.3053 1.4942 4.8129 0.9186
2.250 0.5117 0.0196 4.3285 0.2967 1.5351
3.175 0.3395 -0.2951 2.1773 -4.1803 1.0185
4.500 0.0000 -0.2016 -5.3121 -9.2528 0.0000

Tramo 3: L = 4.00 , p = 2.00

Es una viga apoyada en sus dos extremos de modo que puede tratarse como una viga
común sin influencia de la elasticidad del suelo. Existen procedimientos numéricos que permiten
analizar este tipo de vigas en forma similar a la empleada en los tramos 1 y 2, pero se elige el
método más simple habida cuenta de que se conocen rigideces, cargas y rotaciones, aplicando el
procedimiento común con la convención de signos de resistencia de materiales. Los valores de
interés: momentos y cortantes en los cuartos de la luz de muestran en el siguiente cuadro.

Características en los Cuartos de la Luz

x m V
0.00 -5.312 4.992
1.00 -1.320 2.992
2.00 0.672 0.992
3.00 0.664 -1.008
4.00 -1.344 -3.008

1.7 Programa FUNDAVIGA

Se presenta un programa de aplicación práctica a modelos simples: módulo de balasto

variable de un tramo a otro, juntas en los apoyos de las columnas, sin apoyos fijos, y tramos de
sección constante, con cargas concentradas en las juntas y cargas uniformemente repartidas en los
tramos. Con esta idea se intenta eludir el laborioso cálculo discreto de las soluciones particulares
expuestas en la sección 1.5 para cargas parciales poco importantes. En los casos de cargas
parciales importantes o reacciones de vigas secundarias es posible sustituirlas por cargas
uniformes equivalentes ó aumentar el número de juntas para obtener tramos que satisfagan los
requerimientos expuestos.

Fundaciones Elásticas 17
1.7.1 Codificación de Datos

Cada línea de datos debe estar precedida de la palabra DATA. La última línea del
programa, incluyendo los datos, debe ser Fin:END. Se prevé miembros con sección transversal
arbitraria a fin de incluir el procesamiento de fundaciones imaginarias para estudio ó
investigación. En el programa anexo se mantienen los datos del Ejemplo Nº 4 como muestra de
codificación.

1. Datos Generales

DATA A$ , Ec , Nj , Ntipo , Nq
A$ : identificación del trabajo hasta 50 caracteres alfanuméricos
Ec : módulo de elasticidad longitudinal ( tf/cm2)
Nj = Nº de juntas ; Ntipo = Nº de tipos de viga ; Nq = Nº de casos de carga

2. Datos de Geometría (ver ANEXO A)

Para cada tipo de viga:

Sección rectangular: DATA 1 , Bo , Ho = base, altura (m)
Sección T invertida: 2 , Bo , Ho, Bf , Hf = (base, altura) nervio, ala (m)
Sección arbitraria: 3 , Io , Bj = Inercia (dcm4) , ancho de base (m)
en cualquier caso: Bsto = módulo de balasto (kgf/cm3)

Para cada viga : DATA J , J1 , J2 , L , Typ

J = Nº de la viga ; J1 , J2 = Nº de las juntas de la viga.
L = luz (m) ; Typ = número del tipo de viga.

3. Datos de Cargas

Para cada caso de Carga:

DATA Njq = Nº de juntas con cargas
Si Njq > 0 , Para cada junta cargada: J , Q1 , Mx
J = Nº de la junta, Q1 = carga vertical (tf), Mx = momento (tf.m)

DATA Ptr = Nº de vigas con carga uniforme.

Si Ptr > 0 , Para cada viga cargada: J , Pt
J = Nº del tramo, Pt = carga uniforme (tf/m)

4. Archivo de Resultados

Los resultados de todas las vigas se archivan en el arreglo Rest (5, 6*Nº de vigas). Las 5
filas corresponden a las fracciones de la luz: 0, L/4, L/2, 3L/4 cuando L>2.0 m; si L<=2.0 m las
secciones corresponden a: 0, L/2, L, y, en cada fila en grupos de 6, a la abscisa, flecha, giro,
momento, corte, presión, correspondientes a cada sección.

18 Fundaciones Elásticas
100 ! FUNDAVI GA Viga s de Fu nda c ió n E lást ica.
110 !- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -- - -- - -- - - -
120 OP T I ON BAS E 1
130 DIM A$[50], L$[130], Tit1$[35], Tit2$[100], Tit3$[85], Tit4$$[120], Tit5$[70], DIM D(50)
140 DIM F(50,50), K(50,50), Mb(50,8), Mtyp(50.8), Pjta(20,2), Ptra(20), U(60), X(20.3), W(5,150)
150 L$=RP T $( “- “, 130)

160 RAD ! Dec lar a Rad ia nes

170 RE S T ORE E je mp lo ! E sco ge el pr o ble ma
180 RE AD A$, E c, Nj, Nt ipo , Nq
190 Nvt = Nj - 1
200 Nd j = 2 * Nj

210 RE DI M F( Nd j, Nd j) , K( Nd j, Nd j) , Mt yp( Nvt , 8), Mb( Nvt , 8) , D( Ndj) , U( Nd j)

220 RE DI M M f( 2*N j, 2) , Vf( 2* N j, 2) , P jt a( N j, 2) , Pt r a ( Nj) , W( 5, 6 * Nvt )
230 MAT K=( 0)
240 M AT Mt yp=( 0)

250 FOR I =1 T O Nt ipo

260 B f=0
270 H f=0
280 RE AD T ipo
290 ON T ipo GOT O S _1, S _2, S _3
300 ! S ecc ió n Rect angu lar
310 S _1: RE AD Bo , Ho
320 A= Bo *Ho
330 I o =Bo *Ho ^3/ 12
340 B j= Bo
350 GOT O S _4
360 ! S ecc ió n T inver t ida
370 S _2: RE AD Bo , Ho , Bf, H f ! Ner vio , Ala ( m)
380 Bh=B f- Bo
390 A= Bo *Ho +Bh*H f
400 M=( Bo *Ho ^2+Bh*H f^2) / 2
410 I o =( Bo *Ho ^3+Bh*H f^3) / 3 - M^2/ A
420 B j= B f
430 GOT O S _4
440 ! S ecc ió n Ar bit r ar ia
450 S _3: RE AD I ner c ia, B j ! I ner c ia ( dc m4) , Base ( m)
460 I o =I ner c ia*1. E - 4 ! I ner c ia ( m4)
470 S - 4: RE AD Bst o ! Bala st o ( kgf/ c m3)
480 Mt yp( I , 1) =Bo
490 Mt yp( I , 2) =Ho
500 Mt yp( I , 3) =B f
510 Mt yp( I , 4) =H f
520 Mt yp( I , 5) =I o
530 Mt yp( I , 6) =B j
540 Mt yp( I , 7) =Bst o
550 NE XT I

560 FOR I =1 T O Nvt

570 RE AD J, J1, J2, L, T yp ! N. V iga, Junt as, Luz, N. T ipo
580 I o =Mt yp( T yp, 5)
590 B j= Mt yp( T yp, 6)
600 Bst o =Mt yp( T yp, 7)

Fundaciones Elásticas 19
610 E 1=10*E c *Io ! R ig ide z/ 1000
620 L b=( ( 4*E 1) / ( Bj*B st o) ) ^0. 25 ! La mbda
630 Mb( J, 1) =J1
640 Mb( J, 2) =J2
650 Mb( J, 3) =L
660 Mb( J, 4) =T yp
670 Mb( J, 5) =E 1
680 Mb( J, 6) =Lb
690 Mb( J, 7) =Bst o
700 NE XT I

710 FOR I =1 T O Nvt

720 L i=M b( I , 3)
730 L b=M b( I , 6)
740 E 1=Mb( I , 5)
750 GOS UB S _10
760 Ki=2*E 1/ ( L b*( K6 - K7) )
770 Ko =Ki*( K5 - K8)
780 A= Ki*( K3 - K2)
790 B= Ki*( K6 + K7) / L b
800 Bo =2*Ki* K4/ L b
810 T =2*Ki*( K5 + K8) / L b^2
820 T o =2*Ki*( K2 + K3) / L b^2

830 I 2=2*Mb( I , 1)
840 I 1=I 2 - 1
850 J2=2*M b( I , 2)
860 J1=J2 - 1

870 K( I 1, I 1) =K( I 1 , I 1) +T ! T r iá ng u lo S uper io r de la Mat r iz de R ig idez

880 K( I 1, I 2) =K( I 1, I 2) + B
890 K( I 2, I 2) =K( I 2, I 2) +Ko
900 K( I 1, J1) = - To
910 K( I 1, J2) = Bo
920 K( I 2, J1) = - Bo
930 K( I 2, J2) =A
940 K( J1, J1) = K( J1, J1) +T
950 K( J1, J2) = K( J1 , J2) - B
960 K( J2, J2) = K( J2, J2) + Ko
970 NE XT I

980 FOR I =1 T O Ndj ! Co mp let a la Mat r iz

990 FOR J=I T O Nd j
1000 K( J, I ) =K( I , J)
1010 NE XT J
1020 NE XT I
1030 M AT F= I NV( K)

1040 FOR I = 1 T O Ndj ! M ini m iz a er r or es de r edo ndeo

1050 FOR J=2 T O Nd j
1060 F( I , J) =( F( I , J) +F( J, I ) ) / 2
1070 F( J, I ) =F( I , J)
1080 NE XT J
1090 NE XT I

20 Fundaciones Elásticas
1100 Nq1=0
1110 ! CAS OS DE CARG A
1120 S _5: Nq1=Nq1+1
1130 M AT U=( 0)
1140 M AT P jt a=( 0) ! Nº Junt as C ar gadas
1150 RE AD N jq
1160 I F N jq>0 T HE N
1170 FOR I =1 T O Njq
1180 RE AD J, Q1, Mx ! N. ju nt a , Q , M
1190 P jt a( J, 1) =Q1
1200 P jt a( J, 2) =Mx
1210 NE XT I
1220 E ND I F
1230 M AT Pt r a =( 0)
1240 RE AD P t r ! Nº T r amo s car gado s
1250 I F P t r >0 T HE N
1260 FOR I =1 T O Pt r
1270 RE AD J, P t ! N. Tr a mo , p( unifo r me)
1280 P t r a ( J) =Pt
1290 L i=M b( J, 3)
1300 L b=M b( J, 6 )
1310 G=L i/ L b
1320 Ch=( E XP ( G) +E XP( - G)) / 2
1330 S h=( E XP ( G) - E XP( - G)) / 2
1340 C1=COS ( G)
1350 S 1=S I N( G)
1360 M f( J, 1) =P t *Lb^2*( S h - S 1) / ( 2*( S h+S 1) )
1370 M f( J, 2) = - M f( J, 1)
1380 V f( J, 1) =P t *Lb*( Ch - C1) / ( S h+S 1)
1390 V f( J, 2) =V f( J, 1)
1400 NE XT I
1410 E ND I F

1420 FOR I =1 T O Nj ! Car gas Ge ner a l izada s

1430 I F I =1 T HE N
1440 U( 1) =P jt a( 1, 1) + Vf( 1, 1)
1450 U( 2) =P jt a( 1. 2) +Mf( 1, 1)
1460 E LS E
1470 I F I <N j T HE N
1480 U( 2*I - 1) =P jt a( I, 1) +Vf( I - 1, 2) + Vf( I , 1)
1490 U( 2*I ) =P jt a( I , 2) +M f( I - 1, 2) + Mf( I , 1)
1500 E LS E
1510 U( 2*I - 1) =P jt a( I, 1) + Vf( I - 1, 2)
1520 U( 2*I ) =P jt a( I , 2) +M f( I - 1, 2)
1530 E ND I F
1540 E ND I F
1550 NE XT I

1560 M AT D=F*U
1570 FOR I =1 T O Nvt
1580 Wo =0
1590 J1=M b( I , 1)
1600 J2=M b( I , 2)
1610 L=M b( I , 3)
1620 T yp=M b( I , 4)

Fundaciones Elásticas 21
1630 L b=M b( I , 6)
1640 Bst o =Mb( I , 7)
1650 B j= Mt yp( T yp, 6)
1660 Wo =P t r a (I ) /( Bst o *Bj)

1670 D1=D( 2*J1 - 1) - Wo

1680 D2=D( 2*J1)
1690 D3=D( 2*J2 - 1) - Wo
1700 D4=D( 2*J2)
1710 L i=L
1720 GOS UB S _10
1730 Dn= K6 - K7

1740 A1=D1 ! Co nst ant es de I nt egr ació n

1750 A2=( - D1*( K5 + K8) - D2*L b* K7+D3*( K2+ K3) - D4*L b* K4) / D n
1760 A3=( D1*( K5+ K8) +D2* L b* K6 - D3*( K2+ K3) +D4* L b* K4) / Dn
1770 A4=( - D1*( K6 + K7) - D2*L b*( K5 - K8) +D3 *2*K4+D4* L b*( K2 - K3) ) / D n

1780 Z1=5
1790 Dx=0. 25
1800 I F L=<2. 0 T HE N
1810 Z1=3
1820 Dx=0. 50
1830 E ND I F
1840 S k=2*Mb( I , 5) / Lb^2
1850 J=6*I - 6

1860 FOR Z =1 T O Z1
1870 Li=( Z- 1) *Dx* L
1880 GOS UB S _10
1890 W( Z, J+1) =L i
1900 W( Z, J+2) = Wo +A1* K1+ A2* K2 + A3* K3+A4* K4
1910 W( Z, J+3) =( ( A2+ A3) * K1+( A1+ A4) * K2+ ( A4 - A1) *K3 +( A3 - A2) * K4) / L b
1920 W( Z, J+4) =( - A4* K1 - A3* K2+ A2* K3+ A1* K4) *S k
1930 W( Z, J+5) =( ( A2 - A3) *K1 +( A1 - A4) * K2 +( A1+ A4) * K3+( A2+ A3) * K4) * S k/ L b
1940 W( Z, J+6) = W( Z, J +2) *M b( I , 7)
1950 NE XT Z
1960 NE XT I

1970 GOS UB S _20

1980 I F Nq1=1 T HE N
1990 GOS UB Exi t _1
2000 GOS UB Exi t _2
2010 E ND I F
2020 GOS UB Exi t _3
2030 GOS UB Exi t _4
2040 I F Nq1< Nq T HE N S _5
2050 GOT O Fi n
2060 ! Func io nes
2070 S _10 : G=L i/ L b
2080 Ch=( E XP ( G) +E XP( - G)) / 2
2090 S h=( E XP ( G) - E XP( - G)) / 2
2100 C1=COS ( G)
2110 S 1=S I N( G)
2120 K1= C1*C h

22 Fundaciones Elásticas
2130 K2= C1*S h
2140 K3=S 1*C h
2150 K4=S 1*S h
2160 K5= C h*S h
2170 K6=S h*S h
2180 K7=S 1*S 1
2190 K8= C1*S 1
2200 RE T URN
2210 S _20: ! P ara im pri n i r : L o s pu n t os ( …. .) in di cados en l os t í tu l os , represe n t an
espaci os en bl an co para su com pagi n aci ón con el t ext o .
2220 T it 1$=”Unidad es NO descr it as: m , t f “
2230 T it 2$[ 1, 47] =” . Vig. *. . J1. . J2.T ipo . . .. L. .. *.. T ipo . .. . Bo .. .. .. “
2240 T it 2$[ 48, 100] =”Ho. .. .. . Bf…. . . Hf. . . .. . I nerc ia( m4) . . .. Base. . . Ba last o ”
2250 T it 3$[ 1, 80] =” - “
2260 T it 3$[ 13, 30] =”CARG AS AP LI C AD AS ”
2270 T it 3$[ 41, 85] =”*. .. . . DE FORMACI ONE S …. *…P RE S I ONE S … *”
2280 T it 4$[ 1, 57] =”Junt a…. . . Q. .. .. .. .. M. .. . *. Tra mo . . . pt .. *.. . . .. DZ. . .. .. .. “
2290 T it 4$[ 58, 80] =”RY…. * Junt a…. .P s. . *”
2300 T it 4$[ 8 1, 120] =”*…. .. mm…. Rad*1000. . * ……. kgf/ c m2. *”
2310 T it 5$=”VI G A. . . .. . X…. . Flec ha…Ro t ac ió n…. . Mo me nt o ……. Co rt e….. Pr es ió n”
2320 RE T URN
2330 !- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- -- - - -- - -- - -- -- - -- - - - - -- - -- - -- - -- -- - -- - -- - -- - - - - -- - - -- - -- - -- -- -
2340 Exit_1: PRINT USING “54A, 26A, /, 100A”; A$, “VIGA DE FUNDACION ELASTICA”, L$
2350 P RI NT US I NG “#, 35A, 4A, 3D. D, 15A” ; T it 1$, “E c=”;E c, “ t f / c m2 ”
2360 P RI NT US I NG “18 A, / , 100A”;” Ba last o (kg f / c m3) ”, L$
2370 RE T URN
2380 !- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - - - - -- - -- -- -
2390 Exi t _2 : P RI NT US I NG “100 A, / , 100A” ;T it 2$, L$
2400 FOR I =1 T O Nvt
2410 P RI NT US I NG ”#, 3D, 4A”; I , ” *”
2420 P RI NT US I NG “#, 3( 4D), 4D. 3D, 2A”; Mb( I , 1) , Mb( I , 2), Mb( I , 4), Mb( I , 3), “ *”
2430 I F I >Nt ipo T HE N P RI NT
2440 I F I <=Nt ipo T HE N
2450 T yp=M b( I , 4)
2460 Bo =Mt yp( I , 1)
2470 Ho =Mt yp( I , 2)
2480 Bf= Mt yp( I , 3)
2490 Hf=Mt yp( I , 4)
2500 Io =Mt yp( I , 5)
2510 I f Bo >0 T HE N
2520 IF Bf=0 THEN PRINT USING “#, 4D, 2(4D.3D), 2(6X, A,X),X”;I, Bo, Ho,” -“,”-“
2530 I F Bf>0 P RI NT USI NG “#, 4D, 4(4D. 3D) , X”; I , Bo , Ho , Bf, H f
2540 E LSE
2550 P RI NT US I NG “#, 4D, 4( 6X, A, X), X”; I , “ - “, “- “, “- “, “- “
2560 E ND I F
2570 P RI NT US I NG ”3X, D. 6DE , 2( 6D. 2D) ”; Mt yp( I , 5) , Mt yp( I , 6) , Mt yp( I , 7)
2580 E ND I F
2590 NE XT I
2600 P RI NT US I NG “100 A” ; L$
2610 RE T URN
2620 !- - - - - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - - - -- - -- - -- - -- - - - - -- - - --
2630 Exi t _3 : P RI NT US I NG “/ , 18A, 2D, / , 79A, A”; “C AS O DE CARG A Nº ”, Nq1, L$, “*”
2640 P RI NT US I NG “12X, 17A, 11X, 40A”;T it 3$[ 13, 30], T it 3$[ 41, 85]
2650 P RI NT US I NG “80 A, / , 40X, 50A, /, 79A, A” ; T it 4$[ 1, 80], T it 4$[ 81, 120] , L$, “*”
2660 P ma x=0

Fundaciones Elásticas 23
2670 FOR J=1 T O N j
2680 I F J=1 T HE N Bst o =Mb( J, 7)
2690 I F J >1 AND J< N j T HE N Bst o =MAX( M b( J, 7) , Mb( J+1, 7) )
2700 I F J=N j T HE N Bst o =Mb( J - 1, 7)
2710 J1=2*J- 1
2720 P r es ió n= 0. 1*D( J1) *Bst o
2730 P 1=P jt a( J, 1)
2740 P 2=P jt a( J, 2)
2750 P 3=P t r a( J)
2760 D j1=D( J1)
2770 D j2=D( J1 +1)
2780 P r s= P r es io n
2790 I F J < N j T HE N
2800 PRINT USING “#, 4D,2(7D.2D), 3A, 4D, 4D.2D,3A”; J,P1,P2,” *”, J,P3,” *”
2810 P RI NT US I NG “2( 4D. 5D) , 3A, 4D, 7D.2D, 2A) ” ;D j1, D j2, “ *”, J, Pr s, ” *”
2820 E LS E
2830 P RI NT US I NG “#, 4D, 2( 7D. 2D), 3A, 11X, 2( 3A) ”; J, P 1, P 2, ” *”, ” *”
2840 P RI NT US I NG “2( 4D. 5D), 3A, 4D, 7D.2D, 2A) ” ;D j1, D j2, “ *”, J, Pr s, ” *”
2850 E ND I F
2860 NE XT J
2870 P RI NT US I NG “79 A, A” ; L$, ”*”
2880 RE T URN
2890 !- - - - - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - - - - -- - - --
2900 Exi t _4 : P RI NT US I NG “ / , 70A”;L$
2910 P RI NT US I NG “80 A” ; T it 5$
2920 P RI NT US I NG “70A” ; L$
2930 FOR I =1 T O Nvt
2940 Z1=5
2950 I F Mb( I , 3) <=2. 0 T HE N Z1=3
2960 J=6*I - 6
2970 FOR Z=1 T O Z1
2980 I F Z=1 T HE N
2990 P RI NT US I NG “#, 3D” ; I
3000 E LSE
3010 P RI NT US I NG “# , 3X”
3020 E ND I F
3030 X1= W( Z, J +1)
3040 W1= W( Z, J+2)
3050 G1= W( Z, J+3)
3060 M1= W( Z, J+ 4)
3070 V1= W( Z, J+5)
3080 P r es= W( Z, J+6)
3090 P RI NT US I NG “5D. 2D, 2( 6D. 4D) , 3( 8D.3D) ”;X1, W1, G1, M1;V1 ;P r es
3100 NE XT Z
3110 I F I < Nvt T HE N P RI NT
3120 NE XT I
3130 P RI NT US I NG “70 A” ;L$
3140 RE T URN
!====== ===== ===== ===== ===== === ======= ===== === == ===== ===== === ======= ==
3160 Ejemplo: DATA EJEMPLO Nº 4 [1], 210, 5, 3, 1
3170 DATA 3, 300, 1.80, 2.00, 3, 400, 1.80, 2.50, 3, 450, 1.80, 2.50
3180 DATA 1, 1, 2, 1.50, 1, 2, 2, 3, 6.00, 2, 3, 3, 4, 6.00, 3, 4, 4, 5, 1.80, 1
3190 DATA 3, 2, 65, 20, 3, 85, 35, 4, 75, 25
3200 DATA 4, 1,3.0, 2, 4.0, 3,6.0, 4, 3.0
3210 Fin: END

24 Fundaciones Elásticas
 Ejemplo Nº 4 [1] Aplicación del programa FUNDAVIGA. Viga continua de 4 tramos,
con ancho constante, diferentes secciones y balastos, con cargas en los tramos y en las juntas.

Cargas y Propiedades de la Viga

EJEMPLO Nº 4 [1] VIGA DE FUNDACION ELASTICA

Unidades NO Descritas: m , tf Ec = 210 tf/cm2 Balasto (kg/cm3)
Viga J1 J2 Tipo L * Tipo Bo Ho Bf Hf Inercia (m4) Base Balasto
1 1 2 1 1.50 * 1 1 - - - - 3.00000E-02 1.80 2.00
2 2 3 2 6.00 * 2 - - - - 4.00000E-02 1.80 2.50
3 3 4 3 6.00 * 3 - - - - 4.50000E-02 1.80 2.50
4 4 5 1 1.80 *
CASO DE CARGA Nº 1
CARGAS APLICADAS * DEFORMACIONES * PRESIONES
Junta Q M * Tramo pt * DZ RY * Junta Ps
tf tf.m * tf/m * mm Rad * 1000 * kgf/cm2
1 0.00 0.00 * 1 3.00 * 4.96538 - 0.17806 * 1 0.99
2 65.00 20.00 * 2 4.00 * 4.64913 - 0.30869 * 2 1.16
3 85.00 35.00 * 3 6.00 * 4.20382 0.30268 * 3 1.05
4 75.00 25.00 * 4 3.00 * 5.92742 0.82457 * 4 1.48
5 0.00 0.00 * 6.96976 0.49638 * 5 1.39

Viga X Flecha Rotación Momento Corte Presión

1 0.00 4.9654 - 0.1781 0.000 - 0.000 9.931
0.75 4.8287 - 0.1945 4.138 10.975 9.657
1.50 4.6491 - 0.3087 16.361 21.539 9.298

2 0.00 4.6491 - 0.3087 36.361 - 43.461 11.623

1.50 3.9490 - 0.4847 - 10.854 - 20.303 9.872
3.00 3.4702 - 0.1063 - 26.863 - 1.589 8.676
4.50 3.6521 0.3170 - 16.149 16.085 9.130
6.00 4.2038 0.3027 22.887 36.603 10.510

3 0.00 4.2038 0.3027 57.887 - 48.397 10.510

1.50 4.2274 - 0.1184 0.143 - 28.593 10.569
3.00 4.1895 0.1431 - 28.319 - 9.413 10.474
4.50 4.7680 0.6258 - 27.314 11.404 11.920
6.00 5.9274 0.8246 8.992 38.325 14.819

4 0.00 5.9274 0.8246 33.992 - 36.675 11.855

0.90 6.5136 0.5382 8.729 - 19.152 13.027
1.80 6.9698 0.4964 0.000 - 0.000 13.940

Fundaciones Elásticas 25
 Las cargas generalizadas se encuentran con el mismo procedimiento empleado en el
Ejemplo Nº 2, excepto que hay que añadir las reacciones de la solución particular con signos
compatibles con los de las cargas exteriores aplicadas en las juntas. Las deformaciones
generalizadas se calculan mediante el método de los desplazamientos, a partir de lo cual se
analiza cada tramo independientemente uno de otro con la convención de resistencia de
materiales. Para fijar conceptos se describe a continuación el procesamiento del 2º tramo:

2º Tramo, entre juntas 2 y 3: L = 6.00 m , λ = 2.079 m , L/ λ = 2.886 , p = 4.0 tf/m

β = 2.5 kfg/cm3 , r = 1.80 m , wo = p/(β.r) = 0.8889 mm

I=2 J=3
w10 = wo = 0.8889 0.8889 = w20 = wo
wi = d3 = 4.6491 4.2038 = wj = d5
w10 – wi = -3.7602 -3.3149 = w20 – wj
φ i = d4 = -0.3087 +0.3027 = φ j = d6

con las dos últimas líneas anteriores y el CUADRO Nº 2 se obtienen las constantes de
integración:

A1 = 3.7602 , A2 = -3.7393 , A3 = 2.8319 , A4 = -1.8702

con las cuales, y con el CUADRO Nº 1 se deducen las características de cada tramo mostradas en
la salida del programa.

Conviene hacer referencia a la salida de resultados del programa en lo que se refiere a

cargas y deformaciones de una Junta (i): se presentan separadas la carga (Pi) del momento (Mi) y
el desplazamiento (DZi) de la rotación (RYi), lo que permite una inmediata identificación de
juntas, acciones y resultados. Es evidente que (DZi) corresponde a valores de la deformación
generalizada (d2i-1) y (RYi) a (d2i), según la numeración de casillas en la ecuación (1.17).

26 Fundaciones Elásticas
 2. RETICULOS DE FUNDACION ELASTICA

2.1 Ecuaciones diferenciales de la Torsión

En la Fig. 2.1(a) se muestra un segmento de longitud (dx) perteneciente a un miembro que

descansa en un suelo elástico y que está solicitado por un momento de torsión (mt) con una
rotación (φt). En la Fig. 2.1(b) se muestra la sección transversal del miembro de ancho (r), el
desplazamiento (w) y presión (q) en el extremo de la sección, y la presión total (P) asumiendo
que el suelo es capaz de soportar compresiones y tracciones.

Fig 2.1

De acuerdo con lo asumido en la Sección 1.1 sobre el comportamiento elástico del suelo
se tiene para el segmento elemental:

t .r q.r P.2r
 , q   ..dx , P , dmt  (2.1a)
2 4 3
dmt  .r 3
reemplazando valores:  t (2.1b)
dx 12

Por otra parte, si G = módulo de elasticidad transversal y J = constante de torsión se

dt
tiene, de resistencia de materiales: mt  GJ (2.1c)
dx

12GJ d 2t
cuya 1ª derivada combinada con (2.1b) resulta:  t  0 (2.1d)
 .r 3 dx 2

introduciendo el concepto de longitud elástica de torsión definida como:

1/ 2
 12GJ 
t   3 
(2.2)
 r 

y adoptando la variable adimensional: x = u.λ t ; du/dx = 1/ λ t

d 2t
la ecuación homogénea (2.1d) queda:  t  0 (2.3)
du 2

Fundaciones Elásticas 27
2.2 Vigas de Longitud Intermedia

2.2.1 Solución Homogénea:

Empleando una función del tipo: φt = B.eα u , con B y α constantes se tiene:

d2φt/ du2 =α2.φt que reemplazando en la (2.3) se deduce la ecuación característica:

α2 – 1 = 0 (2.4)

cuyas raíces reales permiten escribir la siguiente solución, usando funciones hiperbólicas:

φt = F1.Ct + F2.St (2.5a)

y de (2.1c) mt = (GJ/λ t ) . (F1.St + F2.Ct) (2.5b)

en las que Ct = cosh (x/ λ t) , St = senh (x/ λ t)

Con valores de u = x /λ t desde 0 hasta L/ λt , se obtienen las rotaciones y momentos en

cualquier punto intermedio de la viga, como se muestra en el CUADRO Nº 7, incluyendo una
posible solución particular.

CUADRO Nº 7 VIGAS DE LONGITUD INTERMEDIA (1.5 < L/λ < 5.0)

(Para cualquier punto interior: u = x / λ t)
Solución
Particular Factor Ct St
φt = φto 1 F1 F2
mt = mto GJ/λ t F 2 F 1

2.2.2 Constantes de Integración

Las constantes de integración F1 , F2 se obtienen de reemplazar en (2.5a,b) los valores de

las funciones hiperbólicas en los extremos de la viga:

así: en el extremo izquierdo (1) con u = 0: φt1 = φt10 +F1 (2.6a)

mt1 = mt10 + (GJ/λt ).F2 (2.6b)

en el extremo derecho(2) con u = L/λ t : φt2 = φt20 +F1.Ct + F2.St (2.6c)

mt2 = mt20 +(GJ/λt ) (F1.St + F2.Ct ) (2.6d)

de (2.6a , 2.6c) se deducen las constantes de integración que se muestran en el Cuadro Nº 8,

incluyendo las posibles soluciones particulares: φt10 , mt10 en el extremo (1) y φt20 , mt20 en el
extremo (2).

28 Fundaciones Elásticas
 CUADRO Nº 8 VIGAS DE LONGITUD INTERMEDIA (1.5 < L/λ < 5.0)
CONSTANTES DE INTEGRACION (u = L / λ t)

φt1 – φt10 φt2 – φt20

F1= 1 0
F2 = -Ct / St 1 / St

De (2.6b) , (2.6d) con el CUADRO Nº 8, se obtiene:

mt1  mt10  kt (t1  t10 )  at (t 2  t 20 ) (2.7a)

mt 2  mt 20  at (t1  t10 )  kt (t 2  t 20 ) (2.7b)

GJ Ct GJ 1
donde las rigideces a torsión, son: kt  ; at  (2.8a,b)
t St t St

2.3 Vigas de Longitud Indefinida

Según lo expuesto en la Sección 1.4 respecto a la propagación de solicitaciones y

deformaciones entre los extremos de una viga de longitud indefinida sobre suelo elástico, deben
eliminarse los términos con (eu) en las funciones hiperbólicas (2.5a,b) de modo que para
cualquier abscisa adimensional (u = x/λ t), incluyendo la solución particular (φot , mot), se obtiene:
t  ot  B1.eu (2.9a,b)
GJ
mt  mot  .B1.e u
t
La constante de integración B1 , depende del tipo de carga.

2.3.1 Momento Torsor Aislado

En la Fig. 2.2(a) se muestra una viga de longitud indefinida con un momento torsor Mt
aplicado en un punto alejado de los extremos y elegido como origen de las abscisas. En las Figs
2.2(b,c) se esquematizan los diagramas de rotaciones y momentos.

Viga con momento torsor aislado

Rotaciones (φt) diagrama simétrico.

Momentos (mt) antimétrico: mt = - Mt/2

Fig 2.2

Fundaciones Elásticas 29
 Para este tipo de carga la solución particular es nula de modo que φot = mot = 0 , y el
momento a la derecha del origen para (u = 0), según (2.9b), es:

Mt GJ
  .B
2 t 1
.M t
y la constante de integración B1, resulta: B1  (2.10)
2.GJ

con lo que las características de la viga en cualquier punto de abscisa adimensional (u = x/λ t) a la
derecha del origen quedan:
 .M
t  t t .eu (2.11a)
2GJ

M t u
mt   .e (2.11b)
2
las características a la izquierda del origen se deducen por simetría o antimetría según el caso.

2.4 Solución Particular en Vigas de Longitud Intermedia

2.4.1 Momentos de Torsión Aislados

El procedimiento es similar al descrito en 1.5.1 para la flexión de vigas con cargas

parciales. En primer lugar se considera una viga solicitada por momentos torsores aislados, sin
restricciones en los extremos. Para todos y cada uno de los momentos se calculan las
características (φto , mto) en las secciones escogidas que se superponen en una solución total. Se
utiliza como solución particular la que corresponde a una viga de longitud indefinida con el
origen en el punto de aplicación del momento y abscisas hacia la derecha utilizando el CUADRO
Nº 7. Los valores de la zona izquierda del origen se determinan por simetría o antimetría, según
el caso. Nuevamente se llama la atención sobre la nomenclatura utilizada para los extremos (φt10 ,
mt10) en el izquierdo (1) y (φt20 , mt20) en el derecho (2).

En segundo lugar, la viga real exige rotaciones nulas en los apoyos, por lo que con los
valores de φt10 , φt20 , haciendo φt1 = mt1 = φt2 = mt2 = 0 porque la solución homogénea es nula,
se calculan las constantes de integración con el CUADRO Nº 8 y las características adicionales
con el CUADRO Nº 7 que permiten determinar los momentos de empotramiento.

Ejemplo Nº 5

En la figura anexa se muestra una viga con dos momentos torsores aislados, aplicados en
las abscisas indicadas, con cinco secciones, incluidas las de los apoyos, escogidas a conveniencia.
En el cuadro anexo (X) representa la abscisa acumulada de cada sección a partir del apoyo
izquierdo, mientras que (x) no mostrada en el cuadro, representa la abscisa de una sección
respecto al origen de la carga en estudio. En la primera columna del cuadro se muestra el
momento (Mt) con el Nº de la sección donde está aplicado.

30 Fundaciones Elásticas
Características de la Viga:
Módulo elasticidad transversal: * G = 84 tf/cm2 ; Luz L= 6.00 m ; Balasto β = 2.0 kgf/cm3
sección: b = 0.60 m ; h = 0.40 m ; base r = 0.60 m ; * J = 75.1 dc

* según ANEXO A: GJ = 6.310 x 103 tf.m2 , β.r3 = 0.432 x 103 tf/m2 , λt = 13.24 m

Mt3= 40.0 Mt4= 25.0

|------------- |--------------|---------- --------|------------------|
S1 S2 S3 S4 S5

X 0.00 1.50 3.00 4.50 6.00

Mt3 33.4544 37.4678 41.9628 37.4678 33.4544
Mt4 18.6693 20.9090 23.4174 26.2267 23.4174
φto (φt10)52.1237 58.3768 65.3802 63.6945 (φt20)56.8718
φtf 0.000 0.000
Mt3 15.9448 17.8577 20.00/-20.00 -17.8577 -15.9448
Mt4 8.8980 9.9655 11.1610 12.50/-12.50 -11.1610
Mto (mt10)24.8428 27.8232 31.1610/-8.8390 -5.3577/-30.3577 (mt20) -27.1058
Mtf 25.5510 -37.9703

Los números en negrilla representan los valores de cada item en el origen de las abscisas
parciales o sea la sección donde está aplicada la carga respectiva. El momento presenta dos
valores en las secciones donde están aplicados, que corresponden a la izquierda y derecha de la
sección.

Con los valores φt10 , φt20 se obtienen, del Cuadro Nº 8, las constantes de integración:
F1 = -52.1237 , F2 = 1.4860 y del Cuadro Nº 7 los valores finales de φtf y Mtf en los apoyos.

2.4.2 Momento Torsor Uniforme en la Viga

Por la acción de un momento torsor uniforme (mtu) la rotación es constante porque el

momento es constante. La solución particular debe cumplir con la ecuación completa (2.1d) que
se escribe en la siguiente forma:
d 2  .r 3 m
2
   tu  0
dx 12GJ GJ

la solución particular (φto) también es constante, la diferencial se anula y por lo tanto:

m . 2
to  tu t (2.12)
GJ
que, por ser el único caso de carga: φt10 = φt20 = φto mientras φt1= φt2 = mt10 = mt20 = 0

Con las constantes de integración del CUADRO Nº 8 aplicadas al CUADRO Nº 7 se

obtiene el momento de empotramiento del apoyo izquierdo con u = 0 y del apoyo derecho con
u = L/λt , tales que:
C 1
mtf 1  mtf 2  t t .mtu (2.13)
St

Fundaciones Elásticas 31
 2.5 Ecuación matricial de una viga:

En la Fig. 2.3 se muestra una viga (m) entre las juntas i, j con la nomenclatura de
solicitaciones y deformaciones escogidas en el siguiente orden: cargas en el eje Z, momentos en
los planos ZOX , ZOY que se asocian con la secuencia normalmente observada en el cálculo de
edificios, con la convención positiva de los ejes de referencia. Obsérvese que, por ejemplo, la
acción sísmica en la dirección X produce momentos en el plano ZOX cuya representación
vectorial coincide con el eje Y.

Fig. 2.3 Coordenadas Locales (q,d)

Añadiendo la torsión con la nomenclatura descrita, la ecuación (1.17) de la viga queda:

(3i-2) (3i-1) (3i) (3j-2) (3j-1) (3j) * numeración de las casillas de la matriz
q1 t b 0 -to bo 0 d1
q2 b k 0 -bo a 0 d2
q3 = 0 0 kt 0 0 -at * d3
q4 -to -bo 0 t -b 0 d4
q5 bo a 0 -b k 0 d5
q6 0 0 -at 0 0 kt d6
qm = km * dm (2.14)

de la que, por comodidad, se recopila la nomenclatura utilizada hasta ahora.

2EI CS  cs 2EI S 2  s 2 4EI CS  cs 4 EI

k b t  4
 S 2  s2  2 S 2  s2  3 S 2  s2  .r

2EI Cs  Sc 4 EI Ss 4 EI Cs  Sc
a bo  to 
 S 2  s2  2
S  s2
2
 3 S 2  s2

GJ Ct GJ 1 12GJ
kt  at  t 
t St t St  .r 3

s = sen (L/λ) c = cos (L/λ) S = senh(L/λ) C = cosh(L/λ)

S t = senh(L/λ t) Ct = cosh(L/λ t)

32 Fundaciones Elásticas
 2.6 Vigas Esviadas:

Son vigas que no están alineadas con los ejes de referencia. Aunque el problema es
conocido [7] se presenta un resumen de la solución porque la numeración escogida de las
coordenadas locales difiere de la comúnmente utilizada lo que modifica la presentación de
algunas matrices.

En la Fig 2.4 se muestra una viga (m) entre las juntas (i , j) esviada el ángulo α respecto al
eje X, con las coordenada locales (q , d) y las globales (Q , D), cuyas relaciones (Q , q) y (D , d)
tienen igual representación. Se denomina Cx , Cy a los cosenos directores de la inclinación de la
viga con los ejes X,Y.

Fig 2.4 Coordenadas Locales (q) y Globales (Q)

Según el método convencional de proyectar (q) en (Q) se obtienen las siguientes

relaciones:
Qm= Rtm .qm (2.15)

Q1 Rte 0 q1
Qi = Q2 Rtm = qi = q2
Q3 0 Rte q3
Qm = qm =
Q4 1 0 0 q4
Qj = Q5 Rte = 0 Cx Cy qj = q5
Q6 0 -Cy Cx q6

Los desplazamientos (d,D) tienen una representación similar a la (2.15), de modo que:

dm = Rm .Dm (2.16)

reemplazando (2.14) y (2.16) en (2.15) se obtiene:

Qm = Rtm .km .Rm .Dm (2.17)

Fundaciones Elásticas 33
y la matriz de transformación Km = Rtm .km .Rm , expandida resulta:

(3i-2) (3i-1) (3i) (3j-2) (3j-1) (3j)

t b Cx -b Cy -to bo Cx -bo Cy
k Cx2+kt Cy2 (kt – k) Cx Cy -bo Cx a Cx2 – at Cy2 -(a +at ) Cx.Cy
kt Cx2 + k Cy2 bo Cy -(a +at) Cx.Cy a Cy2 – at Cx2
t -b Cx b Cy
Simétrica k Cx2 + kt Cy2 (kt – k) Cx.Cy
k Cy2 + kt Cx2

(2.18)
2.7 Ecuación Matricial de un Retículo

El equilibrio de un retículo solicitado por cargas normales a su plano se representa con la

ecuación: Q = K.D cuya solución es: D = K-1.D , donde:

K = matriz de rigidez del retículo formada con el aporte de la matriz Km de cada

miembro siguiendo la numeración de las casillas correspondientes.
Q = vector de cargas generalizadas ; D = vector de desplazamientos generalizados.

Para cada miembro (m) se determinan los desplazamientos generalizados Dm de las juntas
respectivas y, con (2.16) los desplazamientos locales dm = Rm . Dm , con lo que finaliza el uso
del método de los desplazamientos y se retorna a la convención de signos de resistencia de
materiales procesando cada viga independientemente, aplicando los CUADROS Nº 2 y Nº 1 para
las características a flexión y los CUADROS Nº 8 y Nº 7 para el cálculo de la torsión, cambiando
el signo del momento de torsión en la primera junta de cada viga, de modo que la interpretación
de las características, según la resistencia de materiales, es la que se muestra en la Fig. 2.5

Fig. 2.5 Interpretación de las Características de una Viga

Según la Teoría de Resistencia de Materiales

34 Fundaciones Elásticas
2.8 Momento Máximo de Tramo

Para el cálculo de las características en secciones intermedias de cada viga se eligen los
cuartos de la luz. Esta elección puede ser insuficiente para obtener el verdadero momento
máximo del tramo, -entendiéndose por momento máximo el mayor valor absoluto del momento
calculado-, por lo que habría que adoptar intervalos mas cortos. Si se aplica el criterio según el
cual la ubicación del momento máximo coincide con la del corte nulo, la determinación de la
abscisa correspondiente se puede encontrar por aproximaciones sucesivas (Newton, Regula Falsi,
etc) según las relaciones (1.1a,c)

Para las cargas previstas, el diagrama de cortes se asemeja al triangular mostrado en la

Fig. 2.6 con los signos de la Fig. 2.5. Si se cumple la relación (V2/V1)<0 (Caso a), se encuentra
la abscisa Ln del corte nulo en una primera aproximación, que permite determinar el segmento de
viga (i-j) donde se anula el corte (Caso b), en el cual también debe cumplirse la condición
(Vj/Vi)<0 que fija la abscisa Li para el cálculo del momento máximo de tramo que puede ser
negativo o positivo.

Si (V2/V1) > 0 , el diagrama de momentos no tiene puntos de inflexión y es siempre

creciente o decreciente. Las secciones de cálculo corresponden a los cuartos de la Luz.

Caso (a): (V2/V1)< 0 Caso (b): (Vj/Vi)< 0

L L
Ln  Li  x i 
V
1- 2  V 
4 1- j 
V1  Vi 

Fig. 2.6 Cálculo del Momento Máximo de Tramo

Fundaciones Elásticas 35
2.9 Programa FUNDARETIC

Las limitaciones expuestas en el art. 1.7 son aplicables a este programa, excepto que el
módulo de balasto se considera constante.

2.9.1. Preparación de Datos

La numeración de juntas y miembros de retículos es de libre elección del usuario. Sin

embargo, para facilitar la lectura de resultados y realizar el diseño de secciones comunes a dos
tramos contiguos se recomienda lo siguiente:

-Numere las juntas de izquierda a derecha en cada fila paralela al eje X empezando por la
más baja en el plano del esquema.

-Numere los miembros de izquierda a derecha en cada fila paralela al eje X empezando
por la mas baja.

-Numere los miembros de abajo hacia arriba en cada eje paralelo al Y empezando por el
eje ubicado mas a la izquierda.

-Declare como dirección positiva de cada miembro la que resulte más cómoda para la
lectura de resultados, tomando en cuenta que las abscisas de las secciones se miden a partir de la
primera junta del miembro.

-El peso propio de las vigas puede incluirse automáticamente.

-El momento máximo de tramo se calcula según lo indicado en la sección

En la figura anexa se muestra un ejemplo de la numeración de juntas y miembros, así

como la convención positiva y el orden de declaración de las solicitaciones.

Numeración de Juntas y Miembros Solicitaciones

36 Fundaciones Elásticas
2.9.2 Codificación de Datos

Cada línea de datos debe estar precedida de la palabra DATA. La última línea del
programa, incluyendo los datos, debe ser Fin:END. En el programa anexo se mantienen los datos
de los Ejemplos Nº 6 y Nº 7 como muestra de codificación.

1. Datos Generales

DATA A$ , Ec , Gc , Bsto , Nj , Nvt , Ntipo , Nq

A$ : identificación del trabajo hasta 50 caracteres alfanuméricos
Ec : módulo de elasticidad longitudinal ( tf/cm2)
Gc : módulo de elasticidad transversal ( tf/cm2)
Bsto: módulo de balasto ( kgf/cm3)
Nj = Nº de juntas <= 120 , Nvt = Nº total de viga <= 120
Ntipo = Nº tipos de viga <= 50 , Nq = Nº casos de carga = (arbitrario)

2. Datos de Geometría (ver ANEXO A)

Para cada junta : DATA J , Xo , Yo = Nº de la junta , abscisa , ordenada (m)

Para cada tipo de viga:

Sección rectangular : DATA 1 , Bo , Ho = base, altura (m)
Sección T invertida: 2 , Bo , Ho, Bf , Hf = (base, altura) nervio, ala (m)
Sección arbitraria : 3 , Bj , Io , Jo , A = base (m) , Inercia flexión (dcm4)
Inercia torsión (dcm4) , Area (dcm2)

Para cada viga : DATA J , J1 , J2 , Typ

J = Nº de la viga ; J1 , J2 = Nº de las juntas de la viga
Typ = Nº del tipo de viga

3. Datos de Cargas

Para cada caso de Carga:

DATA Njq = Nº de juntas con cargas
Si Njq > 0 , Para cada junta cargada: J , Pjta1 , Pjta2 , Pjta3
J = Nº de la junta , Pjta1 = carga (tf) , Pjta2 = My , Pjta3 = Mx (tf.m)

DATA Npp = código para cálculo automático de peso propio.

Npp > 0 calcula el peso propio con 2.5 tf/m3

DATA Ntr= Nº de vigas con carga uniforme diferente al peso propio.

Si Ntr > 0 , Para cada viga cargado: J , Pt
J = Nº del tramo , Pt = carga uniforme (tf/m)

Fundaciones Elásticas 37
4. Archivo de Resultados

Los resultados de todas las vigas se archivan en el arreglo Rest (5, 8*Nº de vigas). Las 5
filas corresponden a cuartos de la luz y la sección de momento máximo de tramo cuando L>2.0
m; si L<=2.0 m las secciones corresponden a 0, L/2, L, y, en cada fila en grupos de 8, a la
abscisa, flecha, giro por flexión, momento de flexión, corte, giro por torsión, momento de torsión,
presión, correspondientes a cada sección.

100 ! FUNDARETI C Ret ícu lo s de Fu ndac ió n E lást ica.

110 !- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -- - -- - -- - -
120 OP T I ON BAS E 1
130 DIM A$[50], L$[130], Tit1$[35], Tit2$[100], Tit3$[85], Tit4$[120], Tit5$[100], Tit6$[120]
140 DIM Tit7$[80], D(120), F(120,120), K(120,120), Mb(120,15), Mtyp(50,15), Pjta(120,3)
150 DIM Ptra(120,2), Qz(60,5), Rest(120,10), U(120), Xo(120), Yo(120), INTEGER Jn(120)

160 L$=RP T $( “- “, 130)

170 RAD ! Dec lar a Rad i an es
180 RE S T ORE E je mp lo 6
185 ! RE S T ORE E je mp lo 7 ! E sco ge el pr o ble ma
190 RE AD A$, E c, Gc, Bst o
200 RE AD N j, Nvt , Nt ipo , Nq
210 Nd j = 3 * Nj

220 RE DI M D( Ndj) , F( Nd j, Nd j) , K( Nd j, Nd j ) , Mb( Nvt , 15) , Mt yp( N t ipo , 15)

230 RE DI M P jt a( N j, 3) , Pt r a ( Nvt , 2 ) , Rest ( Nvt *8, 10) , U( Ndj) , Xo( Nj) , Yo ( Nj)
240 MAT K=( 0)
250 M AT Mt yp=( 0)
260 MAT Mb=( 0)
270 M AT Xo =( 0)
280 M AT Yo =( 0)
290 M AT Jn=( 0)

300 FOR I =1 T O Nj
310 RE AD J, Xo ( J) , Yo ( J) ! Coor denadas de la s Ju nt as
320 NE XT I

330 FOR I =1 T O Nt ipo

340 B f=0
350 H f=0
360 RE AD J, T ipo ! Nº T ipo _V iga , Nº T ipo _ secc
370 ON T ipo GOT O S _1, S _2, S _3

380 ! S ecc ió n Rect angu lar

390 S _1: RE AD Bo , Ho
400 A= Bo *Ho
410 I o =Bo *Ho ^3/ 12
420 B j= Bo
430 P p=A*2. 5 ! P eso pro pio (t f/ m)
440 Z i=Bo / Ho
450 I F Z i > 1 T HE N
460 Jt =Bo *Ho ^3
470 Zi=1/ Z i
480 E LS E

38 Fundaciones Elásticas
490 Jt =Ho *Bo ^3
500 E ND I F
510 Jt =Jt *( 1/ 3 – 0. 21*Zi*( 1 - Z i^ 4 / 12) )
520 GOT O S _4

530 ! S ecc ió n T inver t ida

540 S _2: RE AD Bo , Ho , Bf, H f ( Ner vio , Ala )
550 Bh=B f- Bo
560 A= Bo *Ho +Bh*H f
570 M=( Bo *Ho ^2+Bh*H f^2) / 2
580 I o =( Bo *Ho ^3+Bh*H f^3) / 3 - M^2/ A
590 B j= B f
600 P p=A*2. 5
610 H i=Ho - H f
620 Z i=Bo / H i
630 I F Z i >1 T HE N
640 Z i=1/ Z i
650 Jt =( 1 - 0. 63*Z i) *B o *H i^3
660 E LS E
670 Jt =( 1 - 0. 63*Z i) *H i* Bo ^3
680 E ND I F

690 Z i=B f/ H f
700 I F Z i >1 T HE N
710 Z i=1/ Z i
720 Jt =( Jt +( 1 - 0. 63*Zi) *B f*H f^3) / 3
730 E LS E
740 Jt =( Jt +( 1 - 0. 63*Zi) *H f* B f^3) / 3
750 E ND I F
760 GOT O S _4

770 ! S ecc ió n Ar bit r ar ia

780 S _3: RE AD B j, I z, Jo , A ! Bj( m) , Iz , J0 ( dc m4) , A ( dcm2)
790 I o =I z*1. E- 4 ! I. F le xió n ( m4)
800 Jt =Jo *1.E - 4 ! J.To r sió n ( m4)
810 P p=A*2. 5*1. E - 2 ! P eso P ro pio (t f/ m)
820 S _4: E 1=10*E c*Io ! Rig ide z. fle xió n / 1000
830 E 2=10*Gc*Jt ! Rig ide z. t o r sió n/ 100 0
840 L b=( 4*E 1/ ( B j*B st o ) ) ^ . 25
850 L bt =( 12*E 2/( Bst o *Bj ^ 3) ) ^. 50
860 Mt yp( J, 1) =Bo
870 Mt yp( J, 2) =Ho
880 Mt yp( J, 3) =B f
890 Mt yp( J, 4) =H f
900 Mt yp( J, 5) = E 1
910 Mt yp( J, 6) =B j
920 Mt yp( J, 7) =L b
930 Mt yp ( J, 8) =E 2
940 Mt yp ( J, 9) =L bt
950 Mt yp( J, 10) =P p
960 Mt yp( J, 11) =I o
970 Mt yp( J, 12) =Jt
980 NE XT I

990 FOR I =1 T O Nvt

Fundaciones Elásticas 39
1000 RE AD J, J1, J2, T yp ! N. V iga, Junt as, N. T ipo . S ecc
1010 Dx=Xo ( J2) - Xo ( J1)
1020 D y=Yo ( J2) - Yo ( J1)
1030 L=S QR( D x^2+ DY^2)
1040 Cx=D x/ L
1050 C y=D y/ L

1060 Mb( J, 1) =J1

1070 Mb( J, 2) =J2
1080 Mb( J, 3) =L
1090 Mb( J, 4) =T yp
1100 Mb( J, 5) = Mt yp( T yp , 5)
1110 Mb( J, 6) = Mt yp( T yp , 7)
1120 Mb( J, 7) =Bst o
1130 Mb( J, 8) =Mt yp( T yp , 8)
1140 Mb( J, 9) =Mt yp( T yp , 9)
1150 Mb( J, 10) =Mt yp( T yp , 10)
1160 Mb( J, 11) =Cx
1170 Mb( J, 12) =C y
1175 Mb( J, 13) =Mt yp( T yp, 6)
1180 NE XT I

1190 FOR I =1 T O Nvt

1200 L i=M b( I , 3)
1210 E 1=Mb( I , 5)
1220 L b=M b( I , 6)
1230 E 2=Mb( I , 8)
1240 L bt =Mb( I , 9)
1250 Cx=M b( I , 11)
1260 C y=M b( I , 12)
1270 GOS UB S _1 5
1280 Ki=2*E 1/ ( L b*( K6 - K7) )
1290 Ko =Ki*( K5 - K8)
1300 A= Ki*( K3 - K2)
1310 B= Ki*( K6 + K7) / L b
1320 Bo =2*Ki* K4/ L b
1330 T =2*Ki*( K5 + K8) / L b^2
1340 T o =2*Ki*( K2 + K3) / L b^2
1350 Kt =E 2*Cht / ( Lbt *S ht )
1360 At =Kt / Cht

1370 I 1=3*Mb( I , 1) - 2
1380 I 2=I 1+1
1390 I 3=I 2+1
1400 J1=3*M b( I , 2) - 2
1410 J2=J1+1
1420 J3=J2+1

1430 K( I 1, I 1) =K( I 1, I 1) +T ! T r iá ngu l o S uper io r de la Mat r iz de R ig idez

1440 K( I 1, I 2) =K( I 1, I 2) + B*Cx
1450 K( I 1, I 3) = K( I 1, I 3)- B*C y
1460 K( I 1, J1) = K( I 1, J1) - To
1470 K( I 1, J2) = K( I 1, J2) +Bo *Cx
1480 K( I 1, J3) = K( I 1, J3) - Bo *C y
1490 K( I 2, I 2) =K( I 2, I 2) +Ko *Cx^2+ Kt *C y^2

40 Fundaciones Elásticas
1500 K( I 2, I 3) = K( I 2, I 3) +( Kt - Ko ) *Cx*C y
1510 K( I 2, J1) = K( I 2, J1) - Bo *Cx
1520 K( I 2, J2) = K( I 2, J2) +( A* C x^2 - At *C y^2)
1530 K( I 2, J3) = K( I 2, J3) - ( A+ At ) *Cx*C y
1540 K( I 3, I 3) = K( I 3, I 3) +Ko *Cy^2+ Kt *Cx^2
1550 K( I 3, J1) = K( I 3, J1) +Bo *C y
1560 K( I 3, J2) = K( I 3, J2) - ( A+ At ) *Cx*C y
1570 K( I 3, J3) = K( I 3, J3) +( A* C y^2 - At *Cx^2)
1580 K( J1, J1) = K( J1, J1) +T
1590 K( J1, J2) = K( J1, J2) - B*C x
1600 K( J1, J3) = K( J1, J3) +B* C y
1610 K( J2, J2) = K( J2, J2) + Ko *C x^2+ Kt *C y^2
1620 K( J2, J3) = K( J2, J3) +( Kt - Ko ) *Cx*C y
1630 K( J3, J3) = K( J3, J3) + Ko * C y^2+ Kt *C x^2
1640 NE XT I

1650 FOR I =1 T O Ndj ! Co mp let a la Mat r iz

1660 FOR J=I +1 T O Nd j
1670 I F K( J, I ) =0 T HE N
1680 K( J, I ) =K( I , J)
1690 E LS E
1700 K( I , J) =K( J, I )
1710 E ND I F
1720 NE XT J
1730 NE XT I

1740 M AT F= I NV( K)
1750 FOR I = 1 T O Nd j ! M ini m iz a er ro r es de r edo ndeo
1760 FOR J=2 T O Nd j
1770 F( I , J) =( F( I , J) +F( J, I ) ) / 2
1780 F( J, I ) =F( I , J)
1790 NE XT J
1800 NE XT I

1820 Nq1=0 ! CAS OS DE CARG A

1830 S _5 : Nq1=Nq1+1
1840 MAT U=( 0)
1850 M AT P jt a=( 0)
1860 RE AD N jq ! Nº Junt as car gadas
1870 I F N jq>0 T HE N
1880 FOR I =1 T O Njq
1890 RE AD J, P jt a( J, 1) , P jt a( J, 2) , P jt a( J, 3) ! N. ju nt a , Q , M y , Mx
1900 U( 3*J- 2) =P jt a( J, 1)
1910 U( 3*J- 1) =P jt a( J, 2)
1920 U( 3*J) =P jt a( J, 3)
1930 NE XT I
1940 E ND I F
1950 M AT Pt r a=( 0)
1960 RE AD Npp ! Có digo de peso pro pio aut o mát ico
1970 I F Npp>0 T HE N
1980 FOR I =1 T O Nvt
2090 P t r a( I, 1) =Mb( I , 10)
2000 NE XT I
2010 E ND I F

Fundaciones Elásticas 41
2020 RE AD N t r ! Nº de T r amo s Car gado s
2030 I F Nt r >0 T HE N
2040 FOR I =1 T O Nt r
2050 RE AD J, P t r a( J, 2) ! N. t r amo , p. unifo r me ( t f/ m)
2060 NE XT I
2070 E ND I F

2080 FOR I =1 T O Nvt

2090 Pt =Pt r a( I, 1) +Pt r a( I, 2) ! Acu mu la car ga s en e l t r a mo
2100 I F P t >0 T HE N
2110 J1=Mb( I , 1)
2120 J2=Mb( I , 2)
2130 Lb=M b( I , 6)
2140 Cx=M b( I , 11)
2150 Cy=M b( I , 12)
2160 Gi= M b( I , 3) / Lb
2170 Ch=( E XP ( G i) +E XP ( - G i) ) / 2
2180 S h=( E XP ( G i) - E XP ( - Gi) ) / 2
2190 C1=COS ( G i)
2200 S 1=SI N ( Gi)
2210 M1=Pt *Lb^2*( S h- S 1) / ( 2*( S h+S 1) )
2220 M2=- M1
2230 V1=Pt *Lb*( C h- C1) / ( S h+S 1)
2240 V2=V 1

2250 U( 3*J1 - 2) =U( 3*J1 - 2) +V1

2260 U( 3*J1 - 1) =U( 3*J1 - 1) +M1*Cx
2270 U( 3*J1) =U( 3*J1) - M1*Cy
2280 U( 3*J2 - 2) =U( 3*J2 - 2) + V2
2290 U( 3*J2 - 1) =U( 3*J2 - 1) +M2*Cx
2300 U( 3*J2) =U( 3*J2) - M2*Cy
2310 E ND I F
2320 NE XT I

2330 MAT D=F*U

2340 MAT Qz=( 0

2350 FOR I =1 T O Nvt
2360 J1=M b( I , 1)
2370 J2=M b( I , 2)
2380 L=M b( I , 3)
2390 T yp=M b( I , 4)
2400 E 1=Mb( I , 5)
2410 L b=M b( I , 6)
2420 Bst o =Mb( I , 7)
2430 E 2=Mb( I , 8)
2440 L bt =Mb( I , 9)
2450 Cx=M b( I , 11)
2460 C y=M b( I , 12)
2470 B j= M b( I , 13)
2480 Wo =0
2490 P t =Pt r a( I, 1) +Pt r a( I, 2)
2500 I F P t >0 T HE N Wo =Pt / ( Bj*B st o) ! Wo = mm
2510 D1=D( 3*J1 - 2) - Wo
2520 D2=D( 3*J1 - 1) *Cx- D( 3*J1) *C y

42 Fundaciones Elásticas
2530 D3=D( 3*J2 - 2) - Wo
2540 D4=D( 3 *J2 - 1) *Cx- D( 3*J2) *C y
2550 Dt 1=D( 3*J1 - 1) *Cy+D( 3*J1) *C x
2560 Dt 2=D( 3*J2 - 1) *Cy+D( 3*J2) *C x
2580 L i=L
2590 GOS UB S _1 5

2600 Dn= K6 - K7
2610 A1=D1 ! Co nst ant e s de I nt egr ació n
2620 A2=( - D1*( K5 + K8) - D2*L b* K7+D3*( K2+ K3) - D4*L b* K4) / D n
2630 A3=( D1*( K5+ K8) +D2* L b* K6 - D3*( K2+ K3) +D4* L b* K4) / Dn
2640 A4=( - D1*( K6 + K7) - D2*L b*( K5 - K8) +D3 *2*K4+D4* L b*( K2 - K3) ) / D n
2650 At 1=Dt 1
2660 At 2=Dt 2/ S ht - Dt 1* Cht / S ht
2670 S k=2*Mb( I , 5) / Lb^2
2680 S t =Mb( I , 8) / Lbt
2690 Z1=3
2700 Dx=0. 50
2710 I F L>2. 0 T HE N
2720 Z1=5
2730 Dx=0. 25
- -- - -- - -- - M(- ) má xi mo de t r amo - -- - -- - -- -
2740 FOR L i=0 T O L S TE P L
2750 GOS UB S _15
2760 V2=( ( A2 - A3) * K1+( A1 - A4) * K2+( A1+ A4 ) *K3+( A2+ A3) * K4) *S k/ L b
2770 I F L i=0 T HE N V1=V2
2780 NE XT L i
2790 I F V2/ V1 <0 T HE N
2800 L n=L/ ( 1 - V2/ V1)
2810 I i=1+I NT ( 4*Ln/ L)

2820 S _ 10 : J i=I i+1

2830 L1=L*( I i- 1) / 4
2840 L2=L1+ L/ 4
2850 FOR L i= L1 T O L2 S T EP L/ 4
2860 GOS UB S _1 5
2870 V2=( ( A2 - A3) * K1+( A1 - A4) * K2+( A1+ A4) * K3+( A2+ A3) * K4) *S k/ L b
2880 I F Li=L1 T HE N V1=V2
2890 NE XT Li
2900 I F V2/ V1>0 T HE N
2910 I F ABS ( V1 ) >ABS ( V2) T HE N I i= I i+1
2920 I F ABS ( V1 ) <ABS ( V2) T HE N I i=I i- 1
2930 GOT O S _10
2940 E ND I F
2950 L n=L1+ L/ ( 4*( 1 - V2/ V1) )
2960 I F J i<=3 T HE N Jn( I ) =5*I - 5+Ji
2970 I F J i>3 T HE N Jn( I ) =5*I - 5+I i
2980 E ND I F
2990 E ND I F

3000 FOR Z =1 T O Z1
3010 J=5*I- 5+Z
3020 Li=( Z- 1) *D x*L
3030 I F J=J n( I ) T HE N Li= L n
3040 GOS UB S _1 5

Fundaciones Elásticas 43
3050 Rest ( J, 1) =Li
3060 Rest ( J, 2) =Wo +A1*K1 + A2* K2+ A3* K3+ A4* K4
3070 Rest ( J, 3) =( ( A2+ A3) *K1 +( A1+ A4) *K2+ ( A4 - A1) * K3+( A3 - A2) * K4) / L b
3080 Rest ( J, 4) =( - A4*K1 - A3* K2+ A2* K3 +A1* K4) *S k
3090 Rest ( J, 5) =( ( A2 - A3) * K1+( A1 - A4) * K2+( A1+ A4) * K3+( A2+ A3) * K4) *S k/ L b
3100 Rest ( J, 6) =At 1*Cht +At 2*S ht
3110 Rest ( J, 7) =- St *( At 1*S ht +At 2*Cht )
3120 Rest ( J, 8) =Rest ( J, 2) *Bst o

3130 I F Z=1 T HE N ! S uma de Car gas/ Ju nt a

3140 Qz( J1, 1) =Qz( J1, 1) - Rest ( J, 5)
3150 Qz( J1, 2) =Qz( J1, 2) +Rest ( J, 4) *Cx+Rest ( J, 7) *C y
3160 Qz( J1, 3 ) =Qz( J1, 3) - Rest ( J, 4) *Cy+Re st ( J,7) *Cx
3170 E LS E
3180 I F Z=Z1 T HE N
3190 Qz( J2, 1) =Qz( J2, 1) + Rest ( J, 5)
3200 Qz( J2, 2) =Qz( J2, 2) - ( +Rest ( J, 4) *Cx+ Rest ( J, 7) *C y)
3210 Qz( J2, 3) =Qz( J2, 3) - ( - Rest ( J, 4) *Cy+ Rest ( J, 7) *Cx)
3220 E ND I F
3230 E ND I F
3240 NE XT Z
3250 N E XT I

3260 GOS UB S _20

3270 I F Nq1=1 T HE N
3280 GOS UB Exi t _1
3290 GOS UB Exi t _2
3300 E ND I F
3301 GOS UB Exi t _ 3
3310 I NP UT “I MP RE S I ON COMP LE T A: S / N”, S $
3320 I F S $=”S ” T HE N
3330 GOS UB Exi t _4
3340 E LS E
3350 GOS U B Exi t _5
3360 E ND I F
3370 I F Nq1< Nq T HE N S _5
3380 GOT O Fi n

3390 ! Func io nes

3400 S _15 : G=L i/ L b
3410 Ch=( E XP ( G) +E XP( - G)) / 2
3420 S h=( E XP ( G) - E XP( - G)) / 2
3430 C1=COS ( G)
3440 S 1=S I N( G)
3450 K1= C1*C h
3460 K2= C1*S h
3470 K3=S 1*C h
3480 K4=S 1*S h
3490 K5= C h*S h
3500 K6=S h*S h
3510 K7=S 1*S 1
3520 K8= C1*S 1
3530 Gt =L i/ L bt
3540 Cht =( E XP ( Gt ) +E xp( - Gt )) / 2
3550 S ht =( E XP ( Gt ) -E XP ( - Gt ) ) / 2

44 Fundaciones Elásticas
3560 RE T URN

3570 S _20: ! P ara im pri n i r : L os pu n t os ( …. .) in di cados en l os t í tu l os , represe n t an

espaci os en bl an co para su com pagi n aci ón con el t ext o .
3580 Tit1$=”FUNDACION ELASTICA: Balasto= “
3590 Tit2$=”Unidades: m , tf , desplazamientos = mm , rotaciones = rad*1000”
3600 Tit3$=”JUNTA…….X…….Y..*.VIGA….J1….J2..Tipo……..L……….CX……….CY”
3610 Tit4$[1, 88]=”V-TIPO…Bo…..Ho…..Bf.....Hf.*..Base.*.....Io.(m4)...Lambda…..Jp.(m4)..Lambda_T..*“
3620 Tit5$[1,100]=”Junta….....Q.......MY.......MX.*.Viga.....pt.*.........DZ.........RY.........RX.*.Junta..Presión.*”
3625 Tit6$[1,120]=” “
3630 Tit6$[1, 49]=”VIGA...J1…J2……..X....DEFLEXIO....G.FLEXION”
3640 Tit6$[54, 111]=”M.FLEXION……CORTE….G.TORSION….M.TORSION….PRESION”
3650 Tit7$[1, 47]=”VIGA…J1…J2.*….X……Mf…….V……Mt.*”
3660 RETURN
3670 !- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- - - - -- - -- -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - - - - -- - - -- - -- - -- -
3680 Exit_1: PRINT USING “58A, 29A, Z.2D, 10A, /, 100A”; A$, Tit1$, Bsto, “ kgf/cm3”,L$
3690 PRINT USING “70A, 2(4A, 3D.2D), 7A, /, 100A ”; Tit2$, “E =”, Ec, “G =“,Gc, “ tf /cm2”,L$
3700 RETURN
3710 !- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - - - - -- - - -- - -- - -- -
3720 Exi t _2: N ma x= M AX( N j, N vt )
3730 P RI NT US I NG “80 A, / , 80A” ;T it 3 $, L$
3740 FOR I =1 T O Nma x
3750 I F I <=N j T HE N
3760 P RI NT US I NG “#, 5D, 2( 5D. 2D), 3A”; I , Xo( I ), Yo( I ), ”. . *”
3770 E LS E
3780 P RI NT US I NG “#, 23X”;””
3790 E ND I F
3800 I F I <= Nvt T HE N
3810 J1=M b( I , 1)
3820 J2=M b( I , 2)
3830 L=M b( I , 3)
3840 T ipo =M b( I , 4)
3850 Cx=M b( I , 11)
3860 C y=M b( I , 12)
3870 PRINT USING “5D, 3(6D), 6D.2D, 2(4X, SZ.5D)”; I, J1, J2, Tipo, L, Cx, Cy
3880 E LS E
3890 P RI NT
3900 E ND I F
3910 NE XT I
3920 P RI NT US I NG “80 A” ; L$
3930 P RI NT US I NG “/ , 15A, 3( / , 88A) ”; ”T I P OS DE VI GAS ”, L$, T it 4$, L$
3940 FOR I =1 T O Nt ipo
3950 Bo =Mt yp( I , 1)
3960 Ho =Mt yp( I , 2)
3970 B f= Mt yp( I , 3)
3980 H f=Mt yp( I , 4)
3990 Base= Mt yp( I , 6)
4000 L b=Mt yp( I , 7)
4010 L bt =Mt yp( I , 9)
4020 I o =Mt yp( I , 11)
4030 Jp=Mt yp( I , 12)
4040 I f Bo >0 T HE N
4050 IF Bf=0 THEN PRINT USING “#, 4D, 2(4D.2D), 2(5X, A,X),X,A”;I, Bo, Ho,”-“,”-“,”*”
4060 I F Bf>0 P RI NT US I NG “#, 4D, 4(4D. 2 D) , 2A” ; I , Bo , Ho , Bf, H f, “ *”
4070 E LS E

Fundaciones Elásticas 45
4080 P RI NT US I NG “#, 4D, 4( 6X, A) , 2A” ; I , “- “, “- “, “- “, “- “, “ *”
4090 E ND I F
4100 PRINT USING ”3D.2D, 2A, 2(2X,D.5DE, 4D.3D, X), 3A; Base, “ *”, Io, Lb, Jp, Lbt, “ *”
4110 NE XT I
4120 P RI NT US I NG “88 A” ; L$
4130 RE T URN
4140 ! - - - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - --- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - - - - -- - -- - -- - - -
4150 Exit_3 : PRINT USING “/, 5A, 2D,2X,10A, /,100A”; “CASO”,Nq1,”DE CARGAº”,L$
4160 P RI NT US I NG “#, 16X, 16A, 15X, 4A”; “C ARG AS AP LI C AD AS ”, “* “
4170 PRINT USING “49A”; DESPLAZAMIENTOS Y ROTACIONES...*…PRESIONES….*”
4180 P RI NT US I NG “100 A, / , 100A” ; T it 5 $, L$
4190 FOR J=1 T O N j
4200 J1=3*J- 2
4210 P t =Pt r a( J, 1) +Pt r a( J, 2)
4220 P s=D( J1) * Bst o
4230 Q=U( J1)
4240 M y=U( J1+1)
4250 Mx= U( J1+2)
4260 D1=D ( J1)
4270 D2=D ( J1+1)
4280 D3=D ( J1+2)
4290 PRINT USING “#, 5D, 3(6D.2D), 2A, 5D, 4D.2D, 2A”; J, Q, My, Mx, “ *” , J, Pt,” *”
4300 P RI NT US I NG “3( 5D. 5D) , 2A, 6D, 6D . 2D, 2A) ”; D 1, D2, D3, “ *”, J, P s, ” *”
4310 NE XT J
4320 P RI NT US I NG “100 A” ; L$,
4330 RE T URN
4340 ! - - - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - --- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -
4350 Exit_4 : PRINT USING “ /, 5A, 2D, 2X, 10A, /, 110A”; ”CASO”, Nq1, ”DE CARGA:”, L$
4360 P RI NT US I NG “11 0 A, / , 110A” ; T it 6$, L$
4370 FOR I =1 T O Nvt
4380 L=M b( I , 3)
4390 FOR J=5*I - 4 T O 5*I
4400 I F J=5*I - 4 T HE N
4410 P RI NT US I NG “#, 4D, 5D, 5D”; I , Mb( I , 1) , Mb( I , 2)
4420 E LS E
4430 P RI NT US I NG “#, 14X”
4440 E ND I F
4450 I F L<=2. 0 AND J >5*I - 2 T HE N S _25
4460 X1=Rest ( J, 1)
4470 W1= Rest ( J, 2)
4480 R f= Rest ( J, 3)
4490 M f= Rest ( J, 4)
4500 V f=Re st ( J, 5)
4510 Rt =Rest ( J, 6)
4520 Mt =Rest ( J, 7)
4530 P s=Rest ( J, 8)
4540 P RI NT US I NG Format 1 : X1, W1, R f, Mf, V f, Rt , Mt ,P s
4545 Fo rmat 1: 6D. 2 D, 10D. 2D, 2( 7D. 5D, 9D. 3D, 7D. 3D)
4550 S _25 : NE XT J
4560 I F I < Nvt T HE N P RI NT
4570 I F I =Nvt T HE N P RI NT US I NG “110 A;L$
4580 NE XT I
4590 GOS UB Ch ek
4600 RE T URN
4610 ! -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - - -- - -- - -- -

46 Fundaciones Elásticas
4620 Exi t _5 : P RI NT US I NG “/ , 5A, 2D, 2X, 10A, / , 100A” ; “C AS O, Nq1, “DE CARG A ”, L$
4630 P RI NT US I NG “#, 15X, 2A, 7X, 15A, 7X, A, 2X”; ”. * ”, ”AP OYO I ZQUI E RD O”, ”*”
4640 PRINT USING “2(7X, 15A, 7X, A, 2X)”;” TRAMO CENTRAL”, “*”, “ APOYO DERECHO ”, ”*”
4650 P RI NT US I NG “48 A, 2( 32A) , /, 111A”; T I T 7$[ 1, 48 ] , TI T 7$[ 17, 47], T IT 7 $[ 17, 47] , L$
4660 FOR I =1 T O Nvt
4670 L=M b( I , 3)
4680 J=3*I - 2
4690 I F L >2. 0 T HE N J=5*I - 4
4700 P RI NT US I NG “#, 4D, 5D, 5D, X”;I , Mb( I , 1) , Mb( I , 2) , ” *”
4710 GOS UB S _30
4720 J=3*I - 1
4730 I F L >2. 0 T HE N J=Jn( I )
4740 GOS UB S _30
4750 J=3*I
4760 I F L >2. 0 T HE N J=5*I
4770 GOS UB S _30
4780 P RI NT
4790 I F I =Nvt T HE N P RI NT US I NG “111 A” ; L$
4800 NE XT I
4810 GOS UB Ch ek
4820 RE T URN
4830 S _30 : L i=Re st ( J, 1)
4840 M i=Re st ( J, 4)
4850 V i= Rest ( J, 5)
4860 T i=Re st ( J, 7)
4870 P RI NT US I NG “#, 3D. 2D, 3( 5D. 2D) , 2A”; L i, M i, V i, T i, ” *”
4880 RE T URN
4890 C h ek: P RI NT US I NG”/ , 30A, / , 35A”; ”S UM A DE CARG AS E N LAS JUNT AS ”, L $
4900 P RI NT US I NG “5 A, 3( 8X, 2A) , / , 35A”; ”Ju nt a” , ”P z”, ”M y”, ”M x”, L$
4910 FOR I =1 T O Nj
4920 P RI NT US I NG “5D, 3( 7D. 2D ) ”;I , Qz( I, 1) ,Qz( I , 2), Qz( I, 3)
4930 NE XT I
4940 P RI NT US I NG “35 A” ;L$
4950 RE T URN
4960 ! - - - -- - -- - -- - -- - -- - -- -- - -- - -- - -- - -- - -- - -- -- - --- - - -- - -- - -- -- - -- - -- - -- - -- - -- - - - -- - -- - -- - -- - -- - -- - - -
4970 Ej emp lo 6: DAT A E JE MP LO Nº 6 RETI CULO ( E sviado ) , 217. 37, 86. 948, 1
4980 DAT A 4, 4, 1, 2
4990 DATA 1. 2.993, 0, 2, 4.993, 2.993, 3, 0, 2.0, 4,2.0,4.993 !coordenadas de juntas
4995 DATA 1, 3, .7, 72.92, 228.87, 28.0 ! tipos de vigas, propiedades
5000 DAT A 1, 1, 2, 1, 2, 3, 4, 1, 3, 1, 3, 1, 4, 2, 4, 1 ! dat o s de vigas
5010 ! ! CAS O 1 DE CARG A
5020 DAT A 4, 1, 40. 0, 7. 491, - 2. 211, 2, 40. 0, - 0. 824, - 7. 767 ! car gas en ju nt as
5030 DAT A 3, 40. 0, 0. 824 , 7. 767 4, 40. 0, - 7. 491, 2. 211
5040 DAT A 0, 0 ! pp aut o mát ico , t r amo s co n p. unifo r me
5050 ! ¡ C AS O 2 DE CARG A
5060 DAT A 0, 0 ! ju nt as co n car ga, pp. aut o mát ico
5070 DAT A 2, 1, 2. 0, 2, 2. 0 ! Nº del t r amo , car ga. unifo r me ( p)
5080 ! ( Dat o s del Ejempl o7 e n la pág. 51 )
5200 Fi n : E ND

Fundaciones Elásticas 47
Ejemplo Nº 6. Retículo Esviado [3]

Se analiza el retículo ortogonal estudiado en la Ref [3], rotado un cierto ángulo respecto a
los ejes X-Y como se muestra en la figura anexa. El Caso(1) de cargas aplicadas en las juntas se
refiere a los momentos M1 = 6.00 , M2 = 5.00 tf.m previstos en el retículo original y señalados
con líneas punteadas en el modelo esviado, cuyas proyecciones según los ejes X-Y, de acuerdo
con la ecuación (2.15) , tienen los valores mostrados en el cuadro de cargas. El Caso (2) se refiere
a cargas uniformes en las vigas (1) y (2). Para comparación de resultados se declaran vigas con
ancho de base, inercia y constante de torsión iguales a las originales.

CASO 1 DE CARGA

Junta Pz My Mx
1 40.00 7.4908 - 2.2107
2 40.00 - 0.8237 - 7.7667
3 40.00 0.8237 7.7667
4 40.00 - 7.4908 2.2107

CASO 2 DE CARGA

Viga p(tf/m)
1 2.00
2 2.00

EJEMPLO Nº 6. RETICULO (Esviado) FUNDACION ELASTICA (Balasto = 1.0 kgf/cm3)

Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 t/cm2

Junta X Y Viga J1 J2 Tipo L Cx Cy

1 2.99 0.00 1 1 2 1 3.60 +0.55560 +0.83145
2 4.99 2.99 2 3 4 1 3.60 +0.55560 +0.83145
3 0.00 2.00 3 1 3 1 3.60 - 0.83145 +0.55560
4 2.00 4.99 4 2 4 1 3.60 - 0.83145 +0.55560

TIPOS DE VIGAS
V-Tipo Bo Ho Bf Hf * Base * Io (m4) Lambda Jp (m4) Lambda_t
1 1 - - - - * 0.70 * 7.29200E-03 3.085 2.28870E-02 26.386

48 Fundaciones Elásticas
EJEMPLO Nº 6. RETICULO (Esviado) FUNDACION ELASTICA (Balasto = 1.0 kgf/cm3)
Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 t/cm2

CASO 1 DE CARGA
CARGAS APLICADAS DEFORMACIONES PRESIONES
Junta Q MY MX Viga Pt DZ RY RX Junta Presión
1 40.00 7.49 - 2.21 1 0.00 16.76803 0.56069 1.75790 1 16.77
2 40.00 - 0.82 - 7.77 2 0.00 16.90552 1.87030 - 0.58306 2 16.91
3 40.00 0.82 7.77 3 0.00 16.90552 - 1.87030 0.58306 3 16.91
4 40.00 - 7.49 2.21 4 0.00 16.76803 - 0.56069 - 1.75790 4 16.77

Viga J1 J2 X Flecha G.Flexión M.Flexión Corte G.Torsión M.Torsión Presión

1 1 2 0.00 16.77 - 1.15009 3.999 - 22.177 1.44287 1.241 16.77
0.90 15.78 - 0.89913 - 11.304 - 11.937 1.38757 1.205 15.78
2.00 15.35 0.18960 - 17.855 - 0.018 1.32210 1.162 15.35
2.70 15.75 0.93829 - 15.230 7.565 1.28176 1.136 15.75
3.60 16.91 1.52392 - 3.862 17.823 1.23112 1.104 16.91
2 3 4 0.00 16.91 - 1.52392 - 3.862 - 17.823 - 1.23112 1.104 16.91
0.90 15.75 - 0.93829 - 15.230 - 7.565 - 1.28176 1.136 15.75
1.60 15.35 - 0.18960 - 17.855 0.018 - 1.32210 1.162 15.35
2.70 15.78 0.89913 - 11.304 11.937 - 1.38757 1.205 15.78
3.60 16.77 1.15009 3.999 22.177 - 1.44287 1.241 16.77
3 1 3 0.00 16.77 - 1.44287 - 3.759 - 17.823 - 1.15009 2.001 16.77
0.90 15.68 - 0.86243 - 15.160 - 7.629 - 1.24127 2.032 15.68
1.61 15.33 - 0.10733 - 17.838 0.014 - 1.31383 2.058 15.33
2.70 15.85 0.97498 - 11.374 11.873 - 1.42808 2.101 15.85
3.60 16.91 1.23112 3.897 22.177 - 1.52392 2.138 16.91
4 2 4 0.00 16.91 - 1.23112 3.897 - 22.177 1.52392 2.138 16.91
0.90 15.85 - 0.97498 - 11.374 - 11.873 1.42808 2.101 15.85
1.99 15.33 0.10733 - 17.838 - 0.014 1.31383 2.058 15.33
2.70 15.68 0.86243 - 15.160 7.629 1.24127 2.032 15.68
3.60 16.77 1.44287 - 3.759 17.823 1.15009 2.001 16.77

COMPROBACION:
SUMA DE CARGAS
Junta Pz My Mx
1 40.00 7.49 - 2.21
2 40.00 - 0.82 - 7.77
3 40.00 0.82 7.77
4 40.00 - 7.49 2.21

Fundaciones Elásticas 49
 EJEMPLO Nº 6. RETICULO (Esviado) FUNDACION ELASTICA (Balasto = 1.0 kgf/cm3)
Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 tf/cm2

CASO 2 DE CARGA
CARGAS APLICADAS DEFORMACIONES PRESIONES
Junta Q MY MX Viga Pt DZ RY RX Junta Presión
1 3.56 1.18 - 1.77 1 2.00 1.42857 0.15735 - 0.03129 1 1.43
2 3.56 - 1.18 1.77 2 2.00 1.42857 0.03129 0.15735 2 1.43
3 3.56 1.18 - 1.77 3 0.00 1.42857 - 0.03129 - 0.15735 3 1.43
4 3.56 - 1.18 1.77 4 0.00 1.42857 - 0.15735 0.03129 4 1.43

Viga J1 J2 X Flecha G.Flexión M.Flexión Corte G.Torsión M.Torsión Presión

1 1 2 0.00 1.43 0.11344 - 0.006 1.697 0.11344 0.006 1.43
0.90 1.52 0.07795 1.126 0.827 0.11325 0.003 1.52
1.80 1.56 - 0.00000 1.496 - 0.000 0.11318 0.000 1.56
2.70 1.52 - 0.07795 1.126 - 0.827 0.11325 - 0.003 1.52
3.60 1.43 - 0.11344 - 0.006 - 1.697 0.11344 - 0.006 1.43
2 3 4 0.00 1.43 0.11344 - 0.006 1.697 - 0.11344 - 0.006 1.43
0.90 1.52 0.07795 1.126 0.827 - 0.11325 - 0.003 1.52
1.80 1.56 - 0.00000 1.496 - 0.000 - 0.11318 0.000 1.56
2.70 1.52 - 0.07795 1.126 - 0.827 - 0.11325 0.003 1.52
3.60 1.43 - 0.11344 - 0.006 - 1.697 - 0.11344 0.006 1.43
3 1 3 0.00 1.43 - 0.11344 0.006 - 1.697 0.11344 0.006 1.43
0.90 1.34 - 0.07795 - 1.126 - 0.827 0.11325 0.003 1.34
1.80 1.30 0.00000 - 1.496 0.000 0.11318 - 0.000 1.30
2.70 1.34 0.07795 - 1.126 0.827 0.11325 - 0.003 1.34
3.60 1.43 0.11344 0.006 1.697 0.11344 - 0.006 1.43
4 2 4 0.00 1.43 - 0.11344 0.006 - 1.697 - 0.11344 - 0.006 1.43
0.90 1.34 - 0.07795 - 1.126 - 0.827 - 0.11325 - 0.003 1.34
1.80 1.30 0.00000 - 1.496 0.000 - 0.11318 0.000 1.30
2.70 1.34 0.07795 - 1.126 0.827 - 0.11325 0.003 1.34
3.60 1.43 0.11344 0.006 1.697 - 0.11344 0.006 1.43

COMPROBACION:
SUMA DE CARGAS
Junta Pz My Mx
1 0.00 0.00 0.00
2 0.00 0.00 0.00
3 0.00 0.00 0.00
4 0.00 0.00 0.00

50 Fundaciones Elásticas
 Ejemplo Nº 7.- Declaración de datos:

Retículo con vigas de sección rectangular y sección T invertida, con dos casos de carga:
(1) cargas aplicadas en las juntas, (2) cargas debidas al peso propio. La numeración de la DATA
muestra que los datos de este ejemplo pueden declararse a continuación de los del Ejemplo Nº 6,
incluídos en el programa. Para procesar el Ejemplo Nº 7 se debe eliminar el símbolo (!) en la
instrucción 185 y colocarlo en la instrucción 180.

5080 Ejemplo7: DATA EJEMPLO Nº 7, 238.75, 99.48, 1.50

5090 DATA 5, 5, 2, 2
5100 DATA 1, 10, 0, 2, -1.50, 3.0, 3, 0, 3.0, 4, 6.0, 3.0, 5, 10.0, 5.0 ¡!
5110 DATA 1, 1,0.30, 0.60, 2, 2, 0.30, 0.60, 0.80, 0.30
5120 DATA 1, 2, 3, 1, 2, 3, 4, 1, 3, 1, 5, 1, 4, 1, 4, 2, 5,4, 5, 2
5130 ¡ CASO 1 DE CARGA
5140 DATA 4, 1,20.0, 6.0, -12.0, 3, 18.0, 8.0, 0, 4, 50.0, 10.0, 20.0, 5,15.0, -4.0, -3.0
5150 DATA 0, 0
5160 ¡ CASO 2 DE CARGA
5170 DATA 0, 1, 0 ! N.Juntas cargadas , pp automático , N.Tramos cargados
5200 Fin: END

EJEMPLO Nº 7 FUNDACION ELASTICA (Balasto = 1.50 kgf/cm3)

Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 tf/cm2

Junta X Y Viga J1 J2 Tipo L Cx Cy

1 10.00 0.00 1 2 3 1 1.50 +1.00000 +0.00000
2 -1.50 3.00 2 3 4 1 6.00 +1.00000 +0.00000
3 0.00 3.00 3 1 5 1 5.00 +0.00000 +1.00000
4 6.00 3.00 4 4 1 2 5.00 +0.80000 - 0.60000
5 10.00 5.00 5 4 5 2 4.47 +0.89443 +0.44721

TIPOS DE VIGAS
V-Tipo Bo Ho Bf Hf * Base * Io (m4) Lambda Jp (m4) Lambda_t
1 1 0.30 0.60 - - * 0.30 * 5.40000E-03 3.272 3.70786E-03 33.059
11 2 0.30 0.60 0.80 0.30 * 0.80 * 8.36591E-03 2.856 6.49800E-03 10.050

Fundaciones Elásticas 51
 EJEMPLO Nº 7. FUNDACION ELASTICA (Balasto = 1.50 kgf/cm3)
Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 tf/cm2

CASO 1 DE CARGA
CARGAS APLICADAS DEFORMACIONES PRESIONES
Junta Q MY MX Viga Pt DZ RY RX Junta Presión
1 20.00 6.00 -12.00 1 0.00 6.16012 1.11212 - 0.10095 1 9.24
2 0.00 0.00 0.00 2 0.00 9.36315 0.16158 1.36702 2 14.04
3 18.00 8.00 0.00 3 0.00 9.53624 - 0.02340 1.36843 3 14.30
4 50.00 10.00 20.00 4 0.00 7.01287 - 0.87251 1.40235 4 10.52
5 15.00 - 4.00 - 3.00 5 0.00 4.54769 0.14790 - 0.20825 5 6.82

Viga J1 J2 X Flecha G.Flexión M.Flexión Corte G.Torsión M.Torsión Presión

1 2 3 0.00 9.36 0.16158 0.000 - 0.00 1.36702 - 0.000 14.05
0.75 9.48 0.13853 1.190 3.180 1.36737 - 0.003 14.22
1.50 9.54 - 0.02340 4.779 6.393 1.36843 -0.007 14.30
2 3 4 0.00 9.54 - 0.02340 12.779 - 11.607 1.36843 - 0.007 14.30
1.50 8.82 - 0.68262 0.119 - 5.357 1.37265 - 0.014 13.23
2.93 7.99 - 0.41832 - 3.626 0.013 1.37929 - 0.020 11.98
4.50 7.61 - 0.15520 0.760 5.513 1.38960 - 0.028 11.41
6.00 7.01 - 0.87251 12.822 10.507 1.40235 - 0.035 10.52
3 1 5 0.00 6.16 0.10095 12.102 - 9.228 1.11212 0.718 9.24
1.25 5.76 - 0.58283 2.712 - 5.835 0.86955 0.714 8.65
2.50 5.00 - 0.55547 - 2.646 - 2.810 0.62821 0.711 7.50
3.83 4.51 - 0.15699 - 4.477 0.001 0.37301 0.708 6.76
5.00 4.55 0.20825 - 3.088 2.371 0.14790 0.707 6.82
4 4 1 0.00 7.01 0.14340 28.087 - 25.222 1.64539 3.228 10.52
1.25 6.46 - 0.76407 3.054 - 14.976 1.03227 3.121 9.69
3.42 5.38 0.04709 - 12.622 0.013 0.00182 3.050 8.07
3.75 5.43 0.25541 - 12.260 2.170 - 0.15524 3.051 8.15
5.00 6.16 0.82912 - 4.287 10.772 - 0.74804 3.088 9.24
5 4 5 0.00 7.01 - 1.40754 1.811 - 14.271 0.86410 1.498 10.52
1.12 5.52 - 1.15501 - 9.274 - 5.894 0.60976 1.446 8.28
2.10 4.64 - 0.60484 - 12.062 0.061 0.39204 1.414 6.97
3.35 4.33 0.06002 - 7.789 6.694 0.12068 1.394 6.50
4.47 4.55 0.22542 2.984 12.629 - 0.12012 1.394 6.82

COMPROBACION:
SUMA DE CARGAS
Junta Pz My Mx
1 20.00 6.00 - 12.00
2 0.00 0.00 - 0.00
3 18.00 8.00 0.00
4 50.00 10.00 20.00
5 15.00 - 4.00 - 3.00

52 Fundaciones Elásticas
EJEMPLO Nº 7. FUNDACION ELASTICA (Balasto = 1.50 kgf/cm3)
Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 tf/cm2

CASO 2 DE CARGA
CARGAS APLICADAS DEFORMACIONES PRESIONES
Junta Q MY MX Viga Pt DZ RY RX Junta Presión
1 3.05 - 1.29 - 1.87 1 0.45 0.77656 0.03950 - 0.01837 1 1.16
2 0.34 0.08 0.00 2 0.45 1.06278 - 0.03225 - 0.00220 2 1.59
3 1.61 1.17 0.00 3 0.45 1.01402 - 0.03324 - 0.00220 3 1.52
4 5.02 1.22 0.38 4 0.82 0.76536 - 0.02763 - 0.00225 4 1.15
5 2.88 - 1.18 1.49 5 0.82 0.77441 0.03267 0.02037 5 1.16

Viga J1 J2 X Flecha G.Flexión M.Flexión Corte G.Torsión M.Torsión Presión

1 2 3 0.00 1.06 - 0.03225 0.000 0.000 - 0.00220 - 0.000 1.59
0.75 1.04 - 0.03238 0.007 0.017 - 0.00220 0.000 1.56
1.50 1.01 - 0.03324 0.024 0.026 - 0.00220 0.000 1.52
2 3 4 0.00 1.01 - 0.03324 0.024 0.026 - 0.00220 0.000 1.52
1.50 0.96 - 0.03828 0.061 0.018 - 0.00221 0.000 1.44
2.07 0.94 - 0.04118 0.068 0.005 - 0.00221 0.000 1.41
4.50 0.83 - 0.04776 - 0.049 - 0.122 - 0.00223 0.000 1.24
6.00 0.77 - 0.02763 - 0.332 - 0.262 - 0.00225 0.000 1.15
3 1 5 0.00 0.78 0.01837 - 0.094 0.231 0.03950 0.005 1.17
1.25 0.80 0.01596 0.119 0.112 0.03771 0.005 1.20
2.54 0.81 - 0.00075 0.191 0.000 0.03592 0.005 1.22
3.75 0.80 - 0.01668 0.128 - 0.105 0.03430 0.005 1.20
5.00 0.77 - 0.02037 - 0.077 - 0.225 0.03267 0.005 1.16
4 4 1 0.00 0.77 - 0.02346 - 0.124 - 0.142 0.01478 0.068 1.15
1.25 0.74 - 0.01153 - 0.237 - 0.044 0.00165 0.068 1.11
1.97 0.74 - 0.00260 - 0.252 0.001 - 0.00589 0.068 1.11
3.75 0.75 0.01696 - 0.153 0.117 - 0.02473 0.069 1.13
5.00 0.78 0.02058 0.061 0.231 - 0.03839 0.072 1.17
5 4 5 0.00 0.77 - 0.02371 - 0.180 - 0.119 - 0.01437 - 0.068 1.15
1.12 0.75 - 0.01088 - 0.261 - 0.030 - 0.00267 - 0.067 1.12
1.57 0.74 - 0.00485 - 0.268 0.001 0.00206 - 0.067 1.11
3.35 0.75 0.01590 - 0.161 0.123 0.02079 - 0.069 1.13
4.47 0.77 0.02011 0.030 0.225 0.03283 - 0.071 1.16

COMPROBACION:
SUMA DE CARGAS
Junta Pz My Mx
1 - 0.00 0.00 0.00
2 - 0.00 - 0.00 0.00
3 0.00 - 0.00 - 0.00
4 - 0.00 0.00 0.00
5 0.00 0.00 - 0.00

Fundaciones Elásticas 53
2.9.3. Salida Resumida

Las máximas deformaciones y presiones del suelo en vigas apoyadas en suelo elástico se
producen en las juntas donde se aplican las cargas y constan en el cuadro de salida de Cargas,
Deformaciones y Presiones como se muestra en los ejemplos Nº 6 y Nº 7. Estos valores sirven
para comprobar que los asentamientos y presiones cumplen con las recomendaciones pautadas en
el estudio de suelos correspondiente y que, en este aspecto, la estructura propuesta es válida.

Tomando en cuenta que los puntos críticos son los extremos y el de momento negativo
máximo absoluto en el tramo, se procede a comprobar si la sección de la viga en geometría,
calidad del concreto y cantidad de refuerzo son, en este otro aspecto, aceptables.

En el programa FUNDARETIC se incluye una salida resumida que el lector puede usarla
a discreción para aceptar o modificar la estructura. Imprime en una sola línea los momentos y
cortes de las tres secciones descritas, como se muestra en el Ejemplo Nº 8. Además, se prevé que
el usuario pueda imprimir todos los valores para los cuartos de la luz.

Ejemplo Nº 8.

Retículo del Ejemplo Nº 7 con el Caso Nº 1 de carga.

EJEMPLO Nº 7. (Impresión Resumida) FUNDACION ELASTICA (Balasto = 1.50 kgf/cm3)

Unidades: m , tf , desplazamientos = mm , rotaciones = rad * 1000 , E = 217.37 , G = 86.95 t/cm2

CASO 1 DE CARGA
CARGAS APLICADAS DEFORMACIONES PRESIONES
Junta Q MY MX Viga Pt DZ RY RX Junta Presión
1 20.00 6.00 -12.00 1 0.00 6.16012 1.11212 - 0.10095 1 9.24
2 0.00 0.00 0.00 2 0.00 9.36315 0.16158 1.36702 2 14.04
3 18.00 8.00 0.00 3 0.00 9.53624 - 0.02340 1.36843 3 14.30
4 50.00 10.00 20.00 4 0.00 7.01287 - 0.87251 1.40235 4 10.52
5 15.00 - 4.00 - 3.00 5 0.00 4.54769 0.14790 - 0.20825 5 6.82

APOYO IZQUIERDO TRAMO CENTRAL APOYO DERECHO

Viga J1 J2 X Mf V Mt X Mf V Mt X Mf V T
1 2 3 0.00 0.00 - 0.00 - 0.00 0.75 1.19 3.18 - 0.00 1.50 4.78 6.39 - 0.01
2 3 4 0.00 12.78 -11.61 - 0.01 2.93 - 3.63 0.01 - 0.02 6.00 12.82 10.51 - 0.03
3 1 5 0.00 12.10 - 9.23 0.72 3.83 - 4.48 0.00 0.71 5.00 - 3.09 2.37 0.71
4 4 1 0.00 28.09 - 25.22 3.23 3.42 - 12.62 - 0.01 3.05 5.00 - 4.29 10.77 3.09
5 4 5 0.00 1.81 -14.27 1.50 2.10 - 12.06 0.06 1.41 4.47 2.98 12.63 1.39

54 Fundaciones Elásticas
3. ANEXO A
PROPIEDADES DEL CONCRETO [6]

f ´c = resistencia especificada del concreto (kgf/cm2)

Ec = módulo de elasticidad longitudinal Ec = 15.1 f¨c (tf/cm2)
Ec
Gc = módulo de elasticidad transversal Gc = (tf/cm2)
2 1+ 
γ = módulo de Poisson = 0.20

PROPIEDADES DEL MODULO DE BALASTO (suelos arenosos) [5]

β30 = módulo de balasto de ensayo de suelos con placa de 30x30 cm , en (kgf/cm3)

β = módulo de balasto de análisis (kgf/cm3)
2
 0.30 + B 
Fundación cuadrada de lados (B x B) : β = β30   , B = ancho en metros
 2B 
L (0.5  m)
Fundación rectangular ( B x L) : m = , β = β30  0.667β30
B 1.5m

En cualquier caso: β (kgf/cm3) = 1.22 qadm , qadm = esfuerzo admisible (kgf/cm2)

con FS =3 , asentamiento 2.54 cm.

PROPIEDADES DE LAS SECCIONES [10 , 6]

Sección Rectangular: [10]

A=b.h ; I = b . h3 / 12

b 1 b 
4
1
J = b3 h βt ; Si: (b / h) ≤ 1. : β t   0.21 1    
3 h  12  h  
 
L
h=3  3β / Ec  L/u 
4
Elegido u=  1.0

Sección T invertida: [6]
1 2
c = bf - bo ; A = bo . ho + c . hf ; m= (h o .bo  h f2 .c)
2
2
1 3 3 m
I= (h o .bo  h f .c) -
3 A
hi = ho – hf : Si: (bo / hi) ≤ 1.
3

J o = b o .h i 1- 0.63 b o / h i 
3

Si: (bf / hf ) ≤ 1. J f = b f .h f 1- 0.63 b f / h f 
3
Jo  Jf
J=
1
 
para la constante de Torsión ( J ) si (b* / h * )  1.0 , alternar las dimensiones: (h * / b* )

Fundaciones Elásticas 55
4. ANEXO B

CONVENCION DE SIGNOS DE LA ESTRUCTURA y DE LA FUNDACION

Salvo estudios de inter-acción Suelo-Estructura los análisis de la estructura del edificio y

de la respectiva fundación se hacen independientemente uno del otro; en cuyo caso es importante
confirmar la compatibilidad de signos entre fundación y estructura.

Los tipos de análisis comúnmente utilizados en el cálculo de la estructura, son:

(B.1) Estructura Espacial cuyos momentos se representan con vectores aplicados al eje
al cual pertenecen.

(B.2) Estructura con Pórticos Planos cuyos momentos están en el plano al que
pertenece el pórtico y cuya representación vectorial está en el plano normal.

B.1 Estructura Espacial: B.2 Estructura con Pórticos Planos: (Px , Py)

qz qz

mx . Px my
Py
my mx

Convención Positiva Convención Positiva

En todo caso, los resultados del análisis de la estructura, acciones o reacciones, deben
ajustarse a la nomenclatura y convención de signos de las acciones en la fundación, según lo
indicado en la Sección 2.9.1 – pág. 36 – del texto. Las cargas incluyen momentos de flexión y
carga axial en planos verticales normales al de la fundación. Se excluyen torsiones y cargas
axiales en el plano de la fundación.

CASOS DE COMBINACION DE SOLICITACIONES

(FONDONORMA 1753-2006) (COVENIN 1756-2001)

NOMENCLATURA PARA FUNDACIONES

CP = Carga Permanente 1) 1.4 CP
CV = Carga Variable 2) 1.2 CP + 1.6 CV
SH = Carga por sismo horizontal 3) 1.1 CP + 1.0 CV ± (SH + SV)
SV = Carga por sismo vertical 4) 0.9 CP ± (SH + SV)
SV = CP(0.2 αφβ Ao)
(SH) y (SV) con igual signo αφβ Ao = parámetros sísmicos

. .

56 Fundaciones Elásticas
5. ANEXO C

OBSERVACIONES AL ESTUDIO DE SUELOS

Con relativa frecuencia, los estudios de suelos recomiendan el uso de fundaciones rígidas,
fijando el módulo de balasto de referencia, varias especificaciones constructivas y escasa
información sobre los requisitos que debe cumplir la fundación para ser considerada rígida. Con
este Anexo se intenta despejar las dudas, comprender mejor el texto y estimular el uso de los
programas.
En la ec(1.5) se define como longitud elástica: λ= 4 4EI / β.r en metros, a diferencia de
la expresión inversa de Hetenyi que resulta en metros-1 .

En el art. 1.2 se define como viga de longitud corta a aquella cuyo valor de la relación
(L/ ) es:  /4 según Hetenyi;  1.50 según Segovia;  1.75 según W. C. Teng, en cuyo caso
la fundación se asimila a una pieza rígida, simplificando el cálculo estático en casos de
fundaciones de un tramo o de retículos con limitaciones de geometría y cargas.

Para estimar la altura de una fundación rectangular, se elige un valor de u = L/, y se

aplica la fórmula del Anexo A: h= 3  3β/E  (L/u)4 que es independiente del ancho de la
fundación. Ejemplo: determinar la altura de una viga de luz L = 9.60 m, de concreto con una
resistencia f ´c = 250 kgf/cm2, un módulo de balasto de diseño  = 1500 tf/m3 y una longitud
elástica  = L/1.40; Entonces E=15.1 250 *104  2.39*106 tf/m2 ; u =(L/) = 1.40 y
h= 3 (3*1500/2.39*106 )(9.60/1.40) 4 =1.61 m.; el ancho de la fundación se determina con una
presión promedio de la carga total / área de apoyo. Estas dimensiones no garantizan una pieza
rígida pero sirven para hacer un análisis inicial como viga de longitud intermedia, según el
Art.1.3, para obtener resultados de asentamientos y presiones que ayuden a tomar decisiones que
mejoren su comportamiento.

De lo expuesto, se concluye que la clasificación de vigas según el Art. 1.2 establece la

diferencia entre el cálculo aproximado de vigas de longitud corta o larga, y el procedimiento más
exacto del Art. 1.3 para fundaciones elásticas de longitud corta, intermedia o larga y, por lo
tanto, aplicable al uso de los programas.

Respecto al módulo de balasto de referencia, los estudios de suelos recurren a tablas de

valores, más o menos estandarizados según el tipo de suelo, que reemplazan a las pruebas de
carga con placas de 30 x 30 cm. El balasto de análisis para el cálculo de la estructura real se
obtiene por extrapolación del balasto de referencia mediante las fórmulas de Therzaghi, o los
criterios de J. Bowles , brevemente resumidos en el Anexo A.

Además, el estudio de suelos debe incluir el cálculo de la presión admisible (qa) con el
factor de seguridad (FS) utilizado y el desplazamiento diferencial (dw) previsto, que permitan
comprobar la validez de las dimensiones de la fundación, utilizadas en el análisis definitivo.

Fundaciones Elásticas 57
6. BIBLIOGRAFIA
[1] Segovia G. A., 1970, Correspondencia Personal

[2] Segovia G. A., 1974, El Método de Rigideces Sucesivas, I Congreso Venezolano de

Sismología e Ingeniería Sísmica, AVIE – FUNVISIS.

[3] Aguiar F. R., 2005, Análisis Estático de Estructuras, ESPE, 1ª edición, 154 págs.

[4] Hetenyi M., 1979, Beams on Elastic Foundation, The University of Michigan Press,
Ann Arbor, eleventh printing, 255 págs.

[5] Bowles J., 1977, Foundation Analysis and Design , Mc Graw Hill, N.Y,
Second edition, 492 págs.

[6] A.C.I., 2005, Requisitos del Reglamento para Concreto Estructural, (ACI 318SR-05),
American Concrete Institute, Farmington Hills, MI 48333, 490 págs.

[7] Weaver W. and Gere J., 1980, Matrix Analysis of Framed Structures, Van Nostrand
Reinhold Company, second edition, 492 págs.

[8] TRANSERA www.htbasic.com

[9] Arias Alban G., 1984, Análisis Sísmico Estático de Edificios, con Aplicación a
Microcomputadores Ediciones Vega, Caracas, 275 págs.

[10] Timoshenko S. and Gere J. , 1961, Theory of Elastic Stability, Mc Graw Hil, N.Y,
Second edition, 541 págs.

[11] Reti Andrés, 1951, Las Fundaciones Continuas y Cruzadas de Concreto Armado
Ministerio de Obras Públicas de Venezuela
Dirección de Edificaciones Médico Asistenciales.
Caracas.

[12] Fortoul P. Celso, 1973, Análisis de Estructuras Planas Solicitadas por Cargas
Normales a su Plano. Notas del Curso: Uso del Computador en
Estática y Dinámica en la Maestría de la USB 1973-1974.
Sartenejas.

[13] Teng Wayne C.., 1975, Mat Foundations, Foundation Engineering Handbook
H.F. Winterkorn, F. Hasai-Yang , Editores, Van Nostrand
Reinhold Company, 751 págs.

58 Fundaciones Elásticas