（1）施加约束及重力荷载。

通过路径Main Menu＞Solution＞Define Loads＞Apply＞Structural＞Inertia＞Gravity＞Global施加重力ACELY=9.8。

通过路径Main Menu＞Solution＞Define Loads＞Apply＞Structural＞Displacement＞On Areas，拾取模型底面，施加UY=0约束；拾取两侧面，施加UX=0约束；拾取前后两立面，施加UZ=0约束。

具体命令流如下。

/SOLU ！进入求解模块

NSEL，ALL ！选择所有节点

ACEL，0，9.8，0， ！施加重力加速度

ASEL，S，，，21，22 ！选择底面

ASEL，A，，，28，35，7

DA，ALL，UY， ！约束Y方向自由度

ASEL，S，，，10，18，8 ！选择侧面

ASEL，A，，，31，34，3

DA，ALL，UX， ！约束X方向自由度

ASEL，S，，，1，8，7 ！选择地铁站台附近前后面

ASEL，A，，，16，26，10

ASEL，A，，，27，33，6

ASEL，A，，，36，45，9

DA，ALL，UZ， ！约束Z方向自由度

FINISH ！退出求解模块

完成约束后，得到模型如图19-40所示。

（2）读入加速度时程。

读入加速度时程前，要先将地震波时程文件（参见配书光盘）放到指定文件夹中，以便调用。存放位置可在调用时自行设定，无特别规定。调入命令如下。

ALLSEL

∗DIM，SHUIPING，ARRAY，2，1000 ！定义数据矩阵

∗CREATE，ANSUITMP ！生成文件

∗VREAD，SHUIPING（1，1），SHUIPING，TXT， ！读入时程文件(www.daowen.com)

E：\ANSYS17.0\EXAMPLE\C12\C12.1，1000，

（E9.3，E11.3） ！读入完毕

∗END！输入数据

/INPUT，ANSUITMP

图19-40 施加约束后的模型

（3）设置求解选项。

在GUI界面中选择Main Menu＞Solution＞Analysis Type＞New Analysis命令，在弹出的对话框中设置分析类型为Transient。在GUI界面中选择Main Menu＞Solution＞Analysis Type＞So′ln Controls＞Transient命令，在弹出的对话框中设定Full Transient Options为Stepped loading。

调入地震波加速度时程文件，以文件内数据为每一步的加速度，荷载的时间间隔为0.02秒，共50步。循环加载方式采用APDL的循环语句。具体命令流如下。

/SOLU ！进入求解模块

ANTYPE，4 ！分析类型为Transient

TRNOPT，FULL ！设定分析项

LUMPM，0

∗DO，T，1，50，1 ！循环读入地震波加速度时程文件数据

TIME，0.02∗T ！每一步间隔0.02

KBC，1 ！荷载设置

NSUB，1 ！子步数1

ACEL，SHUIPING（2，T） ！读入加速度

ALPHAD，0.069 ！质量矩阵

BETAD，0.035 ！刚度

SOLVE ！求解

∗ENDDO ！结束循环

FINISH

求解结束。