有关地震波的输入问题

tangaoming

我在输入地震波时，有一个命令始终不明白，我输入地震波的明尼命令如下

1. *do,t,1,1000,10
   
2. time,0.1*t
   
3. kbc,0
   
4. nsubst,10
   
5. acel,dizhenbo(1,t)
   
6. nsel,all
   
7. solve
   
8. *enddo

复制代码

"nsel,all"这个命令是什么意思，它表示地震波是怎么输入进去的

yanyujing

nsel,all
应该是选择所有的节点

yanyujing

!借助于ANSYS动力分析基本功能和APDL语言以及有关的地震资料编写的SEISMIC程序

1. NT=100
   
2. DT=0.02
   
3. *DIM,ac,,NT
   
4. /INPUT,FF,txt
   
5. /SOLU
   
6. NSUBST,1,,,1
   
7. OUTRES,ALL,1
   
8. ANTYPE,TRANS
   
9. *DO,I,1.NT
   
10. ACEL,0.01*ac(I)*1.3,0,0 !ac(I)的值由地震资料获得
    
11. TIME,I*DT
    
12. SOLVE
    
13. *ENDDO

复制代码

AaronSpark

NSEL，Type，Item，Comp，VMIN，VMAX，VINC，KABS（节点的选择）

选择方式Type=S（选择一组节点为Active的节点）、R（在现有的Active节点中，再选择某些节点为Active节点）、A（再选择某些节点，加入现有Active节点中）、U（在现有Active节点中，排除某些节点）、ALL（全部）；

资料卷标Item=NODE（用节点号码选取）、LOC（用节点坐标选取）；
Comp=无（Item=NODE）、X（Y、Z）（表示节点X（Y、Z）坐标为准，当Item=LOC）。

tangaoming

那它是选择所有节点，难道地震波是对所有节点进行输入的吗？
不是很明白？

AaronSpark

对于地震波的输入，可以把荷载记录做成文件，利用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. 
   
7. *creat,ff
   
8. *vread,aceX(1,1),acex,txt,,1
   
9. (e16.6)
   
10. *vread,aceX(1,0),ACETT,,,1
    
11. (e17.6)
    
12. ACEX(0,1)=1
    
13. *end
    
14. /input,ff
    
15. 
    
16. *creat,ff
    
17. *vread,aceY(1,1),acey,txt,,1
    
18. (e16.6)
    
19. *vread,aceY(1,0),ACETT,,,1
    
20. (e17.6)
    
21. ACEY(0,1)=1
    
22. *end
    
23. /input,ff
    
24. 
    
25. *creat,ff
    
26. *vread,aceZ(1,1),acez,txt,,1
    
27. (e16.6)
    
28. *vread,aceZ(1,0),ACETT,,,1
    
29. (e17.6)
    
30. ACEZ(0,1)=1
    
31. *end
    
32. /input,ff

复制代码

！地震波时程记录分成了3个文件，每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件，这样程序就简单多了。大家可以试看看修改一下。
下面是计算部分语句：

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

复制代码

AaronSpark

nsel并不是说加在所有的点上，而是为了求解而选择所有的点

jasmine927

地震波不是输入到各节点上么？那么如何考虑行波效应？

dreamskyweb

7楼讲的正确！！！