FRP混凝土用ANSYS模拟的命令流文件.pdfVIP

FRP混凝土用ANSYS模拟的命令流文件

Combin39单元学习小结(本文来源于清溪长河)很有用，转过来慢慢学了

(2011-07-2510:01:48)转载标签:ansyscombine39单元杂谈分类:专业学习

这两天在帮助一位同学解决论文中的问题，其间涉及到组合结构不同材料的构

件的连接。原本，她想做的是二者分离的动态过程。要做到这一点，貌似要用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!CFRPet,4,solid45!垫块

et,5,combin39!弹簧单元(沿切向Z向)et,6,combin39!弹簧单元(沿法

向X向)et,7,combin39!弹簧单元(沿法向Y向)

keyopt,1,1,1!不考虑大变形

keyopt,1,5,1!混凝土的线性解输出控制:给出每个积分点的解

ke