声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 10139|回复: 22

[前后处理] 简谐载荷怎么加

[复制链接]
发表于 2007-6-8 17:10 | 显示全部楼层 |阅读模式

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

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

x
我在做一个瞬态动力学分析,要加一个基础振动(或者简谐载荷),问了很多同学, 都做的静力学分析都不太清楚,
不知道ANSYS里面有没有直接可以定义的,还是要用函数来定义,希望做过的朋友能帮忙一下, 先谢谢了。
回复
分享到:

使用道具 举报

发表于 2007-6-8 19:17 | 显示全部楼层
不是用函数,而是分载荷步和子步计算
在一时间末,加对应荷载,分几个子步计算
然后,再取一个时间末,,加对应荷载,分几个子步计算
直至加一个周期就行了
有点笨,谁有简单的方法?探讨探讨

[ 本帖最后由 rodge 于 2007-6-9 21:24 编辑 ]

评分

1

查看全部评分

 楼主| 发表于 2007-6-8 19:48 | 显示全部楼层
呵呵,那对应的载荷是不是要自己先计算好?那简谐载荷一个周期要加好多个载荷步啊,关键好象没什么人做这个方面的分析,问了好多人都不太清楚,大家来帮忙讨论讨论啊,先谢谢大家
发表于 2007-6-8 20:24 | 显示全部楼层

看看ANSYS帮助文件里有没有

 楼主| 发表于 2007-6-8 20:36 | 显示全部楼层
我大概找了一下, 没有发现相关的说明,关键比较冷门,呵呵
发表于 2007-6-9 14:15 | 显示全部楼层
你这个是谐响应分析吧,可以直接加简谐荷载的
 楼主| 发表于 2007-6-9 14:45 | 显示全部楼层
不是谐响应分析,我这个是非线性的,所以不能谐响应分析的,只能瞬态动力分析,要加谐载荷
发表于 2007-6-11 08:50 | 显示全部楼层
LS-DYNA来处理吧
发表于 2007-6-11 08:51 | 显示全部楼层
施加的周期性载荷可以做得拿DYNA
发表于 2007-6-11 18:30 | 显示全部楼层
contact_T=0.6e-3                             !周期
dt=0.3e-4
re_n=contact_T/dt
!dt=0.001                              !设置循环时间间隔
*do,j,1,re_n                          !执行循环的次数是N次
kbc,1                               !荷载变化方式为渐变
outres,all,all                      !把所有子步的所有结果输入到结果文件
nsubst,1,,,1                        !设置子步
time,dt*j                           !设置荷载步
fk,1000,fx,12000*sin(pi*j*dt/contact_T)               !对应加一正弦荷载,

评分

1

查看全部评分

 楼主| 发表于 2007-6-13 11:38 | 显示全部楼层
恩,我试试,谢谢了,
发表于 2007-6-17 11:55 | 显示全部楼层

回复 #1 hailin 的帖子

我这有个和你类似的例子,希望对你有帮助:

梁的正弦加载的瞬态分析

模型描述:
1)悬臂梁顶端施加X向正弦变化的载荷,进行瞬态动力求解;
2)单元为beam188;
3)/post26中能提取末端的位移结果。

模型如下图,已经通过模态、谐响应计算(如下图),但瞬态计算总是无法得到理想的结果。

问题:
我认为,梁顶端的位移输出为一正弦曲线,但ANSYS计算结果的时间-位移历程总是无法得到正弦曲线,不知是何原因?我觉得是力载荷的加载有问题,主要出现在激励频率的确定上,但总是调试不出来。请高手指点。

命令流:
/FILNAME,BTSanlysis
/UNITS,SI
/TITLE,Bond Tool计算

D=3e-3!直径
L=15e-3 !长度
S=3.141593*(D/2)*(D/2)!截面积
I=3.141593*D*D*D*D/64

/PREP7!建模
ET,1,BEAM188!梁单元 3,188,189
MP,EX,1,620.53E9!弹性模量
MP,DENS,1,14880
MP,PRXY,1,0.29
!MP,DAMP,1,0.0085

secnum,1
sectype,1,beam,CSOLID
secdata,D/2,10,3!

N,1,0,0
N,21,0,-L
FILL,1,21
E,1,2
EGEN,20,1,1,1,1

/PNUM,NODE,1!图示节点及编号
/PNUM,ELEM,1!图示单元及编号
EPLOT!显示单元
FINISH

/SOLU
ANTYPE,2!模态分析
EQSLV,SPAR
MXPAND,100, , ,1
LUMPM,0
PSTRES,0
MODOPT,LANB,100,5000,100000, ,OFF

OUTPR,BASIC,ALL
d,1,uy,,,,,uz,rotx,roty,rotz
SOLVE
FINISH
/POST1

!做谐响应分析
/SOLU
ANTYPE,HARMIC! HARMONIC RESPONSE ANALYSIS
HROUT,OFF! Print results as amplitudes and phase angles
HROPT,full

OUTPR,BASIC,1
harfrq,5000,100000! FREQUENCY RANGE FROM 0 TO 7.5 HZ
NSUBST,200! 30 INTERVALS WITHIN FREQ. RANGE
KBC,1! STEP BOUNDARY CONDITION

d,1,uy,,,,,uz,rotx,roty,rotz!very important!!
!D,1,UX,10
F,1,FX,200

ALLS
SOLVE
FINISH

/POST26
NSOL,2,21,U,X,UX! STORE UX DISPLACEMENT OF NODE 21作为变!量2
!ESOL,3,21,,SMISC,1,F1! STORE 摩擦力FORCE F1 OF ELEMENT 1 AS VARIABLE 3
                         !由K1引起的弹性力F1.
!PRVAR,2,3! PRINT VARIABLES 2 AND 3
/GRID,1! TURN GRID ON
/AXLAB,Y,DISP! Y AXIS LABEL AS DISP
/GTHK,CURVE,2! CURVE LINES THICKNESS RATIO OF 2
PLVAR,2

!做瞬态分析
/SOLu
antype,trans
trnopt,full
outpr,all,all
outres,all,all
!alphad,0.014
!betad,0.000092

autots,on
kbc,1
!tintp,,0.25,0.5,0.5
d,1,uy,,,,,uz,rotx,roty,rotz

*do,i,1,2000,1
t=i/100000
time,t
nsubst,50
!施加约束条件
Pressure=10*sin(2*3.14159*91743*t) !如何设置激励力的频率?f为其中的FX方 向振动的固有频率?
F,1,FX,Pressure

solve
lswr,i
*enddo
fini

/POST26
NSOL,2,21,U,X,UX! STORE UX DISPLACEMENT OF NODE 21作为变!量2
!ESOL,3,21,,SMISC,1,F1! STORE 摩擦力FORCE F1 OF ELEMENT 1 AS VARIABLE 3
                         !由K1引起的弹性力F1.
!PRVAR,2,3! PRINT VARIABLES 2 AND 3
/GRID,1! TURN GRID ON
/AXLAB,Y,DISP! Y AXIS LABEL AS DISP
/GTHK,CURVE,2! CURVE LINES THICKNESS RATIO OF 2
PLVAR,2

评分

1

查看全部评分

 楼主| 发表于 2007-6-17 22:56 | 显示全部楼层
恩, 我参考参考, 最近正在看LS-DYNA的例子,准备用DYNA来做了
顺便问下,有人做LS-DYNA的吗,我的QQ540630,希望能请教一些这个方面的问题

[ 本帖最后由 xuruikl 于 2007-6-17 23:16 编辑 ]
发表于 2007-8-22 11:21 | 显示全部楼层
原帖由 cyh3278 于 2007-6-11 18:30 发表
contact_T=0.6e-3                             !周期
dt=0.3e-4
re_n=contact_T/dt
!dt=0.001                              !设置循环时间间隔
*do,j,1,re_n                          !执行循环的次数是 ...

是否应该改为KBC=0
发表于 2008-1-2 14:29 | 显示全部楼层
可以先将载荷保存成一个txt文件(两列:时间  载荷)
然后在parameters_array parameters_define/edit_add_调入txt文件
最后,像施加一般载荷一样,define loads_
注意 选existing table
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 09:12 , Processed in 0.076065 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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