简谐载荷怎么加

hailin

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

mdj219

不是用函数,而是分载荷步和子步计算
在一时间末,加对应荷载，分几个子步计算
然后,再取一个时间末,,加对应荷载，分几个子步计算
直至加一个周期就行了
有点笨，谁有简单的方法？探讨探讨

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

hailin

呵呵，那对应的载荷是不是要自己先计算好？那简谐载荷一个周期要加好多个载荷步啊，关键好象没什么人做这个方面的分析，问了好多人都不太清楚，大家来帮忙讨论讨论啊，先谢谢大家

hao1982

hailin

我大概找了一下， 没有发现相关的说明，关键比较冷门，呵呵

sysh320

你这个是谐响应分析吧，可以直接加简谐荷载的

hailin

不是谐响应分析，我这个是非线性的，所以不能谐响应分析的，只能瞬态动力分析，要加谐载荷

Realplayer

LS-DYNA来处理吧

Realplayer

施加的周期性载荷可以做得拿DYNA

cyh3278

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)               !对应加一正弦荷载，

hailin

恩，我试试，谢谢了，

jiangxuchu

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

梁的正弦加载的瞬态分析

模型描述：
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

hailin

恩， 我参考参考， 最近正在看LS-DYNA的例子，准备用DYNA来做了
顺便问下，有人做LS-DYNA的吗，我的QQ540630，希望能请教一些这个方面的问题

[ 本帖最后由 xuruikl 于 2007-6-17 23:16 编辑 ]

16443

原帖由 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

bonzer

可以先将载荷保存成一个txt文件（两列：时间  载荷）
然后在parameters_array parameters_define/edit_add_调入txt文件
最后，像施加一般载荷一样，define loads_
注意 选existing table