声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4062|回复: 10

[前后处理] [分享]地震波的输入和求解

[复制链接]
发表于 2006-5-16 16:55 | 显示全部楼层 |阅读模式

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

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

x
!对于地震波的输入,可以把荷载记录做成文件,利用apdl的读取功能读入倒数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。有用到的朋友自己作一下把。
  1. fini
  2. /config,nres,1000
  3. *dim,aceX,TABLE,3000,1
  4. *dim,aceY,TABLE,3000,1
  5. *dim,aceZ,TABLE,3000,1

  6. *creat,ff
  7. *vread,aceX(1,1),acex,txt,,1
  8. (e16.6)
  9. *vread,aceX(1,0),ACETT,,,1
  10. (e17.6)
  11. ACEX(0,1)=1
  12. *end
  13. /input,ff

  14. *creat,ff
  15. *vread,aceY(1,1),acey,txt,,1
  16. (e16.6)
  17. *vread,aceY(1,0),ACETT,,,1
  18. (e17.6)
  19. ACEY(0,1)=1
  20. *end
  21. /input,ff

  22. *creat,ff
  23. *vread,aceZ(1,1),acez,txt,,1
  24. (e16.6)
  25. *vread,aceZ(1,0),ACETT,,,1
  26. (e17.6)
  27. ACEZ(0,1)=1
  28. *end
  29. /input,ff  
复制代码
!地震波时程记录分成了3个文件,每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件,这样程序就简单多了。大家可以试看看修改一下。
下面是计算部分语句:
  1. /SOLU

  2. ANTYPE,trans
  3. !求解其自己选了
  4. TM_START=0.01
  5. TM_END=15.00
  6. TM_INCR=0.01
  7. *DO,TM,TM_START,TM_END,TM_INCR
  8. TIME,tm
  9. alpha,
  10. BETAD,
  11. ACEL,acex(tm),acey(tm),acez(tm)
  12. SOLVE
  13. *ENDDO
  14. fini
复制代码

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2006-5-26 09:48 | 显示全部楼层
在地震输入中怎样考虑行波效应呢?
发表于 2006-5-26 12:10 | 显示全部楼层
我正在用ANSYS做抗震分析,有用!
发表于 2006-5-27 06:15 | 显示全部楼层
1)是正弦波的地震波反应
  1. /solu
  2. *dim,baseace1,array,600 ! 定义正弦波的“基础加速度”, 500时间点
  3. *do,i,1,600
  4. baseace1(i)=0.1*sin(i*2) !
  5. *enddo
  6. NSUBST,1, , ,1 !1个子步
  7. OUTRES,ALL,1 !输出每个子步的结果
  8. ANTYPE,TRANS

  9. CSYS,1
  10. NSEL,S,LOC,X,13000,13000
  11. CSYS,0
  12. *do,i,1,600
  13. f,all,fx,3e4*baseace1(i)
  14. f,all,fz,3e4*baseace1(i)
  15. TIME,i*0.02
  16. solve
  17. *enddo
复制代码
2)做z向行波分析
  1. /nerr,,99999999
  2. fini
  3. /solu
  4. *dim,baseace1,array,300 ! 定义正弦波的“基础加速度”, 500时间点
  5. *dim,baseace2,array,300
  6. *do,i,1,300
  7. baseace1(i)=0.1*sin(i*0.2)
  8. baseace2(i)=0.1*sin(i*0.2-31)
  9. *enddo

  10. ANTYPE,TRANS !时程分析
  11. CSYS,1
  12. NSEL,S,LOC,X,13000,13000
  13. CSYS,0
  14. *do,i,1,300
  15. f,all,fx,3e4*baseace1(i)
  16. f,all,fz,3e4*baseace1(i) ! F=m*a
  17. TIME,i*0.02
  18. *enddo

  19. NSUBST,300, , ,1 !500个子步
  20. OUTRES,ALL,1 !输出每个子步的结果
  21. LSWRITE

  22. !在Z方向加行波效应
  23. ANTYPE,TRANS
  24. CSYS,1
  25. NSEL,S,LOC,X,13000,13000
  26. CSYS,0
  27. *do,j,1,300
  28. f,all,fz,3e4*baseace2(j)
  29. f,all,fx,0
  30. TIME,(j+301)*0.02
  31. *enddo

  32. NSUBST,300, , ,1 !500个子步
  33. OUTRES,ALL,1 !输出每个子步的结果
  34. LSWRITE

  35. LSSOLVE,1,2
复制代码
发表于 2006-6-14 18:24 | 显示全部楼层
感谢阿!
发表于 2008-1-13 16:52 | 显示全部楼层
我也要做地震分析 谢谢拉!
发表于 2008-1-15 08:58 | 显示全部楼层
谢谢!
发表于 2008-11-12 18:33 | 显示全部楼层
如何获得位移地震波
发表于 2008-11-12 18:36 | 显示全部楼层
有没有位移地震波啊
发表于 2008-11-14 10:15 | 显示全部楼层
我也正在考虑行波效应,谢谢了!
回复 支持 1 反对 0

使用道具 举报

发表于 2012-4-22 14:07 | 显示全部楼层
谢谢大侠的代码!!!!!!!!!!!!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 11:49 , Processed in 0.069701 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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