声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2470|回复: 10

[前后处理] ANSYS动力学分析的几个入门例子

[复制链接]
发表于 2006-9-2 18:03 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
问题一:悬臂梁受重力作用发生大变形,求其固有频率。

基本过程:
1、建模
2、静力分析
  NLGEOM,ON
  PSTRES,ON
3、求静力解
4、开始新的求解:modal
  PSTRES,ON
  UPCOORD,1,ON 修正坐标
  PSOLVE...
5、扩展模态解
6、察看结果

/PREP7
ET,1,BEAM189 !使用beam189梁单元
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,210e9
MPDATA,PRXY,1,,0.3
MPDATA,DENS,1,,7850
SECTYPE, 1, BEAM, RECT, secA, 0 !定义梁截面secA
SECOFFSET, CENT
SECDATA,0.005,0.01,0,0,0,0,0,0,0,0
K, ,,,, !建模与分网
K, ,2,,,
K, ,2,1,,
LSTR, 1, 2
LATT,1, ,1, , 3, ,1
LESIZE,1, , ,20, , , , ,1
LMESH, 1
FINISH

/SOL !静力大变形求解
ANTYPE,0
NLGEOM,1
PSTRES,ON !计及预应力效果
DK,1, , , ,0,ALL, , , , , ,
ACEL,0,9.8,0, !只考虑重力作用
TIME,1
AUTOTS,1
NSUBST,20, , ,1
KBC,0
SOLVE
FINISH

/SOLUTION
ANTYPE,2 !进行模态求解
MSAVE,0
MODOPT,LANB,10
MXPAND,10, , ,0 !取前十阶模态
PSTRES,1 !打开预应力效应
MODOPT,LANB,10,0,0, ,OFF
UPCOORD,1,ON !修正坐标以得到正确的应力
PSOLVE,TRIANG !三角化矩阵
PSOLVE,EIGLANB !提取特征值和特征向量
FINISH
/SOLU
EXPASS,1 !扩展模态解
PSOLVE,EIGEXP
FINISH

/POST1
SET,LIST !观察结果
FINISH

评分

1

查看全部评分

回复
分享到:

使用道具 举报

 楼主| 发表于 2006-9-2 18:03 | 显示全部楼层
问题二:循环对称结构模态分析
这是ANSYS HELP里的例子,但那个命令流似乎有些问题,下面是整理过的命令流。

基本过程:
1、建模
2、define cyclic
3、定义约束
4、求模态解
5、展开并察看结果

r1=5 !建模
r2=10
d1=2
nsect=24
alpha_deg=360/nsect
alpha_rad=2*acos(-1)/nsect

/prep7
csys,1
k,1,0,0,0
k,2,0,0,d1
k,3,r1,0,0
k,4,r1,0,d1
l,3,4
arotat,1,,,,,,1,2,alpha_deg/2
k,7,r2,0,0
k,8,r2,0,d1
l,7,8
arotat,5,,,,,,1,2,alpha_deg/2
arotat,2,,,,,,1,2,alpha_deg/2
arotat,6,,,,,,1,2,alpha_deg/2
a,5,6,10,9
mshkey,1
et,1,181
r,1,0.20
r,2,0.1
mp,ex,1,10e6
mp,prxy,1,0.3
mp,dens,1,1e-4
esize,0.5
asel,,,,1,4
aatt,,1
asel,,,,5
aatt,,2
allsel

CYCLIC, , , ,'CYCLIC' !定义基本扇区
AMESH,all !分网
FINISH

/SOL !求模态解
ANTYPE,2
MODOPT,LANB,10
MXPAND,10, , ,0
PSTRES,0
MODOPT,LANB,10,0,0, ,OFF
DK,5, , , ,0,UZ, , , , , ,
SOLVE
FINISH

/POST1
SET,LIST
SET, , ,1, , , ,4,
/CYCEXPAND, ,ON
PLNS,U,SUM !观察扩展模态
FINISH
 楼主| 发表于 2006-9-2 18:04 | 显示全部楼层
问题三:三角平台受谐波载荷作用,求结构响应
谐波载荷为作用在平台上面一点的集中力,幅值为10N,频率范围5Hz~70Hz

基本过程:
1、建模
2. 求模态解
3、用模态叠加法作谐分析
4. 观察结果

/PREP7 !建模
ET,1,BEAM189
ET,2,SHELL93
R,1,0.01, , , , , ,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE, 1, BEAM, RECT, secA, 0
SECOFFSET, CENT
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0
K, ,-0.5,,,
K, ,0.5,,,
K, ,,,1,
K, ,,1,1,
K, ,-0.5,1,,
K, ,0.5,1,,
K, ,,,0.5,
A,4,5,6
LSTR, 1, 5
LSTR, 3, 4
LSTR, 2, 6
LSEL,S,LINE,,4,6
LATT,1,1,1, , 7, ,1
LSEL, , , ,ALL
LESIZE,ALL, , ,10, , , , ,1
LMESH,4,6
TYPE, 2
MSHAPE,0,2D
MSHKEY,1
AMESH,1
FINISH

/SOL !为了使用模态叠加法谐分析
ANTYPE,2
MODOPT,LANB,20 !先取结构前20阶模态
MXPAND,20, , ,0
MODOPT,LANB,20,0,0, ,OFF
KSEL,S,KP,,1,3
DK,ALL, , , ,0,ALL, , , , , ,
SOLVE
FINISH

!/POST1 !最好事先看一下模态结果
!SET,LIST !以了解模态频率范围
!FINISH

/SOL !模态叠加法谐响应分析
ANTYPE,3
HROPT,MSUP
HROUT,ON
HROPT,MSUP,20, ,
HROUT,ON,OFF,0
F,177,FY,-10, !施加10N的力在节点177的负Y方向
HARFRQ,5,70, !载荷的强制频率范围
NSUBST,200,
KBC,1
ALPHAD,5,
SOLVE
FINISH

/POST26 !察看位移响应
FILE,,rfrq
NUMVAR,20
NSOL,2,139,U,Y,
PLVAR,2
FINISH
 楼主| 发表于 2006-9-2 18:04 | 显示全部楼层
问题四:三角平台受一地震谱激励,求应力分布和支反力

基本过程:
1、建模
2、求模态解
3、求谱解
4、扩展模态
5、模态合并
6、观察结果

/PREP7 !建模
ET,1,BEAM189
ET,2,SHELL93
R,1,0.01, , , , , ,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE, 1, BEAM, RECT, secA, 0
SECOFFSET, CENT
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0

K, ,-0.5,,,
K, ,0.5,,,
K, ,,,1,
K, ,,1,1,
K, ,-0.5,1,,
K, ,0.5,1,,
K, ,,,0.5,
A,4,5,6
LSTR, 1, 5
LSTR, 3, 4
LSTR, 2, 6
LSEL,S,LINE,,4,6
LATT,1,1,1, , 7, ,1
LSEL, , , ,ALL
LESIZE,ALL, , ,10, , , , ,1
LMESH,4,6
TYPE, 2
MSHAPE,0,2D
MSHKEY,1
AMESH,1
FINISH

/SOL !取前十阶模态
ANTYPE,2
MODOPT,LANB,10
KSEL,S,KP,,1,3
DK,ALL, , , ,0,ALL, , , , , ,
SOLVE
FINISH

/SOL !谱分析
ANTYPE,8
SPOPT,SPRS,10,1
SVTYP,3
SED,0,1,0, !给出激励方向
FREQ,0.25,1.34,6.73,14.6,28.9,0,0,0,0 !激励谱
SV,0,0.00073,0.00016,0.00034,0.00034,0.00052,
SOLVE
FINISH

/SOL !扩展模态
ANTYPE,2
EXPASS,1
MXPAND,10,0,0,1,0.001,
SOLVE
FINISH

/SOL !模态合并
ANTYPE,8
SRSS,0.001,DISP
SOLVE

/POST1 !观察结果
SET,LIST
/INPUT,,mcom
PRRSOL,F
FINISH
 楼主| 发表于 2006-9-2 18:04 | 显示全部楼层
问题五:三角平台受时程载荷作用,求应力分布和变形过程

基本过程:
1、建模
2. 施加随时间变化载荷,定义载荷步
3、求解
4. 在POST1和POST26中观察结果

/PREP7 !建模
ET,1,BEAM189
ET,2,SHELL93
R,1,0.01, , , , , ,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE, 1, BEAM, RECT, secA, 0
SECOFFSET, CENT
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0
K, ,-0.5,,,
K, ,0.5,,,
K, ,,,1,
K, ,,1,1,
K, ,-0.5,1,,
K, ,0.5,1,,
K, ,,,0.5,
A,4,5,6
LSTR, 1, 5
LSTR, 3, 4
LSTR, 2, 6
LSEL,S,LINE,,4,6
LATT,1,1,1, , 7, ,1
LSEL, , , ,ALL
LESIZE,ALL, , ,10, , , , ,1
LMESH,4,6
TYPE, 2
MSHAPE,0,2D
MSHKEY,1
AMESH,1
FINISH

/SOL
ANTYPE,4
TRNOPT,FULL !完全法瞬态分析
KSEL,S,KP,,1,3
DK,ALL, , , ,0,ALL, , , , , ,
OUTRES,ALL,1
ALPHAD,5, !α阻尼

TIME,2 !定义载荷曲线
AUTOTS,1
NSUBST,50, , ,1
KBC,0
SFA,ALL,1,PRES,500
LSWRITE,1,

TIME,3
LSWRITE,2,

TIME,4
SFA,ALL,1,PRES,150
KBC,1
LSWRITE,3,

TIME,5
SFA,ALL,1,PRES,
LSWRITE,4,

LSSOLVE,1,4,1, !求解
FINISH
  
/POST26 !观察变形随时间的变化
NSOL,2,177,U,Y,uyy
PLVAR,2
FINISH

评分

1

查看全部评分

发表于 2006-9-2 20:34 | 显示全部楼层
好贴,好东西,支持
发表于 2006-9-3 14:30 | 显示全部楼层
楼主能搜集这么多同类的例子,确实值得表扬,谢谢拿来共享
本贴就暂不关闭主题,希望大家能把同类的例子都粘过来
发表于 2013-4-15 21:58 | 显示全部楼层
表示是好东东,和ansys书里介绍的操作差不多,只是用APDL写出来了
发表于 2013-4-17 15:02 | 显示全部楼层
发表于 2013-10-21 09:25 | 显示全部楼层
最近在做Blocked Lanczos法的模态分析   各位大大有这方面的例子么
发表于 2013-10-29 21:52 | 显示全部楼层




您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 05:06 , Processed in 0.070433 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表