From 7b360c7d0d7bae3fdcda6d3422bafae2bdf5854e Mon Sep 17 00:00:00 2001
From: Ernesto Oquelis <ernesto.oquelis@gmail.com>
Date: Sun, 23 Sep 2018 19:19:31 +0200
Subject: [PATCH] update
---
TEST_2.py | 8 ++++----
fem_2d.py | 11 +++++++----
2 files changed, 11 insertions(+), 8 deletions(-)
diff --git a/TEST_2.py b/TEST_2.py
index 7b97f7c..cbe7d16 100644
--- a/TEST_2.py
+++ b/TEST_2.py
@@ -19,12 +19,12 @@ STR = Structure2DPlaneStress(mesh_elements, mesh_points,
poisson_ratio = 0.3,
thickness = 0.1)
-NODES = STR.nodes
+STRUCT = list()
+STRUCT.append(STR)
+NODES = STR.nodes
ELEMENTS = STR.elements
-STR.nodes[2].force_2d = Force2D([1,1])
+STR.nodes[2].force_2d = Force2D([1,1])
STR.solve()
-A = list()
-A.append(STR)
\ No newline at end of file
diff --git a/fem_2d.py b/fem_2d.py
index 814674c..40981db 100644
--- a/fem_2d.py
+++ b/fem_2d.py
@@ -176,11 +176,14 @@ class Structure2DPlaneStress:
node_ids = self.elements[element_id].node_ids
k_e = self.elements[element_id].get_k_mat(self.thickness,
self.c_mat)
- # First [0,0] component of k_e
-
+ for (i, j) in ((0, 0), (0, 1), (1, 0), (1, 1)):
+ k_global[node_ids[i], node_ids[j]] = (k_global[node_ids[i],
+ node_ids[j]] +
+ k_e[i, j])
+
+
+
-
-
'''
--
GitLab