我正在研究处理有限元分析(FEA)的多个节点的计算应变/变形的脚本。对于不同的正输入扭矩有多种解决方案。我想推断这些数据并模拟负扭矩的结果。虽然这样做,原始数据会发生变化。
由于我们没有直接更改任何原始值,我认为必须通过推断有限元分析的函数来参考。我一直在尝试copy.deepcopy
,但想通过多个线程,这不会复制类结构。在other threads it was recommended to inherit,但我正在努力适应我的情况。
以下代码位于包含对象半径相同的所有节点的类中。所有节点都在按其角度排序的列表self._nodes
中。每个节点都具有所有扭矩水平的应变。
class RadialNodeContainer:
def __init__(self, radius):
self._r = radius
self._nodes = []
def compute_negatives_radial_container(self): # ... see below
class Node:
def __init__(self, node_id, x, y, z,
strain_0nm, strain_100nm, strain_204nm, strain_369nm):
self._coordinate = Coordinate(x, y, z)
self._torque_levels = [strain_0nm, strain_100nm,
strain_204nm, strain_369nm]
def get_all_strains_complete(self):
return copy.deepcopy(self._torque_levels)
class Strain:
def __init__(self, torque_nm, exx, exy, eyy):
self._torque = torque_nm
self._exx = exx
self._exy = exy
self._eyy = eyy
导致不需要的原始数据更改的功能:
def compute_negatives_radial_container(self):
points_per_360deg = len(self._nodes)
jj = points_per_360deg
corresponding_node_strains = None
for ii in range(points_per_360deg):
jj -= 1
# Mistake is expected below here
corresponding_node_strains = copy.deepcopy(
self._nodes[jj].get_all_strains_complete())
for kk in range(len(corresponding_node_strains)):
torque = corresponding_node_strains[kk].get_torque_level()
if torque != 0:
exx, exy, eyy = corresponding_node_strains[kk].get_raw_strains()
calculated_negative_strain = Strain(torque_nm=-torque,
exx=exx,
exy=-exy,
eyy=eyy)
self._nodes[ii].add_torque_level(calculated_negative_strain)
我想到了创建一个应变元素列表的deepcopy
(Node -> self._torque_level
)。最初这个列表看起来像[Strain(0Nm), Strain(100Nm), ...]
。但我无法理解我必须适应的代码的哪些部分以允许传递类实例的副本。
详细说明。如果我有时间,我会扩展这个,这就是OP正在寻找的东西。
我还建议使用numpy作为你的数组而不是列表,因为mumpy要快得多。
class Strain:
def __init__(self, torque_nm, exx, exy, eyy):
self._torque = torque_nm
self._exx = exx
self._exy = exy
self._eyy = eyy
self._neg = False
def setNeg(self, neg):
self._neg = neg
@propery
def torque(self):
return _torque if not self._neg else -self.torque
@propery
def exx(self):
return _exx
@propery
def exy(self):
return _exy if not self._neg else -self._exy
@propery
def eyy(self):
return _eyy