资源描述
土结动力分析中人工边界还需要考虑其重力作用吗我用的是刘晶波的三维一致粘弹性人工边界。根据他的理论，不考虑人工边界的重力。设置人工边界的密度为一个接近零的值。
很多的话题在讨论动力相互作用问题，也包括相互作用时地震作用的输入方法，还有人工边界的设置和使用等问题。但是恰恰没有讨论（1）使用人工边界时地震波的输入方法，（2）以及相互作用时地震波动的输入方法。而实际上，相互作用和设置了人工边界后地震波的输入方法对能否得到正确的结果是重要的。
关于第1个问题，刘晶波老师参编的抗震工程学等论著中都有关于设置人工边界时地震波的输入方法，如将地震波动化作用于人工边界上的等效荷载实现地震波的输入。对于粘性边界和透射边界等其他形式的人工边界，也都有相应的地震波动输入方法。此时，若依然采用自由边界时加与输入加速度相应的惯性力的方法作为地震波输入方法，是否正确？
第2个问题，在考虑相互作用时，吴世明老师的土动力学和倪汉根老师的大坝抗震特性与抗震计算中有运动方程的详细推导过程，其中也包括了地震波的输入方法，分2类，一是输入结构和基础交界面上的地震波，这时首先要进行自由场反应分析，确定交界面上的加速度。二是在基岩输入加速度，由于一般测量的是地表自由场的加速度，因此要进行地震波反演，首先确定基岩的加速度。
不管使用那种方法，可以看出，地震波的输入都是一个关键的问题，若不能正确的输入，分析结果可能就是错误的。因此，有了1楼的问题，这里再次提出：
1、当采用远置边界时，能否采用等效荷载输入方式。
2、当采用人工边界时，能否采用施加惯性力的方式完成地震波的输入。
-
“若依然采用自由边界时加与输入加速度相应的惯性力的方法作为地震波输入方法，是否正确？”
这种思路是地上结构抗震时程分析思路，主要用于分析结构在一定加速度下的响应，由于地下结构受到岩土介质的约束作用，导致地下结构的自振特性响应不显著，因此以上做法应该是 不对的。所以才需要施加对应的粘性、透射等边界，进行阻尼和能量传递特性的模拟。
-
采用各种人工边界时，可以减小地基的计算范围，此时，应当取多大合适，或者如何根据基础的特征尺寸选取计算范围？虽然可以通过选取不同的计算范围，经过试算确定，但当模型很大时，花费时间可能过多，很难操作。有什么经验方法吗？
另外，计算范围和基础特征尺寸的关系如何，比如，当基础特征尺寸为20m左右，可以取5倍基础特征尺寸作为计算范围，如果基础的特征尺寸为200m时，是否也需要取5倍，怎样大致确定这个关系？
学习了一天,发现用connection来模拟阻尼器很好，不知道上面几位有何进步，来讨论讨论吧
-
abq中有dashpot1 dashhpot2 dashpota element，: n$ S! X# a4 H* ^5 s1 s9 q
包括非线性阻尼。
-
在边界节点附近生成一个坐标有一点差别的节点,然后由这些节点和边界节点形成Jointc单元。在这里可以定义spring单元和Dashpot单元特性。
! S" B0 a2 _+ Q3 m$ F8 h& H/ n2 jJointc单元是Spring单元和Dashpot单元的并联形式。参见Jointc单元说明
在边界附近生成一排有坐标差的节点，是在边界的外侧吧？这些节点和边界节点联结成jointc单元的意思是不是给他们jointc单元类型？, o2 h! A7 a( J
能否提供个实例，inp最好不过了。
按照帮助中给的阻尼系数和弹簧系数设定粘弹性边界后,动力分析的结果和ansys里面的完全不搭边,不知道原因何在?
/ j( E3 g) g; [帮助中说阻尼系数是材料密度与波速的乘积 弹簧系数不知道怎么确定? 是否是取岩基的刚度?
关于参数的选取可以参考刘晶波的相关文献
还有人不知道如果生成粘弹性边界吗？在此给出两种解决方案。- W$ k$ Q! o! j# }
1.通过cae操作，在已划分好的网格上的要加的节点上创建springer和dashpot单元。具体操作：（1）划分网格后，如图1所示，依次点开part-1...后，双击“springs/Dashpot”,出现图中的对话框，选择接地弹簧阻尼器（connect points to ground），再选择continue，在画好的网格上选择要施加粘弹性边界的节点（可以选择单个节点，也可以选多个节点，此次操作我选择2个节点）之后点击鼠标中间（或done），出现图2所示对话框.（2）在该对话框degree of freedom对应的后面选择弹簧或（和）阻尼器的施加方向（共有1、2、3三平动个自由度方向）。（3）在spring stiffenss 和dashpot coefficient后面填入对应的系数（弹簧刚度和阻尼系数的计算公式可参考刘晶波，杜修力的论文）。最后单击“OK”，操作成功。（4）说明：对于同一个节点可以生产多个自由度方向的弹簧和阻尼器单元，只要重复上述操作，选择不同的自由度方向即可。0 t( s3 s( \2 u9 O; d: _* k
2.通过修改inp文件实现。4 ^; V' Q: n2 x6 q q
*Element, type=Dashpot1, elset=D21 **定义接地阻尼单元类型，将一系列节点创建为该类型单元，并将该一系列单元定义一个单元集合D21.
2 h8 I/ `4 S/ W M; N" b10528, 10 **10528为节点号为10生产的阻尼器单元的单元编号
0 i6 P$ b- ^" J0 I( S7 E% a10529, 11
; q t/ Z/ r6 \* a5 B10530, 132 o& F5 k) U+ S" q! f
10531, 14 m3 R4 Z1 ^- `: `# r$ f
*Element, type=Spring1, elset=S21 **同上，创建弹簧单元，并定义单元集合S21% j6 d8 z- W: P% G& e
10502, 10; r+ I- q) M* ?5 [& u
10503, 117 b3 E. X) ]; N4 r
10504, 131 A$ x3 d( {1 x
10505, 14
) S* X. q) `: G! T) |( B3 _**' `. k$ x. ]! _" x7 _# ~; h& L/ F
*Dashpot, elset=D21 **赋给自由度方向和阻尼系数" H7 m: ?! {1 N5 ]+ V, Y
1 **阻尼器自由度方向为1（广义坐标的1）5 O4 q+ q9 E1 b4 M* v- z8 z
216E6. **阻尼系数 大小为216E6.
0 D+ C0 `" J- j/ W, d, E*Spring, elset=S21 **同上，赋给弹簧自由度和刚度系数
y/ w$ y5 f' H: u1
7 q" Z1 B, V# x6 [119E7.
^: g3 N6 s' [2 a: C**
7 p; T4 Y+ ?- n! h6 N; v; o- Q2 ^需要说明的是以上操作适用于aba/standard，对于aba/explicit本人没有尝试，好像帮助中有说明，很简单，懒得看了。有用到的可以去查看。
sd.jpg (74.94 KB, 下载次数: 59)
2008-8-20 16:51:01 上传
下载次数: 59
sd-02.jpg.jpg (72.26 KB, 下载次数: 34)
2008-8-20 16:51:01 上传
下载次数: 34
刘晶波 1998年：结构_地基动力相互作用问题分析的一种直接方法
: ~# S9 y, R) d8 J0 I可以用无限元做动力边界。而且收敛性更好。呵呵。
-
图是我从一篇文章上截取下来的，请帮看看！这个节点支配的面积到底该怎么去算！！？
33.bmp (185.89 KB, 下载次数: 16)
2006-3-16 15:02:43 上传
下载次数: 16
人工, 粘弹性
相关帖子
· • 瞬态分析能不能用人工粘度？
· • 关于人工边界
· • 人工边界条件的声散射问题
· • 求助：谁有地震计算中使用的人工波??
· • 人工肌肉的运动
· • 请问用COMBIN14模拟人工边界的问题
· • [请教]ansys如何实现人工阻尼边界的功能？？
· • 土——地下结构地震作用时人工边界的选取
· • 哪位有人工假体的模型！
· • 人工粘性系数的设定对计算结果有多大影响
回复
使用道具 举报
· 变色卡
·
Patrick_823357 当前离线
最后登录
2011-9-8
在线时间
61 小时
技术积分
0
仿真币
5
注册时间
2005-7-23
阅读权限
5
帖子
207
精华
0
积分
0
UID
124652
·
Re:土的粘弹性人工边界
-
如果我的有限元模型如下图所示（厚度 可以看作单位 1），切向 法向是不是可以就像图上那样设置两个有一定夹角的两个 一维弹簧？？？
YsModel001.jpg (34.82 KB, 下载次数: 5)
2006-3-16 15:11:25 上传
下载次数: 5
Patrick_823357 当前离线
最后登录
2011-9-8
在线时间
61 小时
技术积分
0
仿真币
5
注册时间
2005-7-23
阅读权限
5
帖子
207
精华
0
积分
0
UID
124652
Re:土的粘弹性人工边界
根据波动推倒出来的
Untitled-2 copy.jpg (62.43 KB, 下载次数: 3)
2006-3-16 15:19:44 上传
下载次数: 3
-
用COMBIN14对平面问题进行人工边界模拟，似乎并没有得到预想的结果。用简单截断边界法时，当地基的计算范围取到大于31倍结构尺寸时，计算结果将不再发生变化。但当施加了COMBIN14的人工边界后，地基计算范围也要取到31倍结构尺寸以后，计算结果才不会发生变化，而且其计算结果与简单截断法结果一样。所不同的是，在31倍以前，两种边界的计算结果是不一样（地基计算范围相同的情况下）。换一种说法就是，COMBIN14边界在31倍以后是没有作用的，其计算结果与自由边界相同，亦即当地基的计算范围大于结构尺寸31倍以后，加与不加COMBIN14人工边界都是一样的；但在31倍以前，COMBIN14人工边界是有影响的，但它仅对计算结果的大小有影响（计算范围不同，其计算结果也不同），并没有达到预想的结果。预想的结果是，加了COMBIN14后，可减小地基计算范围来达到简单截断法自由边界的较大计算范围的计算精度。举个例子就是，如果用自由边界地基计算范围要大于31倍结构尺寸后，计算结果才不发生变化，那么采用COMBIN14人工边界后，其地基计算范围应当可能是大于10倍结构尺寸后就不发生变化，而且其计算结果与自由边界地基计算范围大于31倍结构尺寸后的结果相同。不知我的想法对不对？我计算结果的问题又出在哪里？希望在这方面有经验的人指点一下。谢谢。
-
我用DYNA计算地震下边坡的反应
1 v8 v; S) Z6 y5 \7 @5 I* M5 ` 四边我加的是无反射人工边界
% u& O3 W7 ]+ c+ p- E7 t 但计算结果感觉其结果更接近自由边界 3 Q2 q! F4 y7 A* K6 ]7 s) X: {
9 I3 @9 F3 b. p+ S# `! H: _
DYNA的无反射边界到底效果怎么样 / I- X$ x* {( j
可以不可以在同样一组节点上既设无反射边界 ( L' @2 m: L' W- y' S7 t
又设其它固定约束啊
-
同感同感，我计算中的结果也是这样，整个模型总是朝一个方向跑。 / t7 A/ T! B! r4 S( A! A# `
不反射边界是不能加其他约束的，否则不反射边界失效。这是ansys帮助上叙述的。
5 l6 M( g5 ]: Q8 H) m7 Q 希望和楼上多交流！
-
有一定位移是应该的 / p3 v/ z5 x8 I6 V8 p* R3 U2 k
但是得出的结论就有点问题了 ) V, G) J" v* x& L
边坡的地震放大系数应该离边缘越远越小
* n% I2 J* o. _0 y 但是加了无反射边界后 先是减小，然后就变大 了
' ~6 C" Y8 {$ l% X 感觉那边就没加什么边界 5 `9 A1 |- j1 |) S+ B1 d
也当自由边界，按边坡边缘了
不知道楼上的兄弟解决了问题没？或者有什么想法嘛?
5 b3 j9 s* Q* C$ k) } 我想了好久也没有解决的办法，本来想把不反射边界加在里面一层单元，然后最外面加不反射边界，可是不可以，不反射边界只能加在最外面。 - [& c" V+ H2 M" u X- p. _9 C
不知道位移的影响会不会对最后的结果造成致命的影响？你有试验对比嘛？
是不是我们用的方法不对啊 2 w7 } t% L! S+ V2 f3 f
或者这种无反射边界只能用在爆炸这种冲击荷载下 而不能用在地震这种连续动力作用？ ) d. z9 K! {0 t E, S
但应该不会啊 - f. E( f( Q6 D& u
这个人工边界是个很古老的人工边界了 2 {# Y0 G( D- l2 {. F. G: K
不应该不能用在地震啊 ! x w* l v* J" X9 D
3 U# p" c& L8 K' k) x' R& o 实在不行，我只能把边界尽量取远了
你看 他曾经用这个做个地震
2 [( v; v- p: _3 |5 [9 _: N6 w 但是讲得太不详细了 * T: o$ F, N% U/ ^
他还有本书 8 X. l6 ^( C# D+ d
等我去借来看看
0 _& V. ^; s `8 ]) I0 U 但我记得也讲得不清楚
) l C) i& V" s( ^% T 只有结论 没有过程
/ _% i2 t: @3 h! N7 w8 p 你看能不能找他的博士论文看看
唉，真是没办法，感觉不反射边界加后，计算结果就不对，因为边界的运动实在是离谱呀，我加y方向的加速度，本来理论上分析不应该又x方向的位移呀，可是计算出来的结果竟然是x方向的位移很大，而且不对称，这是不可能正确的呀。
8 r. N, e% ^5 k# ~- V* m- v 真不知道不反射边界到底怎么了？还是我们用的不对？ , B, H) ]/ }3 X9 Y( K5 `
希望大家努力一下呀
可以用自由边界代替阻尼边界，两者差别很小
这样的边界条件是很难加上的。
) y. K9 I4 c. _建议拓展模型用子模型方法来传递载荷和位移模拟粘弹性边界。
展开阅读全文
浙ICP备2021020529号-1 | 浙B2-20240490 
