TRABAJO FINAL

INGENIERIA GEOTECNICA

JONATHAN QUISPE, MARIA PAULA MENDOZA Y WILIAMS

CABALLEO
CICLO 2018-02
Profesor: Odilio
Pregunta N° 1
Utilizando el criterio de Meyerhof se le solicita hallar las dimensiones óptimas y de
diseño por resistencia (considerando incrementos de 0.1m) para una zapata rígida
rectangular con relación L/B=1.3 y con cargas Qv, Qh y M (para la carga vertical final
deberá considerar un edificio de 6 pisos tomando en cuenta el plano de cimentación del
Anexo 1). Las dimensiones del edificio son similares, sin embargo, el proyecto se ha
planteado en 3 lugares distintos. La profundidad de cimentación es Df y el nivel freático
se encuentra a NF m de profundidad, el suelo presenta una densidad relativa de denso a
muy denso.

El suelo por debajo del nivel de cimentación tiene un peso específico saturado de 𝛾2𝑠𝑎𝑡 , un
ángulo de fricción igual a 𝜙 y un valor de cohesión de “c”, el suelo que se encuentra por
encima del nivel de cimentación tiene un peso específico natural de 𝛾1𝑛𝑎𝑡 y un peso
específico saturado de 𝛾1𝑠𝑎𝑡 . Se conoce además que el suelo de apoyo tiene un
coeficiente de Poisson µ y a partir de ensayos MASW se ha determinado un valor
promedio de velocidad de onda de corte de Vs.

Por otro lado se conoce que la zapata Z-? que se ubica a una distancia L de la zapata en
análisis presenta un asentamiento individual de 𝛿2. El diseño contempla un edificio
donde no se permite la presencia de grietas. (1) Resultados de ensayos químicos en el
terreno de apoyo.

Grupo A Grupo B Grupo C

Zapata Ejes 3 y F Ejes 2 y F Ejes 2 y
G
Qv Ver Ver Ver
Anexo 1 Anexo 1 Anexo 1
Qh 90 ton 80 ton 92.5 ton
M 48 ton-m 45.5 ton- 58 ton-m
m
Df 1.2 m 1.3 m 1.2 m
NF 0.7 m 0.6 m 0.5 m
𝛾1𝑛𝑎𝑡 16.5 16.5 16.6
kN/m3 kN/m3 kN/m3
𝛾1𝑠𝑎𝑡 18.5 18.3 18.6
kN/m3 kN/m3 kN/m3
𝛾2𝑠𝑎𝑡 19.5 19.2 19.4
kN/m3 kN/m3 kN/m3
𝜙 32° 31° 32°
c 30 kPa 31 kPa 29 kPa
Vs 313 m/s 323 m/s 319 m/s
µ 0.28 0.27 0.28
𝛿2 2.0 cm 1.95 cm 2.03 cm
Cloruro 4 7930 887 6 642
s (1) p.p.m. p.p.m. p.p.m.
Sulfato 8 142 453 1 440
s (1) p.p.m. p.p.m. p.p.m.
SST (1) 6 9150 15 100 17100
p.p.m. p.p.m. p.p.m.

NOTA: 𝛿2 corresponde al asentamiento producido en la zapata adyacente con la que se

evalúa la distorsión angular. El peralte de la zapata para todos los grupos es de 0.70m.

Se pide:
a) Las dimensiones óptimas por resistencia de la zapata. Se deberán colocar por lo
menos 3 procesos de iteración, todos los cálculos deberán ser presentados de
manera clara y ordenada.
 b) Las dimensiones de diseño por resistencia de la zapata, considerar incrementos de
0.1m.
c) Verificar si la zapata cumple con los criterios de asentamiento individual y distorsión
angular. Comentar.
d) Verificación por volteo y deslizamiento. Comentar.
e) Dibujar acotado y a escala un esquema de la zapata diseñada en elevación y en
planta. Considerar una zapata centrada con una columna de área de 0.4 x 0.4 m
f) Evaluar los resultados de los ensayos químicos de acuerdo al ACI-318 y E.060 e
indicar, de ser el caso, los niveles de daño a los que se encontrarían expuestas
las cimentaciones, asimismo indicar y sustentar las medidas de protección a
implementar

SOLUCION:

M =45.5ton . m

Qv=175.74 ton

Qh=80 ton

SUELO 1

γ 1 nat =1.65 g/c m3

0.6 m

N .F
Df =1.3 m

γ 1 sat =1.83 g/c m3

0.7 m
 SUELO 2 ∅=31 °
γ 2 sat =1.92 g/c m3
c=0.31kg /c m2
μ=0.27
Vs=323 m/s

Datos:
Zapata Eje 2 y F
Qv=¿
Qh=80 ton
M =45.5ton . m
Df =1.3 m
NF =0.6 m
γ 1 nat =16.5 kN /m3 =1.65 g /c m 3
γ 1 sat =18.3 kN /m3 =1.83 g /c m3
γ 2 sat =19.2 kN /m3 =1.92 g/c m3
∅=31 °
c=31 kPa=0.31 kg/c m2
Vs=323 m/s
δ 2=1.95 cm
μ=0.27
g=10 m/ s2

Para calcular el Qv de la zapata (Eje 2 y F) se considera una edificación de 6 pisos con una
carga de 1 ton/m2 y la suma de cada cuarta parte de área.
 Calculo de Área Total ( AT ) .

A 1 A2 A3 A 4
AT = + + +
4 4 4 4

A1=BxH =4.60 m x 5.20 m=23.92m 2

A2=BxH =5.50 m x 5.20 m=28.60 m2

A3 =BxH =4.60 m x 6.40 m=29.44 m2
A 4=BxH =5.50 m x 6.40 m=35.20 m2

23.92 28.60 29.44 35.20

AT = + + + =29.29 m2
4 4 4 4

Caculo de Carga Vertical (Qv).

1 ton
Qv= x Numero de pisos x AT
m2
1 ton
Qv= 2
x 6 x 29.29m 2=175.74 ton
m
 Qv=175.74 ton

Se considera falla general ya que el suelo tiene una densidad relativa de denso a muy denso.
∅=31 °
c=0.31kg /c m2

Para el cálculo de β

175.74
tan α=
80
Qv= 175.74 ton
β α =65.82 °
α α + β = 90°

β = 24.18°
Qh= 80 ton

Ahora, hallamos la excentricidad ( e x ) .

M
ex=
QV
45.5ton .m
ex=
175.74 ton
e x =0.259 m

Ahora, calculamos los factores de carga ( N q , N c , N γ ¿ en función de ɸ

( ∅2 +45 ° )
N q =e π . tan ∅ . tan 2
31
N q =e π . tan 31 . tan ( +45 ° )
2
2
N q =20..631

N c =cot ∅ .( N q −1)

N c =cot 31.(20.631−1)
N c =32.671

N γ =( N q−1 ) . tan ⁡( 1.4 ∅)

N γ =( 20.631−1 ) . tan ⁡(1.4 x 31)

 N γ =18.564

Como el Nivel Freatico(N.F) está entre 0 y la profundidad de cimentación(D f) hay una influencia
del N.F en q.

q=γ 1 nat∗D1 + γ 1 ∑ ¿∗D ¿

2

q=1.65 g/cm3∗0.6 m∗0.1+(1.83−1.00) g/cm3∗0.7 m∗0.1

q=0.157 kg /cm 2

Para hallar las dimensiones optimas y de diseño se usará el criterio de Meyerhof.

Como no tenemos las dimensiones de la cimentación se realizará iteraciones.

PRIMERA ITERACION:

B=2.00 m B´ =B−2e x B´ =2.00−2(0.259) B´ =1.482 m

L=1.3 B L=1.3( 2) L=2.6 m L ´=2.600 m
e x =0.259 m
QV 175.74 ton 2
q trans= ´ ´
= =4.561 kg /cm
B L (1.482m)(2.600 m)

Factores por forma:

B 1.482
∗N q ∗20.631
L 2.600
f cs=1+ =1+ =1.360
Nc 32.671
B 1.482
f qs =1+ ∗tan ∅=1+ ∗tan31=¿ 1.342¿
L 2.600

0.4∗B 0.4∗1.482
f ɣs =1− =1− =0.772
L 2.600

Factores por profundidad:

1−f qd 1−1.184
f cd =f qd− =1.184− =1.193
N c∗tan ∅ 32.671∗tan 31

f qd=1+2 tan ∅∗¿ ¿=1.184

f ɣd =1
 Factores por inclinación de carga:

β 2 24.18 2
f ci=f qi =(1− ) =(1− ) =0.535
90 90
β 2 24.18 2
f ɣi=(1− ) =(1− ) =0.048
∅ 31

Después, de calcular los factores ahora hallaremos la capacidad de carga ultima(q ¿¿ u) ¿.

1
q u=c . N c . f cs . f cd . f ci + q . N q . f qs . f qd . f qi + . B . γ 2 . N γ . f γ . f γ d . f γ i
2

1
q u=( 0.31 ) . ( 32.671 ) . (1.360 ) . ( 1.193 ) . ( 0.535 )+ ( 0.157 ) . ( 20.631 ) . ( 1.342 ) . (1.184 ) . ( 0.535 ) + . ( 1.482 ) . ( 0.92 ) .(0.1).(18.564)
2

q u=11.589 kg /cm2

Para el cálculo de carga admisible (q adm ¿ se reduce la capacidad de carga ultima(q ¿¿ u) ¿ con un
factor de seguridad (F.S) que en condiciones estáticas es tres.

qu
q adm=
F .S
qu
q adm =
3

11.589 kg/cm2
q adm=
3

q u=3.863 kg /cm2
Finalmente, la carga transmitida(q trans) debe ser menor que la carga admisible ( q adm ) .

q trans ≤ q adm

4.561 kg /cm2 ≤ 3.863 kg/cm2

Como no se cumple la condición se tendrá que cambiar la geometría de la zapata.

SEGUNDA ITERACION:

B=2.30 m B´ =B−2e x B´ =2.30−2(0.259) B´ =1.782 m

L ´=2.990 m
 L=1.3 B L=1.3(2.3) L=2.990 m

e x =0.259 m
QV 175.74 ton 2
q trans= ´ ´
= =3.298 kg /cm
B L (1.782m)(2.990 m)

Factores por forma:

B 1.782
∗N q ∗20.631
L 2.990 =1.376
f cs=1+ =1+
Nc 32.671 B 1.782
f qs=1+ ∗tan ∅=1+ ∗tan 31=1.358
L 2.990

0.4∗B 0.4∗1.782
f ɣs =1− =1− = 0.762
L 2.990

Factores por profundidad:

1−f qd 1−1.160
f cd =f qd− =1.160− =1.168
N c∗tan ∅ 32.671∗tan 31
f qd=1+2 tan ∅∗¿ ¿

f ɣd =1

Factores por inclinación de carga

β 2 24.18 2
f ci=f qi =(1− ) (1−
= ) =0.535
90 90
β 2 24.18 2
f ɣi=(1− ) =(1− ) =0.048
∅ 31

Después, de calcular los factores ahora hallaremos la capacidad de carga ultima(q ¿¿ u) ¿.

1
q u=c . N c . f cs . f cd . f ci + q . N q . f qs . f qd . f qi + . B . γ 2 . N γ . f γ . f γ d . f γ i
2

1
q u=( 0.31 ) . ( 32.671 ) . (1.376 ) . ( 1.168 ) . ( 0.535 )+ ( 0.157 ) . ( 20.631 ) . (1.358 ) . ( 1.160 ) . ( 0.535 )+ . ( 1.782 ) . ( 0.92 ) .(0.1).(18.564)
2

q u=11.492 kg /cm2
 Para el cálculo de carga admisible (q adm ¿ se reduce la capacidad de carga ultima(q ¿¿ u) ¿ con un
factor de seguridad (F.S) que en condiciones estáticas es tres.

qu
q adm=
F .S
qu
q adm =
3

11.492 kg/cm 2
q adm =
3

q adm=3.831 kg/cm 2

Finalmente, la carga transmitida(q trans) debe ser menor que la carga admisible ( q adm ) .

q trans ≤ q adm

3.298 kg /cm2 ≤ 3.831kg /cm2

Se cumple la condición, pero se desea que la condición se ha precisa.

TERCERA CONDICION:

B=2.160 m B´ =B−2e x B´ =2.160−2(0.259) B´ =1.642 m

L=1.3 B L=1.3( 2.160) L=2.808 m L ´=2.808 m

e x =0.259 m
QV 175.74 ton 2
q trans= ´ ´ = =3.812 kg/cm
B L (1.642m)(2.808 m)

Factores por forma

B 1.642
∗N q ∗20.631
L 2.808 =1.369
f cs=1+ =1+
Nc 32.671 B 1.642
f qs=1+ ∗tan ∅=1+ ∗tan31=1.351
L 2.808
0.4∗B 0.4∗1.642
f ɣs =1− =1− =0.766
L 2.808

Factores por profundidad

 1−f qd 1−1.170
f cd =f qd− =1.170− =1.179
N c∗tan ∅ 32.671∗tan31
f qd=1+2 tan ∅∗¿ ¿

f ɣd =1

Factores por inclinación de carga

β 2 24.18 2
f ci=f qi =(1− ) =(1− ) =0.535
90 90
β 2 24.18 2
f ɣi=(1− ) =(1− ) =0.048
∅ 31

Después, de calcular los factores ahora hallaremos la capacidad de carga ultima(q ¿¿ u) ¿.

1
q u=c . N c . f cs . f cd . f ci + q . N q . f qs . f qd . f qi+ . B . γ 2 . N γ . f γ . f γ d . f γ i
2
1
q u=( 0.31 ) . ( 32.671 ) . (1.369 ) . ( 1.179 ) . ( 0.535 )+ ( 0.157 ) . ( 20.631 ) . ( 1.351 ) . ( 1.170 ) . ( 0.535 ) + . (1.642 ) . ( 0.92 ) .(0.1) .(18.564). ( 0.7
2

q u=11.534 kg/cm 2

Para el cálculo de carga admisible (q adm ¿ se reduce la capacidad de carga ultima(q ¿¿ u) ¿ con un
factor de seguridad (F.S) que en condiciones estáticas es tres.

qu
q adm=
F .S
qu
q adm =
3
11.589 kg/cm2
q adm=
3

q adm =3.845 kg /cm 2

Finalmente, la carga transmitida(q trans) debe ser menor que la carga admisible ( q adm ) .

q trans ≤ q adm

3.812 kg /cm2 ≤3.845 kg /cm2

 Se cumple la condición, por lo tanto, la geometría de la zapata está bien.

a)
Dimensiones optimas por resistencia

B=2.160 m

L=2.808 m

b)
Dimensiones de diseño considerando incrementos de 0.1m.

B=2.200 m

L=2.900 m

c)
Verificación por servicio: Asentamiento individual y distorsión angular
Asentamiento individual:
Se considera las dimensiones de diseño
B=2.200 m
B´ =B−2e x
B´ =2.200−2(0.259)
B´ =1.682 m L ´=2.900 m
e x =0.259 m

μ=0.27

Para calcular la carga transmitida ¿ ¿) se considera las dimensiones efectivas.

175.74 ton
q trans= q trans=3.603 kg /cm2
(1.682 m)(2.900 m)

El módulo de elasticidad del suelo se obtiene del gráfico de módulo de corte y deformación.

γ 1.92
( )
Gmax = 2 sat V S2 =
g 10
x (323)2=20031168.00 kg /m2
Gs=10 % Gmax=10 % ¿) = 2003116.80 kg /m 2
Es=2 ( 1+ μ ) Gs=2 ( 1+0.27 ) (2003116.80 kg/m2 )=5087916.67 kg /m 2
El factor de influencia (IF) se obtiene de la interpolación de la tabla 3.

Tabla 3 de factor de influencia

Según, los datos del problema la zapata es rígida.

RECTANGULAR(L/B) RIGIDA
1.500 1.060
1.318 X
1.000 0.820

1.5−1.318 1.06− X
=
1.5−1.0 1.06−0.82

X =0.973

El factor de influencia (IF) es 0.973.

Finalmente, reemplazamos y resolvemos la ecuación del asentamiento instantáneo (∆ H ¿.

qtrans∗B´ ( 1−μ2 )
∆ H= ∗IF
Es
3.603∗1.682 ( 1−0.272 )
∆ H= ∗0.973
5087916.67
∆ H =1.074 cm

-El asentamiento individual (∆ H ¿ de una zapata debe ser como máximo a una pulgada.
Aplicando la fórmula delasentamiento individual ( ∆ H ) , nos resultó el ( ∆ H ) menor a una pulgada.
Entonces, nuestra zapata cumple con el criterio de asentamiento individual.
∆ H ≤ 2.54 cm
 1.074 cm ≤2.54 cm

Distorsión angular:
Para este cálculo se considera el asentamiento instanteo de otra zapata ( δ 2).

δ 1=1.074 cm
δ 2=1.95 cm
L=4.60 m

Ahora, reemplazamos los datos en esta ecuación.

δ 2−δ 1
α= ≤1 /500
L
1.950−1.074
α=
4.6∗100

α =0.00190345

- Según la Norma Técnica De Edificación E.050 Suelos y Cimentaciones, la distorsión angular (

α) debe ser menor igual a 1/500 que describe el limite seguro para edificios en los que no se
permiten grietas.

Aplicando la fórmula de la distorsión angular, nos resultó (α) menor a 1/500. Entonces, nuestra
zapata cumple con el criterio de distorsión angular.
0.00190 ≤ 0.00200
En conclusión, nuestra zapata cumple con los criterios de asentamiento individual y distorsión
angular.

d)

-Verificación por volteo:

Datos:

B 2.200 m
L 2.900 m
Qv 175.74 ton
M 45.5 ton.m
QH 80 ton
Df 1.3 m
WZAP 19.9056 ton
ex 0.259 m
d1 0.841 m

Usaremos la fórmula:

W ZAP∗B
QV ∗d 1 +
2
FS voleto= ≥ 1.5
M +Q H∗D f

19.9056∗2.2
175.74∗0.841+
2
FS voleto= ≥ 1.5
45.5+80∗1.3

FS voleto=1.135< 1.5

No cumple con el criterio de Factor de Seguridad por volteo.

- Verificación por deslizamiento:

Datos:

Qv 175.74 ton 175740 kg

φ 31 ° 0.541 rad
δ 21.8297 °    
 c 0.31 kg/cm2    
Ca 0.186 kg/cm2    
Qh 80 ton 80000 kg

Usaremos la fórmula:

FSdeslizamiento=
∑ F resistentes ≥ 1.5
∑ F actuantes
2
δ =tan−1 ( tan ∅) C a=0.6 c
3
2 C a=0.6∗0.31
δ =tan−1 ( tan 31 °)
3
C a=0.186
δ =21.8297

∑ F resistentes =¿ Qv ¿ tan δ+ c a B∗L¿

∑ F resistentes =¿ 175740∗tan21.8297+ 0.186∗2.2∗100∗2.9∗100 ¿

∑ F resistentes =82263.499 kg

∑ F actuantes=80000 kg

82263.499
FSdeslizamiento = ≥1.5
80000

FSdeslizamiento =1.0286<1.5

No cumple con el criterio de factor de seguridad de deslizamiento.

e)Plano (anexo)

f)
Exposición a sulfatos
La mayor parte de los procesos de destrucción del concreto causados por la formación de sales
solubles son debidos a la acción agresiva de los sulfatos. El ataque se presenta, cuando a
través del agua, concentraciones relativamente altas de sulfatos entran en contacto con los
compuestos hidratados de la pasta de cemento. Este contacto hace que se produzca una
reacción química que genera expansión en la pasta y crea una presión capaz de romperla y
finalmente desintegrar el concreto.
Según la Norma E.060 Concreto Armado, para un ataque de sulfatos de 453 p.p.m que se
encuentra entre el rango de 150 ≤ SO4< 1500 se considera una exposición de sulfatos
moderada. Un concreto expuesto a sulfatos moderada debe cumplir estos requisitos mínimos:
- Tipo de cemento: II, IP(MS), IS(MS), P(MS), I(PM)(MS), I(SM)(MS)
- Relación máxima agua – material cementante: 0.50

- f´C mínimo (MPa): 28

Otra medida de protección a la exposición moderada de sulfatos es usar un material sintético
para proteger la cimentación, este deberá ser geomembrana, geomanta o geotextil impregnado
con asfalto.
Pregunta N° 2
Diseñar una zapata de concreto armado de un edificio destinado a oficinas. La
profundidad de cimentación será Df y la distorsión angular permisible 𝛼=1/500, sobre la
zapata se apoya una columna de 0.40x0.40m.

Las cargas aplicadas a la zapara son Qv, Qh y M (Incluir en la carga vertical, además de
Qv, el peso propio de la zapata, el peso de la porción de columna y el peso de suelo que
se apoya sobre la zapata). El terreno sobre el cual se apoya la cimentación está
conformado por:

1.1 – H1 Arcilla inorgánica de plasticidad media, medianamente

compacta, húmeda, gris claro. e=0.78 𝛾1 𝛾1𝑠𝑎𝑡 qu=1.39
kg/cm 2
Cc=0.24 Cr=0.023 OCR=2.7

H1 – H2 Arcilla inorgánica de plasticidad media, medianamente compacta, saturada,

color marrón oscuro.
Gs=2.69 𝛾 = 𝛾2𝑠𝑎𝑡 g/cm3 Cc=0.195 Cr=0.02 OCR=?
H2 - ¿? Grava muy densa

Se conoce además, a partir de ensayos de corte directo que el ángulo de fricción del
terreno de apoyo es 𝜙 y la cohesión es “c” (Para el cálculo de las dimensiones de la
zapata considerar el criterio de Meyerhof y la consistencia del terreno).

NOTA: El peralte de la zapata para todos los grupos es de 0.70m.

Grupo Grupo Grupo MASW - 1 (Grupo A) MASW - 2 (Grupo B) MASW - 3 (Grupo
A B C C)
Prof (m) Vs Prof (m) Vs Prof Vs
Qv 250 275 280 (m/s (m/s) (m) (m/s)
ton ton ton )

Qh 70 ton 72 ton 80 ton 0.0 -

1.07
293.
4
0.0 - 1.07 240.
4
0.0 -
1.07
380

M 50 55 58 1.07 - 279. 1.07 - 182. 1.07 - 150

ton-m ton-m ton-m 2.31 5 2.31 4 2.31
2.31 - 328. 2.31 - 190. 2.31 - 210
Df 1.3 m 1.3 m 1.3 m 3.71 4 3.71 4 3.71
L/B 1.25 1.25 1.25 3.71 - 457. 3.71 - 479. 3.71 - 270
5.28 1 5.28 3 5.28
H1 1.3 m 1.3 m 1.3 m 5.28 - 532. 5.28 - 269. 5.28 - 245.6
H2 6.5 m 6.5 m 6.5 m 7.01 4 7.01 7 7.01
NF 0.5 m 0.6 m 0.5 m 7.01 - 367. 7.01 - 239. 7.01 - 222.7
𝛾1𝑛𝑎𝑡 16.5 16.3 16.7 kN/m3
8.90 4 8.90 1 8.90
8.90 - 509. 8.90 - 277. 8.90 - 320.9
kN/m3 kN/m3 10.96 8 10.96 4 10.96
𝛾1𝑠𝑎𝑡 18.0 18.6 18.6 10.96 - 598. 10.96 - 389. 10.96 - 403.4
13.19 9 13.19 7 13.19
kN/m3 kN/m3 kN/m3 13.19 - 582. 13.19 - 457. 13.19 - 452.4
𝛾2𝑠𝑎𝑡 19.7 19.8 19.7 15.58 4 15.58 2 15.58

kN/m3 kN/m3 kN/m3

𝜙 32 31 32

Se pide:
a) Las dimensiones óptimas de la cimentación y las dimensiones de diseño (considerando
incrementos de 0.1m)
b) La capacidad de carga última y la capacidad de carga admisible (en condiciones
estáticas según E.050 Suelos y Cimentaciones)
c) Cálculo del asentamiento por consolidación primaria debajo del centro de la
zapata (considerar 4 subestratos, se sugiere subdividir la zapata).
d) Comentar sobre el orden de magnitud alcanzado para el asentamiento de la zapata
 diseñada.
e) Indicar cuál sería el máximo asentamiento de una zapata adyacente
(asentamiento mayor a la zapata del problema) que se encuentra a una distancia
L para cumplir con la distorsión angular máxima.
f) A partir de los resultados de los ensayos MASW indicar la caracterización sísmica
del terreno tomando en cuenta los criterios de la norma E.030 Diseño
Sismorresistente y el International Building Code. Asimismo indicar un resumen
de los parámetros para el diseño sísmico para el primer caso.

 1era iteración detallada

Se tienen los siguientes datos
Df 1.3
NF 0.6
H2 5.2
H2/4 1.3
Peralte 0.7
Excentricidad 0.176979332
B asumido 4
B' 3.646041337
L 5
L' 5
L/B 1.25
Area de columna
(m2) 0.16
Altura de columna 0.6 Fórmulas empleadas:
Df/B 0.325 M
Ex= , B' =B−2 Ex
Qv
ɣnat (g/cm3) 1.63 ∅ mod=tan −1 ¿ ¿
ɣsat1 (g/cm3) 1.86 2
ɣsum1 Cmod= C
3
(g/cm3) 0.86
Q=0.1∗( Nf ∗γnat+ Peralte∗γsum1 )
ɣsat2 (g/cm3) 1.98
Qv
ɣsum2 Betha=cos−1 ( )
Qresultante
(g/cm3) 0.98
ɣconcreto
(tn/m3) 2.4
ɣagua (g/cm3) 1

Ángulos  
β 12.3722
θ 31
θmod 21.82974
Valores en
kg/cm2  
c 0.3
cmod 0.2
q 0.158

 Hallando carga vertical total

Qvtotal=qv +qsuelo nat .+q zapata+q columna=310.771 Tn
qv=275
q suelo nat =( γnat )∗( ( B∗L∗NF )−( Area col .∗Alt . col ) ) =19.404
qzapata= ( γconc )∗Peralte∗B∗L=33.6
qcolumna=( γconc )∗( Area col )∗( Alt col )=0.2304
 Hallando Carga resultante
Se tiene Qv=328.23 ,Qh=72 ( Tn ) v
Qresultante=√(328.23)2+(72)2=336.04 Tn
 En resumen

Qvtotal 328.23392
qv 275
qsuelo nat 19.40352
qzapata 33.6
qcolumna 0.2304
Qh 72
M 55
Qresultante (Qh y Qv) 336.03795

 Hallando Nq, Nc y Nγ por el método de Meyerhof

Nq 7.687345162
Nc 16.69441702
Nɣ 3.948844143

 Hallando factores de corrección por Forma (s)

 Por forma (s), utilizamos B´
B 3.665
∗N q ∗7.687
L 5.000
f cs=1+ =1+ =1.338
B Nc 16.694
3.665
f qs =1+ ∗tan ∅=1+ ∗tan 21.830=¿ 1.294 ¿
L 5.000
0.4∗B 0.4∗3.665
f ɣs =1− =1− =0.707
L 5.000

 Por Profundidad (d), Df/B<1, θmod>0, utilizamos B

 1−f qd 1−1.184
f cd =f qd− =1.184− =1.118
N c∗tan ∅ 32.671∗tan 31
f qd=1+2 tan ∅∗¿ ¿=1.103

f ɣd =1

 Por inclinación (i)

β 2 12.37 2
f ci=f qi =(1− ) =(1− ) =0.744
90 90
β 2 12.37 2
f ɣi=(1− ) =(1− ) =0.188
∅ 21.830

 Hallando qu

1
q u=c . N c . f cs . f cd . f ci + q . N q .. f qs . f qd . f qi + . B . γ 2 . N γ . f γ . f γ d . f γ i
2

1
q u=0.3 .(16.694).(1.338).(1.118 ).(0.744)+ 0.158 .(7.687).(1.294).(1.103) .(0.744)+ .(3.665).( 0.98).(3.94
2

q u=5.098 kg /cm2

 Hallando Qtransmitido
Qv 328.234 Tn
Qtrans= = =17.912 2
B ∗L' 3.66∗5
'
m
 Hallando Qadmisible
qu 5.098 kg Tn
Qadm= = =1.699 2 =16.993 2
FS 3 cm m
 Comparando Qtransmitido y Qadmisible
Se debe cumplir que Qtrans ≤ Qadm- 17.912≠ 16.993 NO CUMPLE
 Al no cumplirse la condición, debemos iterar nuevamente, cambiando el valor de B y así
poder llegar a la igualdad requerida. En este caso, al ser el Qtrans mayor al Qadm,
debemos aumentarle el valor a B para que el valor del Qtransmitido disminuya, de tal
forma que el Qadm también lo hará, pero en menor amplitud y ambos valores de
acercarán más.
Haremos e número de iteraciones que sean necesarias para llegar a un valor exacto
entre Qtrans y Qadm.
En este caso, realizamos 4 iteraciones y a continuación se muestra el cuadro resumen
de los resultados.
  1era 2da iteración 3era 4ta iteración
 iteración iteración
B 4 4.1 4.13 4.11545
3.64604133 3.748079970 3.77869636 3.78347160
B' 7 2 6 7
L=L' 5 5.125 5.1625 5.1443125
16.9928161 17.0288164
Qadm 3 17.02018542 6 17.02461
17.9124300 16.9177224
Qtrans 7 17.13843958 9 17.02461

 Resultados detallados de última iteración

qu (kg/cm2) 5.107382007
Qtrans 17.02413462
Qadm Qadm 17.0246
(kg/cm2) 1.702460669 (tn/m2) 1
El valor de B que se obtuvo fue el siguiente
B 4.11545
Por lo que llegamos a los siguientes resultados óptimos (valores exactos) y se diseñó
(aproximación al 0.1)

Respuestas metros
Bóptimo 4.11545
Lóptimo 5.1443125
Bdiseño 4.2
Ldiseño 5.2

 Para hallar el asentamiento por consolidación primaria necesitamos los esfuerzos

efectivos en A, B, C y D.

σ ´ A =q+ γ 2∑ ¿∗H =0.158+0.98∗0.65∗0.1=0.2217 kg/ cm2 ¿

A

σ ´ B=q +γ 2∑ ¿∗H =0.158+0.98∗( 1.3+ 0.65)∗0.1=0.3491 kg/cm2 ¿

B

σ ´ C =q+ γ 2 ∑ ¿∗H =0.158+0.98∗( 1.3+1.3 +0.65) ∗0.1 =0.4765kg/ cm2 ¿

C

σ ´ D=q+ γ 2 ∑ ¿∗H =0.158+ 0.98∗( 1.3 +1.3+1.3+0.65) ∗0.1=0.6039 kg/ cm2 ¿

D

Presiones efectivas kg/cm2 Tn/m2

Ꝺ'A 0.2217 2.217
Ꝺ'B 0.3491 3.491
Ꝺ'C 0.4765 4.765
Ꝺ'D 0.6039 6.039
 Con la Teoría de Boussinesq calculamos las presiones transmitidas a las profundidades
A, B, C y D bajo el centro de la zapata.
Como nos piden los incrementos de presiones al centro de la zapata, calcularemos con
un B´ y L´.

B=4.2 m
L¿ 5.2 m

B 4.2
B´ =b= = =2.1m
2 2
L 5.2
L´ =a= = =2.6 m
2 2

a 2.6
M= = =
b 2.1
1.2380
Para el cálculo de los incrementos de presiones efectivas en A, B, C y D.

Para ∆ σ A :

z 0.65
N= = =0.3095
b 2.1

A=M∗N=1.238∗0.3095=0.3832

2 2 2 2
B=1+ M ∗N =1+ ( 1.2380 ) ∗( 0.3095 ) =2.6287

N 2∗B−M 2=( 0.3095 )2∗2.6287−( 1.2380 )2=−1.2810

N 2∗B−M 2< 0
−1.2810<0

2 A √B
W =arctg( )
N ∗B−M 2
2

W ´ =W + π
2A
∗B+ N 2
q 2
N ∗B+ M 2
2 A√B
σ A = trans ( + arctg( 2 ))
4π √B N ∗B−M 2
 2∗0.3832
∗2.6287+ 0.30952
q 2
0.3095 ∗2.6287+1.2380 2
2∗0.3832 √ 2.6287
σ A = trans ( +arctg ( )+ π)
4π √2.6287 0.30952∗2.6287−1.23802

∆ σ A=4.1903 ton /m 2

Para ∆ σ B :

z 1.95
N= = =0.9286
b 2.1

A=M∗N=1.238∗0.9286=1.1497

B=1+ M 2∗N 2 =1+ ( 1.2380 )2∗( 0.9286 )2=3.3951

N 2∗B−M 2=( 0.9286 )2∗3.3951−( 1.2380 )2=1.3945

N 2∗B−M 2< 0

1.3945>0

2A
∗B+ N 2
q trans N ∗B + M 2
2
2 A √B
σ B= ( +arctg ( 2 ))
4π √B N ∗B−M 2

1.7024
σ B= ¿
4π

∆ σ B =3.3109ton /m2

Para ∆ σ C:
 z 3.25
N= = =1.5476
b 2.1

A=M∗N=1.238∗1.5476=1.9161

B=1+ M 2∗N 2 =1+ ( 1.2380 )2∗( 1.5476 )2=4.9280

N 2∗B−M 2=( 1.5476 )2∗4.9280−( 1.2380 )2=10.2703

N 2∗B−M 2< 0

10.2703>0
2A
∗B+ N 2
q 2
N ∗B+ M 2
2 A √B
σ C = trans ( + arctg( 2 ))
4π √B N ∗B−M 2
1.7024
σC = ¿
4π

∆ σ C =2.2214 ton /m2

Para ∆ σ D :

z 4.55
N= = =2.1667
b 2.1

A=M∗N=1.238∗2.1667=2.6825

2 2 2 2
B=1+ M ∗N =1+ ( 1.2380 ) ∗( 2.1667 ) =7.2273
2 2 2 2
N ∗B−M =( 2.1667 ) ∗7.2273−( 1.2380 ) =32.3953

N 2∗B−M 2< 0

32.3953>0
2A
∗B+ N 2
qtrans N ∗B+ M 2
2
2 A √B
σ D= ( +arctg ( 2 ))
4π √B N ∗B−M 2
qtrans
σ D= ¿
4π
 ∆ σ D =1.4764 ton /m2

  A B C D
0.92857142 2.16666666
N 0.30952381 9 1.547619048 7
1.14965986 2.68253968
A 0.383219955 4 1.916099773 3
3.39512471 7.22732426
B 2.628684807 7 4.928004535 3
1.39454914 32.3953924
n^2*B-M^2 -1.2810387 6 10.27030565 2
n^2*B-M^2 <0 >0 >0 >0

Ahora, multiplicamos cuatro veces el incremento de presiones ya que nuestros cálculos fueron
para la cuarta parte de la zapata.

Incrementos de
presión Tn/m2 x4 Tn/m2
4.19027316 16.761092 16.761092
ΔꝹA 6 7 7
3.31085128 13.243405 13.243405
ΔꝹB 1 1 1
2.22138523 8.8855409 8.8855409
ΔꝹC 5 4 4
1.47639940 5.9055976 5.9055976
ΔꝹD 7 3 3

El incremento de presiones más los esfuerzos efectivos en A, B, C y D.

∆ σ A +σ ´ A=18.9781 ton /m 2

∆ σ B + σ ´ B=16.7344 ton/m 2

∆ σ C +σ ´ C =13.6501 ton /m2

∆ σ D + σ ´ D=11.9446 ton/m2

Como el OCR del suelo dos (OCR=1.5) es mayor a uno, el suelo es una arcilla preconsolidada.
 OCR>1 → Preconsolidado

Ahora, calculamos el esfuerzo de preconsolidacion de A, B, C y D.

σ ´P
OCR=
σ ´V

σ ´P
OCR= =1.5
σ ´V

OCR=σ ´ PA =1.5∗σ ´ A =3.3255 ton/m2

2
OCR=σ ´ PB=1.5∗σ ´ B=5.2365 ton/m

OCR=σ ´ PC =1.5∗σ ´ C =7.1475 ton /m2

OCR=σ ´ PD=1.5∗σ ´ D=9.0585ton/m2

OCR 1.5
OCR=1.5>1
- PC   Ꝺ'P
Ꝺ'A+ΔꝹA 18.97809266 3.3255
Ꝺ'B+ΔꝹB 16.73440512 5.2365
Ꝺ'C+ΔꝹC 13.65054094 7.1475
Ꝺ'D+ΔꝹD 11.94459763 9.0585
Conclusión Ꝺ'+ΔꝹ>Ꝺ'P

Como el usaremos esta fórmula para hallar los asentamientos por

σ ´V + ∆ σ > σ ´ P
consolidación primaria en A, B, C, y D.

σ ´V + ∆ σ > σ ´ P

Hi
δ P= ¿
1+ e0

C C =1.95 e 0=0.7744
 C r=0.02 H i=1.3

Hi
δ PA = ¿
1+e 0

Hi
δ PB = ¿
1+e 0

Hi
δ PC = ¿
1+e 0

Hi
δ PD = ¿
1+ e0

El asentamiento por consolidación primaria total es la suma de los asentamientos por

consolidación primaria en A, B, C y D.

δ T =δ PA +δ PB +δ PC +δ PD = 25.4958 cm

Calculo del asentamiento instantáneo.

Asentamiento individual:
Se considera las dimensiones de diseño
B=4.2 m
B´ =B−2e x
B´ =4.2−2(0.1775)¿
B´ =3.845 m L ´=5.200 m
e x =0.1775 m

μ=0.33

Para calcular la carga transmitida ¿ ¿) se considera las dimensiones efectivas.

q trans =1.702 kg/cm 2

El módulo de elasticidad del suelo se obtiene del gráfico de módulo de corte y deformación.

γ 1.98
( )
Gmax = 2 sat V S2 =
g 9.81
x(389.7508)2=11779.456 ton/m2
Gs=10 % G max=10 % ¿) = 1177.9456 ton/m2

ton 2
(
Es=2 ( 1+ μ ) Gs=2 ( 1+0.33 ) 1177.9456
m )
=3133.3352 ton/m 2
 30
Vs=
1.07 1.24 1.4 1.57 1.73 1.89 2.06 2.23 2.39 2.55 2.72 2.89 3.05 3.21
+ + + + + + + + + + + + +
240.4 182.4 190.4 479.3 269.7 239.1 277.4 389.7 457.2 632.1 596.7 611 642 742

Vs=389.7508 m/ s

Hallando el If para hallar el asentamiento instantáneo

Consideramos una zapata rígida, hallando asentamiento en el centro, tal como lo indica el
enunciado del ejercicio.
L/B (diseño)=5.1/4.1=1.24
Interpolamos en tabla

−( 1.06−0.82 )∗( 1.5−1.24 )

x= + 1.06
1.5−1
x=0.9352=If

Aplicando la fórmula de asentamiento instantáneo.

qtrans∗B´ ( 1−μ2 )
∆ H= ∗IF
Es

1.702∗3.8450 ( 1−0.332 )
∆ H= ∗0.9352
3133.3352

∆ H =0.0017 cm
El asentamiento total es la suma de asentamiento por consolidación primaria y del
asentamiento instantáneo.
 ∑ Asentamientos=δ I +δ P

∑ Asentamientos=25.4958+0.0017
∑ Asentamientos=25.4975 cm

Distorsión angular:
Para calcular el asentamiento mayo adyacente se usa esta fórmula:

δ 1=0.254975 m
δ 2−δ 1
α= ≤1 /500
L L=5.90 m

δ 2−0.254975
α= ≤ 1/500
5.9

δ 2=0.266775 m

 Comentario del orden de magnitud alcanzado para el asentamiento de la zapata

diseñada.
El asentamiento total alcanzado para la zapata diseñada es de 25,498cm, lo cual significa
que la zapata analizada, al llegar al 100% de su asentamiento después de haberse
producido el asentamiento inicial al haberle ejercido cargas al suelo que la contiene, es el
mencionado anteriormente.
Con este valor se concluye que, con las dimensiones de diseño halladas, el asentamiento
generado no supera los valores máximos de asentamiento que podrían provocar daños
mayores al suelo. Por otro lado, con el asentamiento hallado no se verá afectada la
estructura que soporte dicha zapata, lo cual quiere decir que esta no sufrirá daños (grietas,
que no causen daño alguno, e inclinaciones de terreno) ni fallas.

2f
Resultados de ensayo MASW (caracterización sísmica del terreno)
Antes de proceder al análisis de las velocidades de ondas S, consideramos importante resaltar
que la forma de la curva de dispersión o diferencia entre velocidades según profundidad
generada por el ensayo presentado, está fuertemente relacionada a la variación de la rigidez
con la profundidad (hasta los primeros 30m), por lo tanto el método MASW tiene la capacidad
de identificar estratos más blandos entre estratos rígidos.
Según el promedio armónico determinado (fórmula en Excel) dentro de la profundidad permitida
por el ensayo MASW(30m), se observa que su valor es de 389.75 m/s, lo cual indica, según la
tabla incluida en la norma E.030 Diseño Sismo resistente y el International Building Code, que el
suelo en análisis pertenece al tipo S2. El suelo de tipo S2 es, según la norma mencionada, un
“suelo muy denso o roca blanda”, es decir contiene, mayor cantidad de partículas gruesas que
finas, lo cual significa que podría ser más permeable.
Por otro lado, analizando el promedio armónico hallado (fórmula en Excel) dentro del estrato
que se encuentra por debajo de la zapata en análisis (241,58 m/s), concluimos que esta
“porción” de suelo pertenece, según la tabla incluida en la norma, al tipo S2 también, pues su
promedio armónico se encuentra dentro del rango indicado en ella. Por lo tanto, cuenta con las
mismas características mencionadas en el párrafo anterior.
A continuación, se muestran tabla y definición de suelo de tipo S2, extraídos de norma E.030
Analizando la zona en donde se encuentra el suelo hallado y considerando que nuestro suelo
se encuentra en lima, obtuvimos según la norma mencionada, que este se encuentra en la zona
4, con un factor z de 0.45. Dicho factor representa “la aceleración máxima horizontal en suelo
rígido con una probabilidad de 10 % de ser excedida en 50 años.”

Finalmente, como resumen de los parámetros hallados según las condiciones del suelo en
análisis, podemos determinar los parámetros de sitio (S, Tp, Tl)
 S Tp Tl
1.05 0.6 2.0

Pregunta N° 3
En el estribo derecho del puente mostrado en el Anexo 2, se requiere transmitir al
terreno un total de 650 ton (incluye peso propio del encepado y la mitad del metrado de
cargas del puente), se le solicita, utilizando el método general de Meyerh, calcular el
número mínimo de pilotes de L metros de longitud (medidos desde el nivel de fondo del
encepado NFE) para soportar dicha carga,

DATOS DEL PROBLEMA

GRUPO B
L(m) 15
δ 0.8∅
Ca Interpolando
diapositiva La interpolación se mostrará en las
páginas siguientes
K 0.8Ko
Eficiencia 0.65
Diámetro(m) 0.75
NFE(m) 3.5
Prof(m) c(Kpa) ∅ (°) g
γ( )
cm3
0.00-1.20 0 29 1.6
1.20-2.40 0 30 1.7
2.40-3.45 6 30 1.8
3.45-4.40 6 32 1.85
4.40-5.00 6 31 1.85
5.00-8.50 28 0 1.9
8.50-11.50 12 33 1.9
11.50-13.40 8 31 1.95
13.40-16.00 10 33 2.2
16.00-20.00 15 35 2.2
 Prof(m)
3.5-4.40
4.40-5.00

8.50-11.50
11.50-13.40
13.40-16.00
*Primero hallamos la carga lateral Qs

 En el caso de suelos granulares se plantea lo siguiente:

 En el caso de suelos
cohesivos se plantea lo siguiente:

 Como el NFE está a 3.5m, entonces se comienza a analizar a partir de esa cota hasta los
16.00 m
 Los estratos con suelo granular se presentan a continuación.
  El perímetro será igual a 2∗π∗r=2.3562
 A continuación se mostrarán los ∆ Lpara cada estrato

∆L(m)
0.9
0.6
3.5
3
1.9
2.6

Ahora hallamos el K para cada estrato

 K=0.8∗Ko

 Ko=1−seno (∅ )

K0

∅ (°)

32 K
31 0.470080736 0.376064589
0 0.484961925 0.38796954
33 - -
31 0.455360965 0.364288772
33 0.484961925 0.38796954
0.455360965 0.364288772

Ahora hallamos el σ v para cada estrato

  σ v =γ∗Profundidad desde el 0.00m hastala mitad del estrato a analizarse∗0.1

 Se le multiplica por un factor de 0.1 para que las unidades estén en kg/cm2

g Profundidad desde el σ v (kg/cm2)

γ( )
cm3 0.00m hasta la mitad
de cada estrato
1.85 3.95 0.73075
1.85 4.7 0.8695
1.9 6.75 1.2825
1.9 10 1.9
1.95 12.45 2.42775
2.2 14.7 3.234

El siguiente paso es hallar el tan ( δ ) para cada estrato

 Se mostrará a continuación los valores de δ para cada estrato

 δ =0.8∗∅

 Luego se procederá a hallar la tan(δ )

δ
δ

25.6
24.8
0
26.4
24.8
26.4

tan

0.479119721
Ahora, ya que se tiene el 0.46206487 perímetro, ∆ L y f podemos hallar el Qs ,
mediante una sumatoria de los - valores calculados por cada estrato.
0.496404311
 Como se había 0.46206487 mencionado antes
0.496404311 ∑ Qs=∑ ( perimetro)∗∆ L∗f
  f =k∗σ v∗tanδ

 La sumatoria por cada estrato se muestra a continuación:

∑ Qs( ton)

Tipo de
suelo
2.792093602
granular
2.203604088
Resultado de la granular
25.729704
෍ ܳ‫ݏ‬ሺ‫݊݋‬݊ cohesivo
‫ݐ‬ ሻ
sumatoria :
24.28669045
para los estratos congranular
19.48364896
suelos granulares
=110.322
granular
35.82670316
granular

Hay un único estrato el cual presenta suelo cohesivo, para su cálculo de Qs será de la siguiente
manera:

 ∑ Qs=∑ ( perimetro)∗∆ L∗S

 S=Ca

 El valor de Ca, lo hallamos interpolando teniendo como referencia a Cu=0.375 con un

Ca= 0.35 y Cu=0.28 (dato) y “X” que es el valor que deseamos obtener

0.28−0.375 x−0.35
=
0.28−0.75 x−0.5

x=0.3120
  El perímetro es 2.3562 m
 El ∆ L correspondiente a ese estrato es 3.5m
 Por lo tanto, el valor de Qs para este estrato de suelo cohesivo será igual a:

2.3562∗3.5∗0.3120=25.729704
Finalmente, como ya se tienen los valores para los suelos granulares y cohesivos, se procederá
a sumar para obtener el ∑ Qs

 La ∑ Qs para suelos granulares es 84.59

 La ∑ Qs para suelos cohesivos es 25.729704

∑ Qs(ton)
 Sumando estos valores obtendremos la sumatoria de carga lateral Qs

2.792093602
2.203604088
25.729704
24.28669045
19.48364896
35.82670316
110.322
Sumatoria de todos los estratos

*Segundo, hallamos la carga por punta “Qp”

 Según Meyerhof Qp se calcula de la siguiente manera:

Q p= A p∗q p= A p (c '∗N c +q' ∗N q )

 En el caso de arenas (c=0 Kpa):

Q p= A p∗q p= A p ( q'∗N q )
  En el caso de arcillas (ø=0°):

Q p= A p∗9∗Cu

 Para hallar el Qp, se debe analizar el ultimo estrato, es decir en este caso de 16.00
-18.50 m, sin embargo, necesitamos saber el q’ (kg/cm2) de cada estrato, ya que se
necesita los acumulados.

 q’=γ ( cmg )∗0.1∗∆ L ( m )+ q ' anterior

3

 Diámetro=0.75m
 Ap(m2) =0.4417 m2
 Nq*=50
 Qp=q’*Nq*
g
γ( )
cm3
∆L(m) ∆L(m) q'(kg/cm2)
Prof(m)
1.2 1.6 1.2 0.192
0.00-1.20
1.2 1.7 1.2 0.396
1.20-2.40
1.05 1.8 1.05 0.585
2.40-3.45
0.95 1.85 0.95 0.76075
3.45-4.40
0.6 1.85 0.6 0.87175
4.40-5.00
3.5 1.9 3.5 1.53675
5.00-8.50
3 1.9 3 2.10675
8.50-11.50
1.9 1.95 1.9 2.47725
11.50-13.40
1.6 2.2 2.6 3.04925
13.40-16.00
2.5 2.2 2.5 3.59925
16.00-18.50

 Por lo tanto, Qp = Ap*qp = 795.0.49 (ton)

 Aquí se presenta la gráfica para el cálculo del Nq* para el último estrato, tiene un ∅=35 °
 La relación longitud/ancho = 15/0.75= 20
Como se ve en la
gráfica le
correspondería
un Nq*= 50 para
un ángulo de
fricción 35°

 Pero,
como ya tenemos nuestro Qp = 795.0499 (ton), este no cumpliría con esta restricción el
cual dice 𝑸𝒑 = 𝑨𝒑. 𝒒𝒑 = 𝑨𝒑. (𝒒 ′. 𝑵𝒒 ∗) ≤ 𝑨𝒑. 𝒒𝟏
𝒒𝟏 (𝒌𝑵 𝒎𝟐) = 𝟓𝟎. 𝑵𝒒 ∗. 𝒕𝒂𝒏∅

 El valor de la restricción es = 773.356 KN=78.833(ton), por ende, tomo este valor “78.833
ton” como mi Qp.
 El siguiente paso para hallar el número de pilotes sería calcular el Qu
 Qu= Qp+∑Qs=78.833+110.322=189.155
 Ahora se halla El Qadm, el cual tienes 2 formas de hallarlos, pero se escoge el menor
1. Qadm= Qu/2.5=75.662
2. Qadm=(∑Qs/1.5) +(Qp/3) =99.826
 Si colocamos un número de pilotes= 14, si cumpliría con el Qgrupo, el cual es:
Qgrupo=Eficiencia*Qu*Número de pilotes=, y con este valor si supera los 650 ton a las
que se desea transmitir.
*OBSERVACIÓN
 Si consideramos que hay presencia de un río debajo del puente, entonces tomaríamos
nuestro σ con un γsumergidopara cada estrato.
 La forma de calcular ∑Qs y Qp no variará

 σ v =0.1∗γ −1∗Prof desdeel 0.0 mhasta la mitad del estrato a analizar

 Eficiencia=0.65

g
γ( ) σ v (kg/cm2)
cm 3
Prof(m)
1.85 0.33575
3.5-4.40
1.85 0.3995
4.40-5.00
1.9 0.6075
5.00-8.50
1.9 0.9
8.50-11.50
1.95 1.18275
11.50-13.40
2.2 1.764
13.40-16.00

 Para hallar los valores de Ko, K, tanδ y f se emplea las mismas mecánicas de lo antes
calculado sin considerar presencia de un río.

∑ Qs( ton)=f∗perímetro∗∆ L

∅ (° )
δ K0
Ca(kg/c
m2) K
0.376064
25.6 0.470080736 589
32
0.387969
24.8 0.484961925 54
31
0 0.3120 - -
0
0.364288
26.4 0.455360965 772
33
 0.387969
24.8 0.484961925 54
31
0.364288
26.4 0.455360965 772
33
0.341138
28 0.426423564 851
35

δ
tan
f(kg/cm2)
∆L(m) ∑ Qs(ton)
0.479119721 0.060495422
0.9
1.282853817
0.46206487 0.071617204
1.012466743
0.6
25.729704
- -
3.5
11.50422179
0.496404311 0.162751065
9.492034111
3
19.54183809 0.46206487 0.212028157
1.9
68.56311
0.496404311 0.318992088
2.6
∑Qs =
68.56311 Ton
Por otro lado, para la carga por punta Qp:
 Siguiendo la misma mecánica, para hallar el Qp anterior, el q’ se halla desde el nivel 0.0
hasta el ∆ L para cada estrato *γ sum*0.1 , este q’ es acumulado.
 A continuación, se muestran los q’ de cada estrato, el más importante es el último estrato
para hallar el Qp.
 El Nq es el mismo = 50 y la gráfica se mostró anteriormente
 qp(kg/cm2) =87.4625
 Qp(ton)= Ap*qp*10 = 386.397

g
γ( )
cm 3
∆L(m) q'(kg/cm2)
Prof(m)
1.6 1.2 0.072
0.00-1.20
1.7 1.2 0.156
1.20-2.40
1.8 1.05 0.24
2.40-3.45
 1.85 0.95 0.32075
3.45-4.40
1.85 0.6 0.37175
4.40-5.00
1.9 3.5 0.68675
5.00-8.50
1.9 3 0.95675
8.50-11.50
1.95 1.9 1.13725
11.50-13.40
2.2 2.6 1.44925
13.40-16.00
2.2 2.5 1.74925
16.00-18.50

 Pero, como ya tenemos nuestro Qp = 386.397 (ton), este no cumpliría con esta
restricción el cual dice 𝑸𝒑 = 𝑨𝒑. 𝒒𝒑 = 𝑨𝒑. (𝒒 ′. 𝑵𝒒 ∗) ≤ 𝑨𝒑. 𝒒𝟏
𝒒𝟏 (𝒌𝑵 𝒎𝟐) = 𝟓𝟎. 𝑵𝒒 ∗. 𝒕𝒂𝒏∅
 El valor de esta restricción es 773.35KN=78.833 ton
 El siguiente paso sería calcular Qu = Qp+∑Qs =147.39
 Ahora se halla El Qadm, el cual tienes 2 formas de hallarlos, pero se escoge el menor
1.-Qadm= Qu/2.5=58.958
2.-Qadm=(∑Qs/1.5) +(Qp/3) =71.98

 Si colocamos un número de pilotes= 17, si cumpliría con el Qgrupo, el cual es:

Qgrupo=Eficiencia*Qu*Número de pilotes=651.49, y con este valor si supera los 650 ton
a las que se desea transmitir.

Pregunta 4
 Cálculos de las primeras dos filas de la tabla
1ra FILA
Profundidad: 1
Profundidad real: 1.3
CF (%)=20
γnat =1.65
γsum=0.82
Sigma prima=( Prof ∗γnat + ( Profreal−Prof )∗γsum ) x 0.1=0.1896
Sigma total=( Prof 1 m∗γnat∗0.1 )+ ( Prof real 1.3− prof 1 m )∗γsum∗0.1=0.2196
rd=1−0.00765 xZ=1−0.00765 x 1.3=0.990055
1 1 /2
Cn=( ' ) ; ≤1.7
σ
Cr: De acuerdo a la longitud de barra <3m corresponde 0.75
Cb= 1
N1(60)=Nspt∗Cn∗Cr∗Cb=6.375
amax σ
CSR=0.65∗ (g )( )
∗ ' ∗rd =0.327958852
σ
190
Alfa=e (1.76− ( CF ))=3.614668
2

CF 1.5
Betha=(0.99+
1000 ( )
)=1.07944

N160cs= Alfa+ Betha∗N 1(60)=10.4961

Factor de corrección Mw= 0.967
1 ( N 1 ) 60 cs 50 1
CRR sin corregir= + + − =0.117518534
34−( N 1 ) 60 cs 135 (10∗( N 1 ) 60 cs + 45) 200
2

CRR corregido=Crr sin corregir∗Factor de corrección Mw=0.113640423

FS=CRRcorregido /CSR=0.3465082
amax
a/g=( )∗FS=0.1524636
g
PL=¿=0.6283293

2da fila
Profundidad=2
Profundidad real=2.3
CF (%)=20
Sigma prima= ( Profxγnat+ ( Profreal−Prof )∗γsum ) x 0.1=¿0.2716
Sigma total=σ T 1+ ( Profreal 2.3−Prof real 1.3 ) x 1.82 x 0.1=0.4016
Rd=1−0.00765 xZ=1−0.00765 x 1.3=0.982405
1 1 /2
Cn=( ' ) ; ≤1.7
σ
Cr: De acuerdo a la longitud de barra <3m corresponde 0.75
Cb=1
N1(60) =Nspt∗Cn∗Cr∗Cb= 5.1
amax σ
CSR=0.65∗ ( g )( )
∗ ' ∗rd =0.361658149
σ
190
Alfa= e (1.76− ( CF ))=3.614668
2

CF 1.5
Betha=(0.99+
1000 ( )
)=1.0794427

N160cs= Alfa+ Betha∗N 1(60)=9.119825825

Factor de corrección Mw= 0.967
1 ( N 1 ) 60 cs 50 1
CRR sin corregir= + + − =0.105442333
34−( N 1 ) 60 cs 135 (10∗( N 1 ) 60 cs + 45) 200
2

CRR corregido=Crr sin corregir∗Factor de corrección Mw=0.101962736

FS=CRRcorregido /CSR=0.2819313
amax
a/g==( )∗FS= 0.1240498
g
1
PL= =0.6328684
¿¿

 Cálculo del PGA (empleando escalimetro 1:200)

PGA más
próximo 0.44
 Ubicación del ensayo SPT (11.01° Latitud, 77.62° Longitud)
 Cálculos del factor de corrección para el CRR corregido.

0.96
7
La magnitud del sismo es de 7.8, por lo que interpolamos en la tabla mostrada mediante los
valores entre M=8-1/2 Y M=7-1/2, mediante la siguiente expresión:
7.5−7.8 1−x
=
7.5−8.5 1−0.89
x=0.967

Pregunta 5
 Datos a emplear:
1.-Diámetro(m)=1.1
2-NFC (m) = 8
No se toma en consideración el último estrato: 30.00-48.00

Prof(m) g
∅ (° ) γ( )
cm3
 28 1.55
0.00-3.00
26 1.5
3.00-8.00
34 1.9
8.00-12.00
36 2.0
12.00-30.00

Comenzaremos hallando ∑Qs:

σ v (kg/cm2) = 0.1*Prof desde el 0.0 hasta la mitad del estrato*γ

σ v (kg/cm2)

0.2325
0.825
1.9
4.2

Ahora hallamos el K para cada estrato:

 K=0.855∗Ko

 Ko=1−seno (∅ )

K0

K
0.530528 0.453601
44 81
0.561628 0.480192
85 67
0.440807 0.376890
1 07
0.412214 0.352443
75 61
Ahora hallaremos f (kg/cm2) ) = K*σ v*tanδ

 δ =0.9∗∅ (considero concreto prefabricado/arena)

f(kg/cm2)
0.049626
85
0.171433
29
0.423495
11
0.939403
57

Ahora hallaremos ∑Qs (ton)= Perímetro* ∆ L*f*10

∑ Qs( ton)
5.1449423
3
29.621545
6
58.539761 Ahora hallaremos la carga por punta Qp (ton)
9
584.34222  Se considera ∑Qs(ton)=642.881 para el q’(kg/cm2) desde el nivel
2 0.0 hasta la profundidad donde acaba cada
642.881 estrato, este q’ es acumulado.
 Se analiza el último estrato, el cual es el más importante
 q’= γ∗∆ L∗0.1
 Diámetro=1.1 m
 Ap(m2) = 0.9503
 Nq*=95 (según la gráfica adjuntada), ángulo de fricción (último estrato) = 38°
 qp(kg/cm2) =q’ *Nq = 905.825
 Qp(ton)=Ap*qp*10=8608.34
 ∆ L( m)
g
q’(kg/cm2) γ( )
Prof(m) cm 3
1.55 3 0.465  Pero,
0.00-3.00 como ya
1.5 5 1.215 tenemos
3.00-8.00
1.9 4 1.975 nuestro
8.00-12.00 Qp =
2.0 18 5.575 86048.34
12.00-30.00
2.2 18 9.535 8608.342825 (ton), este
30.00-48.00 no
Qp=86048.34 Ton

cumpliría con esta restricción el

cual dice 𝑸𝒑 = 𝑨𝒑. 𝒒𝒑 = 𝑨𝒑. (𝒒 ′.
Nq*=95 �
� Cálculo del Nq* mediante la gráfica, se
� utilizó una relación 20 y un angulo de
� fricción 38°.
Corresponde un Nq =95

∗) ≤ 𝑨𝒑. 𝒒𝟏
𝒒𝟏 (𝒌𝑵 𝒎𝟐) = 𝟓𝟎. 𝑵𝒒 ∗. 𝒕𝒂𝒏∅
 El valor de la restricción es = 3526.782 KN/m2=359.508 Ton

Prueba de carga bidireccional en un pilote de prueba construido con el mismo diámetro de

diseño y a la misma profundidad de cimentación proyectada con una longitud de 40m y los
cálculos correspondientes fueron los siguientes:
Ap(m2) =0.9503
Perímetro=3.455
Recordar: Qs =f*perímetro* Δl

Prof(m Def.Pil(mm) Fric.fust(kPa) Qp(tn) f(kg/cm2) Δl(m Qs(tn)

) )
0-4.5 0.2 560.6 0.11 4.5 17.1059
5
4.5-12 0.9 0.48 7.5 124.407
12-20 2.7 0.98 8 270.93
20-32 4.5 1.35 12 559.83
32-40 5900 8
Sumatoria Qs972.275

5900*0.9503*10^-2*10

Respecto para el cálculo de Qadm se debe tener estos valores

Qu=Qp+∑Qs=1331.783 ton
Qp=359.508 ton
Qs=972.275 ton
Qdm1=(Qs/1.5) +(Qp/3) =532.71
Qadm2=(Qs/1.5) +(Qp/3) =768.02
Se escoge el Qadm1 (el menor)

 Como el Qadm1 es menor a mi Qp (tn) =560.65 sí cumple y este tiene más jerarquía que
el ∑
Pregunta N° 6
Un depósito rectangular de 55 x 12 m se apoya superficialmente sobre un terreno de
40m de espesor con un peso específico de 𝛾𝑛𝑎𝑡, un módulo de deformación E’ y un
coeficiente de Poisson de n. La cohesión c’ y el ángulo de rozamiento f. El nivel freático
se encuentra a NF m bajo la superficie.
Gru Gru Gru
po A po B po C
NF (m) 4.2 3.6 3.2
h (m) 15 16 17
E' 300 270 280
(kg/cm
2
)
 0.25 0.29 0.27
 2.05 2.05 2.1
sat(ton/
m2 )
 1.9 1.95 1.98
(ton/m2
)
' (°) 31 30 32
C' 19 20 18
(kPa)

El depósito está lleno de agua hasta una altura de h metros sobre

el nivel de cimentación. Aplicando el método de los elementos

finitos se le solicita:

a) Los asentamientos que se producen inmediatamente debajo de la

cimentación, considerando que esta fuera completamente flexible (caso
hipotético).
b) Los asentamientos que se producen inmediatamente debajo de la cimentación
considerando una losa de 45 cm de espesor (Econcreto=250000 kg/cm2).
c) Identificar los incrementos de presión a 1m, 2m, 4m, 6m y 8m debajo del centro del
nivel de cimentación.
d) Utilizando la formulación propuesta por Boussinesq identificar los incrementos
de presión a las mismas profundidades indicadas en el ítem c) y luego
comparar y comentar.

Deberá explicar claramente el procedimiento utilizado en cada etapa de cálculo.

SOLUCION

d) Utilizando la fórmula de Boussinesq determinaremos los incrementos de presiones en las

profundidades de 1m, 2m, 4m, 6m y 8m.
 B=12 m
L¿ 55 m

B 12
B´ =b= = =6.0 m
2 2
L 55
L´ =a= = =27.5 m
2 2

a 27.5
M= = = 4.583
b 6

2A
∗B+ N 2
q 2
N ∗B+ M 2
2 A√B
∆ σ= trans ( + arctg( 2 ))
4π √B N ∗B−M 2

Para ∆ σ A (z= 1m):

z 1
N= = =0.167
b 6

A=M∗N=4.583∗0.167=0.7639

2 2 2 2
B=1+ M ∗N =1+ ( 4.583 ) ∗ ( 0.167 ) =22.395

2 2 2 2
N ∗B−M =( 0.167 ) ∗22.395−( 4.583 ) =−20.395

N 2∗B−M 2< 0
−20.395<0

2 A √B
W =arctg( )
N ∗B−M 2
2

W ´ =W + π
2A
∗B+ N 2
q trans N 2∗B+ M 2 2 A√B
σ A= ( + arctg( 2 + π ))
4π √B N ∗B−M 2
 2∗0.3832
∗2.6287+ 0.30952
q trans 0.3095 ∗2.6287+1.23802
2
2∗0.3832 √ 2.6287
σ A= ( +arctg ( )+ π)
4π √2.6287 0.30952∗2.6287−1.23802
∆ σ A=39.924 kN /m2

Siguiendo este procedimiento calculamos ∆ σ B ( z=2 m ) , ∆ σ C ( z=4 m ) , ∆ σ D ( z=6 m ) y ∆ σ C ( z=8 m ) .

Después, multiplicamos por cuatro ya que la Teoría de Bousinessq solo se aplica para el
incremento de esfuerzos en una esquina y para el ejercicio estamos calculamos al centro, por
eso se dividió en cuatro áreas para aplicar Bousinnesq.
En resumen:

Increment
os de
presión KN/m2 x4 KN/m2
39.92375 159.6950 159.6950
ΔꝹA 42 17 17
39.44457 157.7783 157.7783
ΔꝹB 81 12 12
36.76754 147.0701 147.0701
ΔꝹC 33 73 73
32.69332 130.7733 130.7733
ΔꝹD 54 02 02
28.52226 114.0890 114.0890
ΔꝹE 34 54 54
PREGUNTA 7
 Se tienen los siguientes datos:

betha 12
X1(m) 0.8 '
X2(m) 1.05 δ∅ '1
∅2
X3(m) 2.6  Hallando Ea y Ep
X4(m) 0.55
H2(m) 6.7
H3(m) 0.7
H4(m) 1.5
Yc(ton/m3
) 2.4
Y1(g/cm3) 1.85
Y2(g/cm3) 1.98

31
  29
c'1(Kpa) 0
c'2(Kpa) 33

17
Ca 0.85c'

 Se obtiene Ka y Kp mediante los siguientes datos

Ka 0.338
fi 31
w 0
behta 12
H1 0.55264706
 Kp 4.477
fi 29
w tan-1(x2-x4/h2) = 4.268
betha 0

 Reemplazando en las fórmulas de Ea y Ep obtenemos

kg kg
Ea=γ 1∗( H 1+ H 2+ H 3 )∗Ka=0.4973 2 ,
Ep=γ 2∗ ( H 4 )∗Kp=1.329 2
cm cm
 Hallando Distancia desconocida de cateto en triángulo formado por Ea
x=sin ( 90−δ )∗( H 1+ H 2+ H 3 )=7.605

 Hallando el Ea y Ep resultante
Ea∗x∗10 ton Ep∗H 4∗10 ton
Earesultante= =18.62 2 , Ep resultante= =9.973 2
2 m 2 m
 Con los datos anteriores y con ayuda del grafico obtenemos la siguiente tabla:
F
horizonta Mom. Mom.Desestabilizant
N° nombre F vertical l Brazo Estabilizante e
1 Ea horiz   17.806 2.650   47.202
2 Ea vertical 5.444   4.45 24.226  
3 E pasivo   9.973 0.5 4.98  
4 peso 1 8.844   1.575 13.9293  
5 peso 2 7.476   2.225 16.6341  
1.1333333
6 peso 3 4.02   3 4.556  
7 peso 4 32.227   3.15 101.51505  
3.5833333
8 peso 5 1.329   3 4.762  
∑ 59.340 170.609 47.202

 Verificación por volteo

∑ M . estabilizantes =3.614>2 , CUMPLE

∑ M . desestabilizantes
 Verificación por deslizamiento

∑ FH actuante ( ton )=Ea horizontal(no consideramos el Ep) =17.806

∑ FH resistente ( ton )=¿ ¿ ∑Fv.tg∅+cb (pasivo)= 47.577 , ∑Fv.tgδ+ca.b (activo)= 30.624
Escogemos el menor = 30.624Tn
Hallando FS (deslizamiento)

FS=
∑ FH resistente ( ton ) = 30.624 =1.71>1.5 ,CUMPLE
∑ FH actuante ( ton) 17.806
 Verificación por excentricidad

∑ Moresultante=∑ Mestabilizante−∑ Mdesestabilizante=123.406

∑ Fvertical =59.340
∑ Mresultante =Eo=2.080
∑ Fvertical
b X 1+ X 2+ X 3
Ec= −Eo= −Eo=0.145
2 2
b
Verificando que Ec< ,=0.145<0.741CUMPLE
6
 Cálculo de presiones por punta y talón

∑ Fvertical ∗ 1+ 6 e tn
qmáx=
B ( B )=15.948
m

∑ Fvertical ∗ 1− 6 e tn
qmin=
B ( B )=10.721
m

 Verificación por resistencia

 ∅ 2 π∗tan ∅
N q =tan ⁡(45+ ) ∗e =16.443
2
N c =( N q−1 ) cot ∅=27.860
N γ =2 ( N q +1 )∗tan ∅
Ψ 2
(
F γi = 1−
∅ )=0.18
F γd=1
1−F qd
F cd =Fqd − =1.106
N c∗tan ∅

Ψ 2
(
F ci=F qi = 1−
90) =0.663

D
( )
F qd=1+2 tan ∅ ( 1−sen ∅ )2∗
B'
=1.099

Ψ =tan −1 ( Pa∗cosα
∑V )
=16.703

H4
Ahora calculamos q(Tn/m2) = =2.97
γ2
qu=

qu(tn/m2) =117.342

qu
FS= =7.358
qmax

Ahora se resolverá el problema considerando que se nivelará el terreno que corresponde al

talud superior
 Los datos para resolverlos se mostrarán a
continuación

betha 0
X1(m) 0.8
X2(m) 1.05
X3(m) 2.6
X4(m) 0.55
H2(m) 6.7
 H3(m) 0.7
H4(m) 1.5
Yc(ton/m3) 2.4
Y1(g/cm3) 1.85
Y2(g/cm3) 1.98
∅' 2
1

31
29
c'1(Kpa) 0
c'2(Kpa) 33
δ 17

Ca 0.85c'
 Las fórmulas para hallar el Ka y Kp son las siguientes:

Ka 0.288
fi 31
w 0
behta 0
H1 0.000

Kp 4.477
fi 29
tan-1(x2-
w x4/h2) 4.268
betha 0

Ahora
se hallará
Ea y Ep en
kg/cm2
Ea = γ∗z∗ka =0.394
Ep=γ∗z∗kp=1.329
Ea resultante(ton/m2) =Ea∗0.5∗10∗( H 2+ H 3 ) =14.588
Ep resultante (ton/m2) =Ep*0.5*H4*10=9.972
 La imagen que se muestra a continuación es el bosquejo de los diagramas que se presentan
debido a la sobrecarga de 20KN/m2

Luego se procederá a hacer un cuadro de cálculos para hallar sumatoria de fuerzas, verticales,
brazos y momentos estabilizantes o desestabilizantes
F F Mom. Mom.Desestabilizant
N° nombre vertical horizontal Brazo Estabilizante e
14.58806
1 Ea horiz   4 2.467   35.98875389
2 Ea vertical   0      
9.972517
3 E pasivo   5 0.5 4.98625875  
4 peso 1 8.844   1.575 13.9293  
5 peso 2 7.476   2.225 16.6341  
6 peso 3 4.02   1.133 4.55466  
7 peso 4 32.227   3.583 115.469341  
8 sobrecarga   14.8 3.7   54.76

∑ 52.567     155.5736598 90.74875389

Ahora se procederá a hacer la verificación por volteo:

FS= Mestab./M.Desestab= 1.714 NO CUMPLE PORQUE NO ES MAYOR A 2
Ahora se procederá a hacer la verificación por deslizamiento:
 Primero hallaremos sumatoria FH actuantes(ton) y sumatoria FH resistente(ton)
∑FH actuante(ton)=Ea horiz+F horiz de la sobrecarga =29.388
∑FH resistente se eligirá el menor valor reemplazado en las fórmulas

∑Fv.tg∅ +cb=43.82
Se escoge el menor
∑Fv.tgδ +ca.b=28.553
Ahora calculamos el FS deslizamiento
mediante:
∑FH actuante(ton)/ ∑FH resistente(ton)=0.971 POR LO TANTO NO CUMPLE YA QUE NO ES
MAYOR A 1.5
Ahora se procederá a hacer la verificación por excentricidad:
Se mostrará previamente ∑Mo resultante, ∑Fv , eo y ec
 ∑Mo resultante =∑M.estab-∑M.desestab=64.824
 ∑Fvertical=52.567 ec no es menor a b/6
por lo tanto no cumple
 eo=∑Fvertical/∑Mo resultante=1.233
la verificación por
b excentricidad
 e c = −e 0=0.991
2
b
 =0.741
6

Como conclusión se podría decir que si se somete a una sobrecarga 20Kn/m2 no

cumplirá con ninguna de las verificaciones mostradas

Pregunta N° 8

Para la evaluación de la estabilidad de un talud se le pide lo siguiente:

a) Mostrar paso a paso la deducción de la ecuación general asociada al factor de
seguridad aplicando el método de Fellenius.
b) Tomando en cuenta la figura mostrada, determinar el factor de seguridad estático
 y pseudoestático para la superficie de falla mostrada aplicando el método de
Fellenius. Se deberán realizar los cálculos de manera tabulada utilizando como
mínimo 20 dovelas.
c) Utilizando programas de cálculo evaluar la estabilidad del talud mostrado
mediante la técnica de equilibrio límite, llevar a cabo la evaluación utilizando los
métodos de Spencer, Fellenius, Bishop y Morgestern Price. Deberá evaluar la
misma superficie de falla del paso anterior. Comentar los resultados, asimismo
comparar las diferencias que encontró con el cálculo manual realizado en el paso
anterior.
d) Utilizando programas de cálculo identificar el mínimo factor de seguridad
mediante la técnica de equilibrio límite y también mediante la técnica de
elementos finitos. Comentar los resultados, respecto a los cálculos anteriores.
e) Deberá realizar los cálculos en los pasos b), c) y d) para condiciones estáticas y
para condiciones pseudoestáticas.

Deberá indicar y explicar claramente el procedimiento seguido en cada paso

y comentar los resultados. NOTA: La línea de color verde representa el nivel
freático

Coordenadas geográficas - Grupo B: -14.50° LAT; -75.50° LONG

8ª. El método de fellenius es un proceso más rápido y simplificado que los
demás, pues este considera de igual magnitud las fuerzas de tangenciales y
normales, ya que es complejo determinarlas mediante tediosos procesos. En
otras palabras, se asume cero la sumatoria de las fuerzas mencionadas, lo que
genera que estas sean eliminadas.
Teniendo lo anteriormente mencionado en cuenta, se determina la ecuación de
equilibrio:
Nr=Wn∗cos α
Sin embargo, se debe considerar la presencia de agua, lo cual genera:
Nr=( Wn+Ws )∗cos α
Además, el agua ejerce empuje en la dovela, lo cual es transmitido a ella como
una fuerza negativa y se relaciona también con el volumen que ocupa. De esta
forma, se llega a la siguiente expresión:
Nr=( Wn+Ws )∗cos ∝−(9.81∗b∗zw) , donde b es el acho de la dovela y Zw la
altura que ocupa el agua en ella.
Teniendo en cuenta la fórmula de resistencia del suelo, el esfuerzo cortante
sería el siguiente:
τ =( c +σ '∗tan ∅ )∗ΔL
Consideramos también, el esfuerzo normal donde incluimos a Nr definido
anteriormente:
Nr
σ' =
ΔL
Combinando las dos últimas ecuaciones planteadas llegamos a lo siguiente:
τ =C∗ΔL+ Nr∗tan ∅
A continuación, realizamos el cálculo de la dovela:
F=c∗b+ tan ∅∗( w∗cos α −u∗b)
Descomponemos el peso, y de esta manera hallamos las fuerzas actuantes:
F=W ∗sin α
Finalmente, llegamos a la conclusión de que el método de fellenius satisface el
equilibrio de momentos, además de la relación entre las fuerzas actuantes y
dividiéndolas, obtenemos el FS.
F . S .=R∗∑ ¿ ¿ ¿

8b. Desarrollo detallado del excel

Para aplicar el método de fellenius, debemos en primer lugar graficar (según
coordenadas y distancias), el talud a analizar en autocad. Acto seguido, ubicar
el punto C (especificado en el enunciado del ejercicio) y con él, graficar
mediante una circunferencia la superficie de falla y finalmente dividirla entre el
número de dovelas que se requieran (en nuestro caso 20 de 1.2m de ancho
cada una). Finalmente, trazamos el nivel freático que dividirá el talud y la
superficie de falla en dos secciones.
Con el gráfico ya hecho procedemos a hallar, mediante el programa autocad
también, las áreas, tanto saturadas como naturales, de cada dovela creada (se
combinan las rectas y curvas que componen cada dovela para poder hallar su
área mediante el comando que lleva su nombre). Luego, se trazan rectas
horizontales desde el centro de cada curva de las dovelas, para después poder
trazar una tangente a las mismas y así determinar el ángulo entre ambas

rectas.

*Observación: las tangentes del centro de la curva y las rectas horizontales de

cada dovela, no están completamente pegadas a ellas, pues es un error del
autocad separarlas automáticamente. Sin embargo, está confirmado que el
ángulo que se forma entre ellas es el real, pues se ha comparado con otros.
A continuación, se procede a calcular mediante las áreas de las dovelas, los
volúmenes de las mismas, multiplicándolas por su espesor (1m) y
posteriormente por su densidad (dato del problema), ya sea saturada o natural,
para poder hallar su peso en Kn.
Hasta el momento se tienen los siguientes datos: Radio de superficie de falla,
Ángulo alpha, espesor de dovela (B) y pesos de dovelas (W). Sin embargo, el
ejercicio nos brinda datos extras, entre ellos la cohesión del suelo en Kn/m2 y
el ángulo de fricción.
Acto seguido, determinamos la presión de poros (u) en cada porción de dovela
incluida dentro del nivel freático, multiplicando la gravedad (9.81m/s2) por los
volúmenes saturados de las mismas.
Por último, hallamos el momento resistente, el cual será el numerador de la
división posteriormente realizada para determinar el facto de seguridad
estático. Reemplazamos datos en la siguiente fórmula: c*b+tan(Φ)*(W*cos(α)-

u*b/cos(alpha)+cb/cos(alpha)). Dicha fórmula se aplica en cada dovela,

considerando el peso total de la misma como la suma del peso saturado más el
natural, para luego sumar los 20 resultados y multiplicarlos por el radio de la
superficie de falla.
Para determinar el denominador de la fórmula del factor de seguridad estático,
el cual lleva el nombre de momento actuante, se realiza la siguiente operación
en cada dovela y cada resultado se suma para luego multiplicarlo nuevamente
por el radio de la superficie de falla: W*sen(α).
Para hallar el F.S. pseudoestático, se divide el momento resistente hallado
anteriormente, entre la suma del momento actuante y el mismo multiplicado por
la fuerza sísmica actuante: W*sen(α) + W*sen(α)*(PGA/2). En este caso,
nuestro PGA será 0.48.

A continuación, se presentarán los cálculos detallados de la primera dovela,

pues estos serán los mismos a lo largo de las 20 analizadas. Sin embargo, sí
estarán especificados todos los datos hallados en cada una de ellas mediante
tablas.
 Áreas de dovelas extraídas de dibujo en autocad
N° Área sumergida Área total
Dovela (m2) Área NO sum. (m2) (m2)
1 0.1032 0.3455 0.4487
2 0.2701 1.0366 1.3067
3 0.3766 1.7277 2.1043
4 0.4211 2.4188 2.8399
5 0.4016 3.1099 3.5115
6 0.3157 3.801 4.1167
7 0.1602 4.4921 4.6523
8 0.0073 5.1075 5.1148
9 0 5.4996 5.4996
10 0 5.8016 5.8016
11 0 6.0139 6.0139
12 0 6.1284 6.1284
13 0 6.1348 6.1348
14 0 6.0197 6.0197
15 0 5.7659 5.7659
16 0 5.3498 5.3498
17 0 4.7382 4.7382
 18 0 3.8816 3.8816
19 0 2.7 2.7
20 0 1.0474 1.0478
Hallamos el área total sumando el área sumergida más el área natural.

 Ángulos alpha hallados también en autocad

Alpha (α)
8.775
11.079
13.402
15.747
18.119
20.524
22.986
25.457
27.998
30.601
33.276
36.037
38.897
41.879
45.008
48.319
51.862
55.714
59.997
64.946

 Alturas de dovelas según autocad

altura(z)
0.376
1.091
1.755
2.368
2.928
3.432
3.878
4.263
4.584
4.835
5.012
5.107
 5.112
5.015
4.802
4.454
3.941
3.223
2.23
0.836

Altura de cada dovela desde el centro de la superficie de falla hasta el talud.

 El ancho de la dovela es igual a la distancia entre los puntos de

intersección entre la superficie de falla y el talud, dividida entre la
cantidad de dovelas analizadas (20).
24 m
=1.2m
20
 Se sabe que el espesor de cada dovela es 1m, por lo que se multiplica el
área de la misma por este y así hallamos su volumen.

Vtotal=0.4487 m∗1. m=0.4487 m2

Sin embargo, se requiere hallar el volumen sumergido y natural multiplicando,
de igual manera, el área respectiva por el espesor de la dovela.

Vsumergido=0.1032 m∗1 m=0.1032 m 3

Vnatural=0.3455 m∗1 m=0.3455 m3
 Hallando el peso específico sumergido
19.8 Kn 9.81 Kn 9.99 Kn
γsumergido=γsaturado−γh 2 o= − =
m3 m3 m3
  Multiplicando el peso específico (sumergido o natural) por el volumen
(sumergido o natural), hallamos los pesos naturales y sumergidos, los
cuales se sumarán para hallar el peso total.

m3∗9.99 kn
Wsumergido=0.12384 =1.2136 Kn
m3
0.4146∗18.4 kn
Wnatural= =7.62864 Kn
m3
Wtotal=1.2136+ 7.62864=8.8658 Kn
 La presión de poros se determina multiplicando el área sumergida por el
peso específico del agua

m 3∗9.81 Kn
u=0.1032 =1.012392
m3
 El enunciado del problema nos brinda la cohesión y el ángulo de fricción.
c=14.5 KPa
∅=33 °
 El radio del círculo de falla se halla en autocad, definiendo la distancia
del pto. C hacia alguno de los extremos del talud.
R=30.3031m

*Observación: el radio izquierdo mide 30.3031m y el de la derecha 30.3033m.

Esto se debe a que realmente la superficie de falla no llega a la esquina
superior derecha del talud, por lo que la distancia desde ahí al punto C, será
mayor. Por cuestiones de exactitud, decidimos tomar como radio de superficie
a 30.3031m.
  Momento resistente

u∗b
( (
Mresistente=R∗∑ c∗b+ tan ∅∗ Wtotal∗cos ∝−
cos ∝) )
+c∗b/cos ∝ ¿ ¿

1.012 Kn∗1.2 m
Mresistente=30.3031m∗
m2 (
14.5 Kn
(
∗1.2 m+ tan 33∗ 7.388 Kn∗cos 8.775−
m2
cos 8.775
+
c∗b
))
cos 8.775
=3

 Momento actuante
Mactuante=R∗∑ Wtotal∗sin ∝

Mactuante=30.3031 m∗7.388 Kn∗sin 8.775=30.3031∗1.1271

 ∑ Momento resistente=
u∗b c∗b
( (
R∗∑ c∗b+ tan ∅∗ Wtotal∗cos ∝− )
+
cos ∝ cos ∝ )
¿ ¿=30.3031m∗1569.92

 ∑ Momento actuante=¿ R∗∑ Wtotal∗sin ∝=30.3031 m∗1024.7855¿

 Factor de seguridad estático
Mresistente
FS= =1.53194
Mactuante
 Factor de seguridad pseudoestático
Mresistente 1569.92
FS= = =1.2354
R∗∑ Wtotal∗sin∝+ R∗∑ W ∗sen (α )∗(PGA / 2) 1024.785503+245.9485

*Comentarios de resultados al final de la pregunta 8c.

8c.
Datos del suelo en análisis:
 C (Kpa) ∅ kn kn
ρnat ( ) ρsat ( )
m3 m3
14.5 33° 18.4 19.8

Talud dibujado en slide con características del suelo detalladas

 Para un search grid de coordenadas –(35;95), se obtuvieron los

siguientes factores de seguridad según cada método (condiciones
estáticas, sin NF):
Se observa que los factores de seguridad mínimos pertenecen al color naranja,
lo cual significa que cualquier punto diferente se alejará de ese tono, es decir
tendrá factores de seguridad mayores.
- Se analizaron 22500 superficies a lo largo del talud mostrado y a
continuación se muestra el gráfico hallado
Fuerza sísmica:
Hallando el PGA en mapa de distribución de isoaceleraciones para un 10% de
excedencia de 50 años, este es 0.48.
 Valor de PGA

 Ingresamos el PGA en slide como 0.24=0.48/2. Además, 0.24<0.25

(cumple). Al determinar el FS obtenemos lo siguiente. (Condición
pseudoestática, sin NF)
 FS<1, Talud no
construible,
falla asegurada

Al ingresar el PGA en slide como 0.16=0.48/3. Además, 0.16<0.25 (cumple). Al

determinar el FS obtenemos lo siguiente.

Cuando FS=1 – fuerza exactamente igual a resistencia del suelo, suelo a punto
de fallar.

Nivel freático:
NOTA: Sabemos que en clase se dijo que al ingresar solo una recta como nivel
freático sería inconsistente, pues por un lado abría agua y por otro no. Sin
embargo, nosotros decidimos copiar el NF de la figura del enunciado del
trabajo tal cual se mostraba ahí para poder obtener resultados más cercanos a
la realidad.
Se observa que la presencia de agua (nivel freático) en el talud, disminuye el
FS, esto se debe a que el agua afecta principalmente a la presión de poros y al
esfuerzo cortante. Este, al estar en función de la cohesión y del ángulo de
fricción, se incorpora el efecto reductor del agua en ellos también, lo cual
significa que se utilizarán los valores totales de dichos factores (Ctotal y ∅total),
que son menores a los efectivos.
Se observa también que el impacto del agua sobre la superficie de falla es
pequeño. Si abarcase mayor área de la superficie de falla, el FS disminuiría
aún más.

Área de impacto
de agua en
superficie de falla

 Condición pseudoestática (con NF)

En el caso de utilizar el PGA como 0.16, el FS que genera el NF en cada uno
de los métodos es el siguiente:
En el caso de utilizar el PGA como 0.48, el FS que genera el NF en cada uno
de los métodos es el siguiente:
Observamos que, en todos los casos, el FS<1, lo cual indica que el talud ya
falló y por tal razón no es posible implementarlo.

Análisis extras:
 FS estático con NF
 FS pseudoestático sin NF
  FS en punto c
Al utilizar la opción de hallar un FS en cualquier parte del talud, ya sea dentro o
fuera de él, el programa solo permite hallar factores de seguridad en
coordenadas específicas del mismo, por lo que no hemos podido hallar el factor
de seguridad exacto. Determinamos un FS aproximado, este se encuentra
entre los puntos (-4.400;31.367), (-4.400;32.133), (-3.800;31.367) y (-
3.800;32.133). Dichos puntos conforman un rectángulo de base 0.6m y altura
0.766m. El punto C tiene las coordenadas (-4.022;31.535).
FS C condición estática con NF

Condición estática
Coordenada FS
Método s FS mínimo
(- 1.53
Fellenius 3.8;32.133) 7 1.335
(- 1.57
  4.4;32.133) 1  
(- 1.57
  4.4;31.367) 5  
(-
  3.8;31.367) 1.54  
(- 1.58
Bishop simpl. 3.8;32.133) 7 1.42
(- 1.61
  4.4;32.133) 6  
(- 1.62
  4.4;31.367) 1  
(- 1.59
  3.8;31.367) 2  
Morgenstern- (- 1.58
price 3.8;32.133) 2 1.416
(-
  4.4;32.133) 1.61  
 (- 1.61
  4.4;31.367) 5  
(- 1.58
  3.8;31.367) 7  
(- 1.58
Spencer 3.8;32.133) 3 1.418
(- 1.60
  4.4;32.133) 9  
(- 1.61
  4.4;31.367) 4  
(- 1.58
  3.8;31.367) 7  

FS C condición pseudoestática con NF

Condición pseudoestática
Coordenada FS
Método s FS mínimo
1.05
Fellenius (-3.8;32.133) 3 0.88
1.07
  (-4.4;32.133) 9  
1.08
  (-4.4;31.367) 2  
1.05
  (-3.8;31.367) 5  
1.09
Bishop simpl. (-3.8;32.133) 6 0.961
1.11
  (-4.4;32.133) 8  
1.12
  (-4.4;31.367) 2  
1.09
  (-3.8;31.367) 9  
1.09
Morgenstern-price (-3.8;32.133) 4 0.968
  (-4.4;32.133) 1.12  
1.12
  (-4.4;31.367) 4  
1.09
  (-3.8;31.367) 7  
1.09
Spencer (-3.8;32.133) 7 0.971
1.11
  (-4.4;32.133) 8  
1.12
  (-4.4;31.367) 2  
 1.10
  (-3.8;31.367) 1  

Resumen detallado de resultados:

Condición estática Método  

Con Nivel Freático
(1) Fellenius 1.335
Bishop
  simplificado 1.42
  Spencer 1.418
Morgenstern-
  Price 1.416
Sin Nivel Freático (2) Fellenius 1.388
Bishop
  simplificado 1.461
  Spencer 1.455
Morgenstern-
  Price 1.456
Condición  PGA=0.4
pseudoestática Método 8
Con Nivel Freático
(3) Fellenius 0.88
Bishop
  simplificado 0.962
  Spencer 0.971
Morgenstern-
    Price 0.968
Sin Nivel Freático (4) Fellenius 0.928
Bishop
  simplificado 0.987
  Spencer 0.987
Morgenstern-
    Price 0.986

Comentarios:
En caso
1. F<MP<S<B
2. F<S<MP<B
3. F<B<MP<S
4. F<B=S<MP
En todos los casos se cumple que el FS en el método de fellenius es el menor,
lo cual indica que el talud se encuentra en las condiciones más desfavorables y
por lo tanto, más próximo a fallar. El método de fellenius es el más simple, ya
que no considera las fuerzas normales y de fricción (fuerzas entre dovelas, se
asume que suman cero, ya que son muy parecidas y contrarias) de cada
dovela de la superficie de falla del talud, por lo que sus resultados serán
menores y menos precisos. Sin embargo, tiene mayor precisión de poros.

Por otro lado, el FS mínimo hallado mediante el método de Bishop en la

condición estática es el mayor, pues este método utiliza todas las fuerzas
generadas en la dovela, lo cual reduce el porcentaje de error para el resultado
final y por lo tanto, aumenta el valor del FS. No obstante, en la condición
pseudoestática el FS mínimo es segundo menor. Esto indica que el método
cumple con una muy buena precisión al momento de calcular FS’s mínimos en

condiciones estáticas, más no en condiciones pseudoestáticas, pues los FS se

reducen bastante.
El método de Spencer considera también todas las fuerzas de la dovela. No
obstante, direcciona las fuerzas de conexión a lo largo de las superficies de
división de cada dovela con un ángulo tetha de manera que sean paralelas a la
superficie de la dovela en análisis (utiliza misma inclinación para cada dovela),
satisfaciendo la estática en su totalidad. Por esta razón, se cumple que el FS
mínimo es el mayor en la condición pseudoestática y el segundo mayor en
condición estática. Se observa que dicho método logra una mayor precisión de
resultados en condiciones pseudoestáticas, a diferencia del método de Bishop.
Finalmente, el método de morgenstern-price tiene la particularidad de dividir la
masa en franjas infinitesimales, a quienes se les aplica las ecuaciones de
equilibrio, lo cual genera también precisión. Además, se cumple que es
equivalente con el método de Spencer (sin ser más preciso que él) y por eso
sus valores de FS también son próximos y el valor de su FS mínimo en la
condición pseudoestática, es mayor al de la condición estática, igual que en el
método de Spencer.
Se observa también, una clara disminución del FS al agregarle fuerzas
sísmicas al talud, lo cual es lógico porque evidentemente su resistencia no se
verá afectada negativamente después de un sismo (dependiendo de su
magnitud e intensidad).
Además, se observa que en todos los casos pseudoestáticos, el FS del talud se
encuentra por debajo de la unidad. Esto quiere decir, como ya se mencionó en
reiteradas oportunidades, que el talud fallará y no soportará las cargas que se
le apliquen. Por lo tanto, el talud en análisis no podrá ser implementado en la
obra y se tendrá que rediseñar o de ser el caso, reubicar.
Lo mismo sucede con la presencia de agua, dentro de la superficie de falla del
talud. Esta también genera una disminución del FS (efecto explicado a detalle
anteriormente).
Según norma, existen distintos valores del FS mínimo según cada caso, estos
son los siguientes:

Asumimos un talud permanente

En la condición estática, todos los FS mínimos de los 4 métodos utilizados
pertenecen al rango de la norma AASHTO LRFD. Sin embargo, el valor de las
normas NAVFAC y CE.020, no coincide con ningún FS hallado, pues el más
cercano es el de Bishop simplificado, el cual es 1.42.
En la condición pseudoestática, ninguno de los valores hallados pertenece al
rango presentado por la norma NAVFAC, pues todos son menores a 1. Esto,
como ya se mencionó anteriormente, significa que el talud ha fallado y por lo
tanto no podía ser utilizado.
Comentarios acerca de comparación de resultados en slide y a mano:
 Método de Fellenius Slide (aprox. Según Excel
datos en coordenadas
cercanas a c)
Condición estática 1.55575 1.53195
Condición 1.06725 1.23544
pseudoestática

Se observa una diferencia de 0.0238 entre el FS en la condición estática

hallado en slide y en Excel. Consideramos que esta diferencia se debe a los
errores automáticos que autocad posee, pues esta herramienta es muy útil
pero también puede jugarnos malas pasadas. Un claro ejemplo es el caso de
las tangentes y rectas horizontales trazadas en el centro de la superficie de
falla hallado en cada dovela (error especificado anteriormente). Es posible que
sí exista algún error milésimo en las alturas, pues al darle el máximo zoom a la
intersección de la altura trazada con la recta tangente, se observa que la altura
sobrepasa ligeramente la superficie de falla. Sin embargo, el programa no nos
permite medir la altura real por errores automáticos del mismo.
Finalmente, llegamos a la conclusión de que otro factor que imposibilita la
igualdad exacta de ambos factores de seguridad, son los decimales no
utilizados al hallar las áreas de las dovelas y sus alturas en autocad. Si bien es
cierto, este programa si cuenta con la opción de aumentar decimales, pero en
muchos casos el propio programa redondea automáticamente los resultados
(este es el ejemplo de las áreas), sin tener opción a aumentar más decimales.
Por otro lado, al hallar distancias de alturas, sí podemos aumentar decimales,
pero en los casos donde estos son infinitos, es por razones obvias, imposible
copiarlos tal cual AutoCAD los indica. Buscamos en internet y confirmamos que
no existe ninguna opción que brinde AutoCAD de copiar valores infinitos.

Diferencia estática= 0.0238

Diferencia pseudoestática= 0.16819

8d.
En primer lugar, definimos las dimensiones sobre las cuales se hará todo el
análisis, donde pueda entrar el talud completo. Acto seguido, dibujamos el talud
con sus coordenadas respectivas.

A continuación, definimos los parámetros del suelo del talud (cohesión y Angulo
de fricción indicados en el enunciado del ejercicio). Elegimos un valor de 30
000kn/m2 para el módulo de elasticidad, de 0.24 para el módulo de poisson,
extrayendo de internet el dato para arena limosa y de 1° para el ángulo de
reducción.
Definimos el peso específico del agua con un valor de 9.81kn/m3, colocándole
también las restricciones al talud
Cabe resaltar, que debimos regeneramos las mallas para llegar a una mayor
exactitud.
A continuación, ingresamos el nivel freático con sus respectivas coordenadas.
Finalmente, obtenemos el resultado final del desplazamiento del talud al
momento de fallar:

Presión de poros
Resultados finales
Factor de seguridad mínimo=Msf= 1.322
El FS mínimo hallado difiere del FS hallado mediante el programa slide con el
método de Morgenstern-price, pues este tiene el valor de 1.415.
Consideramos que la diferencia entre ambos se debe a la carencia de los
datos, como ya se mencionó anteriormente, del módulo de elasticidad y del
módulo de poisson, los cuales tuvimos que asumir, sustentando sus valores en
la web.