请问在瞬态分析中如何施加几个频率不同的正弦载荷？相位该如何加进去？

zhangjianzaijiu · 发表于 2009-7-7 16:43

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

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

x

请问下面的命令可以么？如何改变正弦力的相位呢？
/SOLU
ANTYPE,TRANS
TRNOPT,MSUP
ntime=0.0001
kbc,0
DELTIM,1e-4
!OUTPR,none
OUTRES,nsol,1
alphad,2.9768-3
betad,1.01897e-4
/PREP7
CSYS,1
NROTAT,all
FINISH
/SOL
LSWRITE
*do,i,1,50
myf=-1100*sin(2*3.14*22.751*i)
f,650,fx,myf
myf=-10*sin(2*3.14*36.387*i)
f,652,fx,myf
myf=-1*sin(2*3.14*65.375*i)
f,423,fx,myf1
myf=-2000*sin(2*3.14*111.53*i)
f,425,fx,myf
myf=-2*sin(2*3.14*140.31*i)
f,428,fx,myf
kbc,1
SOLVE
time,1
f,652,fx,0,1
f,650,fx,0,1
f,420,fx,0,1
f,423,fx,0,1
f,425,fx,0,1
f,428,fx,0,1
!lswrite
SOLVE
*enddo

jxxansys · 发表于 2009-7-7 19:23

你现在施加的载荷在[0，1]内为正弦规律变化，但time<0以及time>1时为零或其它定值，所以是非周期函数，谈不上频率、相位等。
是不是可以用谐响应分析分别分析不同频率下的结果，然后用载荷工况进行叠加

zhangjianzaijiu · 发表于 2009-7-8 09:52

我想的是1秒时加载为零是因为观测时间仅到1秒，否则运算时间会太长，我试一下用谐响应分析试一下

zhangjianzaijiu · 发表于 2009-7-8 11:30

不知道大家谁有在同一个结构中加不同频率的载荷而且初相位也不同的命令流文件？如果有就感激不尽了

zhangjianzaijiu · 发表于 2009-7-12 21:22

不知道下面这种加法对不对？
但是为什么循环次数会影响结果呢？
/SOLU
ANTYPE,TRANS
TRNOPT,MSUP

DELTIM,1e-4
OUTRES,nsol,1
kbc,1

alphad,2.9768-3
betad,1.01897e-4
LSWRITE





F1 = -4286.4125977
F2 = 4491.6909180
F3 = 3357.7712402
F4 = 698.9247437
F5 = -3582.6000977
F6 = -3787.8789062
L1 = 538.0000000
L2 = 588.0000000
L3 = 468.0000000
L4 = 418.0000000
L5 = 620.0000000
L6 = 433.0000000
P1 = -2.2674365
P2 = -1.3457143
P3 = 0.0675930
P4 = -1.2596869
P5 = 3.0662622
P6 = -0.5837573
P3 = 0.7926810
P4 = -0.1167515
P5 = 2.3288846
P6 = -0.7803914




*do,i,0,1,0.005
time,i
myf1=F1*sin(2*3.14*22.751*i+P1)
f,L1,fx,myf1
myf2=F2*sin(2*3.14*36.387*i+P2)
f,L2,fx,myf2
myf3=F3*sin(2*3.14*65.375*i+P3)
f,L3,fx,myf3
myf4=F4*sin(2*3.14*80.541*i+P4)
f,L4,fx,myf4
myf5=F5*sin(2*3.14*111.53*i+P5)
f,L5,fx,myf5
myf6=F6*sin(2*3.14*140.31*i+P6)
f,L6,fx,myf6

!DELTIM,1e-4
!OUTPR,none
OUTRES,nsol,1
kbc,1
SOLVE
*enddo
fini

jxxansys · 发表于 2009-7-12 21:50

!供参考
/CLEAR
/FILNAME, EXAMPLE

/PREP7
ET,1,MASS21
ET,2,COMBIN14
R,1,1
R,2,10000,63
K,1,0,0,0
K,2,1,0,0
LSTR,1,2
LESIZE,1,,,1
LATT, 1, 2, 2
LMESH,1
E,2
FINISH

/SOLU
ANTYPE, HARMIC
D,1,ALL
D,2,UY
D,2,UZ
D,2,ROTX
D,2,ROTY
D,2,ROTZ
HARFRQ,50,50    !频率1
NSUBST,1
KBC,1
F,2,FX,2000 !频率1实部、虚部
SOLVE
HARFRQ,90,90    !频率2
F,2,FX,1000 ,1000  !频率2实部、虚部
SOLVE
HARFRQ,100,100    !频率3
F,2,FX,1000 ,-200  !频率3实部、虚部
SOLVE
SAVE
FINISH

/POST26
NSOL,2,2,U,X,DispX
PLVAR,2
PLCPLX,1
PLVAR,2
FINISH
!然后在POST1中查询各频率对应的结果，即使用手工处理、叠加一下也可以
!用ANSYS怎样叠加本人没有研究，有高手请指教！

[ 本帖最后由 jxxansys 于 2009-7-12 22:12 编辑 ]

zhouxman · 发表于 2009-7-13 09:23

恩高手就是高手受益不少的

happysky_24 · 发表于 2009-7-18 23:01

你好，小弟有点疑问哈，你那个分成几种工况算然后叠加跟同时作用几个方向的力是不能等效的吧，叠加应该是在线性才拥有的特性啊。本人觉得是不是可以参考一下地震波施加三个方向的加速度的情况，如下：
*dim,ACCEXYZ,TABLE,2000,3                   !01行
*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行
(3e16.6)                                     !03行
*vread,ACCEXYZ(1,0),time,txt                !04行
(e16.6)                                     !05行
ACCEXYZ(0,1)=1                               !06行
ACCEXYZ(0,2)=2                               !07行，同上
ACCEXYZ(0,3)=3                               !08行，同上
finish
/SOLU
ANTYPE,trans
btime=0.01             !定义计算起始时间
etime=15.00             !定义计算结束时间
dtime=0.01             !定义计算时间步长
*DO,itime,btime,etime,dtime
time,itime
AUTOTS,0
NSUBST,1, , ,1
KBC,1
acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3)    !施加三个方向的地震加速度
SOLVE
*ENDDO
希望对大家有所启示

jxxansys · 发表于 2009-7-19 00:30

确实叠加应该是在线性才拥有的特性，如果用瞬态分析代替谐响应分析，实际是对周期信号进行了加窗处理，据信号处理理论应该是可以的（本人道听途说），但是必须分析由此所产生的误差是否满足需要

账号		自动登录	找回密码
密码			我要加入

[结构分析] 请问在瞬态分析中如何施加几个频率不同的正弦载荷？相位该如何加进去？

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

评分

回复 6楼 jxxansys 的帖子

评分

评分