matlab 有限單元法

平面8節(jié)點(diǎn)等參元完整程序

module Elem_Rect8 ! 八節(jié)點(diǎn)等參元

implicit none

integer (kind(1)),parameter ::ikind=(kind(1))

integer (kind(1)),parameter ::rkind=(kind(0.d0))

type :: typ_Kcol

real(rkind),pointer :: Row(:)

end type typ_Kcol

type :: typ_GValue !總體控制變量

integer(ikind) :: NNode, NElem, NLoad, NMat, NSupport

integer(ikind) :: NGlbDOF !整體自由度總數(shù)

integer(ikind) :: NGENS, NodeDOF,ElemNodeNo

integer(ikind) :: NInt

end type typ_GValue

type Typ_Node !定義節(jié)點(diǎn)類型

real(rkind) :: coord(2) !節(jié)點(diǎn)坐標(biāo)

integer(ikind) :: GDOF(2) !整體自由度編碼

real(rkind) :: DISP(2) !節(jié)點(diǎn)位移

real(rkind) :: dDISP(2) !節(jié)點(diǎn)位移增量

real(rkind) :: dForce(2) !節(jié)點(diǎn)不平衡力

end type typ_Node

!=============================================================================

Type typ_IntPoint !定義積分點(diǎn)參數(shù)

real(rkind) :: EPS(3) !應(yīng)變

real(rkind) :: SIG(3) !應(yīng)力

real(rkind) :: D(3,3) !本構(gòu)矩陣

real(rkind) :: B(3,16) !幾何矩陣

real(rkind) :: DETJ !雅克比行列式

end type Typ_IntPoint

type Typ_Rect8 !定義實(shí)體單元

integer(ikind) :: NodeNo(8) !節(jié)點(diǎn)編號

real(rkind) :: E !彈性模量

real(rkind) :: u !泊松比

real(rkind) :: t !單元厚度

real(rkind) :: EK(16,16) !單元剛度矩陣

type(typ_intpoint) :: IntP(9) !積分點(diǎn)

總結(jié)

以上是生活随笔為你收集整理的matlab八节点六面体程序,平面8节点等参元完整程序的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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