matlab 有限單元法

平面8節點等參元完整程序

module Elem_Rect8 ! 八節點等參元

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 !整體自由度總數

integer(ikind) :: NGENS, NodeDOF,ElemNodeNo

integer(ikind) :: NInt

end type typ_GValue

type Typ_Node !定義節點類型

real(rkind) :: coord(2) !節點坐標

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

real(rkind) :: DISP(2) !節點位移

real(rkind) :: dDISP(2) !節點位移增量

real(rkind) :: dForce(2) !節點不平衡力

end type typ_Node

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

Type typ_IntPoint !定義積分點參數

real(rkind) :: EPS(3) !應變

real(rkind) :: SIG(3) !應力

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

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

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

end type Typ_IntPoint

type Typ_Rect8 !定義實體單元

integer(ikind) :: NodeNo(8) !節點編號

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

real(rkind) :: u !泊松比

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

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

type(typ_intpoint) :: IntP(9) !積分點

總結

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

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。