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

FRP混凝土用ANSYS模拟的命令流文件 Combin39单元学习小结(本文来源于清溪长河) 很有用，转过来慢慢学了(2011-07-25 10: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 ! CFRP et,4,solid45 ! 垫块 et,5,combin39 ! 弹簧单元 （沿切向Z向） et,6,combin39 ! 弹簧单元 （沿法向X向） et,7,combin39 ! 弹簧单元 （沿法向Y向） keyopt,1,1,1 ! 不考虑大变形 keyopt,1,5,1 !混凝土的线性解输出控制：给出每个积分点的解 keyopt,1,6,3 !混凝土的线非性解输出控制：给出每个积分点的解 keyopt,1,7,1 ! 开裂后考虑应力松弛 keyopt,3,1,2 ! 作用在拉力状态下的刚度，压力状态下崩溃 keyopt,5,3,3 !沿Z 向平动 keyopt,5,1,0 ! 卸载路径与加载路径相同 keyopt,5,2,1 !不能承受压力 keyopt,5,4,0