施加法向值为0的速度为什么对结果有影响?
陈育民书上的例子。new
conf dyn
gen zone brick size 10 5 10
mod elas
mod null range x=0,5 z=5,10
fix z range x=-.1 .1 z=.1 10.1
fix z range x=9.9,10.1 z=.1 10.1
fix y range y=-.1 .1
fix y range y=4.9 5.1
prop bulk 2e8 shear 1e8 ;设置土体参数
prop bulk 4e9 shear 2e9 range x=5,6 z=5,10 ;设置墙体参数(土体参数的20倍)
ini dens 2000 ;设置密度
def setup ;动荷载中的变量赋值
freq = 1.0
omega = 2.0 * pi * freq
old_time = clock
end
setup ;执行变量赋值
def wave ;定义动荷载函数
wave = sin(omega * dytime) ;定义动荷载变量
end
apply xvel = 1 hist wave range z=-.1 .1 ;施加动荷载
apply zvel = 0 range z=-.1 .1
hist gp xvel 5,2,0
hist gp xvel 5,2,10
hist gp zvel 5,2,10
hist dytime
def tim ;估算程序运行的时间
tim = 0.01 * (clock - old_time)
end
set dyn multi on ;设置动态多步
solve age 1.0
new
config dyn
gen zone brick size 1,1,50
model elas
prop shear 1e7 bulk 2e7
ini dens 1000
def setup
omega = 2.0 * pi * freq
pulse = 1.0 / freq
end
set freq=4.0
setup
def wave
if dytime > pulse
wave = 0.0
else
wave = 0.5 * (1.0 - cos(omega * dytime))
endif
end
range name bottom z=-.1 .1
fix z range z=.5 55 ;将上部网格都施加数值向约束
apply dquiet squiet range bottom
apply sxz -2e5 hist wave syz 0.0 szz 0.0 range bottom ;-2e5的系数来源于 的值
apply nvel 0 plane norm 0,0,1 range bottom
hist gp xvel 0,0,0
hist gp xvel 0,0,25
hist gp xvel 0,0,50
hist dytime
hist wave
plot create hhh
plot add hist 1 2 3 vs 4
plot show
solve age 2
我试了下。红字的语句加与不加结果是有差别的。
我想知道为什么要加上红字的语句?都是0值。为什么加与不加结果还不一样? 这个很简单,因为它要考虑的是S波对模型的德影响,不考虑P波,所以不用考虑竖向地震波的影响。 这个很简单,因为它要考虑的是S波对模型的德影响,不考虑P波,所以不用考虑竖向地震波的影响。
cloud328 发表于 2009-8-25 19:14 http://forum.simwe.com/images/common/back.gif
谢谢你的指点。我还是有点困惑。希望您能再说明一下。
apply sxz -2e5 hist wave syz 0.0 szz 0.0 range bottom
apply nvel 0 plane norm 0,0,1 range bottom
载荷只是加了水平X方向的S波。没有加P波啊?你说的P波哪来的?怎么产生的?
谢谢你的指点。我还是有点困惑。希望您能再说明一下。
apply sxz -2e5 hist wave syz 0.0 szz 0.0 range bottom
apply nvel 0 plane norm 0,0,1 range bottom
载荷只是加了水平X方向的S波。没有加P波啊?你说的 ...
mihang 发表于 2009-8-26 08:27 http://forum.simwe.com/images/common/back.gif
这个你看下书就知道,地震时,一般水平地震荷载破坏力度大,施加在模型底部就叫S波,竖向地震荷载就一般就叫P波,他这个就是不考虑竖向地震荷载。所以为0.
页:
[1]