(word完整版)opensees总结2,推荐文档
1、瑞利阻尼
在OPENSEES 中，结构采用瑞利（Rayleigh ）阻尼，即阻尼矩阵的大小与结构的质量矩阵，刚度矩阵都相关，瑞利阻尼的计算公式如下，阻尼与刚度质量的关系如下图所示。
[][][]01c a m a k =+
式中，ξ为阻尼比，a 0为质量相关系数，a 1为刚度相关系数，[c]为阻尼矩阵，[m]为质量矩阵，[k]为刚度矩阵；ωm 、ωn 为结构两个主振型的圆频率，由于 OPENSEES 能够直接求解振型的特征值，那么特征值与圆频率的关系：=λω。 命令流的解读如下：
set xDamp 0.05 ;————设置阻尼比为0.05
set nEigenI 1;————主振型1为第1振型
set nEigenJ 2;————主振型2为第2振型
set lambdaN [eigen [expr $nEigenJ]];————求解两阶振型即可
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];————提取第1阶特征值 set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];————提取第2阶特征值 set omegaI [expr pow ($lambdaI,0.5)];——— —从特征值求圆频率
set omegaJ [expr pow($lambdaJ,0.5)];————从特征值求圆频率
set alphaM [expr $xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
————alphaM 为a 0，即质量相关系数；
set betaKcurr [expr 2.*$xDam p/($omegaI+$omegaJ)]; ———betaKcurr 为a1，即刚度相关系数；
rayleigh $alphaM $betaKcurr 0 0———定义瑞利阻尼，只需要填写a0、a1，其它值为0。
2、OPENSEES 质量源
OPENSEES基本上采用结点质量源的形式，与大部分有限元程序一样，也就是每个结点有6 个自由度，6 个自由度上都有广义质量，平动方向UX，UY，UZ 称为质量，而转动方向RX，RY，RZ称为转动惯量（除了刚体，基本上结点没有转动惯量）。普通结点只有UX，UY，UZ的平动质量，且三个值是相等的。puts "mass"
mass $NUM $MUX $MUY $MUZ $MRX $MRY $MRZ
其中，$NUM为结点编号，$MUX $MUY $MUZ $MRX $MRY $MRZ代表各个自由度的质量，以单位制（N，mm）的规定，质量的单位应该为ton（吨）。普通结点的质量定义如下：
mass $NUM $M $M $M 0 0 0
3、模态分析
模态分析的命令流与普通静力分析的命令流最大的区别在于记录与分析设置。记录命令如下：
puts "recorder"
recorder Node -file eigen1_node0.out -time -nodeRange 1 28 -dof 1 2 3 "eigen 1"
该命令流用于输出振型位移，即振型形状。
其中，eigen1_node0.out 为输出振型位移的文件名，-dof 1 2 3 代输出的自由度；"eigen1"代表输出的为“第1 振型”。
模态分析设置的命令流：
set numModes 12
set lambda [eigen $numModes]
set period "Periods.txt"
set Periods [open $period "w"]
puts $Periods " $lambda"
close $Periods
record
其中，set lambda [eigen $numModes] ，代表计算n 阶振型，将特征值计算结果存为lambda 数组；
set period "Periods.txt" ，定义输出文本文件的名字为“Periods.txt ”；用于存计算的特征值。
set Periods [open $period "w"]，代表打开文本文件进行记录。
puts $Periods " $lambda"，代表记录特征值数据至文本当中。
close $Periods，代表记录完成，关闭文本文件。
record 代表记录命令。
4、单向地震波数据设置的命令流如下：
s et IDloadTag 1001; ————地震波工况号为1001
set iGMfile "GM1X.txt"; ————地震波数据文件名为GM1X.txt，需要与TCL 文件放在同一个目录下。
set iGMdirection "1"; ————地震波方向为x 方向，即系统1 方向
set iGMfact "0.138"; ————地震波峰值放大系数为0.138mm/s2
set dt 0.02; ————地震波时间间隔为0.02s
foreach GMdirection $iGMdirection GMfile$iGMfile GMfact $iGMfact {
incr IDloadTag;
set GMfatt [expr 1*$GMfact];
set AccelSeries "Series -dt $dt -file Path $iGMfile -factor $GMfatt";
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries;
} ————以上为多维地震波的输入标准样式，可用于三向地震波，也可以用于单向，数组如果只有一个值，如算例所示，只作为单向波。
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries;
为荷载工况的命令流，用于输入地震波数据，$IDloadTag 为地震波工况号，$GMdirection 为地震波的方向，-accel $AccelSeries，为地震波的其它参数，包括文件名，时间间隔，峰值放大倍数等。
5、弹性时程分析命令流如下：
constraints Transformation;————用于动力时程分析
numberer Plain;————普通的排位方法
system UmfPack;————采用Umfpack 自由度排列
test EnergyIncr 1.0e-4 200;————采用能量收敛准则
algorithm Newton————采用牛顿迭代法
integrator Newmark 0.5 0.25————采用Newmark法对时间进行离散analysis Transient————采用时程分析
analyze 1000 0.02————分析步数为1000，时间间隔为0.02s，即分析20s。