风载是有经验公式的，比如
笛卡儿坐标系：q=(-0.1033*(Y-2.802)^2+6.3136*(Y-2.802)+269.5626)*ABS(sin(w*TIME))，式中Y(/m)为高度、w(/rad)为激励角速度、q(/Pa)为风压。
地震也需要原始数据，就是加速度列表，包括水平加速度和竖直加速度

一。

地震载荷用惯性载荷的施加方法，风载荷按梯度面载荷的方法施加。采用实体单元。

做塔的详细分析和大型球罐必须考虑这两个载荷，帮助文件里这方面的介绍。

ANSYS地震分析算例
 
 
! 一个四层弹簧－质点模型的地震ANSYS分析
! Example of seismic analysis in ANSYS
! 包括频率分析，谱分析和时程分析
! With model & frequency analysis, spectrum analysis and transient analysis
! 作者：陆新征，清华大学土木工程系
! Author: Lu Xinzheng Dept. Civil Engrg. of Tsinghua University

[Money=20]

FINI
/CLEAR
/UNITS,SI
/PREP7
DAMPRATIO=0.02 !振型阻尼系数都为0.02
ET,1,BEAM188
KEYOPT,1,7,1
KEYOPT,1,8,1
ET,2,MASS21,,,4
MP,EX,1,210E9
MP,NUXY,1,0.3
MP,DENS,1,7800
MPDATA,DAMP,1,,DAMPRATIO
SECTYPE,3,BEAM,I
SECDATA,0.2,0.2,0.5,0.01,0.01,0.008
R,1,1.6E2,1.6E2
R,2,1.2E2,1.2E2
K,1,0,0,0
K,2,0,3,0
K,3,0,6,0
K,4,0,9,0
K,5,0,12,0
K,6,0,0,10000
L,1,2
L,2,3
L,3,4
L,4,5
LSEL,,,,1,4
LATT,1,0,1,,6,,3
ALLS
LESIZE,ALL, , ,3 , , , , ,1
LMESH,ALL
NUMMRG,ALL
NUMCMP,ALL
TYPE,2
REAL,1
E,2
TYPE,2
REAL,2
E,8
E,14
E,20
*SET,NT,1001
*SET,DT,0.02
*DIM,AC,,NT
*VREAD,AC(1),RECORD,TXT
(F8.3)
/SOLU
D,1,ALL,
ANTYPE,0
SOLVE
FINISH
/SOLU !模态分析
ANTYPE,2
MODOPT,SUBSP,8
MXPAND,8, , ,1
SOLVE
FINI
!得到自振频率1
*GET,FREQ1,MODE,1,FREQ
/ESHAPE,1.0
/SOLU
ANTYPE,SPECTR !谱分析
SPOPT,SPRS
MDAMP,1,DAMPRATIO, , , , , ,
SED,1
SVTYPE,2,1
FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667
SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329514922,0.484352764,0.598723486,0.786897371
FREQ,2.222,10,100000
SV,,1.034212766,1.034212766,0.3528
!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45
SRSS,0.01,DISP
SOLVE
FINISH
/POST1
/INP,,MCOM
NSORT,U,SUM
*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形

!开始地震计算
FINISH
/CONFIG,NRES,20000
/SOLU
ANTYPE,TRANS
TRNOPT,FULL
ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926
BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)

*DO,I,1,1001
ACEL,AC(I),0,0
TIME,I*0.02
OUTRES,ALL,ALL
SOLVE
*ENDDO
FINISH

/POST26
FILE
NSOL,2,20,U,X,
PLVAR,2, , , , , , , , , ,
STORE,MERGE
XVAR,1
PLVAR,2,
*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移

数据文件

record.txt

-0.0014
-0.0108
-0.0101
-0.0088
-0.0095
-0.012
-0.0142
-0.0128
-0.011
-0.0085
-0.0085
-0.0131
-0.0176
-0.0194
-0.0162
-0.0144
-0.0108
-0.0082
-0.0042
-0.0066
-0.0131
-0.019
-0.0196
-0.0066
0.003
0.0141
-0.0049
-0.0128
-0.0144
-0.0203
-0.026
-0.0325
-0.0306
-0.0172
-0.0197
-0.0163
-0.0164
-0.0067
0.0025
0.015
0.0236
0.0252
0.0336
0.0463
0.0492
0.0419
0.0359
0.0271
0.0235
0.0339
0.0412
0.053
0.0639
0.0732
0.0652
0.0599
0.04
0.04
0.0063
-0.0515
-0.0787
-0.0603
-0.0484
-0.025
-0.0059
0.0134
0.0308
0.0499
0.071
0.0995
0.1219
0.1529
0.1449
0.1155
0.0935
0.0892
0.0926
0.0839
0.0901
0.0993
0.1209
0.0328
-0.1475
-0.2066
-0.1989
-0.2034
-0.1816
-0.1725
-0.1752
-0.1753
-0.1805
-0.163
-0.1347
-0.1087
-0.0782
-0.0429
-0.0017
0.036
0.0785
0.1164
0.1598
0.196
0.2412
0.2729
0.3036
0.32
0.3417
0.2821
0.2324
-0.1198
-0.2373
-0.164
-0.1865
-0.1095
-0.0753
-0.0173
0.0113
0.0533
0.0895
0.1186
0.1757
0.0576
-0.2631
-0.1547
-0.1729
-0.1012
-0.0579
0.0237
-0.067
-0.198
-0.1641
-0.1685
-0.1481
-0.1231
-0.1001
-0.0751
-0.0523
-0.0271
-0.0044
0.0188
-0.0095
-0.0433
-0.0838
-0.0951
-0.0716
-0.0599
-0.0334
-0.0108
0.0185
0.042
0.0673
-0.0097
-0.0372
-0.004
0.0011
0.0344
0.0565
0.0883
0.113
0.1363
0.0219
0.0241
0.0683
0.0689
0.1318
0.1353
0.204
-0.0931
-0.1308
-0.0692
-0.0546
0.0072
0.0675
-0.1067
-0.1488
-0.1071
-0.1162
-0.0762
-0.0559
-0.0215
-0.0126
-0.0674
-0.0324
-0.0337
-0.0109
0.0017
0.0299
0.0488
0.0608
0.0222
-0.0032
-0.0245
0.0077
0.0211
0.0568
0.0826
0.1206
0.1478
0.1737
0.0421
0.0029
0.0259
0.0293
-0.0055
-0.0147
0.0143
0.0206
0.0499
0.0645
0.0957
0.1128
0.1447
0.1629
0.1945
0.1856
0.1984
0.1769
0.125
-0.1207
-0.0542
-0.0384
-0.0311
-0.1118
-0.1661
-0.2464
-0.2025
-0.1835
-0.1317
-0.096
-0.0325
0.0154
0.0816
0.1319
0.1818
-0.0058
-0.0169
0.0285
0.0447
0.0983
0.1424
0.1853
0.2456
0.1685
-0.138
-0.0999
-0.1089
-0.0907
-0.0469
-0.125
-0.2111
-0.1617
-0.1692
-0.1306
-0.1111
-0.0773
-0.051
-0.0544
-0.12
-0.1209
-0.1158
-0.1145
-0.0717
-0.0546
0.0064
-0.0804
-0.1634
-0.0859
-0.0961
-0.0396
-0.0147
0.0319
0.0648
0.0876
0.0472
0.0198
-0.0027
0.0292
0.0445
0.0785
0.1033
0.1352
0.1606
0.1861
0.1281
0.064
0.0204
0.0314
0.0373
0.0496
0.0235
-0.0084
-0.0168
-0.0113
-0.0229
-0.0248
-0.0157
-0.0069
0.0147
0.0379
0.0579
0.0255
-0.0041
-0.0428
-0.0133
0.0095
0.023
-0.0129
-0.005
0.008
0.021
0.038
0.051
0.0157
-0.0032
-0.0111
0.0005
0.0076
0.0035
-0.0095
-0.0036
-0.0016
0.0038
0.0085
-0.0056
-0.0304
-0.0421
-0.0244
-0.0236
-0.0177
-0.0129
-0.0018
0.0203
-0.0108
-0.0091
-0.0034
-0.0106
-0.0111
-0.0099
-0.0002
0.0073
0.0235
0.0355
0.0705
0.0779
0.0184
-0.0263
-0.0124
-0.0042
0.0159
0.0048
-0.0219
-0.0467
-0.0428
-0.0216
-0.0043
0.0159
0.032
0.0419
0.0123
-0.016
-0.0204
-0.0082
-0.0206
-0.0137
-0.0055
0.0053
0.0134
0.0266
0.0232
0.0079
-0.0008
0.02
0.0435
0.0492
0.0191
0.0092
-0.0022
-0.0021
0.0052
0.0093
0.0255
0.0368
0.0525
0.0541
0.0425
0.0398
0.0559
0.0756
0.0365
0.0411
0.0098
-0.0204
-0.0249
-0.0405
-0.0413
-0.0471
-0.0433
-0.0458
-0.0057
0.0178
-0.0208
-0.0492
-0.053
-0.0362
-0.0405
-0.0308
-0.0316
-0.0265
-0.0265
-0.0269
-0.0345
-0.0309
-0.0217
-0.0078
0.0087
0.0281
0.031
0.0358
0.0341
0.0358
0.0287
0.0305
0.0112
0.0214
0.0136
0.0384
-0.0861
-0.1349
-0.1342
-0.1354
-0.1193
-0.1042
-0.0829
-0.0651
-0.0444
-0.0258
-0.006
-0.0091
-0.0182
-0.0147
0.0085
0.0163
0.005
0.0264
0.0582
0.0867
0.12
0.1695
0.1111
-0.11
-0.0366
-0.0445
-0.0236
-0.096
-0.0656
-0.0597
-0.067
-0.0552
-0.0027
0.0378
0.1072
0.1669
0.0947
0.0408
0.0667
0.0132
-0.0095
-0.052
-0.0827
-0.1152
-0.115
-0.0803
-0.0369
0.0029
0.0545
0.1178
0.161
-0.027
0.0034
-0.0056
0.002
0.0146
0.0537
0.0798
-0.0205
-0.059
-0.0169
-0.0175
-0.0028
0.0074
0.0382
0.0567
0.0753
0.0801
0.0592
0.0304
0.0023
0.0064
-0.0406
-0.0451
-0.0079
0.0168
0.0567
0.0093
-0.0055
0.0044
-0.0123
-0.0282
-0.0437
-0.0352
-0.0255
-0.0111
0.0205
0.0519
0.0854
0.1144
0.0733
0.0237
-0.0368
-0.0271
-0.0217
-0.0873
-0.0973
-0.0589
-0.0336
0.0077
0.0259
0.0508
0.0361
0.0081
-0.0056
-0.0209
-0.0317
-0.0238
-0.0376
-0.055
-0.0722
-0.0803
-0.0523
-0.034
-0.0011
0.0065
-0.0037
-0.0005
-0.0168
-0.041
-0.008
0.0079
0.0374
0.0615
0.0665
0.0254
-0.0057
-0.0474
-0.0356
-0.0243
-0.0048
0.0126
0.0379
0.0241
-0.0227
-0.0428
-0.0679
-0.0661
-0.059
-0.0513
-0.0408
-0.0309
-0.0266
-0.0541
-0.0628
-0.0908
-0.1107
-0.0881
-0.077
-0.0582
-0.0473
-0.0333
-0.0199
0.002
0.0211
0.0432
0.0613
0.0767
0.0933
0.1066
0.113
0.1187
0.1247
0.1334
0.1594
0.1797
0.2037
0.1236
0.0442
-0.014
-0.0666
-0.0555
-0.0693
-0.0984
-0.1246
-0.1179
-0.105
-0.092
-0.0743
-0.0809
-0.085
-0.086
-0.0863
-0.0873
-0.0868
-0.0885
-0.0537
0.0052
0.0215
0.0245
0.058
0.0314
0.0236
0.0485
0.0589
0.0525
0.0355
0.0197
0.0199
0.0492
0.0343
0.0288
0.0432
0.0239
0.0088
0.0077
-0.0148
-0.0077
-0.0019
0.0075
0.0044
-0.0145
-0.0316
-0.0241
-0.0028
0.0182
0.0426
0.0439
0.0512
0.0466
0.0479
0.0193
0.0222
0.0274
0.0393
0.0504
0.0577
0.0588
0.0822
0.0797
0.0949
0.0345
0.0045
-0.0123
-0.0347
-0.0426
-0.0416
-0.0275
-0.027
0.0074
0.0428
-0.0231
-0.0387
-0.0083
0.0139
0.0445
0.0027
-0.0697
-0.0796
-0.0251
-0.0135
0.0079
-0.0115
-0.0251
-0.0333
-0.0269
-0.0301
-0.02
-0.0067
-0.0038
0.0105
0.0296
0.0344
0.0957
0.0898
0.0179
-0.0362
-0.0994
-0.0807
-0.0744
-0.0539
-0.033
-0.0128
0.0031
0.0148
0.0508
-0.0022
-0.0489
-0.0358
-0.0691
-0.0516
-0.0371
0.0088
0.0632
0.0841
0.1276
0.1388
0.1193
0.0751
0.0225
-0.0088
-0.0227
0.0074
0.0181
0.0544
0.0399
0.0045
-0.0082
-0.0185
-0.002
0.0006
-0.0117
-0.021
-0.0303
-0.0512
-0.0727
-0.0579
-0.0266
-0.0178
0.004
0.0098
0.0137
0.0221
0.0437
0.0091
-0.0548
-0.0555
-0.0243
-0.0081
0.025
0.041
0.0182
-0.0027
-0.0243
-0.0015
0.0247
0.0482
0.0783
0.0622
0.0331
-0.0014
-0.0195
-0.0247
-0.0212
-0.011
0.005
0.0241
-0.0034
-0.0216
-0.0471
-0.0363
-0.0195
-0.0018
0.017
-0.008
0.0005
0.023
0.0374
0.0601
0.0516
0.0432
0.0344
0.0505
0.0653
0.0683
0.0172
-0.017
-0.0527
-0.0664
-0.0387
-0.0222
-0.0033
0.0119
-0.0128
-0.0351
-0.0514
-0.0335
-0.0218
-0.0012
0.0142
0.007
-0.0063
-0.012
-0.0322
-0.0346
-0.0091
0.0073
0.0309
0.0472
0.0603
0.0576
0.033
-0.0073
-0.0777
-0.0608
-0.0438
-0.0209
0.0031
0.035
0.0293
0.0121
0.0338
0.0317
0.0254
0.0206
0.0198
0.0174
0.0021
-0.0144
-0.0343
-0.0339
-0.0145
-0.0028
0.017
-0.0096
-0.0255
-0.0279
-0.0388
-0.0242
-0.0215
-0.0182
-0.0174
-0.0038
-0.0027
-0.0185
-0.0123
0.0087
0.0343
0.0695
0.091
0.0853
0.076
0.0513
0.0186
0.0015
-0.019
-0.0151
-0.0073
0.0021
0.0129
0.0215
0.0024
-0.0124
-0.0329
-0.0519
-0.0708
-0.0579
-0.0462
-0.0307
-0.0145
-0.0009
-0.018
-0.0318
-0.0465
-0.0391
-0.0345
-0.0316
-0.0435
-0.0491
-0.0475
-0.042
-0.0361
-0.0277
-0.0258
-0.0139
-0.0068
0.0507
0.0722
0.0878
0.0782
0.0765
0.0439
0.008
0.0013
-0.0126
-0.0015
0.003
0.0104
0.0104
0.0193
0.0205
0.0074
-0.0056
-0.0072
0.007
0.0106
0.0147
-0.0009
-0.0159
-0.0187
-0.0007
0.0155
0.0105
-0.0115
-0.0302
-0.0309
-0.0095
-0.0058
0.0004
0.002
0.005
0.0057
0.0097
0.0134
0.0177
0.0218
0.0261
0.0302
0.0346
0.0386
0.0474
0.0393
0.0238
0.0115
-0.0079
-0.0124
0.0054
0.0027
-0.025
-0.0566
-0.063
-0.0591
-0.0413
-0.0068
0.0272
0.0277
-0.0021
-0.006
-0.011
-0.0221
-0.0416
-0.0519
-0.0222
0.003
0.0079
0.0139
0.0171
0.0253
0.0323
0.0391
0.0164
-0.0136
-0.0323
-0.0291
-0.0287
-0.0304
-0.0339
-0.0245
-0.0076
0.0125
0.0376
0.0402
0.0245
0.0156
-0.004
-0.0153
-0.0289
-0.0316
-0.0111
0.0094
0.0335
0.0576
0.0424
0.0143
-0.0007
-0.0135
-0.027
-0.0341
-0.0357
-0.0396
-0.0402
-0.0488
-0.048
-0.0406
-0.0407
-0.0351
-0.0187
-0.0057
0.0044
-0.0019
-0.0072
-0.0169
-0.0115
0.0126

一三层 钢框架的地震反映谱分析，请大家笑纳！！

CODE:[Copy to clipboard]!钢框架的振型分解反应谱分析
/prep7
et,1,beam188
mptemp,1
mpdata,ex,1,,2.06e11
mpdata,dens,1,,7850
mpdata,prxy,1,,0.2
mpdata,ex,2,,25.5e15
mpdata,dens,2,,1088e2
mpdata,prxy,2,,0.2
sectype,1,beam,hrec,colum
secoffset,cent
secdata,0.5,0.5,0.02,0.02,0.02,0.02
sectype,2,beam,i,beam
secoffst,user,0,0.5
secdata,0.4,0.4,0.5,0.015,0.015,0.015
k,1,-4.5
k,2,-4.5,4
k,3,-4.5,8.0
k,4,-4.5,12 
k,5,4.5,12 
k,6,4.5,8  
k,7,4.5,4
k,8,4.5

l,1,2
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,7,8
l,3,6
l,2,7
l,7,8
nummrg,all
numcmp,all
lesize,all,,,20
lsel,s,loc,x,-4.5
lsel,a,loc,x,4.5
k,15,0,15
latt,1,1,1,,15,,1
lmesh,all,all
allsel,all
lsel,s,loc,y,4
lsel,a,loc,y,8
lsel,a,loc,y,12
latt,2,1,1,,15,,2
lmesh,all,all
allsel
nsel,s,loc,y,0
d,all,all
allsel,all
/eshape,1
/replot
/solu
antype,modal
modopt,lanb,10,0,0,,off
mxpand,10,,,yes
acel,,9.8
solve
finish
/solu
antype,spectr
spopt,sprs,,yes
dmprat,0.02
svtyp,2,1
sed,1
!!ii类建筑场地，6度设防，多遇地震（近震）
!曲线的下降段衰减系数r=0.95
!amax=0.11
!阻尼比=0.02
!特征周期Tg=0.34s
!直线段的下降段斜率调整系数r1=0.02375
!阻尼调整系数r2=1.319
!定义反应谱曲线
freq,0.420168067,0.452488688,0.49016078,0.534759358,0.588235294,0.653594771,0.735294118,0.840336134,0.980392157
freq,1.176470588,1.470588253,1.960784314,2.941176471,10,1000000
sv,0.02,0.296731616,0.301172866,0.305614116,0.310055366,0.314496616,0.347604662,0.388759029,0.441339541,0.510942806
sv,0.02,0.607567413,0.751032939,0.987076437,1.451,1.451,0.495
solve
finish
/solu
antype,spectr
srss,0.15,disp
solve
finish
/post1
allsel,all
sumtype,prin
/input,,mcom
plnsol,u,sum
etable,mx,smisc,4,17
etable,my,smisc,2,15
etable,bk,smisc,28,30
pletab,mx
pletab,my
pletab,bk
lcwrite,1,'spectrum','l01',''
finish
/solu
antype,static,

solve
fini
/post1
set,1
lcdef,1,1,1
lcfile,2,'spectrum','l01',''
lcfact,2,1.3
sumtype,prin
lcase,1
lcoper,add,2,,2
esel,u,type,,1
plnsol,u,sum

ansys之——地震波的输入和求解                                      

!对于地震波的输入，可以把荷载记录做成文件，利用apdl的读取功能读入倒数据库中。下面的例子是自

己编的一个小文件。修改一下可以更简洁。有用到的朋友自己作一下把。
fini
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*creat,ff
*vread,aceX(1,1),acex,txt,,1
(e16.6)
*vread,aceX(1,0),ACETT,,,1
(e17.6)
ACEX(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceY(1,1),acey,txt,,1
(e16.6)
*vread,aceY(1,0),ACETT,,,1
(e17.6)
ACEY(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceZ(1,1),acez,txt,,1
(e16.6)
*vread,aceZ(1,0),ACETT,,,1
(e17.6)
ACEZ(0,1)=1
*end
/input,ff
！地震波时程记录分成了3个文件，每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件，这样程序就简单多了。大家可以试看看修改一下。
下面是计算部分语句：
/SOLU
ANTYPE,trans
！求解其自己选了
TM_START=0.01
TM_END=15.00
TM_INCR=0.01
*DO,TM,TM_START,TM_END,TM_INCR
TIME,tm
alpha,
BETAD,
ACEL,acex(tm),acey(tm),acez(tm)
SOLVE
*ENDDO
fini

ansys中地震载荷的施加

在许多有限元分析软件中，对于施加地震载荷，多是考虑采用施加体加速度的方法，也就是施加一个跟激励方向想法的惯性力，在实际的地震中是否是这样可能很难说，因为地震是通过地面的运动传递到结构上去的，结构上肯定有一个滞后效应。在纯固体的计算中可能问题还不是很大。
可是在流固耦合结构的计算中，比如大型油罐的抗震计算中，如果给油也施加一个惯性激励力，那就与实际地震差别很大了，毕竟油是由于油罐的运动带动的，而直接加在油上的惯性力会给油罐较大的冲击作用。
  现在有点越来越迷惑了，不知道直接施加惯性力是否正确

目前在做一座桥梁结构的反应谱分析，学习到一些也有一些问题：
1、一般结构用单点响应谱分析，个人感觉单点响应谱分析与单向反应谱作用不是一个概念，感觉单点响应谱分析应该与一致激励相一直。
2、ansys里根据sed定义反应谱方向，但是如何同时假如两个方向的反应谱呢，我先后定义sed,1,0,0和sed，0，0，1,并先后定义两个反应谱的曲线，但是最后solve和模态组合以后，得到的只是后一个反应谱作用下的结构响应，那又如何来同时定义两个方向的呢?
3、隔震结构的响应谱曲线的确定，叠层橡胶隔震规程里根据建筑抗震设计规范并对某些系数做了调整，但是总体还是跟建筑设计规范基本一致的，但是隔震结构的阻尼比明显比非隔震结构的阻尼比大很多，但是这种阻尼比的增加并不是应为结构整体的阻尼比的增加，仅仅是隔震层的阻尼比造成了结构在某个方向的阻尼比的增加，那么也应用规范上的那个反应谱曲线合适吗，我计算出来的结构隔震与非隔震结构的反应谱曲线差一倍左右，感觉不是很对？

 对于问题2,响应谱分析是线性的,两个方向可先单独求解,然后组合结果.

二。