各向同性双管混凝土柱承载力计算方法及系统

本发明属于土木工程领域，具体涉及一种各向同性双管混凝土柱承载力计算方法及系统。

技术背景

相比于普通单钢管混凝土，双钢管混凝土柱具有更优异的抗火性能。双钢管混凝土还可以避免厚壁钢管的采用，而厚壁钢管的供应和加工较普通钢管困难。双层钢管的截面形式还有助于避免混凝土干缩对结构受力性能造成的不利影响。例如，《钢管混凝土结构技术规程》(cecs28：2012)规定，当钢管直径大于2000mm时，应采取同心双层或多层的截面形式，以避免混凝土干缩的影响。双钢管混凝土结构的截面形式还可被用于加固既有cfst构件。

迄今为止，已有的承载力理论分析大多基于极限平衡法，该方法简单将各组成部分分别计算承载力叠加，并乘以折减系数，以考虑各组成部分峰值应变不同步对于承载力的影响。然而，这种不同步对于承载力的影响，应当和外套内外钢管径厚比、内部和夹层混凝土强度等因素有关，但已有的折减系数并未考虑这些因素的影响，而取为某个设定的固定值。这就导致对于一些试件而言，计算承载力低估了实际的承载力，计算结果偏于保守；而对于另外一些试件而言，计算值小于真实承载力；无法准确考虑各组成部分达到各自峰值应变不同步的影响；因而，无法得到准确计算的结果，容易造成结构安全隐患。

技术实现要素：

本发明是为了解决上述问题而进行的，目的在于提供一种各向同性双管混凝土柱承载力计算方法及系统。本发明为了实现上述目的，采用了以下方案：

<方法>

本发明提供一种各向同性双管混凝土柱承载力计算方法，其特征在于，包括以下步骤：

步骤1.获取双管混凝土柱的计算用参数，包括：内部核心混凝土和外部夹层混凝土的泊松比、弹性模量，内部管材和外部管材的泊松比、弹性模量、内径、壁厚；

步骤2.建立双管混凝土柱各部分之间受力和变形关系

(1)内部核心混凝土的受力和变形分析：

(1-1)内部核心混凝土的竖向应力-应变关系

峰值荷载前阶段，内部核心混凝土竖向应力σz,ic表达式为：

峰值荷载后阶段，内部核心混凝土竖向应力σz,ic表达式为：

式中，eic为内部核心混凝土的弹性模量，εz为纵向应变，f′ic,co＝f′ic+4.1p1，p1为内部钢管和内部核心混凝土之间的挤压应力，fc′为混凝土强度，β＝12.16p1/f′ic-3.49，残余应力fic,re：a＝795.7-3.291f′ic，k＝5.79(p1/f′ic)^0.694+1.301，|fic,re|≤0.25|f′ic,co|；

(1-2)内部核心混凝土的环向应变-竖向应变关系

混凝土环向应变εθ,ic包括弹性部分和塑性部分：

弹性部分表达式：

塑性部分表达式：

式中，νc为混凝土的泊松比，εic,z0为混凝土开裂时的应变，表达式如下：

式中，εic为素混凝土达到峰值荷载时的应变；

(2)内部管的受力和变形分析：

内部管环向应力σθ,is表达式：

内部管竖向应力σz,is表达式：

内部管环向应变εθ,is表达式：

式中，ric为内部管的内径，tis为内部管的壁厚，p2为内部管和外部夹层混凝土之间的挤压应力，fis为内部管的屈服强度，es为内部管的弹性模量，νs为内部管的泊松比；

(3)外部夹层混凝土的受力和变形分析：

(3-1)外部夹层混凝土的竖向应力-应变关系

在峰值荷载前，外部夹层混凝土的竖向应力σr,sc表达式为：

式中，是外部夹层混凝土的切线模量；εsc,co为约束状态下外部夹层混凝土受压峰值应力对应的应变：p3为外部夹层混凝土和外部管之间的挤压应力，f′sc为外部夹层混凝土单轴受压时的峰值强度，εsc为外部夹层混凝土单轴受压峰值应力对应的应变；为外部约束状态下夹层混凝土达到峰值应力时的割线模量；f′sc,co为约束状态下外部夹层混凝土受压时的峰值强度，f′sc,co＝f′sc+4.1(p2+p3)/2；

峰值荷载后，外部夹层混凝土竖向应力σr,sc表达式为：

式中，β＝6.08(p2+p3)/f′sc-3.49，残余应力fsc,re：a＝795.7-3.291f′sc，|fsc,re|≤0.25|f′sc,co|；

(3-2)外部夹层混凝土的环向应变-竖向应变关系

外部夹层混凝土的环向膨胀变形εθ,sc包括为弹性部分和塑性部分

弹性部分表达式：

塑性部分表达式：

式中，rsc为外部管的内径，εsc,z0为夹层混凝土开裂时的应变：

(4)外部管的受力和变形分析

外部管环向应力σθ,os表达式：

外部管竖向应力σz,os表达式：

外部管环向应变εθ,os表达式：

(5)变形协调

根据内部混凝土和内部管环向变形协调，得到变形协调方程如下：

εθ,ic＝εθ,is(式19)

根据夹层混凝土和外部管环向变形协调，得到变形协调方程如下：

εθ,sc＝εθ,os(式20)

步骤3.将步骤1获取到的计算用参数带入步骤2的公式1至20中求解，得到内部管和内部核心混凝土之间的挤压应力p1，内部管和外部夹层混凝土之间的挤压应力p2，外部夹层混凝土和外部管之间的挤压应力p3；

步骤4.基于步骤1获取到的计算用参数和步骤3求解得到的结果求解承载力

试件承担的荷载为：

n＝k1(nic+nis+nsc+nos)(式21)

式中，k1为折减系数，nic为内部核心混凝土承担的纵向荷载nic＝aicσz,ic，aic为内部核心混凝土截面面积，nis为内部管承担的纵向荷载nis＝aisσz,is，ais为内部管截面面积，nsc为夹层混凝土承担的纵向荷载nsc＝ascσz,sc，asc为外部夹层混凝土截面面积，nos为外部管承担的纵向荷载，nos＝aosσz,os，aos为外部管截面面积，nos为外部管承担的纵向荷载：nos＝aosσz,os，aos为外部管截面面积；

将n对εz求导，并令求导结果等于0：

根据求导等于0得到的εz,peak，反代回公式21中，求得加载过程中的峰值荷载nu,m：

nu,m＝k1(nic,p+nis,p+nsc,p+nos,p)(式23)

式中，nic,p，nis,p，nsc,p，和nos,p分别为竖向应变等于εz,peak时内部核心混凝土、内部管、外部夹层混凝土和外部管承担的荷载。

<系统>

进一步，本发明还提供一种各向同性双管混凝土柱承载力计算系统，其特征在于，包括：

参数获取模块，获取双管混凝土柱的计算用参数，包括：内部核心混凝土和外部夹层混凝土的泊松比、弹性模量，内部管材和外部管材的泊松比、弹性模量、内径、壁厚；

内部核心混凝土受力变形分析模块，基于下式1对峰值荷载前阶段内部核心混凝土的竖向应力应变关系进行分析，基于下式2对峰值荷载后阶段内部核心混凝土的竖向应力应变关系进行分析，并且基于下式3至6对内部核心混凝土的环向应变-竖向应变关系进行分析：

式中，σz,ic为内部核心混凝土竖向应力，eic为内部核心混凝土的弹性模量，εz为纵向应变，f′ic,co＝f′ic+4.1p1，p1为内部管和内部核心混凝土之间的挤压应力，fc′为混凝土强度，β＝12.16p1/f′ic-3.49，残余应力fic,re：a＝795.7-3.291f′ic，k＝5.79(p1/f′ic)^0.694+1.301，|fic,re|≤0.25|f′ic,co|，εθ,ic为混凝土环向应变，νc为混凝土的泊松比，εic,z0为混凝土开裂时的应变，εic为素混凝土达到峰值荷载时的应变；

内部管受力变形分析模块，基于下式7对内部管的环向应力σθ,is进行分析，基于下式8对内部管的竖向应力σz,is进行分析，并且基于下式9对内部管的环向应变εθ,is进行分析：

式中，ric为内部管的内径，tis为内部管的壁厚，p2为内部管和外部夹层混凝土之间的挤压应力，fis为内部管的屈服强度，es为内部管的弹性模量，νs为内部管的泊松比；

外部夹层混凝土受力变形分析模块，基于下式10对峰值荷载前阶段外部夹层混凝土的竖向应力应变关系进行分析，基于下式11对峰值荷载后阶段外部夹层混凝土的竖向应力应变关系进行分析，并且基于下式12至15对外部夹层混凝土的环向应变-竖向应变关系进行分析：

式中，σr,sc为外部夹层混凝土的竖向应力，σr,sc为外部夹层混凝土竖向应力，εθ,sc为外部夹层混凝土的环向膨胀变形，是外部夹层混凝土的切线模量；εsc,co为约束状态下外部夹层混凝土受压峰值应力对应的应变：为外部夹层混凝土和外部管之间的挤压应力，f′sc为外部夹层混凝土单轴受压时的峰值强度，εsc为外部夹层混凝土单轴受压峰值应力对应的应变；为外部约束状态下夹层混凝土达到峰值应力时的割线模量；f′sc,co为约束状态下外部夹层混凝土受压时的峰值强度，f′sc,co＝f′sc+4.1(p2+p3)/2；β＝6.08(p2+p3)/f′sc-3.49，残余应力fsc,re：a＝795.7-3.291f′sc，|fsc,re|≤0.25|f′sc,co|；rsc为外部管的内径，εsc,z0为夹层混凝土开裂时的应变；

外部管受力变形分析模块，基于下式16对外部管的环向应力σθ,os进行分析，基于下式17对外部管的竖向应力σz,os进行分析，基于下式18对外部管的环向应变εθ,os进行分析：

挤压应力分析模块，基于下式19和20对挤压应力进行分析：

εθ,ic＝εθ,is(式19)

εθ,sc＝εθ,os(式20)

计算模块，与参数获取模块、内部核心混凝土受力变形分析模块、内部管受力变形分析模块、外部夹层混凝土受力变形分析模块、外部管受力变形分析模块和挤压应力分析模块均通信相连；将获取到的计算用参数带入式1至20中计算得到内部管和内部核心混凝土之间的挤压应力p1，内部管和外部夹层混凝土之间的挤压应力p2，外部夹层混凝土和外部管之间的挤压应力p3；进一步，将计算用参数和计算得到的挤压应力p1、p2、p3带入下式21至23计算得到承载力：

n＝k1(nic+nis+nsc+nos)(式21)

nu,m＝k1(nic,p+nis,p+nsc,p+nos,p)(式23)

式中，k1为折减系数，nic为内部核心混凝土承担的纵向荷载nic＝aicσz,ic，aic为内部核心混凝土截面面积，nis为内部管承担的纵向荷载nis＝aisσz,is，ais为内部管截面面积，nsc为夹层混凝土承担的纵向荷载nsc＝ascσz,sc，asc为外部夹层混凝土截面面积，nos为外部管承担的纵向荷载，nos＝aosσz,os，aos为外部管截面面积，nos为外部管承担的纵向荷载：nos＝aosσz,os，aos为外部管截面面积；nu,m为加载过程中的峰值荷载，nic,p，nis,p，nsc,p，和nos,p分别为竖向应变等于εz,peak时内部核心混凝土、内部管、外部夹层混凝土和外部管承担的荷载。

优选地，本发明提供的各向同性双管混凝土柱承载力计算系统，还可以包括：输入显示模块，与参数获取模块和计算模块通信相连，用于对获取的计算用参数和计算得到的结果进行显示；和控制模块，与参数获取模块、内部核心混凝土受力变形分析模块、内部管受力变形分析模块、外部夹层混凝土受力变形分析模块、外部管受力变形分析模块、挤压应力分析模块、计算模块和输入显示模块均通信相连，控制各模块的运行。

优选地，本发明提供的各向同性双管混凝土柱承载力计算系统，还可以包括：图像形成模块，与参数获取模块、计算模块和控制模块均通信相连，用于根据参数获取模块获取到的计算用参数生成相应的双管混凝土柱的图形，并在图形上对应位置处标示出计算用参数信息和计算模块计算得到的结果信息；其中，输入显示模块还用于对图像形成模块生成的图形和标示信息进行显示。

发明的作用与效果

本发明所提供的各向同性双管混凝土柱承载力计算方法和系统，对各组成部分在弹塑性阶段的受力和变形进行分析，并通过建立变形协调方程，给出试件承载力的计算方法，求解得到各组成部分之间径向挤压应力，进而求得各组成部分承担的荷载，可以更加真实有效地反映材料各组成部分无法同时达到各自峰值应力对于承载力计算结果造成的影响，结果更为客观和准确，这对于预测不同强度材料制作的试件的承载力意义重大。并且，本发明适用范围非常广泛，不仅可以适用于常见的双钢管混凝土(管材种类包括但不限于低碳钢、高强钢、不锈钢等)，对其他采用各向同性管材的双管混凝土结构(管材种类如铝合金或其他材料)，也同样适用；此外，还适用于截面形式类似、内部钢管混凝土没有受损的外套钢管夹层混凝土加固cfst柱。

附图说明

图1为本发明实施例一所涉及的双钢管混凝土柱各组成部分受力分析示意图，其中(a)对应内部核心混凝土，(b)对应内部钢管，(c)对应外部夹层混凝土，(d)对应外部钢管；

图2为本发明实施例一所涉及的内部核心混凝土的竖向应力-应变关系示意图；

图3为本发明实施例一所涉及的外部夹层混凝土的竖向应力-应变关系示意图；

图4为本发明实施例一所涉及的计算结果与试验承载力结果对比图。

具体实施方式

以下结合附图对本发明涉及的各向同性双管混凝土柱承载力计算方法及系统的具体实施方案进行详细地说明。

<实施例一>

本实施例中，各组成部分的受力分析如图1所示。并且，在本实施例中：

(1)内部钢管和外部钢管的泊松比相同，均为νs；内部钢管和外部钢管的弹性模量视为相同，均为es。

(2)在轴压荷载作用下，内部核心混凝土、内部钢管、外部夹层混凝土、外部钢管等各组成部分的纵向应变相同，均为εz。

(3)内部核心混凝土和外部夹层混凝土泊松比均设为vc。

(4)内部钢管与夹层混凝土之间的挤压应力p2和夹层混凝土与外部钢管之间的挤压应力p3相等，p2＝p3。

(5)试件达到承载力时，内外钢管均已进入屈服状态。

关于应力和应变的正负号问题，本方案采用弹塑性力学所规定的符号原则，以指向外法线方向的应力为正，与截面外法线方向相反的方向为负。因此，试件在轴压荷载下的竖向应变εz、竖向应力σz,ic均为负数，而各材料在轴压荷载下的竖向应力峰值大小等于各自竖向受压强度，符号为负。

本实施例所提供的各向同性双管混凝土柱承载力计算方法包括如下内容：

步骤1.获取双管混凝土柱的计算用参数，包括：内部核心混凝土和外部夹层混凝土的泊松比、弹性模量，内部管材和外部管材的泊松比、弹性模量、内径、壁厚；

步骤2.建立双管混凝土柱各部分之间受力和变形关系

(1)内部核心混凝土的受力和变形分析

关于内部核心混凝土在约束作用下的竖向应力-应变关系，采用二阶段模型，如图2所示的模型。

其中，峰值荷载前阶段，内部核心混凝土竖向应力σic,z的表达式如下：

其中，εic,co为约束状态下内部核心混凝土受压峰值应力对应的应变。-f′ic,co＝-f′ic+4.1σr,ic；σr,ic＝-p1；是内部核心混凝土的切线模量；为约束状态下内部核心混凝土达到峰值应力时的割线模量；εz为内部核心混凝土的纵向受压应变；εic为内部核心混凝土单轴受压峰值应力对应的应变，可假设为-0.0022；f′ic为无约束作用下内部核心混凝土单轴受压强度(正值)，-f′ic表示相应的峰值压应力，负号表示受压，以和弹塑性力学符号规定一致；f′ic,co为约束状态下内部核心混凝土受压强度(正值)，-f′ic,co表示相应的峰值压应力，负号表示受压，以和弹塑性力学符号规定一致。

因此，整理得：

其中，f′ic,co＝f′ic+4.1p1，

关于峰值荷载后阶段，混凝土竖向应力取为：

其中，β＝12.16p1/f′ic-3.49。残余应力fic,re取为如下：

其中，a＝795.7-3.291f′ic，k＝5.79(p1/f′ic)^0.694+1.301。|fic,re|≤0.25|f′ic,co|。

除竖向应力-应变关系外，为建立内部混凝土和内部钢管在环向的变形协调条件，还需确定在某一竖向应变下内部核心混凝土的环向应变，即需要得到内部核心混凝土的环向应变-竖向应变关系。核心混凝土和夹层混凝土在轴心受压作用下发生膨胀。影响混凝土环向应变的因素包括竖向应变εz，挤压应力p1、p2、p3，以及混凝土强度fc′。混凝土环向应变包括弹性部分和塑性部分：

其中，弹性部分的表达式参考弹性阶段的分析：

其中，塑性部分的表达式如下：

其中，εic,z0为混凝土开裂时的应变，

其中，εic为素混凝土达到峰值荷载时的应变，假设为-0.0022。

(2)内部钢管的受力和变形分析

内部钢管的平面受力分析如图1所示。建立受力平衡方程可得：

即

在得到内部钢管环向应力的表达式后，根据屈服状态方程，可以求得钢管的竖向应力。忽略径向挤压应力的影响，钢材在双向受力下的屈服状态方程为：

其中fis为内部钢管屈服强度。根据上式可以求得内部钢管竖向应力σz,is的表达式为：

其中，σθ,is表达式如公式(10)所示。

为了建立混凝土和钢管接触界面的环向变形协调关系，还需得到钢管的环向应变-竖向应变关系。当钢材在进入塑性状态之初，由于应力变化较小，可以采用全量理论进行钢管应力的分析，并得到如下关于钢管环向应力的表达式：

其中，对公式(13)整理可得钢管环向应变-竖向应变关系如下：

其中，σθ,is表达式如公式(10)所示。

(3)外部夹层混凝土的受力和变形分析

外部夹层混凝土的竖向应力-应变关系，采用和内部核心混凝土类似的二阶段模型，如图3所示。

在峰值荷载前，外部夹层混凝土的竖向应力应变关系表达式为：

其中，是夹层混凝土的切线模量；

为约束状态下外部夹层混凝土达到峰值应力时的割线模量；

εz为约束混凝土的纵向受压应变；

εsc,co为约束状态下外部夹层混凝土受压峰值应力对应的应变。将约束应力取为p2和p3的平均值，可以得到

f′sc,co为约束状态下外部夹层混凝土受压时的峰值强度，f′sc,co＝f′sc+4.1(p2+p3)/2；

εsc为外部夹层混凝土单轴受压峰值应力对应的应变，根据经验取为-0.0022；f′sc为外部夹层混凝土单轴受压时的峰值强度。

关于峰值荷载后，外部夹层混凝土竖向应力表达式取为：

其中β＝6.08(p2+p3)/f′sc-3.49。残余应力fsc,re取为如下：

其中，a＝795.7-3.291f′sc，规定|fsc,re|≤0.25|f′sc,co|。

关于外部夹层混凝土的环向膨胀变形，可以分解为弹性部分和塑性部分

其中，弹性部分，按照弹性阶段的分析可以求得：

而塑性部分则假设为如下：

其中，εsc,z0为外部夹层混凝土开裂时的应变。

(4)外部钢管应变分析

对于外部钢管，应力状态分析如图1所示，通过建立受力平衡方程可得：

化简得：

在得到外部钢管环向应力的表达式后，根据屈服状态方程，可以求得外部钢管的竖向应力。忽略径向挤压应力对于屈服状态方程的影响，假设为双向受力下的屈服状态方程：

其中，fos为外部钢管屈服强度。根据上式可以求得外部钢管竖向应力σz,os的表达式为：

其中，σθ,os表达式如公式(23)所示。

为了建立外部夹层混凝土和外部钢管接触界面的环向变形协调关系，还需得到外部钢管的环向应变-竖向应变关系。当钢材在进入塑性状态之初，由于应力变化较小，可以采用全量理论进行钢管应力的分析，并整理可以得到如下关于外部钢管环向应力的表达式：

其中，对公式(26)整理可得外部钢管环向应变-竖向应变关系如下：

其中，σθ,os表达式如公式(23)所示。

(5)变形协调及挤压应力求解

基于本文关于p2＝p3的假设，所需求解的未知参数仅为p1和p2。因此，仅需要建立两个变形协调条件即可。

根据内部混凝土和内部钢管环向变形协调，可以得到变形协调方程如下：

εθ,ic＝εθ,is(28)

其中，εθ,ic和εθ,is分别如公式(5)和(14)所示。

根据夹层混凝土和外部钢管环向变形协调，可以得到变形协调方程如下：

εθ,sc＝εθ,os(29)

其中，εθ,sc和εθ,os分别如公式(18)和(27)所示。

步骤3.将步骤1获取到的计算用参数带入步骤2的公式中求解，即可得到内部钢管和内部核心混凝土之间的挤压应力p1，内部钢管和外部夹层混凝土之间的挤压应力p2，外部夹层混凝土和外部钢管之间的挤压应力p3；

步骤4.基于步骤1获取到的计算用参数和步骤3求解得到的结果求解承载力

当纵向应变为εz时，试件承担的荷载为：

n＝k1(nic+nis+nsc+nos)(30)

其中，k1为折减系数。经计算结果与试验结果对比，取为1.19。

nic为内部核心混凝土承担的纵向荷载：

nic＝aicσz,ic(31)

其中，aic为内部核心混凝土截面面积，σz,ic如公式(2)和(3)所示。

nis为内部钢管承担的纵向荷载：

nis＝aisσz,is(32)

其中，ais为内部钢管截面面积，σz,is如公式(12)所示。

nsc为夹层混凝土承担的纵向荷载：

nsc＝ascσz,sc(33)

其中，asc为外部夹层混凝土截面面积，σz,sc如公式(15)和(16)所示。

nos为外部钢管承担的纵向荷载：

nos＝aosσz,os(34)

其中，aos为外部钢管截面面积，σz,os如公式(25)所示。

将n对εz求导，并令求导结果等于0：

根据求导等于0得到的εz,peak，反代回公式(30)中，即可求得加载过程中的峰值荷载nu,m：

nu,m＝k1(nic,p+nis,p+nsc,p+nos,p)(36)

其中，nic,p，nis,p，nsc,p，和nos,p分别为竖向应变等于εz,peak时内部核心混凝土、内部钢管、外部夹层混凝土和外部钢管承担的荷载。值得指出的是，在matlab编程计算过程中，还可以采用设置较小的步长和循环语句，计算不同εz下的结果，计算得到的最大值即为峰值荷载。

下表1给出了试件尺寸和材料强度的信息。图4展示了根据所建立的承载力计算模型得到的计算结果nu,m与试验结果nu,exp的对比结果。计算结果与试验结果比值的平均值为1.02，变异系数为2.7％，验证了计算模型的准确性。

表1

<实施例二>

本实施例二中提供一种各向同性双管混凝土柱承载力计算系统，能够实现上述<实施例一>中所描述的方法。

本实施例二所提供的各向同性双管混凝土柱承载力计算系统，包括参数获取模块、内部核心混凝土受力变形分析模块、内部管受力变形分析模块、外部夹层混凝土受力变形分析模块、外部管受力和变形分析模块、挤压应力分析模块、计算模块、图像形成模块、输入显示模块以及控制模块。

参数获取模块用于获取双管混凝土柱的计算用参数，计算用参数包括：内部核心混凝土和外部夹层混凝土的泊松比、弹性模量，内部管材和外部管材的泊松比、弹性模量、内径、壁厚等。

内部核心混凝土受力变形分析模块，基于下式1对峰值荷载前阶段内部核心混凝土的竖向应力应变关系进行分析，基于下式2对峰值荷载后阶段内部核心混凝土的竖向应力应变关系进行分析，并且基于下式3至6对内部核心混凝土的环向应变-竖向应变关系进行分析：

式中，σz,ic为内部核心混凝土竖向应力，eic为内部核心混凝土的弹性模量，εz为纵向应变，f′ic,co＝f′ic+4.1p1，p1为内部管和内部核心混凝土之间的挤压应力，fc′为混凝土强度，β＝12.16p1/f′ic-3.49，残余应力fic,re：a＝795.7-3.291f′ic，k＝5.79(p1/f′ic)^0.694+1.301，|fic,re|≤0.25|fic′,co|，εθ,ic为混凝土环向应变，νc为混凝土的泊松比，εic,z0为混凝土开裂时的应变，εic为素混凝土达到峰值荷载时的应变；

内部管受力变形分析模块，基于下式7对内部管的环向应力σθ,is进行分析，基于下式8对内部管的竖向应力σz,is进行分析，并且基于下式9对内部管的环向应变εθ,is进行分析：

式中，ric为内部管的内径，tis为内部管的壁厚，p2为内部管和外部夹层混凝土之间的挤压应力，fis为内部管的屈服强度，es为内部管的弹性模量，νs为内部管的泊松比；

外部夹层混凝土受力变形分析模块，基于下式10对峰值荷载前阶段外部夹层混凝土的竖向应力应变关系进行分析，基于下式11对峰值荷载后阶段外部夹层混凝土的竖向应力应变关系进行分析，并且基于下式12至15对外部夹层混凝土的环向应变-竖向应变关系进行分析：

式中，σr,sc为外部夹层混凝土的竖向应力，σr,sc为外部夹层混凝土竖向应力，εθ,sc为外部夹层混凝土的环向膨胀变形，是外部夹层混凝土的切线模量；εsc,co为约束状态下外部夹层混凝土受压峰值应力对应的应变：p3为外部夹层混凝土和外部管之间的挤压应力，f′sc为外部夹层混凝土单轴受压时的峰值强度，εsc为外部夹层混凝土单轴受压峰值应力对应的应变；为外部约束状态下夹层混凝土达到峰值应力时的割线模量；f′sc,co为约束状态下外部夹层混凝土受压时的峰值强度，f′sc,co＝f′sc+4.1(p2+p3)/2；β＝6.08(p2+p3)/f′sc-3.49，残余应力fsc,re：a＝795.7-3.291f′sc，|fsc,re|≤0.25|f′sc,co|；rsc为外部管的内径，εsc,z0为夹层混凝土开裂时的应变；

外部管受力变形分析模块，基于下式16对外部管的环向应力σθ,os进行分析，基于下式17对外部管的竖向应力σz,os进行分析，基于下式18对外部管的环向应变εθ,os进行分析：

挤压应力分析模块，基于下式19和20对挤压应力进行分析：

εθ,ic＝εθ,is(式19)

εθ,sc＝εθ,os(式20)

计算模块，与参数获取模块、内部核心混凝土受力变形分析模块、内部管受力变形分析模块、外部夹层混凝土受力变形分析模块、外部管受力变形分析模块和挤压应力分析模块均通信相连；将获取到的计算用参数带入式1至20中计算得到内部管和内部核心混凝土之间的挤压应力p1，内部管和外部夹层混凝土之间的挤压应力p2，外部夹层混凝土和外部管之间的挤压应力p3；进一步，将计算用参数和计算得到的挤压应力p1、p2、p3带入下式21至23计算得到承载力：

n＝k1(nic+nis+nsc+nos)(式21)

nu,m＝k1(nic,p+nis,p+nsc,p+nos,p)(式23)

式中，k1为折减系数，nic为内部核心混凝土承担的纵向荷载nic＝aicσz,ic，aic为内部核心混凝土截面面积，nis为内部管承担的纵向荷载nis＝aisσz,is，ais为内部管截面面积，nsc为夹层混凝土承担的纵向荷载nsc＝ascσz,sc，asc为外部夹层混凝土截面面积，nos为外部管承担的纵向荷载，nos＝aosσz,os，aos为外部管截面面积，nos为外部管承担的纵向荷载：nos＝aosσz,os，aos为外部管截面面积；nu,m为加载过程中的峰值荷载，nic,p，nis,p，nsc,p，和nos,p分别为竖向应变等于εz,peak时内部核心混凝土、内部管、外部夹层混凝土和外部管承担的荷载。

图像形成模块与参数获取模块、计算模块和控制模块均通信相连，根据参数获取模块获取到的计算用参数生成相应的双管混凝土柱弹性的图形，并在图形上对应位置处标示出计算用参数信息和计算模块计算得到的结果信息。

输入显示模块与参数获取模块和计算模块通信相连，用于让操作员输入指令信息，并基于指令信息对获取的计算用参数和计算得到的结果以及图像形成模块生成的图形和标示信息进行显示。

控制模块与参数获取模块、内部核心混凝土受力变形分析模块、内部管受力变形分析模块、外部夹层混凝土受力变形分析模块、外部管受力和变形分析模块、挤压应力分析模块、计算模块、图像形成模块和输入显示模块均通信相连，控制各模块的运行。

以上实施例仅仅是对本发明技术方案所做的举例说明。本发明所涉及的各向同性双管混凝土柱承载力计算方法及系统并不仅仅限定于在以上实施例中所描述的内容，而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换，都在本发明的权利要求所要求保护的范围内。