ansys combine39单元模拟滑移

Combin39单元学习小结(本文来源于清溪长河) 很有用，转过来慢慢学了 (2011-07-25 10:01:48)转载▼
标签： ansys combine39单元 杂谈 分类： 专业学习
这两天在帮助一位同学解决论文中的问题，其间涉及到组合结构不同材料的构件的连接。原本，她想做的是二者分离的动态过程。要做到这一点，貌似要用ls-dyna，而这我知之甚匮，加上时间非常有限，这个方案被我首先排除。
至于接触分析，在参考了一些资料后，感觉即便用点-点接触来进行模拟，也会造成计算机资源的浪费。而且此分析经过了MTS的实验，显示连接很紧密，甚至直接耦合的话也相去不远。经过慎重考虑和探讨，也放弃了做接触分析的想法。
在陆新征的那篇引用数高达230的文章的启发下，终于确定了分析方案：用Combin39弹簧单元来模拟组合结构之间的相对黏结滑移。
其实，用combin39模拟钢筋-混凝土之间的滑移，已经是比较成熟的选择了，可我以前并不十分清楚这一点_{囧……看到网上有不少和我一样的combin新手，在问一些低级问题(正是由于他们，所以我只需要Google就OK了})，我就把最近学到的一些知识总结如下，欢迎讨论或者批评指正。
1.Combin39是一种包含两个节点的一维的具有非线性性能的单向受力单元，在1-D, 2-D, 或3-D应用中具有纵向的和扭转的能力，单元的每个节点都具有X、Y、Z三个自由度。

2.该单元具有大变形的能力，此时需要打开NLGEOM选项时，并且使KEYOPT(4) = 1 或者3。

3.combin39需要通过一个荷载～位移（F—D）曲线来定义弹簧的受力性质，至于这个曲线，可以通过查询相关文献获得，比如[杨勇. 型钢混凝土粘结滑移基本理论及应用研究. 西安建筑科技大学，2003.]
无需定义材料属性。

4.和Combin40不同的是，Combin39的刚度不是用户手动输入的，而是程序自动根据F-D曲线求的。刚度K=dF/dD；F=切应力(粘结面本构关系)*粘结面积；随着滑移的发展，刚度在不断退化，“确定刚度”，便是F-D曲线之所以重要的症结所在。
5.曲线需要经过原点(0,0),初始斜率为正，其他各段的斜率可以为正或为负;但如果KEYOPT(1) = 1，则曲线最后的斜率必须为正。

6.如果建立2D、3D的弹簧单元，即希望Combin39除了轴向拉压之外还能有切向刚度，则此时弹簧的两个节点不可重合——程序需根据弹簧轴向确定单元力的方向。
7.相对于上一条，更常用的，则是另模型两部分节点重合，然后在重合的节点建立多个(比如三个方向)Combin39单元；换言之，2个节点之间可以同时存在多个没有长度、方向不同的弹簧。
8.Combin39的方向由Keyopt控制：1D单元自由度由KEYOPT(3)控制，2D和3D则靠KEYOPT(4)。

9.Combin39的计算结果，可输出弹簧力、弹簧两端节点位移等数据。

至于实例，我找到了一个“无主命令流”——即被转的太多，也不清楚最初来源。由于和今天讨论的话题契合甚佳，故附带贴出。作者可联系我，我一定给出说明或者做上链接(比较地道~)
【注】此方法中使用eintf命令生成弹簧，长河以为实践效果不佳，还是通过循环用e,i,j比较托底
Finish
/clear,all
/com,structural
/config,nres,50000
/config,fsplit,750
/prep7
!定义单元类型
et,1,solid65 ! 混凝土
et,2,link8 ! 钢筋
et,3,shell41 ! CFRP
et,4,