声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1831|回复: 8

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

[复制链接]
发表于 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
回复
分享到:

使用道具 举报

发表于 2009-7-7 19:23 | 显示全部楼层
你现在施加的载荷在[0,1]内为正弦规律变化,但time<0以及time>1时为零或其它定值,所以是非周期函数,谈不上频率、相位等。
是不是可以用谐响应分析分别分析不同频率下的结果,然后用载荷工况进行叠加

评分

1

查看全部评分

 楼主| 发表于 2009-7-8 09:52 | 显示全部楼层
我想的是1秒时加载为零是因为观测时间仅到1秒,否则运算时间会太长,我试一下用谐响应分析试一下
 楼主| 发表于 2009-7-8 11:30 | 显示全部楼层
不知道大家谁有在同一个结构中加不同频率的载荷而且初相位也不同的命令流文件?如果有就感激不尽了
 楼主| 发表于 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
发表于 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 编辑 ]
发表于 2009-7-13 09:23 | 显示全部楼层
恩 高手就是高手 受益不少 的
发表于 2009-7-18 23:01 | 显示全部楼层

回复 6楼 jxxansys 的帖子

你好,小弟有点疑问哈,你那个分成几种工况算然后叠加跟同时作用几个方向的力是不能等效的吧,叠加应该是在线性才拥有的特性啊。本人觉得是不是可以参考一下地震波施加三个方向的加速度的情况,如下:
*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
希望对大家有所启示

评分

1

查看全部评分

发表于 2009-7-19 00:30 | 显示全部楼层
确实叠加应该是在线性才拥有的特性,如果用瞬态分析代替谐响应分析,实际是对周期信号进行了加窗处理,据信号处理理论应该是可以的(本人道听途说),但是必须分析由此所产生的误差是否满足需要

评分

1

查看全部评分

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

本版积分规则

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

GMT+8, 2024-11-14 19:22 , Processed in 0.074064 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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