参考网址:http://blog.sina.com.cn/s/blog_65936c230100l20i.html
案例1
说明
首先建立两个节点1和2,然后由节点1和节点2生成弹簧单元,再由2节点生成质量单元,然后进行分析。一阶弹簧系统,可以分析出其固有频率。
APDL代码
FINISH
/CLEAR
/PREP7
ET,1,14 !定义弹簧combin14单元
R,1,1 !刚度为1,方向为两个节点的连线方向
ET,2,21 !定义质量单元mass21
R,2,1 !设置实常数,质量为1
!定义两个节点,1、2为弹簧的两个端点,并且2节点作为为质点附着在弹簧的一端,从而形成弹簧振子
N,1
N,2,10 !任意
!单元划分:定义弹簧单元
TYPE,1
REAL,1
E,1,2
!定义质点单元
TYPE,2
REAL,2
E,2
!求解
/SOLVE
ANTYPE,2 !分析类型:模态分析
MODOPT,LANB,2 !取前两阶(实际只有一阶:单自由度)
MXPAND,2, , ,0 !扩展模态
D,1,ALL !将弹簧的另一端点全约束
!求解
SOLVE
FINISH
/POST1
SET,LIST
结果
案例2
说明
建立三个节点1,、2、3,其中节点2和节点3重合,由1、2节点生成弹簧单元,再由3节点生成质量单元,然后将2和3节点进行耦合,然后进行分析。
APDL代码
FINISH
/CLEAR
/PREP7
ET,1,14 !定义弹簧combin14单元
R,1,1 !刚度为1,方向为两个节点的连线方向
ET,2,21 !定义质量单元mass21
R,2,1 !设置实常数,质量为1
!定义三个节点,1、2为弹簧的两个端点,3为质点附着在弹簧的一端,从而形成弹簧振子
N,1
N,2,10 !此处的长度可任意,此处取10
N,3,10 !与节点2的坐标相同
!单元划分:定义弹簧单元
TYPE,1
REAL,1
E,1,2
!定义质点单元
TYPE,2
REAL,2
E,3
!将弹簧的端点2与质点3全耦合
CP,1,ALL,2,3
!求解
/SOLVE
ANTYPE,2 !分析类型:模态分析
MODOPT,LANB,2 !取前两阶(实际只有一阶:单自由度)
MXPAND,2, , ,0 !扩展模态
D,1,ALL !将弹簧的另一端点全约束
!求解
SOLVE
FINISH
/POST1
SET,LIST
结果
结果与上图相同
案例3
说明
beam188和combin14联合使用。首先建立梁单元,坐标为(0,0,0)(0,1,0)。然后再梁的两端分别建立3个节点,用于生成弹簧单元,方法如按照上述例子。弹簧刚度为100,梁的材料为钢,然后进行模态分析。不采用耦合的方法
建立点的顺序为
1(0,0,0)12
2(1,0,0)13
(-0.02,0,0)14
(0,0.02,0)15
(0,0,0.02)16
(1.02,0,0)17
(1,0.02,0)18
(1,0,0.02)19
APDL代码
FINISH
/CLEAR
/PREP7
!建立单元beam188
ET,1,BEAM188
!设置beam的section
SECTYPE, 1, BEAM, CTUBE, BEAM1, 0
SECOFFSET, CENT
SECDATA,0.01,0.015,0,0,0,0,0,0,0,0,0,0
!材料1 钢
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1E11
MPDATA,PRXY,1,,0.269
MPDATA,DENS,1,,7850
!建立弹簧单元
ET,2,14 !定义弹簧combin14单元
R,2,100, , , , , ,
!建立梁单元
FLST,3,2,8
FITEM,3,0,0,0
FITEM,3,1,0,0
K, ,P51X
LSTR, 1, 2
!梁结构网格划分
CM,_Y,LINE
LSEL, , , , 1
CM,_Y1,LINE
CMSEL,S,_Y
CMSEL,S,_Y1
LATT,1,2,1, , , ,1
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
FLST,5,1,4,ORDE,1
FITEM,5,1
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
LMESH, 1
/UI,MESH,OFF
!建立节点,用于生成弹簧
FLST,3,8,8
FITEM,3,0,0,0
FITEM,3,1,0,0
FITEM,3,-0.02,0,0
FITEM,3,0,0.02,0
FITEM,3,0,0,0.02
FITEM,3,1.02,0,0
FITEM,3,1,0.02,0
FITEM,3,1,0,0.02
N, ,P51X
!弹簧单元属性
TYPE, 2
REAL, 2
!生成弹簧单元
E,1,14
E,1,15
E,1,16
E,2,17
E,2,18
E,2,19
!耦合梁的端点和弹簧端点
!CP,1,ALL,1,12
!CP,2,ALL,2,13
!添加固定约束--弹簧端点
FLST,2,6,1,ORDE,2
FITEM,2,14
FITEM,2,-19
/GO
D,P51X, , , , , ,ALL, , , , ,
!求解
FINISH
/SOL
ANTYPE,2
MODOPT,LANB,12!分析前十二阶
/STATUS,SOLU
SOLVE
!后处理
/POST1
SET,LIST!打开模态列表