怎样在三维框架建立Iso2spring单元?

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
luochuizi
Posts: 18
Joined: Fri Jun 07, 2019 6:54 pm

怎样在三维框架建立Iso2spring单元?

Post by luochuizi »

我在框架底部加了Iso2spring单元,用的是0长度截面,支座单元一端固接,一端只约束转角,
约束转角的一端连接框架主题。对框架进行时程分析,运行出错了。部分命令如下:
fix 117 1 1 1 1 1 1;
fix 126 0 0 0 1 1 1;
section Iso2spring 1 1.0e-10 2000 100000 300 720000 610.0 982000000 0.0;#单位mm
element zeroLengthSection 260 117 126 1
et pi 3.14159265;
set xDamp 0.05;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
puts "lambdaI = $lambdaI"
puts "lambdaJ = $lambdaJ"
set omegaI [expr pow(abs($lambdaI),0.5)];
set omegaJ [expr pow(abs($lambdaJ),0.5)];
set alphaM [expr $xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaKcurr [expr 2.*$xDamp/($omegaI+$omegaJ)];
set T1 [expr 2.0*$pi/$omegaI];
set T2 [expr 2.0*$pi/$omegaJ];
puts "T1 = $T1 s";
puts "T2 = $T2 s";
rayleigh $alphaM $betaKcurr 0 0
运行结果出现lambdaI<0的情况:
lambdaI = -4.562712e-012
lambdaJ = -1.232921e-013
T1 = 2941496.5235992153 s
T2 = 17894198.181779195 s
我把0长度单元去掉或者两端固接模型能正常运行,到底出了什么问题?
Post Reply