Ejercicio 3 Ef

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 6

x1 ≔ 4 m x2 ≔ 5 m x3 ≔ 0.2 m x4 ≔ 0.

2 m

y1 ≔ 3 m

m
M ≔ 1800 kg δr ≔ 2.1 g ≔ 9.81 ―
s2
E1 ≔ 200 GPa E2 ≔ 70 GPa

A1 ≔ 800 mm 2 A2 ≔ 1000 mm 2

solución:

calculo de fuerza T

T + ml ⋅ g - M ⋅ g = 0

T = M ⋅ g - δh20 ⋅ Vh20 ⋅ g
M
T = M ⋅ g - δh20 ⋅ ―― ⋅g
δobj

⎛ 1⎞
T = ⎜1 - ―⎟ ⋅ M ⋅ g
⎝ δr ⎠

⎛ 1⎞
T ≔ ⎜1 - ―⎟ ⋅ M ⋅ g = 9249.429 N
⎝ δr ⎠
calculo de longitud de cables

tomando en cuenta que es una elipse

x2 y2
― + ― =1
a2 b2

x1 + x2 + x3 + x4
para a ≔ ―――――― = 4.7 m b ≔ y1 = 3 m
2

despejando y
‾‾‾‾‾‾‾‾‾‾
⎛ x2 ⎞ 2
y ( x) ≔ ⎜1 - ― ⎟⋅b
⎝ a2 ⎠

y ((-2.5 m)) = 2.54 m y ((1.5 m)) = 2.843 m

Condiciones de frontera:

asumiendo que el sistema tendra una relacion proporcional de desplazamientos

Q1 = 0 Q3 = 0
Haciendo relacion de triangulos:

Q2 Q5 Q5 -2
―= ― Q2 - 2 ⋅ ― = 0 B1 ≔ 1 B2 ≔ ―― = -0.5
2 4 4 4

Q4 Q5 3 ⋅ Q5 -3
―= ― Q3 - ――= 0 B1 ≔ 1 B2 ≔ ―― = -0.6
3 5 5 5

numero de grados de libertad:

m1 ≔ 5

Paso 1: Hallamos la matriz de rigidez local de cada elemento finito

Elemento finito 1:

E1 = 200 GPa A1 = 800 mm 2 L1 ≔ y ⎛⎝-2.5 m⎞⎠ = 2.54 m

Cx ≔ 1 Lx ≔ 2

A1 ⋅ E1 ⎡ 1 -1 ⎤ Cx ≔ 1
K1 ≔ ―― ⋅⎢ ⎥
L1 ⎣ -1 1 ⎦ Lx ≔ 2

Remplazando la matriz queda como:

Cx ≔ 1 Lx ≔ 2

⎡ 62.982 -62.982 ⎤ kN Cx ≔ 1
K1 = ⎢ ⎥ ――
⎣ -62.982 62.982 ⎦ mm Lx ≔ 2
Elemento finito 2:

E2 = 70 GPa A2 = 1000 mm 2 L2 ≔ y ⎛⎝1.5 m⎞⎠ = 2.843 m

Cx ≔ 3 Lx ≔ 4

A2 ⋅ E2 ⎡ 1 -1 ⎤ Cx ≔ 3
K2 ≔ ――⋅ ⎢ ⎥
L2 ⎣ -1 1 ⎦ Lx ≔ 4

Remplazando la matriz queda como:

Cx = 3 Lx = 4

⎡ 24.621 -24.621 ⎤ kN Cx = 3
K2 = ⎢ ⎥ ―― Lx = 4
⎣ -24.621 24.621 ⎦ mm
Paso 2: Ensamblamos la matriz

⎡ 62.982 -62.982 0 0 0⎤
⎢ -62.982 62.982 0 0 0⎥
⎢ ⎥ kN
K ≔ k1 + k2 = ⎢ 0 0 24.621 -24.621 0 ⎥ ――
⎢ 0 0 -24.621 24.621 0 ⎥ mm
⎢⎣ 0 0 0 0 0 ⎥⎦

Paso 3: Hallamos la matriz modificada:

Calculo de "C":
N
C ≔ max ⎛⎝K⎞⎠ ⋅ 10 4 = 629823697.562 ――
mm

Condicion mul punto:


-1
B1 ≔ 1 B2 ≔ ―― = -0.5
2
p1 ≔ 2 p2 ≔ 5
⎡ C ⋅ B1 2 C ⋅ B1 ⋅ B2 ⎤ p1 ≔ 2
c1 ≔ ⎢ 2 ⎥
⎣ C ⋅ B1 ⋅ B2 C ⋅ B2 ⎦ p2 ≔ 5

p1 = 2 p2 = 5
⎡ 629823697.562 -314911848.781 ⎤ N p1 = 2
c1 = ⎢ ⎥ ――
⎣ -314911848.781 157455924.391 ⎦ mm p = 5
2

Condicion mul punto:


-3
B1 ≔ 1 B2 ≔ ―― = -0.6
5

p1 ≔ 4 p2 ≔ 5

⎡ C ⋅ B1 2 C ⋅ B1 ⋅ B2 ⎤ p1 ≔ 4
c2 ≔ ⎢ 2 ⎥
p2 ≔ 5
⎣ C ⋅ B1 ⋅ B2 C ⋅ B2 ⎦

p1 = 4 p2 = 5
⎡ 629823697.562 -377894218.537 ⎤ N p1 = 4
c2 = ⎢ ⎥ ―― p2 = 5
⎣ -377894218.537 226736531.122 ⎦ mm

La matriz modificada sera:


⎡ 629886.68 -62.982 0 0 0 ⎤
⎢ -62.982 629886.68 0 0 -314911.849 ⎥
⎢ ⎥ kN
Kmod ≔ K + C1 + C2 = ⎢ 0 0 629848.318 -24.621 0 ⎥ ――
⎢ 0 0 -24.621 629848.318 -377894.219 ⎥ mm
⎢⎣ 0 -314911.849 0 -377894.219 384192.456 ⎥⎦
Paso 4: Determinar los desplazamientos
-1
⎡ Q1 ⎤ ⎛⎡ 629886.68 -62.982 0 0 0 ⎤ ⎞ ⎡ 0 ⎤ ⎡ 0.00001879 ⎤
⎢ Q2 ⎥ ⎜⎢ -62.982 629886.68 0 0 -314911.849 ⎥ ⎟ ⎢ 0 ⎥ ⎢ 0.18794075 ⎥
⎢ ⎥ ⎜⎢ ⎥ kN ⎟ ⎢ ⎥ ⎢ ⎥
⎢ Q3 ⎥ ≔ ⎜⎢ 0 0 629848.318 -24.621 0 ⎥ ―― ⎟ ⋅ ⎢ 0 ⎥ = ⎢ 0.00000882 ⎥ mm
⎢ Q4 ⎥ ⎜⎢ 0 0 -24.621 629848.318 -377894.219 ⎥ mm ⎟ ⎢ 0 ⎥ ⎢ 0.22554264 ⎥
⎢⎣ Q5 ⎥⎦ ⎜⎝⎢⎣ 0 -314911.849 0 -377894.219 384192.456 ⎦⎥ ⎟
⎠ ⎢⎣ T ⎥⎦ ⎢⎣ 0.37591909 ⎥⎦

Q1 = 0.000018792 mm

Q2 = 0.187940753 mm

Q3 = 0.000008817 mm

Q4 = 0.225542638 mm

Q5 = 0.375919091 mm

Paso 5: Determinacion de esfuerzos en cada elemento finito:

Elemento finito 1 :

1 ⎡ Q1 ⎤
σ1 ≔ E1 ⋅ ―⋅ ⎡⎣ -1 1 ⎤⎦ ⋅ ⎢ ⎥ = 14.795 MPa σ1 = 14.795 MPa
L1 ⎣ Q2 ⎦

Elemento finito 2

1 ⎡ Q3 ⎤
σ2 ≔ E2 ⋅ ―⋅ ⎡⎣ -1 1 ⎤⎦ ⋅ ⎢ ⎥ = 5.553 MPa σ2 = 5.553 MPa
L2 ⎣ Q4 ⎦

También podría gustarte