wangguobo4220 发表于 2006-7-21 21:51:20

请教attach命令

两个不同的介质采用实体单元模拟（由于二者之间主要是挤压作用，接触面参数太麻烦），所以不考虑设置接触面，我就想用attach命令，但对这个命令不是很熟，我只尝试用了命令：attach face就默认对整个模型起作用，但我想attach gp 呢？比如一个面上的某个或某些gp，用命令：attach gp id1 face zid fid ，那还要先知道gp号码以及zid、fid，浩麻烦！另外，还想知道用了这个命令后到底起到一个什么样的效果？

[ 本帖最后由 wangguobo4220 于 2006-7-22 18:25 编辑 ]

wangguobo4220 发表于 2006-7-22 09:34:43

建一个箱，如下
gen zone brick p0 -0.4 -0.20.0 p1 -0.2 -0.20.0 p2 -0.4 2.60.0 p3 -0.4 -0.2 1.2 size 114 6   ;left
gen zone brick p03.2 -0.20.0 p13.4 -0.20.0 p23.2 2.60.0 p33.2 -0.2 1.2 size 114 6   ;right
gen zone brick p0 -0.2 -0.20.0 p13.2 -0.20.0 p2 -0.2 0.00.0 p3 -0.2 -0.2 1.2 size 171 6    ;front
gen zone brick p0 -0.22.40.0 p13.22.40.0 p2 -0.2 2.60.0 p3 -0.22.4 1.2 size 171 6;behind
gen zone brick p0 -0.4 -0.2 -0.2 p13.4 -0.2 -0.2 p2 -0.4 2.6 -0.2 p3 -0.4 -0.2 0.0 size 19 14 1 ;bottom

各部分的单元尺寸吻合的很好啊，可用attach face 后有一个gd被attach了，这是怎么回事？

wangguobo4220 发表于 2006-7-22 18:46:44

做了两个简单的例子来理解attach命令，权当抛砖引玉吧！
一、左右两个块体，故意将网格尺寸不一致：
gen zone brick p0 0 0 0 p1 6 0 0 p2 0 6 0 p3 0 0 4 size 3 2 2
gen zone brick p0 6 0 0 p1 7 0 0 p2 6 6 0 p3 6 0 4 size 1 3 4
模型土和attach网格图如下。进行动力分析时发现左面块体右侧面（定义为面1）的各节点均有加速度值，而右边块体的左侧面（定义为面2）上被attach的各节点加速度为0！！！结合attach图，本人的愚见是：可以将面1理解为目标面，attach命令就是把面2上与面1上不重合的各节点“粘”在面1上，使二者成为一体。还有各现象：不管动力分析是应力、速度或加速度输入，被“粘”的点有位移、速度，但就是没有加速度（指加速度一直为0）。

wangguobo4220 发表于 2006-7-22 18:58:10

二

算例二：建上下两块体，网格分为节点一致和节点不一致两中情况。上面块体如同悬臂梁，在上部块体一节点上施加集中力进行静力分析，比较位移云图
1、节点一致：
gen zone brick p0 0.0 0.0 0.0 p1 8.0 0.0 0.0 p2 0.0 6.0 0.0 p3 0.0 0.0 2.0 size 8 6 2
gen zone brick p0 3.0 2.0 2.0 p1 5.0 2.0 2.0 p2 3.0 4.0 2.0 p3 3.0 2.0 12.0 size 2 2 10
2、节点不一致
gen zone brick p0 0.0 0.0 0.0 p1 8.0 0.0 0.0 p2 0.0 6.0 0.0 p3 0.0 0.0 2.0 size 8 6 2
gen zone brick p0 3.5 2.5 2.0 p1 5.5 2.5 2.0 p2 3.5 4.5 2.0 p3 3.5 2.5 12.0 size 2 2 10
给出三种情况下的位移云图：
1、节点一致时的位移云图
2、节点不一致，用attach命令的位移云图
3、节点不一致，不用attach命令的位移云图

比较三个图可很清楚的看到：不用attach命令时，上下块体的位移云图完全分离，位移不连续，而用了attach命令后，节点不一致与节点一致的位移云图相差无几。

这但愿能给用实体单元建模时提供一点启示吧！

lakewater 发表于 2006-7-23 22:26:47

用这种方式理解命令是一种很好的方法，理解深刻！

wangguobo4220 发表于 2006-7-24 10:33:09

attach与接触面

呵呵，谢谢版主了！加了3分啊，惭愧！其实我还是不是很清楚这个命令，还有两个小问题：
1、就是前面提到的，为什么被attach后的节点的节点加速度为零，而位移和速度还存在？
2、如果存在接触面的话，用attach后会不会把两个面“粘”上后，使接触面失去作用？我也做了个例子验证，说明attach命令对接触面有影响。例子如下，同时有几个关于接触面的问题想请教下各位：
config dynamic
set log on

;------create soil and placity board---------------
gen zone brick p00.0 0.0 0.0 p1 3.0 0.0 0.0 p20.0 2.4 0.0 p30.0 0.0 1.0 size 10 8 5;soil
gen zone brick p0 -0.2 0.0 0.0 p1 0.0 0.0 0.0 p2 -0.2 2.4 0.0 p3 -0.2 0.0 1.0 size 18 5;left placity
gen zone brick p03.0 0.0 0.0 p1 3.2 0.0 0.0 p23.0 2.4 0.0 p33.0 0.0 1.0 size 18 5;right placity

gen zone brick p0 -0.2 -1.20.0 p13.2 -1.20.0 p2 -0.2 -1.00.0 p3 -0.2 -1.2 1.0 size 7 1 5            ;front
gen zone brick p0 -0.23.40.0 p13.23.40.0 p2 -0.23.60.0 p3 -0.23.4 1.0 size 7 1 5            ;behind

;------create group-----------
group placity range x -0.2 0.0
group placity range x3.0 3.2
group soil range x 0 3
group box-1 range x -0.2 3.2 y -1.2 -0.2 z 0.0 1.0
group box-2 range x -0.2 3.2 y3.43.6 z 0.0 1.0

;------create interface-------
interface 1 face range y 0.0
interface 2 face range y 2.4
ini y add1.0 range group box-1
ini y add -1.0 range group box-2

plot block group

;----- interface parameters-------
interface 1 prop ks 1e8 kn 1e5 fric 30
interface 2 prop ks 1e8 kn 1e5 fric 30

;-----assign material parameters-----
model elas range group placity
prop bulk 6.88e6 shear 1.475e6
ini dens 15

model elas range group soil
prop bulk 5.2e7 shear 1.1e7 range group soil
ini dens 1760 range group soil

model elas range group box-1
prop bulk 6.88e6 shear 1.475e6 range group box-1
ini dens 1800 range group box-1

model elas range group box-2
prop bulk 5.2e9 shear 1.1e9 range group box-2
ini dens 1800 range group box-2

;--------fix boundary----------
fix x range x -0.2001 -0.1999
fix x range x3.19993.2001
fix y range y -0.2001 -0.1999
fix y range y2.59992.6001
fix z range z -0.00010.0001
set dyna off
set grav 0 0 -10

attach face

hist unbal
hist gp zdis 1.5 1.2 1.0

cy 1000
save static-attach.sav
后面附上用与不用attach命令两情况下接触面的nd、sd、ns和ss云图。前面四个云图是没有attach命令的情况，对比这些图可以很清楚的看出影响了，具体就不写了，如位移降、应力都低很多，说明attach命令会影响接触面的正常工作。

另外，顺着这个例子，想向Lakewater帮助请教两个接触面的问题：这个例子中就是重力作用，因此，两侧的接触面以剪切为主（不用attach命令情况），因此应该ss与sd应该比较大，而且结构、荷载对称，受力和位移也应对成，但从云图上看：1、后面接触面的sd为零！！！后面做动力分析时，也是如此，因此我怀疑接触面是否正常工作了？2、ss只有ns的一半！3、位移和力也不对称。期待各位的指点！

bbs_bbs1000 发表于 2007-11-16 15:02:29

厉害啊 受教了

wangyi2002598 发表于 2007-11-20 13:52:36

好帖啊:D

chchyao 发表于 2007-11-25 12:41:59

非常好 顶起

april_1981 发表于 2007-11-26 13:25:10

受益非浅啊:)
谢谢lz

diwurou 发表于 2008-3-28 12:15:51

如果从别处移来的子网格（移过去后节点一致网格一致）需要用ATTACH粘结嘛？

stevenson06423 发表于 2008-4-1 14:43:13

wangguobo4220,nihao!
我正在做你例子的问题，还是由许多问题不明白，想向你请教，我的qq：405164175，我的邮箱：405164175@qq.com
:handshake

wsf2008 发表于 2008-4-30 00:23:39

我觉得用不用attach关系很大的，尤其是接触面的attach问题。我模拟了两个东西，一个加attach，显示的跟实际情况比较符合，一个不加attach，跟实际情况较吻合，有没有高手给解释一下？

minhrac 发表于 2008-9-5 20:34:05

用attach以后接触面好像不起作用了吧？

risigt 发表于 2010-4-4 18:06:20

这个贴必须顶起，以后好找

张如林 发表于 2010-8-27 14:37:48

受益匪浅！谢谢

yang200909 发表于 2010-10-10 10:12:58

好帖子，谢谢大家了

donglei 发表于 2010-10-11 10:31:28

用attach以后interface 就没用了吗

hnlgdxnyckwl 发表于 2011-5-31 23:49:08

:)我觉得除了网格整数倍关系影响网格连接以外，数据精度应该对网格节点连接有影响，例如：赋值时使用2与2.00000不同精度，从其他软件导入到flac中，由于不同软件精度不同会产生网格节点和面的不连续。。。不知道这样理解对不对？

long_7081 发表于 2011-6-9 13:52:05

楼主下面的那个怎么又加了接触面了！！不是说一attach都连接了么？？？
页: [1] 2
查看完整版本: 请教attach命令