rogerkiki 发表于 2009-4-17 04:26

随机振动风载下的UX(高手帮看下)模拟外墙结构的

模拟外墙结构的
但是最终的某点的X方向位移并不符合随机振动的变化,哪位高手帮我看下
另外问下ANSYS可以看整个shell单元的位移吗,新手...请指教
http://tkkljq.bay.livefilestore.com/y1pmpTWjZgsDPnbLZ0zOrYAarNv8MVWYvZtH5MqUVjhRyBkoJBFTONTu4BOP2JkRQ5P4wxHwxvcDOB1uTDdq7rJ3YWGx4R0PFey/UX.bmp

fini
/clear,all
/filname,muqiang
/title,muqiang
/prep7
ET,1,beam188
MP,EX,1,7.0E11
MP,PRXY,1,0.2
MP,DENS,1,2800
sectype,1,beam,rect
secdata,0.15,0.2

ET,5,beam188
MP,EX,5,3.5E5
MP,PRXY,5,0.3
MP,DENS,5,2500
sectype,2,beam,rect
secdata,0.3,0.3

ET,2,combin14,,,2
R,2,1000000

ET,3,combin14,,,0
R,3,10000000

ET,4,SHELL63
R,4,0.012      !定义玻璃厚度
KEYPOT,4,2,1
KEYPOT,4,6,2
MP,EX,4,7.2E11
MP,PRXY,4,0.2
MP,DENS,4,2560      

k,1,0,0,0
k,2,4,0,0
kgen,11,1,2,,,,4
kgen,4,1,22,,4.5
kgen,2,1,88,,,-0.4
/view,1,1,1,g
/angle,1,270,xm,0
/replot

*do,i,1,87,2
l,i,i+1
*enddo

*do,i,1,19,2
l,i,i+2
l,i+1,i+3
*enddo
*do,i,23,41,2
l,i,i+2
l,i+1,i+3
*enddo
*do,i,45,63,2
l,i,i+2
l,i+1,i+3
*enddo
*do,i,67,85,2
l,i,i+2
l,i+1,i+3
*enddo

*do,i,2,66,2
l,i,i+21
*enddo


*do,i,1,19,2
a,i,i+1,i+3,i+2
*enddo
*do,i,23,41,2
a,i,i+1,i+3,i+2
*enddo
*do,i,45,63,2
a,i,i+1,i+3,i+2
*enddo
*do,i,67,85,2
a,i,i+1,i+3,i+2
*enddo



*do,i,89,175,2
l,i,i+1
*enddo

*do,i,89,107,2
l,i,i+2
l,i+1,i+3
*enddo
*do,i,111,129,2
l,i,i+2
l,i+1,i+3
*enddo
*do,i,133,151,2
l,i,i+2
l,i+1,i+3
*enddo
*do,i,155,173,2
l,i,i+2
l,i+1,i+3
*enddo

*do,i,90,154,2
l,i,i+21
*enddo



*do,i,1,88,1
l,i,i+88
*enddo





lsel,s,line,,1,124,1
latt,1,,1,,,,1      
lesize,all,,,2
lmesh,all         
allsel

lsel,s,line,,125,157,1
latt,2,2   
lesize,all,,,1
lmesh,all         
allsel

lsel,s,line,,158,316,1
latt,5,,5,,,,2      
lesize,all,,,1
lmesh,all         
allsel   

lsel,s,line,,317,404,1
latt,3,3   
lesize,all,,,1
lmesh,all         
allsel

allsel
asel,s,,,1,40,1
aatt,4,4,
amesh,all
allsel
nummrg,all
numcmp,all

nsel,s,loc,y,-0.4
d,all,all
allsel

nsel,s,loc,z,0
nsel,r,loc,y,0
d,all,all
allsel

acel,,,-9.8


save
finish



/solu          !加载求解模块
antype,modal   !指定分析类型为模态(modal)分析
modopt,lanb,20!指定使用子block lanczos法计算前20阶模态
mxpand,20   !设置模态扩展数为20
allsel
solve          !开始模态分析求解
save         !保存求解结果
finish         !退出模态分析求解模块      

finish               !完成上一阶段设置
/config,nres,200000   !设置最大子步数为20000
/config,fsplit,750   

/solu
*dim,sidewind1,table,256,10
                              
*tread,sidewind1,c:\sidewind1,txt   !读入荷载数据
antype,trans      !指定分析类型为瞬态(transient)动力分析
trnopt,full   !指定瞬态动力学分析采用full法


*do,t,0.25,64,0.25
time,t               
nsubst,1               
kbc,0            
i=t/0.25
            
      f,1,fx,sidewind1(i,1)/3
      f,133,fx,sidewind1(i,1)/3
      f,4,fx,sidewind1(i,1)/3+sidewind1(i,2)/3
      f,135,fx,sidewind1(i,2)/3
      f,7,fx,sidewind1(i,2)/3+sidewind1(i,3)/3
      f,137,fx,sidewind1(i,3)/3
      f,10,fx,sidewind1(i,3)/3+sidewind1(i,4)/3
      f,139,fx,sidewind1(i,4)/3
      f,13,fx,sidewind1(i,4)/3+sidewind1(i,5)/3

      f,141,fx,sidewind1(i,5)/3
      f,16,fx,sidewind1(i,5)/3+sidewind1(i,6)/3
      f,143,fx,sidewind1(i,6)/3
      f,19,fx,sidewind1(i,6)/3+sidewind1(i,7)/3
      f,145,fx,sidewind1(i,7)/3
      f,22,fx,sidewind1(i,7)/3+sidewind1(i,8)/3
      f,147,fx,sidewind1(i,8)/3
      f,25,fx,sidewind1(i,8)/3+sidewind1(i,9)/3
      f,149,fx,sidewind1(i,9)/3
      f,28,fx,sidewind1(i,9)/3+sidewind1(i,10)/3
      f,139,fx,sidewind1(i,10)/3
      f,13,fx,sidewind1(i,10)/3      
acel,0,0,-9.8               
allsel,all
solve                     
*enddo
outres,all,all            
finish

[ 本帖最后由 rogerkiki 于 2009-4-17 04:47 编辑 ]

rogerkiki 发表于 2009-4-17 17:41

我不加弹簧振动还是会衰减...

rogerkiki 发表于 2009-4-19 22:06

哪位高手帮帮忙

ab7969 发表于 2009-5-25 12:28

能否提供您风载荷的数据啊 急需!!!感谢!!
页: [1]
查看完整版本: 随机振动风载下的UX(高手帮看下)模拟外墙结构的