-
Notifications
You must be signed in to change notification settings - Fork 0
/
exemplo.py
82 lines (60 loc) · 1.31 KB
/
exemplo.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
# -*- coding: utf-8 -*-
"""
@author: César Eduardo Petersen
@date: 18/08/2021
Utiliza o método dos elementos finitos para resolver um pórtico
esquema da estrutura:
D-------E----G
/| b2 | b5
/ | |b3
b0/ |b1 F
/ | |b4
/ | |
A B C
"""
from pyFFEM.pyFFEM import *
# unidades: kN, m
# define material e seção
mat = Material(2e8, 4e-3, 1.5e-4)
# define os nós
a=5.2
b=3.6
A = Node(0, 0)
B = Node(a, 0)
C = Node(2*a, 0)
D = Node(a, b)
E = Node(2*a, b)
F = Node(2*a, b/2)
G = Node(2.25*a, b)
# define os apoios
A.apoio = [True, True, False]
B.apoio = [True, True, True]
C.apoio = [True, True, False]
# aplica forças nodais
P = 36
H = 26
D += ForcaNodal(0, -P)
G += ForcaNodal(0, -P)
F += ForcaNodal(H, 0)
# define as barras
b0 = Barra(A, D, mat)
b1 = Barra(B, D, mat)
b2 = Barra(D, E, mat)
b3 = Barra(E, F, mat)
b4 = Barra(F, C, mat)
b5 = Barra(E, G, mat)
# aplica forças distribuidas
q = 16
b1 += ForcaDistribuida(q)
b2 += ForcaDistribuida(0, -q)
# monta a estrutura
estrutura = Estrutura([A, B, C, D, E, F, G],
[b0, b1, b2, b3, b4, b5])
# exibe a estrutura
#estrutura.exibir()
# calcula
estrutura.calcular()
estrutura.exibirDeformado(750)
#estrutura.diagramaMomento()
#estrutura.diagramaCortante()
#estrutura.diagramaNormal()