《OpenSEES学习笔记(优选)》由会员分享，可在线阅读，更多相关《OpenSEES学习笔记(优选)（16页珍藏版）》请在金锄头文库上搜索。
1、最新文件仅供参考已改成word文本 方便更改1、利用零长单元模拟阻尼，uniaxialMaterial Elastic 1 6.8098e6;uniaxialMaterial Viscous 2 3.24e5 1;uniaxialMaterial Parallel 3 3 5;element zeroLength 1 $iNode $jNode -mat 3 -dir 1;通常有两种方式：(1) truss element and viscous material.(桁架单元和阻尼材料)(2) force-based beam-column element and Maxwell materi
2、al (基于力的梁柱单元和 Maxwell 材料)。、如何运行 OpenSEES有三种方法可以执行 OpenSees/Tcl 命令：1、interactive 交互式直接将命令输入 Prompt。2、执行文件输入这种方法是最常用的一种，以 source inputfile.tcl 方式执行已写好的外部命令文件。3、Batch 模式即以 Opensees inputFile.tcl 方式在 MS DOS/Unix promt 中运行。二、定义单位和常数在编写一个较大的 Opensees 命令时。最好先定义好单位及常数。在 Opensees 中，编译器 不能自行转换单位。所以一开始就要先定义好。单
3、位定义包括两部分：首先定义基本单位；再定义合成单位。其中基本单位要相互独立。 同时，在定义单位时，既可以按国际公制单位，也可以按私制单位。因些在单位定义文件 中可能是混合的。我个人建议，还是采用国际公制单位较好。像国外常用英制单位。很不 习惯。对于一些常数，如 和 g 等常数要事先定义好。在定义这些单位时所用的命令是“set”。1 / 15# - U nits&C on stants .tc Iset in 1.；# define basic unitsset sec 1.;set kip 1.;set ksi expr $kip/pow($in,2)； set psi expr $ksi 1
4、000.;# define dependent units三、set ft expr 12/$in;set lb expr $kip. 1000;set pcf expr $lbpow($ft,3);expr $kippow($in,2);expr $ksi1000.;set cm expr $in/2.54; set meter expr 100.*$cm;set MPa expr 145*$psi；# define metric unitsset PI expr 2*asin(1.0);生成 Matla命令($sec,2)；# define constantsset U 1.e10；# a
5、 really large number是后处理最常用的工具，通过Tcil脚本语言可以得到Matlab命令文件。同时保证相Matlab同的分析参数。如下例:尸 # script to generate .m file to be read by matlab#MatlabOutput.tcl-set Xframe 1 ；# this parameter would be passed inset fDir Data/；file mkdir $fDir;# create directoryset outFilelD open $fDir DataFrame$Xframe.m 训;# Open o
6、utput file tor writingputs $outFilelD ,Xframe($Xframe)=:$Xframe;M;# trame IDputs$outFilelD uHcol($Xframe) = $Hcol;,R;# column diameterputs $outFilelD uLcol($Xframe) = $Lcol;,R;# column lengthputs $outFilelD uLbeam($Xframe)=$Lbeam；# beam lengthputs $outFilelD ,rHbeam($Xframe):= $Hbeam；R,；# beam depth
7、puts $outFilelD ,rBbeam($Xframe):=$Bbeam；R,；# beam widthputs $outFilelD Weight($Xframe)=$Weight;,r;# superstructure weightclose $outFilelD四、定义 Tcl 命令的方法这种方法是从 Tcl 语言获得的工具，他是一种广义上的函数或者子程序协议。 这个不是很 明白五、读取外部文件外部文件可能是 Tcl 命令或者是数据。有以下三种方法：（1）普通读取文件法；（2）重复计算读取文件法；（3）外部数据文件读 取法。（1）普通读取文件法：source.inputfile.
8、tcl。当外部文件中包括一系列的常用分析命令时， 可用这个方法。（2）重复计算读取文件法：外部文件可能包括一系列重复的计算，则采用以下方法。随时 改变参数。set Hcolumn 66; source analysis.tel set Hcolumn 78;合oiiixE白n自ly各i占.tul也是一种交互式的输入(3) 外部数据文件读取法：这个命令可以用于打开一个数据文件。并按行读取，指定每行 的值为单变量。如果一行中有多个值时，则变量是列阵。而其中每个元素可用“ lindex”命 令提取。# -ReadData.tclif catch open $inFilename r inFilelD
9、 ; # Open the input file and check for errorputs stderr Cannot open $inFilename for reading;# output error statement else foreach line split read SinFilelDn ;# Look at each line in the fileif llength $line = 0 ; # Blank line - do nothingcontinue; else 六、建立模型这项是我们做分析时很关键的一个前提。常规流程包括以下几个方面：（1）定义变量和参数；
10、（2）建立模型和定义节点；（3）建立模型和使用变量定义节点 ；（4）定义材料；（5）定义单元。（1）定义变量和参数：在 Tcl 脚本语言中变量是用来替代数字的。一旦被定义好，则可用 来代替数字。当重新要使用这些数字的时候，只需要用“ $”字符则可以调用。如果不用此 字符。则变量名将会被视为字符串命令将导致错误结果。在定义模型的变量时，通常首先 定义好材料参数，然后单元参数，再是重力（质量）参数，最后是分析参数。 设定好变量 定义的命令文件有很多好处，当需要对材料或尺寸等变量进行修改时，只需在这个文件中 修改即可。无需修改后面的文件。在材料参数中，主要涉及到材料的本构关系变量；单元参数中，主要涉
11、及到单元几何尺寸 及数量；重力参数中，主要是模型质量分布；分析参数中，主要是增量步长及数值计算参 数。（2）建立模型和定义节点：这些模型节点通常是置于输入文件的开头。最后的定义结点质 量的命令会覆盖先前定义的节点质量。在定义这些节点时也可以使用变量来定义，效果是 一样的。（4）定义材料：以下的例子是用于定义钢筋混凝土结构的材料。假定变量都已定义。当这 些命令编制在一个外部文件中时，可以用source命令执行。# MaterialsRC.tclset Concrete Mate rialType inelastic； # options: elastic, inelasticset SteelM
12、aterialType ,rhysteretie；# options: elasticilinearhysteretic# CONCRETEif ConcreteMaterialType =elastic uniaxialMaterial Elastic $IDcone $EcuniaxialMaterial Elastic $IDcover $Ecif $ConcreteMaterialType = inelastic uniaxialMaterial ConcreteOl $IDcore $fc1C $eps1C $fc2C $eps2C; # Core concrete uniaxial
13、Material ConcreteOl $IDcover $fc1U $eps1U $fc2U $eps2U;# Cover concrete# STEELif $SteelMaterialType = elastic uniaxialMaterial Elastic $IDsteel $Esif $SteelMaterialType = bilinear uniaxialMaterial SteelOl $IDsteel $Fy $Es $Bsif $SteelMaterialType = hysteretic uniaxialMaterial Hysteretic $IDsteel $Fy
14、 $epsY $Fy1 $epsY1 $Fu $epsll -$Fy -$epsY -$Fy1 -$epsY1 2555）定义单元：ELE 師 ENTS.tcb# COLUMNSset ColumnType inelastic; set np 5;# source RCcircSection.tcl;# options: rigid elastic inelastic# number ot integration points# proc to define circular fiber section for flexural characteristicsRCcircSection $ID
15、colFlex $riCol $roCol $cover$IDcore $IDcover$IDsteel $NbCol $AbCol $nfCoreR $nfCoreT $nfCoverR $nfCoverT# Define torsional stiffnesssection Aggregator $IDcolSec $IDcolTors T-section $IDcolFlex; # attach torsion and tlexuregeomTransf Linear $IDcolTrans if $ColumnType= rigid set $lyCol expr$lyCol*$lyCol; set $lzCol expr $lzCor$lzCol; element elasticBeamColumn 10 0 1;1 3
