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

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

请问下面的命令可以么？如何改变正弦力的相位呢？
/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

你现在施加的载荷在内为正弦规律变化，但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

回复 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
希望对大家有所启示

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

确实叠加应该是在线性才拥有的特性，如果用瞬态分析代替谐响应分析，实际是对周期信号进行了加窗处理，据信号处理理论应该是可以的（本人道听途说），但是必须分析由此所产生的误差是否满足需要
页: [1]
查看完整版本: 请问在瞬态分析中如何施加几个频率不同的正弦载荷？相位该如何加进去？