單元對象是構(gòu)成整個(gè)結(jié)構(gòu)對象的基本要素，如桿單元，梁單元，板單元，殼單元等等。雖然單元形狀和特性各不相同，但基本特征和功能是相同的。比如都具有一定的幾何形狀，通過節(jié)點(diǎn)與其它單元連接，包含材料信息，在結(jié)構(gòu)分析中各單元皆以單元?jiǎng)偠染仃嚨男问浇M裝成整體結(jié)構(gòu)。各種單元的層次關(guān)系如圖所示

▲單元層次關(guān)系

#?抽象單元類(基類)
class?AbstractElement:
????def?__init__?(self,?id):
????????self.id?=?id

對于一個(gè)桁架單元，其特征有編號、材料類型、橫截面面積、兩個(gè)節(jié)點(diǎn)的信息、單元長度、單元局部坐標(biāo)與總體坐標(biāo)的夾角、單元局部坐標(biāo)與總體坐標(biāo)的轉(zhuǎn)換矩陣、單元局部坐標(biāo)系下的剛度矩陣、單元整體坐標(biāo)系下的剛度矩陣、應(yīng)力矩陣、應(yīng)變矩陣等等。

#?平面桁架單元類
class?TrussElement2D(AbstractElement):
????def?__init__?(self,id,?A,?mat,?node1,?node2):
????????super().__init__(id)
????????self.A?=?A
????????self.mat?=?mat
????????self.node1?=?node1
????????self.node2?=?node2

????def?elemDisp(self):
????????return?np.array(?[?[0.],?[0.],?[0.],?[0.]?]?)
??????????????????????????


????def?elemLength(self):
????????dx?=?self.node2.coord_X?-?self.node1.coord_X
????????dy?=?self.node2.coord_Y?-?self.node1.coord_Y
????????return??math.sqrt(dx*dx?+?dy*dy)

????def?COS(self):
????????dx?=?self.node2.coord_X?-?self.node1.coord_X
????????return?dx?/?self.elemLength()

????def?SIN(self):
????????dy?=?self.node2.coord_Y?-?self.node1.coord_Y
????????return??dy?/?self.elemLength()

????def?Trans(self):
????????S?=?self.SIN()
????????C?=?self.COS()
????????return?np.array([?[C,?S,?0,?0],
??????????????????????????[0,?0,?C,?S]?])

????def?elemStiffnessMatrix(self):
????????EA?=?self.mat.E?*?self.A
????????S?=?self.SIN()
????????C?=?self.COS()
????????ek?=?np.array([?[C*C,??C*S,?-C*C,?-C*S?],?
???????????????????????[C*S,??S*S,?-C*S,?-S*S?],?
???????????????????????[-C*C,?-C*S,?C*C,??C*S],?
???????????????????????[-C*S,?-S*S,?C*S,??S*S]???]?)

????????return?EA/self.elemLength()?*?ek

????def?strainMatrix(self):
????????L??=?self.elemLength()
????????tmp1?=?np.array(?[-1.0,?1.0]?)
????????return?1/L?*?np.array(?[-1.0,?1.0]?)
????????
????def?normalForce(self,?T,?elemDisp):
????????EA?=?self.mat.E?*?self.A
????????b?=?self.strainMatrix()?
????????tmp1?=?EA?*?b
????????
????????tmp2?=?np.dot(tmp1,?T)
????????return?np.dot(tmp2,?elemDisp)

采用從基類按層次繼承來建立單元類的方法，除了重用父類代碼之外，另一個(gè)很重要的優(yōu)點(diǎn)就是可以運(yùn)用多態(tài)。通過指向單元基類對象的指針，調(diào)用各派生單元類的成員函數(shù)。來看一個(gè)例子：從一個(gè)抽象單元類派生兩個(gè)不同的單元類?，F(xiàn)在要獲取這兩個(gè)單元的編號。

#include?
#include?
using?namespace?std;

class?AbstractElement
{
public:

????virtual?~AbstractElement()?=?default;
????virtual?size_t?getID()?const?=?0;???//?純虛函數(shù)?
};


class?Truss?:?public?AbstractElement
{
protected:
????size_t?id{?0?};

public:
????Truss()?=?default;
????Truss(size_t?id_)?:?id{?id_?}?{}
????size_t?getID()?const?override?{?return?id;?}


};


class?Beam?:?public?AbstractElement
{
protected:
????size_t?id{?0?};

public:
????Beam()?=?default;
????Beam(size_t?id_)?:?id{?id_?}?{}
????size_t?getID()?const?override?{?return?id;?}
};

int?main(){
????Truss?t1{?1?};
????Beam??b1{?2?};

????std::vector?elems{?&t1,?&b1?};for?(auto*?el?:?elems)
????{std::cout?<getID()?<std::endl;??//這里體現(xiàn)了多態(tài)
????}return?0;
}

總結(jié)

以上是生活随笔為你收集整理的c++矩阵类_面向对象有限元编程|单元类的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。