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 materia（l 基于力的梁柱单元和Maxwell 材料）。-、如何运行OpenSEES有三种方法可以执行0 pe nSees/Tc l命令：1、interactive 交互式直接将命令输入Prompt。2、执行文件输入这种方法是最常用的一种，以source in putfile.tc l方式执行已写好的外部命令文件。3、Batch 模式即以 Opensees inputFile.tcl 方式在 MS - DOS/Unix promt 中运行。二、定义单位和常数在编写一个较大的Opensees命令时。最好先定义好单位及常数。在Opensees中，编译 器不能自行转换单位。所以一开始就要先定义好。单位定义包括两部分：首先定义基本单位；再定义合成单位。其中基本单位要相互独立。同 时，在定义单位时，既可以按国际公制单位，也可以按私制单位。因些在单位定义文件中可能是混合的。我个人建议，还是采用国际公制单位较好。像国外常用英制单位。很不习惯。对于一些常数，如 和 g 等常数要事先定义好。在定义这些单位时所用的命令是“set”# - 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 1000.;# define dependent unitsset 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 units# define constantsMatlabset U 1.e10；# a really large number是后处理最常用的工具，通过脚本语言可以得到Matlab命令文件。同时保证相同的分析参数。如下例： # script 6 genefsie .mse ff be Feaid bx ma萤 b#Mar-abourpul.lc- sei xfame； # rh parameter wau-d be passem in se=DWa:=e mkdwsmm# creaCDdirecoryoLHFi-e-D【open sfDiwDsiaFrarnesxframemlsoumi-e-D soumi-e-D soumi-e-D soumi-e-D soumi-e-D solHFi-e-D solHFi-e-Dxhame(sxf 畐 me) u SXhame 飞 =HCO_ 筈 rame7SHCo:= = LCssx=amrt7SLCO_.= :Lbeam(sxfame) H SLbeam:? = Hbmamxham H SHbeam:= Ebeam(sxfame) n SBbemmrz :weigssxf恳 me) n sw.g 三飞cowe SQLHFi-e-DEk# Open clsPLJr f=eor ws-5g# 孟me -D# coumn disuneQlr甘 coumn 一 engih# beam -engih# beam depih# beam waih- 韩 supepswuceFe weighrm羞善凹常Tc_血左-H川雷讯二一)哺聲貝畫二2)刪mm棘貝畫二3) 5前|善棘 sfff。(一)msffi棘sw亠dim-dource.inpurfi_e.-tc_。味孚茨W弃召曲庶洲旦吉糠油 滸部33、 凹油滋-邑讯。(2 ) -孚茨总弃凹恶曲庶洲旦刪wsi+w、亘采油F-Hbfff。鬧母 埒艰瞬磬。sei Hcoumn s- souixe Hna-siwn-sei Hco-umn 7 souixetuna-ywiifllc-也是一种交互式的输入（3）外部数据文件读取法：这个命令可以用于打开一个数据文件。并按行读取，指定每行的值为单变量。如果一行中有多个值时，则变量是列阵。而其中每个元素可用lindex”命 令提取。# -ReadData.tclif catch open $inFilename r inFilelD ; # Open the input tile and check for errorputs stderr ,rCannot open $inFilename for reading;# output error statement else foreach line split read 帛inFilelDn ;# Look at each line in the fileif llength $line = 0 ; # Blank line - do nothingcontin ue； else set Xvalues $line;# execute operation on read dataclose $inFilelD; ;# Close the input tile六、建立模型这项是我们做分析时很关键的一个前提。 常规流程包括以下几个方面：（1）定义变量和参数；（2）建立模型和定义节点；（3）建立模型和使用变量定义节点；（4） 定义材料；（5）定义单元。（1）定义变量和参数：在Tcl脚本语言中变量是用来替代数字的。一旦被定义好，则可用 来代替数字。当重新要使用这些数字的时候，只需要用“$”字符则可以调用。如果不用此 字符。则变量名将会被视为字符串命令将导致错误结果。在定义模型的变量时，通常首先定 义好材料参数，然后单元参数，再是重力（质量）参数，最后是分析参数 设定好变量定义 的命令文件有很多好处，当需要对材料或尺寸等变量进行修改时，只需在这个文件中修改即 可。无需修改后面的文件。在材料参数中，主要涉及到材料的本构关系变量；单元参数中，主要涉及到单元几何尺寸及 数量；重力参数中，主要是模型质量分布；分析参数中，主要是增量步长及数值计算参数。（2）建立模型和定义节点：这些模型节点通常是置于输入文件的开头。最后的定义结点质 量的命令会覆盖先前定义的节点质量。在定义这些节点时也可以使用变量来定义，效果是一 样的。（4）定义材料：以下的例子是用于定义钢筋混凝土结构的材料。假定变量都已定义。当这 些命令编制在一个外部文件中时，可以用source命令执行。# MaterialsRC.tclset Concrete Mate rialType inelastic； # options: elastic, inelasticset SteelMaterialType ,rhysteretie;# options: elastic丁tiilinearVh押teretic# CONCRETEif ConcreteMaterialType =elastic uniaxialMaterial Elastic $IDcone $EcuniaxialMaterial Elastic $IDcover $Ecif $ConcreteMaterialType = inelastic uniaxialMaterial ConcreteOl $IDcore $fc1C $eps1C $fc2C $eps2C; # Core concrete uniaxialMaterial 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 $epsY $Fy1 $epsY1 $Fu $epsll -$Fy -$epsY -$Fy1 -$epsY1 -$Fu -$epsU $pinchX $pinchY $damage1 $damage2 $betaMUsteel（5）定义单元：set ColumnType inelastic; set np 5;# source RCcircSection.tcl;# options: rigid elastic inelastic# number ot integration paints# proc to define circular fiber section for flexural characteristics-ELEMENTS.tcl# COLUMNSRCcircSection $IDcolFlex $riCol $roCol $cover$IDcore $IDcover$IDsteel $NbCol $AbCol $nfCoreR $nfCoreT $nfCoverR $nfCoverTuniaxialMaterial Elastic $IDcolTors $GJ; section Aggregator $IDcolSec $IDcolTors# Define torsional stitfnessT -section $IDcolFlex; # attach torsion and flexuregeomTransf Linear $IDcolTrans 0 0 1; if $ColumnType = rigid set $lyCol expr$lyCol*$lyCol;set $lzCol expr $lzCor$lzCol;element elasticBeamColumn