file.delete删除不了文件_opensees 初学-关于输出文件以及输出部分的书写

本文档详细介绍了在OpenSees中进行计算分析和设置加载方式的方法。从文件管理、节点和元素信息的输出,到各种记录器(Recorder)的使用,如节点位移、反力和元素抵抗力等的记录,深入解析了OpenSees的建模后期步骤。通过实例代码,展示了如何指定输出文件、设定不同类型的记录参数,以及针对不同元素类型调用setResponse方法获取特定响应数据。
摘要由CSDN通过智能技术生成

opensees中,总是觉得在前面建模,比较简单,但是一到第二部分,计算分析以及写加载方式的时候,就会出一堆的错误,决定,今天自己好好的研究一下这个部分。

首先,在建模之初,就应该写定好你的文档放置的位置,自己比较喜欢将所有的文档都放在一个一个文件夹里,所以,我就会先命名一个文件夹:代码如下:

if { [ file exists output] == 0 } {

file mkdir output;

}

这几行的意思是如果存在输出文件,建立输出文件的保存路径,即结果保存在名为output的文件夹里。

(file mkdir 创建一个新目录,file delete,删除文件,file copy 复制文件,file exists 如果文件存在返回1,否则返回0)

file mkdir OUTPUT ; # 创建OUTPUT文件夹 (需要注意的是,如果你想在这一行中的后面添加你的注释,一定要加‘;’,否则在opensees中会报错)

建立好文件夹后,写你建的模,此处省略一万字,

开始写输出部分的时候,需要注意的是,opensees 并不保存和输出所有节点,单元,材料等信息,所以只能输出用户只能信息:

如果是需要node信息文件:

写node比较简单, 就用几个例子来解释:

recorder Node - file output/reaction_1.out - time -node 1 -dof 1 2 reaction

recorder 后Node表示的是节点信息,此外还有Elementhe 和 Graphics两大记录方式,(element后面有介绍)

time表示每计算时步都记录

node 1 1号节点

dof 1 2 表示xy方向

reaction 指的是反力

recorder Node - file output/disp_1.out - time -node 4 -dof 1 2disp

time表示每计算时步都记录

node 4 4号节点

dof 1 2 表示xy方向

disp 指的是位移

如果是需要element信息文件:

recorder Element <-file $fileName> <-time> <-ele ($ele1 $ele2 ...)> <-eleRange $startEle $endEle> <-region $regTag> <-ele all> ($arg1 $arg2 ...)

$ele1 $ele2 ...先前定义的element number

$startEle $endEle 这是对于eleRange 来说

$regTag 这个是对于eleRange 来说,这个域中元素你需要哪些被记录(可选)

$arg1 $arg2 ...传递给setResponse()元素方法的参数

setResponse() 的方法主要取决于元素的类型。例如:

Beam-Column Elements (Beam With Hinges Element, Displacement-Based Beam-Column Element, Elastic Beam Column Element, Nonlinear Beam Column Element) : (opensees中大家最常用的元素)

globalForce – 坐标中元素的抵抗力element resisting force in global coordinates

比如:

recorder Element -file ele1global.out -time -ele 1 globalForce

localForce – 局部坐标中元素的抵抗力element resisting force in local coordinates

比如:

recorder Element -file ele1local.out -time -ele 1 localForce

这个意思就是,记录元素1的在局部坐标中的抵抗力,输出文件为:ele1local.out

Sections:

section $secNum –沿着元素长度方向的特定的截面需要相应的数量

$secNum 指一个整合点关于谁的数据需要被输出

force – 截面的力

比如:

recorder Element -file ele1sec1Force.out –time -ele 1 section 1 force

deformation – 截面的变形

比如: recorder Element -file ele1sec1Defo.out –time -ele 1 section 1 deformation

stiffness – 截面的刚度

比如: recorder Element -file ele1sec1Stiff.out –time -ele 1 section 1 stiffness

stressStrain – 记录截面的应力应变响应

example: recorder Element -file ele1sec1StressStrain.out –time -ele 1 section 1 fiber $y $z <$matID> stressStrain

其中: $y:指的是 (fiber截面的局部y坐标)

local y coordinate of fiber to be monitored

$z:指的是 (fiber截面的局部z坐标)

local z coordinate of fiber to be monitored*

$matID:前面定义的材料标签

ZeroLength section element

只有一个截面,所以不需要在recorder命令中去定义截面;

例如:

recorder Element -file Element1.out -time -ele 1 section 1 fiber 0.10 0.10 stressStrain

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
当然可以!下面是一个简单的opensees连续梁-拱组合桥建模示例: ``` # 定义材料属性 uniaxialMaterial Elastic 1 3000 # 弹性材料 uniaxialMaterial Steel01 2 42000 200000 0.01 0.1 # 钢材材料 uniaxialMaterial Concrete01 3 -4.0 0.002 0.1 2.0 # 混凝土材料 # 定义截面属性 section Fiber 1 { patch rect 0.0 0.0 1.0 1.0 3.0 } section Fiber 2 { patch rect 0.0 0.0 1.0 1.0 2.5 } # 定义节点和单元 node 1 0.0 0.0 node 2 5.0 0.0 node 3 10.0 0.0 node 4 15.0 0.0 node 5 20.0 0.0 node 6 25.0 0.0 node 7 30.0 0.0 node 8 35.0 0.0 node 9 40.0 0.0 node 10 45.0 0.0 node 11 50.0 0.0 node 12 0.0 10.0 node 13 5.0 10.0 node 14 10.0 10.0 node 15 15.0 10.0 node 16 20.0 10.0 node 17 25.0 10.0 node 18 30.0 10.0 node 19 35.0 10.0 node 20 40.0 10.0 node 21 45.0 10.0 node 22 50.0 10.0 element elasticBeamColumn 1 1 12 1 1 $secTag # 连续梁1 element elasticBeamColumn 2 12 13 1 2 $secTag # 支座1 element elasticBeamColumn 3 13 14 1 1 $secTag # 连续梁2 element elasticBeamColumn 4 14 15 1 1 $secTag # 连续梁3 element elasticBeamColumn 5 15 16 1 1 $secTag # 连续梁4 element elasticBeamColumn 6 16 17 1 2 $secTag # 拱1 element elasticBeamColumn 7 17 18 1 2 $secTag # 拱2 element elasticBeamColumn 8 18 19 1 2 $secTag # 拱3 element elasticBeamColumn 9 19 20 1 1 $secTag # 连续梁5 element elasticBeamColumn 10 20 21 1 1 $secTag # 连续梁6 element elasticBeamColumn 11 21 22 1 2 $secTag # 支座2 ``` 这个示例定义了一个简单的三跨连续梁-拱组合桥模型,它由10个弹性梁单元和3个弹性拱单元组成,节点上施加了一个垂直向下的荷载。你可以根据自己的需要修改材料、截面和节点位置等参数。注意,这个示例中的支座单元使用了不同的截面属性。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值