ABAQUS添加钢筋混凝土弹簧单元

在inp里修改弹簧单元设置之前，先要做2件事。1、在cae中建立弹簧单元。建模后进入1缸。由田。口，specialfsprings/dashpotsfcreate,在connectivitytype选择connecttwopoints,如下图：分别选择混凝土和钢筋上面各1个点，或者随便选，后面都可以在inp中修改。选好以后按下图修改弹簧参数，1,1,900也是随便的，后面都可以在inp中修改，含义在后面解释。Type:ConnsettaopcincsEpdng/'Ds5h^otPointPairsPoint11zhd-1Verte>c[5]Point11zhd-1Verte>c[5]Point2zhu-1Vertex[4]k:Specifyfixeddireilicn1Direction(StandardOnly)PointDirection(StandardOnly)Point1degreeof ▼Point2degreeoffreedc/□rirntation:(Global)片X2、找出混凝土和钢筋弹簧单元节点处的节点编号建模完成后提交运算，运算出现几步的结果后，进入氏5面15,optionsfcommon一labels—shownodelabels即可显示单元的节点。

若钢筋和混凝土节点重合，想单独查看钢筋或混凝土的节点，tools-displaygroup-create-partinstances,选择钢筋或混凝土，remove,即可隐藏Dkpl即GrcnjpMakeaSelecdcnAllP日rtinst日masElementsNodBESurfacesDisplaygroupsCoordinateAllP日rtinst日masElementsNodBESurfacesDisplaygroupsCoordinatesystemsTiesSheli-Eo-S&lidcouplingsDistnbuhngcoupJin^sKinematiccouplingsRigidBodiesMPC?ItemZHU-1ASSEMBLYDIANP3AN-1JIN-1>nviwportNamefiIter：importb-t+wef-nlinkedviewport?需查看混凝土内部的节点编号，可以用切片功能Wl喻，博珏断鼻-

fflM一，节招飞％直驻nn■图笛人EHShewAllow喻，博珏断鼻-

fflM一，节招飞％直驻nn■图笛人EHShewAllowfoSensitivity:登ViewCutMaMotionofSelecte-dCut以下是修改inp的方法，转载：ABAQUS在inp里添加弹簧单元全过程，作者:天天飞ABAQUS在inp里添加弹簧单元全过程模拟钢筋混凝土粘结滑移，以我现在的水平来看，大概有三个思路，一个是

用cohesive,一个是用Spring2弹簧单元，再有就是Umat了。我先学习的是弹簧单元。ABAQUS里有三种弹簧单元，Springl,Spring2和SpringA。其中，ABAQUS文档的解释如下：SPR.IflG2clcmcnlaIherelalivedisplacemenlacrasia苫阳时GEdemenlis!hediffereruobetweentherwcomporcitold行口日匚ernETtofmespring用firstandthe/hccmponRrttrfdi甲lacemenidthesunng'sMcondnode-wheietdefinedasdescribedbelowandcanbeiniocaldiredions(see''Defningthgdii^tkiriofacl-iunrcfSPRir^jldndSP陶/百2次ftW)Spring2有两个节点，每个节点只有一个自由度。钢筋采用Wire建模，Truss单元，混凝土实体单元。mesh之后会出节点，所以尽量mesh成相同的节点，用钢筋和混凝土重合的节点,在三个方向分别建立弹簧。垂直于钢筋纵向的两个弹簧设置成刚度较大,弹性。沿钢筋方向的弹簧,需要确定粘结滑移本构,定义力和滑移位移的关系。对比一下：Springl是弹性的,只有一个节点，而SpringA的作用是沿两点连线的,每个点有三个自由度。线形的弹簧单元可以在Interaction模块里，找到Special里面添加，非线性的就只能的通过修改inp文件了。内容添加Assembly里，我的习惯是写在EndAssembly前一行。对于添加内容,大概有两部分,一部分是对节点的定义,另一部分是对力-位移曲线的定义。两者是分开的，通过elset识别。elset是Elementset的意思，就是定义一个单元集合，名字根据在ab叫us里添加线性弹簧单元所生成的inp来看，是集合名字后边跟"-spring"。type=Spring2是固定格式，*E1ement,工yp已二Spring2jelset=X-spring这句后边是节点编号,如下图。Wire-1.9,

Wire-1.&

Wire-1.7aWire-1.Wire-1.9,

Wire-1.&

Wire-1.7aWire-1.6.Wire-1.5,Wire-1.4,值得注意的是，Instance名是Part名后边跟编号形成的（如“-1”），如果Part的名字有空格,就需要加引号了。如下图

*E1已in已*E1已in已nt,1,"Steel21"Steel瓦^Steel《"Steel瓦VJSteel6,^SteelPlate-1?\18."CouplingBeam-l^.415Plate-l^.19j"CouplingBeam-1416Pla-te-r;20."CouplingBeam-1\417Plate-l^.21』"CouplingBeam-1^.41SPlate-122,"CouplingBeam-141^Pla-te-r；23,"Cogli口芸Beam-1^420最左边一列的数字是单元的编号，从1开始，顺着排下去就行。中间和右边两列是单元的节点，Wire-1是Instance的名字。查看单元编号和节点编号，在View--PartDisplayOptions里，在mesh一栏分别勾选Shownodelabels和Showelementlabels。在mesh模块的视图里就可以显示了。实体名称和对应的节点编号之间用“.”连接。*Springelset=X-sprin.gL1le+Q7线性弹簧，一般给个大刚度就好，用于定义型钢和混凝土接触面发现法向的接触。+SpringJNonLinear,elset=Y-spring2,2-1770,000,-6,0001-1770,000；-4.000-2040,000；-1.090-3300,000；-0..210-2700,000,0,000d；'ooo,clooo2700.000/0.0003300.000；-0.2102040.000；1.0901770.000；t0001770.000-.6.000非线性弹簧，根据粘结滑移本构得来。关于粘结滑移本构，公式很多，相差也非常大，这里因为规范给出的最好算，还是采取用规范的值。对于型钢混凝土，建议看看杨勇的论文。关于粘结滑移本构，公式很多，相差也非常大，这里因为规范给出的最好算，还是采取用规范的值。对于型钢混凝土，建议看看杨勇的论文。表C3.1混凝土与钢筋间粘结应力-滑移曲域的参数值特征点磨裂(er)峰值（u）残余（0粘结应力（N/xrm3）2.54/t.r相对滑移（mm）0.Q25d004dSr0.55d注：表中d为钢筋直程（mmh人『为混凝土的抗接颐例函蝌帅而蜘4-得到切向的应力之后，再根据网格的划分，确定了每个单元的钢筋的长度之后，可以得到F。Fi=q。"XTldl就是应力乘以面积。*Spring,elset=X-Spring（非线性的Spring2,多一jNonlinear,这个是很多教程都是这么说的，我测试后发现，去掉之后也没有关系）。之后定义弹簧的方向，X轴方向为1,1；Y轴方向为2.2；Z轴方向为3,3。两个数字中，第一个数字的含义是节点1沿着1轴（即X轴）运动，第二个数字的含义是节点2沿着1轴（即X轴）方向运动；F=kX，X就是两点在1轴方向运动的位移差。如果前后数字不同，我理解的，就是向量的加减了。这些定义好之后，inp就可以计算了。我这边有一个完整的inp文件，类似于型钢混凝土的小练习，有需要看这个inp的，在回复中留下邮箱，我会发送到邮箱里。大多数教程讲到这里就没有然后了。后边的一些东西，我感觉还是挺关键的。两点，一个是inp文件做好后，不要导入ab叫us再计算，因为ab叫us识别不了；二是直接在job里提交inp的方法。打开ab叫us，不用新建任何项目，直接到job模块