lica De Chile Pontificia Universidad Cato Escuela de Ingenier a cnica Departamento de Ingenier a Estructural y Geote nica de So lidos ICE
CE 2313 - Meca
Tarea 2
Ayudante : Juan Santib an ez Millacura Profesor : Mat as Hube Ginestar
Problema 1.
a)
Debido a la simetr a de la secci on transversal, conviene calcular los momentos de inercia seg un el sistema de coordenadas referido en el enunciado. Luego, se tiene:
Iz = =
y 2 dA
10 10 10
z 2 +10 z 2 10
y 2 dydz
y3 z 2 +10 |z dz 2 10 10 3 1 10 = (15z 2 + 2000) dz 3 10 = 16666,67 cm4
Iy = =
z 2 dA
10 10 10
z 2 +10 z 2 10
z 2 dydz
=
10
z 2 20 dz
= 13333,33 cm4
Iyz = =
yz dA
10 10 10
z 2 +10
yz dydz
z 2 10
y2 z 2 +10 = z |z dz 2 2 10 10
10
=
10
10z 2 dz
= 6666,67 cm4
b)
El tensor de inercia queda entonces representado por I= 16666,67 6666,67 6666,67 13333,33 cm4
Calculando los valores propios de este tensor se tiene que: 1 = 21871,85 cm4 1 = 37,98 1 . 2 = 8128,15 cm4 2 = 127,98
Utilizando los vectores propios asociados se hallan los siguientes angulos:
En la siguiente gura se aprecia la conguraci on original de la secci on de la viga junto con la rotada seg un
Figura 1: Secci on de la viga Por otro lado, dada la vinculaci on de la viga y la carga presenta en ella, se tiene la siguiente situaci on:
Figura 2: Reacciones de v nculo 2
Adem as, considerando la convenci on habitual de esfuerzos internos (tracci on positiva de la bra inferior) se tiene que V y ( x) = q ( L x) Mz (x) = qL x qL2 qx2 2 2
q = ( L x) 2 2 V z ( x) = M y ( x) = 0
Al reemplazar con los par ametros brindados se obtiene que Vy (0) = 384 kgf = 3797,04 N Mz (0) = 76800 kgf cm = 7534,08 N m
c)
Ahora para calcular las tensiones de tracci on y compresi on utilizamos la siguiente expresi on: x = (Iy Mz Iyz My )y (Iyz Mz + Iz My )z 2 Iy Iz Iyz
Entonces, evaluando la expresi on para z = 10, y = 15 se obtiene la tensi on m axima de tracci on.
T max = 57,60
kg = cm2
5,65 M P a
Tambi en, evaluando la expresi on para z = 10, y = 15 se obtiene la tensi on m axima de compresi on.
C max = 57,60
kg = cm2
5,65 M P a
Para encontrar el eje neutro utilizamos la siguiente expres on: tan = De donde se obtiene que = 26,56 Iyz Mz + Iz My Iy Mz + Iyz My
d)
Para calcular la tensi on de corte promedio, xy es necesario evaluar la siguiente expresi on: xy = V Q(y ) Iz b
en donde V (y ) es la magnitud del esfuerzo de corte, Q(y ) corresponde al momento est atico, Iz momento de inercia y b el ancho de la secci on de estudio. Es conocido que el esfuerzo de corte V en la secc on x = 0, que tiene una magnitud de 3797.04 N y el ancho de la secci on es 20 cm. Se deben calcular los otros par ametros. Para conocer el momento est atico consideremos cm, pero para ser consistentes con nuestro el tr angulo superior. Su centroide (medido desde su base) est a a 10 3 sistema de coordenadas debemos empezar a medir desde el origen, por lo que y = 5 + 10 cm. El a rea de dicho 3 tri angulo tiene magnitud 100 cm2 , entonces Q(y ) = A y = 833,33 cm3 = 8,33 104 m3 Finalmente, como todos los componentes son conocidos, y por la simetr a de la secci on, el valor de la tensi on de corte promedio es xy = 0,094 M P a
e)
Figura 3: x , secci on x = 0
f)
Figura 4: x , plano z = 0
C odigos
q4[:,i]=[Link](x4[:,i]) #Defino los graficos
Tarea [Link]
""" Created on Sun Nov 10 12:15:35 2013 @author: juantik """ import numpy as np import [Link] as plt
p5=[Link](q1[0],q1[1],color=red,label="Rotada") p6=[Link](q2[0],q2[1],color=red) p7=[Link](q3[0],q3[1],color=red) p8=[Link](q4[0],q4[1],color=red) #dibujo todo [Link]([Link]) [Link](bbox_to_anchor=(1.05, 1.05), borderaxespad=0.) [Link](Seccion Original y Rotada) [Link](equal) [Link](on)
I=[Link]([[16666.67,6666.67],[6666.67,13333.33]]) l1=[Link](I)[0][0] #Inercia principal 1 l2=[Link](I)[0][1] #Inercia principal 2 v1=[Link](I)[1][0] #Vector propio 1 a1=np.arctan2(v1[1],v1[0]) #alfa 1 mz=-76800. #momento en torno a z def sigma_x(y,z): #funcion que entrega sigma x a=-I[1,1]*mz b=I[0,1]*mz c=I[1,1]*I[0,0]-I[0,1]**2 d=(a*y+b*z)/c return d #eje neutro psi=np.arctan2((I[0,1]*mz),(I[1,1]*mz))+[Link] #VALORES EXTREMOS (MPA) s1=sigma_x(15.,10.)/10.19716 s2=sigma_x(-5.,10.)/10.19716 s3=sigma_x(5.,-10.)/10.19716 s4=sigma_x(-15.,-10.)/10.19716
v2=[Link](I)[1][1] #vector propio 2 a2=np.arctan2(v2[1],v2[0]) #alfa 2 a1g=[Link](a1+[Link]) #valor en grados de alfa 1 a2g=[Link](a2+[Link]) #valor en grados de alfa 2 print a1g,a2g a3=a2+[Link] #angulo de rotacion (arbitrario) #Grafico de la seccion original z1=[Link](-10.,10.,100) y1=-0.5*z1-10 x1=[Link]([z1,y1]) z2=[Link](10.,10.,100) y2=[Link](-15.,5.,100) x2=[Link]([z2,y2]) z3=[Link](-10.,10.,100) y3=-0.5*z3+10 x3=[Link]([z3,y3]) z4=[Link](-10.,-10.,100) y4=[Link](15.,-5.,100) x4=[Link]([z4,y4])
Tarea 2 - [Link]
# -*- coding: utf-8 -*""" Created on Sun Nov 10 11:56:59 2013 @author: juantik """ import numpy as np import [Link] as plt from mpl_toolkits.mplot3d import Axes3D import Plot2DTriangleGeometry Iy=16666.67 Iz=13333.33 Iyz=6666.67 Mzx0=-76800
#Grafico [Link]() p1=[Link](z1,y1,r--, p2=[Link](z2,y2,r--, p3=[Link](z3,y3,r--, p4=[Link](z4,y4,r--,
xyz=[Link]([Link]) IEN=[Link]([Link]) #Matriz de Rotacion suckads = np.int_(IEN) R=[Link]([[[Link](a3),-[Link](a3)],[[Link](a3), [Link](a3)]]) vals=[Link]([len(xyz),2]) sigma=[Link]([len(xyz),1]) #Defino los vectores rotados a dibujar q1=[Link]((2,100)) for i in range(100): for i in [Link](len(xyz)): q1[:,i]=[Link](x1[:,i]) vals[i,0]= ((-Iy*Mzx0)*(xyz[i,1]-15.)+ (Iyz*Mzx0)*(15.-xyz[i,0]))/(Iy*Iz-Iyz*Iyz) q2=[Link]((2,100)) for i in range(100): [Link]() q2[:,i]=[Link](x2[:,i]) [Link](equal) [Link](xyz[:,0],xyz[:,1],np.int_(IEN)-1,vals[:,0], q3=[Link]((2,100)) levels=[Link](-350,85,0.5)) for i in range(100): Plot2DTriangleGeometry.Plot2DTriangleGeometry(xyz,IEN) q3[:,i]=[Link](x3[:,i]) #Tricontourf(COORD_X, COORD_Y,RELACION_XY, VALORES_COORD,LEVELS=RANGO_DEL_CONTOUR) q4=[Link]((2,100)) [Link](equal) for i in range(100):
color=blue,label="Original") color=blue) color=blue) color=blue)
[Link]() [Link]()
xyz2=[Link]([Link]) IEN2=[Link]([Link]) vals2=[Link]([len(xyz2),1]) for i in [Link](len(xyz2)): #en este caso, z=0 vals2[i,0]=((-Iy*(0.96*400*xyz2[i,0]0.5*0.96*xyz2[i,0]**2-Mzx0))*xyz2[i,1])/(Iy*Iz-Iyz**2) [Link]() [Link](xyz2[:,0],xyz2[:,1],np.int_(IEN2)1,vals2[:,0],levels=[Link](-320,80,0.5)) Plot2DTriangleGeometry.Plot2DTriangleGeometry(xyz2,IEN2) [Link](equal) [Link]() [Link]()