柱的尺寸为600x600x4000。内嵌钢管Φ400×14。钢管外侧混凝土的强度等级为C40，钢管芯内混凝土强度等级为C60。钢管采用Q345B，钢筋为HRB400。配筋率同上。试件通过施加预应力来模拟二次受力（暂时不考虑这个），同时考虑偏心，偏心距为1/400高度加2cm，即30mm。构造图如下图所示。

1.创建部件 （1）C40混凝土柱 把柱中间部分掏空，如下图所示，创建完部件就切割以便划分网格

（2）C60混凝土柱 （3）钢管 （4）两种箍筋和竖向受压筋 （5）栓钉杆和栓钉头 2.材性设置 混凝土选择CDP本构模型，HRB400和Q345B选择理想弹塑性模型 创建截面并赋予截面属性，注意将箍筋和受压筋的截面类型设置为truss类型，同时注意修改之后装配完的钢筋网的截面类型为truss（merge成钢筋网后其截面属性自动设置为beam） 3.装配 将所有部件都创建为实例，注意所有的部件的实例类型均为***非独立***，即网格划分在部件上。将箍筋和受压筋通过阵列的方式merge成钢筋网，同时将栓钉杆和栓钉头merge成栓钉，并创建新的非独立实例。 4.划分网格 划分网格不需要很密。 （1）C40混凝土 （2）钢筋网 （3）钢管：钢管管壁较薄，但至少应该划分为两层，将钢管长度方向种子布置密度小一些，厚度方向和径向的种子按个数布置（为边布种）。 （4）栓钉 由于将栓钉和钢管外表面相互作用设置为tie方式，钢管外表面和C40混凝土相互作用设置为接触，栓钉和C40混凝土相互作用设置为内置区域，如果直接这么做会显示如下的错误信息。

2184 nodes are used more than once as a slave node in *TIE keyword. These constraints are not removed either due to the presence of *CONSTRAINT CONTROLS,NO CHANGES keyword or because removing the constraints might affect the model. The nodes have been identified in node set ErrNodeOverconTieSlave.

原因是栓钉底面的所有节点被重复利用。要解决这个问题，切割栓钉时在栓钉底面往上切10mm的平面，栓钉内置于C40混凝土时不选择这10mm的部分。 划分完网格之后再将栓钉阵列在钢管的外表面 5.分析步 创建一个静力通用的step-1即可，时间长度500，最大增量步数1000，初始增量步0.02，最小1e-5，最大500。场变量输出只勾选需要的即可，可加快计算速度。

6.相互作用 （1）耦合点RP-1、RP-2 偏心位移荷载施加在RP-2上 （2）绑定tie 将栓钉底表面和钢管外表面设置为绑定，并设置1mm的位置公差。（选择从表面要一个一个点栓钉底表面。。小心点） （3）内置 将栓钉和钢筋网内置于C40混凝土内。不赘述。 （4）接触 创建一个“接触”的相互作用属性，摩擦公式为“罚”，摩擦系数设为0.3，其他默认。 再创建“接触”的相互作用，使C40和C60混凝土的内表面和外表面分别和钢管的外表面和内表面分别接触，表面平滑选择自动平滑。混凝土和钢管的表面接触一般选择混凝土表面设为主表面，钢管表面设为从表面。 7.创建set 由于要输出参考点RP-2的荷载位移曲线，因此提交作业前应该设置一个RP-2的历程变量输出H-Output-RP-2。工具–>集–>创建，类型选为几何，选中RP2点。在step-1新建历程输出H-Output-RP-2，作用域选择集，选择刚刚创建的集，输出变量选择位移和作用力。 8.载荷 施加位移荷载，8mm，位移荷载均匀施加。 9.结果 （1）整体变形 （2）C40混凝土 最小主应力达到26.33MPa，即最大受压应力为26.33MPa （3）C60混凝土 最大受压应力为39.74MPa （4）钢管 （5）钢筋网 （6）RP-2荷载位移曲线 由于位移是随时间均布施加，荷载位移曲线即荷载时间曲线。 可见最大承载力达到12067kN。手算结果为12528.4kN。