*Model change 应用技巧进阶

toosky

谢谢版主的解答～

我的理解，这样的话对于土木的模型来说，为了建立追踪单元，需要在CAE中新定义三个梁单元的set：柱set、X方向梁set、Y方向梁set。然后在inp中用*elcopy命令分别复制这三个set，之后在赋予这三个set相应的梁单元截面。总感觉有点麻烦啊

我还有想到另外一种方法，就是不用*elcopy这个关键词。直接在CAE中建立好追踪单元，只需让追踪单元和原模型共节点就行了。这样可以在CAE中给追踪单元赋予截面，后处理中还可观察追中单元的变形。我觉得这样做的原理应该是与inp中*elcopy的原理是相同的。

有时间的话准备这两种方法都尝试下，呵呵

bananaliuchao

toosky 发表于 2011-11-10 09:22
谢谢版主的解答～

我的理解，这样的话对于土木的模型来说，为了建立追踪单元，需要在CAE中新定义三个梁单 ...

其实追踪单元的截面属性不用和原单元一致，随便一点就可以了。对于梁单元，追踪单元其实最重要的作用是或得其orientation，这个很重要。关于cae里面设置追踪单元，这个我还真不知道可行性如何，应该不可以的吧我认为，因为毕竟节点编号是相同的。你可以做个小例子试试，然后把经验发上来以飨坛油~

toosky

bananaliuchao 发表于 2011-11-10 11:02
其实追踪单元的截面属性不用和原单元一致，随便一点就可以了。对于梁单元，追踪单元其实最重要的作用是或 ...

在blueshell的启发和bananaliuchao版主的激励下，做了一个简单的算例供大家探讨，采用在CAE中建立追踪单元，实现施工加载。建模过程简述如下：
1、在CAE中建立了与模型位置完全相同的追踪单元（set none），追踪单元与模型节点位置相同（共节点）。模型材料为混凝土，追踪单元材料为none（轻密度、轻弹模）。
2、模型为一混凝土长柱，下半部分设定为concrete-f1（底层），上半部分设定为concrete-f2（顶层）。
3、先杀死二层混凝土单元，给模型底层混凝土单元施加重力（追踪单元没有施加重力）。
4、激活模型二层混凝土单元，并施加重力。（杀死和激活单元是在inp中用*model change实现的，剩下的操作都在CAE中实现）

inp文件直接计算，因为有*model change关键词，不要导入CAE中计算。
结果显示的时候，消隐掉none的set后，就是混凝土单元的变形和应力了。将变形倍数调大后，结果简述如下：
1、模型初始状态，可以看到二层混凝土单元被杀死。


2、底层混凝土单元重力施加完毕后的状态，可以看到，二层混凝土单元几乎没有变形（可能有微小的变形，我猜测是追踪单元自身刚度造成的），由于底层被压缩，二层的混凝土单元整体向下也有位移。

3、二层混凝土单元施加重力完毕后的状态，可以看到二层混凝土单元也被压缩了。


以上仅仅是一个简单模型通过CAE实现建立追踪单元的步骤，也许其中有不太完善的地方，但计算结果比较符合实际情况。
通过这个小例子，发现一个有趣的现象就是，追踪单元本身未施加重力，但最终也产生了变形，我的猜测是因为它与混凝土单元共节点，混凝土单元变形时，追踪单元被迫跟着一起变形。。。。

模型inp文件如下，感兴趣的同行可以自己算一下结果看看，呵呵

droin2010

要使用“追踪单元”，就得在part里面生成“钝化单元”，然后生成“追踪单元”。所以如果直接在instance内部set“钝化单元”集，就没法生成“追踪单元”。
最后，我只能用放弃在assembly中定义”钝化单元“集的这一方法。

lijingjing230

写的真详细，有心了

vab123

toosky 发表于 2011-11-10 17:10
在blueshell的启发和bananaliuchao版主的激励下，做了一个简单的算例供大家探讨，采用在CAE中建立追踪 ...

要是有两种材料呢 钢筋混凝土的柱子 两边增设翼墙 想用生死单元使墙在第一步不起作用，第二部起作用 怎么设置啊

toosky

vab123 发表于 2012-4-5 19:55
要是有两种材料呢 钢筋混凝土的柱子 两边增设翼墙 想用生死单元使墙在第一步不起作用，第二部起作用 怎么 ...

道理是一样的哦，建议你好好看一下本帖楼主的帖子，model change的原理与材料关系不大，激活的本质是激活单元，材料参数和截面仅仅是附着在单元上的浮云，本质上还是对单元进行激活。不知道我这么解释能不能稍微好理解些，呵呵

hanmumu

学习了，受益匪浅

xwqu525

要是在早点看到就好了，谢谢楼主分享，学习了，做切削可能用得到

cyl568563088206

大家好，请教下钝化单元上施加边界条件的问题，我现在正模拟施工顺序对4层悬挑阳台的影响，其中一个施工顺序是这样的：浇筑2楼楼板和梁————砌筑二楼阳台围墙————浇筑3楼楼板和梁————砌筑三楼阳台围墙————浇筑4楼楼板和梁————砌筑四楼阳台围墙，边界是在各层楼板和梁端部施加固定约束。
现在我用两种方法模拟边界条件：
1、按施工顺序激活相应单元，并施加相应的边界条件，比如在激活三层楼板和梁的分析步里相应建立它的边界条件，结果疑问如下：三层楼板追踪单元由于没有约束，所以会产生一个由于二楼结构影响的初始位移，但是我激活三楼并施加相应约束后，为什么三楼结构的位移就没有了呢？个人感觉应该还有一个追踪单元残留下的位移的？其加载命令流如下：
** BOUNDARY CONDITIONS
**
** Name: bound2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound2, ENCASTRE
** ----------------------------------------------------------------
**
** STEP: beam2
**
*Step, name=beam2
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,remove
beam3,beam4,beam5,beam6,wall2,wall3,wall4,wall5
**
** Name: selfweighbeam2   Type: Gravity
*Dload
beam2, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall2
**
*Step, name=wall2
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall2
**
*model change,remove
Part-combine-1.wallcopy2
**
** Name: selfweighwall2   Type: Gravity
*Dload
wall2, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam3
**
*Step, name=beam3
*Static
1., 1., 1e-05, 1.
**
** BOUNDARY CONDITIONS
**
*model change,add
beam3
**
*model change,remove
Part-combine-1.beamcopy3
**
** Name: bound3 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound3, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam3   Type: Gravity
*Dload
beam3, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall3
**
*Step, name=wall3
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall3
**
*model change,remove
Part-combine-1.wallcopy3
**
** Name: selfweighwall3   Type: Gravity
*Dload
wall3, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam4
**
*Step, name=beam4
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam4
**
*model change,remove
Part-combine-1.beamcopy4
**  
** BOUNDARY CONDITIONS
**
** Name: bound4 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound4, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam4   Type: Gravity
*Dload
beam4, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall4
**
*Step, name=wall4
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall4
**
*model change,remove
Part-combine-1.wallcopy4
**
** Name: selfweighwall4   Type: Gravity
*Dload
wall4, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam5
**
*Step, name=beam5
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam5
**
*model change,remove
Part-combine-1.beamcopy5
**  
** BOUNDARY CONDITIONS
**
** Name: bound5 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound5, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam5   Type: Gravity
*Dload
beam5, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall5
**
*Step, name=wall5
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall5
**
*model change,remove
Part-combine-1.wallcopy5
**
** Name: selfweighwall5   Type: Gravity
*Dload
wall5, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam6
**
*Step, name=beam6
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam6
**
*model change,remove
Part-combine-1.beamcopy6
**  
** BOUNDARY CONDITIONS
**
** Name: bound6 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound6, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam6   Type: Gravity
*Dload
beam6, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
2、由于初始时没有给追踪单元施加约束，激活部分结构会引起追踪单元产生一个初始位移，这与实际结构的位置产生不一致，因此在初始步中我就施加了全部的边界约束（由于钝化单元与追踪单元共节点，所以施加在钝化单元上的边界与施加在追踪单元上的约束应该是一致的），便这样又产生了如下问题：按如上所述的施工顺序进行时，按道理底层的楼板和梁应该是不断增大的（各楼板与梁及围墙等构造截面一样），但我计算结果却不是如此，比如二楼围墙加上后，二楼楼板与梁应力增大了（3.5MPA），当把三楼楼板与梁激活后，二楼的楼板与应力反而减小了一点（3.2MPA），不知道什么原因？相应命令如下：请各位帮忙诊断下，谢谢啦！
** BOUNDARY CONDITIONS
**
** Name: bound2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound2, ENCASTRE
** Name: bound3-1 Type: Displacement/Rotation
*Boundary, fixed
bound3, 2, 2
** Name: bound4-1 Type: Displacement/Rotation
*Boundary
bound4, 2, 2
** Name: bound5-1 Type: Displacement/Rotation
*Boundary
bound5, 2, 2
** Name: bound6-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
bound6,ENCASTRE
** ----------------------------------------------------------------
**
** STEP: beam2
**
*Step, name=beam2
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,remove
beam3,beam4,beam5,beam6,wall2,wall3,wall4,wall5
**
** Name: selfweighbeam2   Type: Gravity
*Dload
beam2, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall2
**
*Step, name=wall2
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall2
**
*model change,remove
Part-combine-1.wallcopy2
**
** Name: selfweighwall2   Type: Gravity
*Dload
wall2, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam3
**
*Step, name=beam3
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam3
**
*model change,remove
Part-combine-1.beamcopy3
**  
** BOUNDARY CONDITIONS
**
** Name: bound3-1 Type: Displacement/Rotation
*Boundary, op=NEW
** Name: bound3-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
bound3, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam3   Type: Gravity
*Dload
beam3, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall3
**
*Step, name=wall3
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall3
**
*model change,remove
Part-combine-1.wallcopy3
**
** Name: selfweighwall3   Type: Gravity
*Dload
wall3, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam4
**
*Step, name=beam4
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam4
**
*model change,remove
Part-combine-1.beamcopy4
**
** BOUNDARY CONDITIONS
**
** Name: bound4-1 Type: Displacement/Rotation
*Boundary, op=NEW
** Name: bound4-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
bound4, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam4   Type: Gravity
*Dload
beam4, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall4
**
*Step, name=wall4
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall4
**
*model change,remove
Part-combine-1.wallcopy4
**
** Name: selfweighwall4   Type: Gravity
*Dload
wall4, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam5
**
*Step, name=beam5
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam5
**
*model change,remove
Part-combine-1.beamcopy5
**
** BOUNDARY CONDITIONS
**
** Name: bound5-1 Type: Displacement/Rotation
*Boundary, op=NEW
** Name: bound5-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
bound5, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam5   Type: Gravity
*Dload
beam5, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: wall5
**
*Step, name=wall5
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
*model change,add
wall5
**
*model change,remove
Part-combine-1.wallcopy5
**
** Name: selfweighwall5   Type: Gravity
*Dload
wall5, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
** Name: bound6-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
** Name: bound6-1 Type: Displacement/Rotation
*Boundary, op=NEW
bound6, 2, 2
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: beam6
**
*Step, name=beam6
*Static
1., 1., 1e-05, 1.
**
*model change,add
beam6
**
*model change,remove
Part-combine-1.beamcopy6
**  
** BOUNDARY CONDITIONS
**
** Name: bound6-1 Type: Displacement/Rotation
*Boundary, op=NEW
** Name: bound6-3 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
bound6, ENCASTRE
**
** LOADS
**
** Name: selfweighbeam6   Type: Gravity
*Dload
beam6, GRAV, 9800., 0., 0., 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
不知道哪种边界条件模拟更好呢？

tumuseng

bananaliuchao 发表于 2011-11-10 11:02
其实追踪单元的截面属性不用和原单元一致，随便一点就可以了。对于梁单元，追踪单元其实最重要的作用是或 ...

请问在基坑开挖支护中若使用elcopy来模拟支护结构（譬如说连续墙、内支撑等）在开挖过程中的施工，那么在支护结构与土体之间还需要添加接触吗？不添加的话感觉不太对，可是两者又采用的是共节点，到底怎么才算对

jaylin628

请问单元钝化是否有应力重分布现象，是等效成实际施工中的哪一类型工序呀？

peterjony

非常不错的内容

peterjony

很不错

mmx

这个技术关于生死单元的变形协调太有用了。牛逼

tumuseng

rocksoul 发表于 2009-4-18 13:31
老兄弄得越来越精细了。，感谢老兄对我的多次帮忙，呵呵
对了，第5步中.在需要激活的阶段激活已 ...

你好，你现在隧道开挖研究的很好了把？能不能留个qq号相互联系一下？

xxyyjjbm

toosky 发表于 2011-11-10 09:22
谢谢版主的解答～

我的理解，这样的话对于土木的模型来说，为了建立追踪单元，需要在CAE中新定义三个梁单 ...

比如在钝化单元和非钝化单元间设置的约束（比如铰接），在激活单元前还是后设置？

peterjony

楼主的资料非常好，我想请教一个问题，如何把收缩徐变考虑进去呢，最近在做，一直没有头绪，恳请赐教

ty443555288

学习学习，模拟焊接好像是能用到焊接

zxmtzxq

顶起！感谢！