钢管混凝土梁柱节点承载力分析算得力和位移太小了

andyqiujie

我正在做钢管混凝土梁柱节点承载力分析，考虑接触，在计算柱子的极限承载力时，老是算得位移和应力很小，一般只有10e-13左右，但是我加的面力已经很大了，开始以为是加弹簧单元的问题，后来我把节点合并了，还是算得那么小，我把命令流贴下，希望各位高手帮我指点一下错误，我算了一个多月了，就没整明白怎么错了，急死了！
/prep7  
et,1,solid65
!*  
R,1, , , , , , ,
RMORE, , , , , , ,  
!*  
!*  
et,2,shell143   !定义钢管
R,2,5, , , , ,  
RMORE, , ,  
!*  
!*  
R,3,11, , , , ,
RMORE, , ,  
!*  
!*  
ET,3,PIPE20 !定义钢筋
R,4,10,4.99,
!*  
!*  
!*  
R,5,6,2.99,
!*
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,3e4
MPDATA,PRXY,1,,0.2  
TB,MELA,1,1,8,  
TBTEMP,0
TBPT,,0.0005,15
TBPT,,0.0010,21.37
TBPT,,0.0015,25.095
TBPT,,0.0020,26
TBPT,,0.0025,25.4496
TBPT,,0.0033,23.404
TB,CONC,1,1,9,  
TBTEMP,0
TBDATA,,0.5,0.9,2.01,-1,,   
TBDATA,,,,,,,   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,2.1E5  
MPDATA,PRXY,2,,0.3  
TB,BKIN,2,1,2,1
TBTEMP,0
TBDATA,,235,,,,,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,3,,2.1E5  
MPDATA,PRXY,3,,0.3  
TB,BKIN,3,1,2,1
TBTEMP,0
TBDATA,,300,,,,,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,4,,2.0E5
MPDATA,PRXY,4,,0.25
TB,BKIN,4,1,2,1
TBTEMP,0
TBDATA,,210,,,,,
ET,4,COMBIN39   !定义钢管与混凝土法向弹簧单元
R,6,-3,-90000,-2.5,-75000,-2,-60000,
RMORE,-1.5,-45000,-1,-30000,-0.5,-15000,
RMORE,-0.3,-9000,-0.1,-3000,0,0,
RMORE,0.1,0.001,0.2,0.001,0.3,0.001,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,   
!*  
ET,5,COMBIN39 !定义钢管与混凝土切向弹簧单元
R,7,0,0,0.4,304,0.8,608,  
RMORE,1.2,912,1.5985,1214.86,1.9182,1166.266,   
RMORE,2,1153.832,2.4,1093.032,2.8,1032.232,
RMORE,3.2,971.432,3.6,910.632,4.2,819.432,
RMORE,4.7955,728.916,5,728.916,,,
RMORE,,,,,,,
RMORE,,,,,   
ET,6,COMBIN39 !定义钢管与混凝土横向弹簧单元
R,8,0,0,0.4,304,0.8,608,  
RMORE,1.2,912,1.5985,1214.86,1.9182,1166.266,   
RMORE,2,1153.832,2.4,1093.032,2.8,1032.232,
RMORE,3.2,971.432,3.6,910.632,4.2,819.432,
RMORE,4.7955,728.916,5,728.916,,,
RMORE,,,,,,,
RMORE,,,,,   
!*  
ET,7,COMBIN39  !定义钢筋与混凝土法向弹簧单元
R,9,0.008,8084.83,0.01,9214.84,0.012,10101.67,   
RMORE,0.015,11075.8,0.018,11737.57,0.02,12056.9,
RMORE,0.022,12306.35,0.026,12658.1,0.03,12859.15,
RMORE,0.035,12847.72,0.04,12291.29,0.05,7359.76,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,  
R,10,0.008,10252.34,0.01,11685.31,0.015,14045.18,   
RMORE,0.018,14884.38,0.02,15289.31,0.022,15605.64,
RMORE,0.026,16051.7,0.03,16306.64,0.035,16292.15,
RMORE,0.04,15586.55,0.05,9332.89,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,  
R,11,0.008,4335.03,0.01,4940.93,0.015,5938.77,   
RMORE,0.018,6293.61,0.02,6464.82,0.022,6598.58,
RMORE,0.026,6787.19,0.03,6894.98,0.035,6888.86,  
RMORE,0.04,6590.51,0.05,3946.25,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,
R,12,0.008,6502.54,0.01,7411.4,0.015,8908.15,   
RMORE,0.018,9440.41,0.02,9697.24,0.022,9897.87,
RMORE,0.026,10180.78,0.03,10342.48,0.035,10333.29,  
RMORE,0.04,9885.76,0.05,5919.38,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,
R,13,0.008,8670.05,0.01,9881.87,0.015,11877.53,   
RMORE,0.018,12587.21,0.02,12929.65,0.022,13197.16,
RMORE,0.026,13574.38,0.03,13789.97,0.035,13777.72,  
RMORE,0.04,13181.01,0.05,7892.51,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,0,0,  
RMORE,0,0,0,0,   
ET,8,COMBIN39  !定义钢筋与混凝土横向弹簧单元
R,14,-0.03,-600,-0.026,-520,-0.022,-440,  
RMORE,-0.02,-400,-0.018,-360,-0.015,-300,
RMORE,-0.01,-200,-0.008,-160,0,0,  
RMORE,0.008,0.001,0.01,0.001,0.03,0.001,
RMORE,,,,,,,
RMORE,,,,,,,
RMORE,,,,,
bopt,numb,off
!*  
KEYOPT,4,1,0 !设置各弹簧单元的keyoption
KEYOPT,4,2,0
KEYOPT,4,3,2
KEYOPT,4,4,0
KEYOPT,4,6,0
!*  
KEYOPT,5,1,0
KEYOPT,5,2,0
KEYOPT,5,3,3
KEYOPT,5,4,0
KEYOPT,5,6,0
!*  
KEYOPT,6,1,0
KEYOPT,6,2,0
KEYOPT,6,3,0
KEYOPT,6,4,0
KEYOPT,6,6,0
!*  
KEYOPT,7,1,0
KEYOPT,7,2,0
KEYOPT,7,3,2
KEYOPT,7,4,0
KEYOPT,7,6,0
!*
KEYOPT,8,1,0
KEYOPT,8,2,0
KEYOPT,8,3,3
KEYOPT,8,4,0
KEYOPT,8,6,0   
K,1,0,0,0,   !建立钢管模型
K,2,0,0,1100,   
K,3,136.5,0,0,  
K,4,136.5,0,1100,  
LSTR,       3,       4  
FLST,2,1,4,ORDE,1   
FITEM,2,1   
FLST,8,2,3  
FITEM,8,1   
FITEM,8,2   
AROTAT,P51X, , , , , ,P51X, ,180, ,
BLOCK,-65,65,0,200,476,516
BLOCK,-65,65,0,200,584,624  
APLOT  
FLST,5,2,5,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
ASEL,S, , ,P51X
ASBV,all,1,,DELETE,DELETE   
ASBV,all,2,,DELETE,DELETE
wpoff, , ,450   !切割钢管模型
asbw,all      
wpoff, , ,26
asbw,all
wpoff, , ,40
asbw,all
wpoff, , ,22
asbw,all
wpoff, , ,22
asbw,all
wpoff, , ,24
asbw,all
wpoff, , ,40
asbw,all                                    
wpoff, , ,26
asbw,all
wprot, , ,90  !旋转工作平面 ,,绕z轴旋转逆时针为正  
wpoff, , ,65
asbw,all
wpoff, , ,25
asbw,all
wpoff, , ,-65   
asbw,all
wpoff, , ,-50   
asbw,all
wpoff, , ,-40   
asbw,all
wpoff, , ,-25   
asbw,all
APLOT   
WPAVE,0,0,0  
wprot, , ,-90  
/VIEW,1,1,1,1  
/ANG,1  
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/USER,  1   
/VIEW,  1,  0.239133333941    ,  0.895762107447    ,  0.374733899534
/ANG,   1,  -9.77832658122  
/REPLO  
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/VIEW,  1,  0.112808728823    ,  0.966372814528    ,  0.231079583786
/ANG,   1,  -15.8625566834  
/REPLO  
/FOC,   1,  -19.5924332489    ,   82.4127661952    ,   500.336093946
/REPLO  
FLST,5,40,5,ORDE,27
FITEM,5,1   
FITEM,5,3   
FITEM,5,13  
FITEM,5,15  
FITEM,5,17  
FITEM,5,-18
FITEM,5,23  
FITEM,5,-25
FITEM,5,27  
FITEM,5,-33
FITEM,5,36  
FITEM,5,-40
FITEM,5,42  
FITEM,5,44  
FITEM,5,-45
FITEM,5,47  
FITEM,5,51  
FITEM,5,53  
FITEM,5,55  
FITEM,5,57  
FITEM,5,-60
FITEM,5,62  
FITEM,5,-63
FITEM,5,65  
FITEM,5,-67
FITEM,5,70  
FITEM,5,-72
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2,       3,   2,       0,   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/VIEW,  1,  0.438347630648    ,  0.883670022155    ,  0.164251778225
/ANG,   1,   25.2712901221  
/REPLO  
/VIEW,  1,  0.117659251263    ,  0.970054785760    ,  0.212485324710
/ANG,   1,  -12.7934857285  
/REPLO  
/VIEW,  1, -0.225776211131    ,  0.967233611042    ,  0.116121678241
/ANG,   1,  -102.391347163  
/REPLO  
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/FOC,   1,  -32.7996522297    ,   68.3635755172    ,   591.679706364
/REPLO  
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/FOC,   1,   4.95206866618    ,   64.5225299170    ,   697.074574049
/REPLO  
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/VIEW,  1, -0.544369262615E-01,  0.998181935969    , -0.258736113363E-01
/ANG,   1,  -154.676889784  
/REPLO  
/FOC,   1,   7.43062351300    ,   62.2889657372    ,   605.690797204
/REPLO  
/FOC,   1,   1.93456278278    ,   58.8742916952    ,   485.519041789
/REPLO  
/FOC,   1,   11.2142932218    ,   56.1739057909    ,   361.816308923
/REPLO  
/REPLO  
/FOC,   1,   15.9244116780    ,   57.7871304850    ,   414.143267060
/REPLO  
/FOC,   1,   15.8653851413    ,   57.7990327265    ,   414.726634628
/REPLO  
/REPLO  
/REPLO  
SAVE
/REPLOT,RESIZE  
/REPLOT,RESIZE  
SAVE
/FOC,   1,   14.8340689529    ,   59.1377175798    ,   468.541801893
/REPLO  
FLST,5,32,5,ORDE,28
FITEM,5,1   
FITEM,5,5   
FITEM,5,7   
FITEM,5,9   
FITEM,5,-11
FITEM,5,13  
FITEM,5,19  
FITEM,5,21  
FITEM,5,-22
FITEM,5,24  
FITEM,5,26  
FITEM,5,28  
FITEM,5,-29
FITEM,5,34  
FITEM,5,-36
FITEM,5,41  
FITEM,5,43  
FITEM,5,46  
FITEM,5,48  
FITEM,5,-50
FITEM,5,52  
FITEM,5,-54
FITEM,5,57  
FITEM,5,61  
FITEM,5,64  
FITEM,5,68  
FITEM,5,-69
FITEM,5,73  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2,       2,   2,       0,    !赋予钢管模型属性
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
SAVE
CYLIND,0,136.5,0,1100,0,180,     !建立管内混凝土模型和混凝土梁
BLOCK,-90,90,0,1336.500,450,650
vsba,2,6     
VDELE,3, , ,1
APLOT   
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,4   
VGLUE,P51X  !粘结管内混凝土和混凝土梁交接面
wpoff, , ,450   !切割管内混凝土模型和混凝土梁
vsbw,all
wpoff, , ,26
vsbw,all
wpoff, , ,40
vsbw,all
wpoff, , ,22
vsbw,all
wpoff, , ,22
vsbw,all
wpoff, , ,24
vsbw,all
wpoff, , ,40
vsbw,all
wpof, , ,26
vsbw,all
KWPAVE,101
wprot, , ,-90   
vsbw,all
wpoff, , ,25
vsbw,all
wpoff, , ,40
vsbw,all
wpoff, , ,25
vsbw,all
wpoff , ,, 25   
vsbw,all
wpoff, , ,40
vsbw,all
wpoff, , ,25
vsbw,all
VPLOT
wprot, , 90,
wpoff, , ,25
vsbw,all
wpoff, , ,75
vsbw,all
wpoff, , ,100   
vsbw,all
*do,i,1,7,1
wpoffs, , ,100  
vsbw,all
*enddo  
wpoff, , ,100   
vsbw,all
wpoff, , ,50
vsbw,all
wpoff, , ,50
vsbw,all
wpoff, ,, 50
vsbw,all
vsel,all
VATT,       1,       1,   1,       0  !赋予混凝土属性
n,1,65,0,625     !开始建梁内钢筋
n,2,25,0,625
n,3,-25,0,625
n,4,-65,0,625
FLST,4,4,1,ORDE,2
FITEM,4,1
FITEM,4,-4  
NGEN,2,4,P51X, , , , ,-148,1,   
VPLOT
FLST,4,8,1,ORDE,2   
FITEM,4,1   
FITEM,4,-8  
NGEN,2,8,P51X, , , ,186.5, ,1,
FLST,4,8,1,ORDE,2
FITEM,4,9   
FITEM,4,-16
NGEN,4,8,P51X, , , ,50, ,1,               
FLST,4,8,1,ORDE,2   
FITEM,4,33  
FITEM,4,-40
NGEN,10,8,P51X, , , ,100, ,1,FLST,4,8,1,ORDE,2   
FLST,4,8,1,ORDE,2   
FITEM,4,105
FITEM,4,-112
NGEN,2,8,P51X, , , ,75, ,1,
type,3                     
real,5
mat,4      
*do,i,9,11,1   
e,i,i+1
*enddo
*do,i,13,15,1
e,i,i+1
*enddo
e,9,13              
e,12,16                                   
!      4   @50  
egen,4,8,1,8,1  
FLST,4,8,2,ORDE,2   
FITEM,4,25  
FITEM,4,-32
EGEN,10,8,P51X, , , , , , , , ,100, ,   
/prep7
type,3
real,4
mat,3
*do,i,1,105,8     
e,i,i+8
*enddo
*do,i,2,106,8
e,i,i+8
*enddo
*do,i,3,107,8
e,i,i+8
*enddo
*do,i,4,108,8
e,i,i+8
*enddo
*do,i,5,109,8
e,i,i+8
*enddo
*do,i,6,110,8
e,i,i+8
*enddo
*do,i,7,111,8
e,i,i+8
*enddo
*do,i,8,112,8
e,i,i+8
*enddo
/prep7
type,3
real,4
mat,3
e,112,120
e,113,117
e,114,118
e,115,119
e,116,120
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
esize,50
lesize,49,,,9  !划分体网格
lesize,56,,,9
lesize,678,,,4
lesize,160,,,4
lesize,90,,,4
lesize,310,,,4
vmesh,all   
mshape,0,3d  
mshkey,1
lesize,49,,,9 !划分面网格，生成与体网格重合的节点，好在上面加弹簧单元
lesize,56,,,9
lesize,678,,,4
lesize,160,,,4
lesize,90,,,4
lesize,310,,,4
amesh,all
save
/prep7     
type,4           !加钢管和混凝土法向弹簧单元
real,6
nsel,s,loc,z,50
nsel,a,loc,z,100
nsel,a,loc,z,150
nsel,a,loc,z,200
nsel,a,loc,z,250
nsel,a,loc,z,300
nsel,a,loc,z,350
nsel,a,loc,z,400
nsel,a,loc,z,450
nsel,a,loc,z,476
nsel,a,loc,z,516
nsel,a,loc,z,538
nsel,a,loc,z,560
nsel,a,loc,z,584
nsel,a,loc,z,624
nsel,a,loc,z,650
nsel,a,loc,z,700
nsel,a,loc,z,750
nsel,a,loc,z,800
nsel,a,loc,z,850
nsel,a,loc,z,900
nsel,a,loc,z,950
nsel,a,loc,z,1000
nsel,a,loc,z,1050
nsel,a,loc,z,1100
nsel,u,,,1,120
eintf
type,5        !加钢管和混凝土切向弹簧单元
real,7
nsel,s,loc,z,50
nsel,a,loc,z,100
nsel,a,loc,z,150
nsel,a,loc,z,200
nsel,a,loc,z,250
nsel,a,loc,z,300
nsel,a,loc,z,350
nsel,a,loc,z,400
nsel,a,loc,z,450
nsel,a,loc,z,476
nsel,a,loc,z,516
nsel,a,loc,z,538
nsel,a,loc,z,560
nsel,a,loc,z,584
nsel,a,loc,z,624
nsel,a,loc,z,650
nsel,a,loc,z,700
nsel,a,loc,z,750
nsel,a,loc,z,800
nsel,a,loc,z,850
nsel,a,loc,z,900
nsel,a,loc,z,950
nsel,a,loc,z,1000
nsel,a,loc,z,1050
nsel,a,loc,z,1100
nsel,u,,,1,120
eintf
type,6         !加钢管和混凝土横向弹簧单元
real,8
nsel,s,loc,z,50
nsel,a,loc,z,100
nsel,a,loc,z,150
nsel,a,loc,z,200
nsel,a,loc,z,250
nsel,a,loc,z,300
nsel,a,loc,z,350
nsel,a,loc,z,400
nsel,a,loc,z,450
nsel,a,loc,z,476
nsel,a,loc,z,516
nsel,a,loc,z,538
nsel,a,loc,z,560
nsel,a,loc,z,584
nsel,a,loc,z,624
nsel,a,loc,z,650
nsel,a,loc,z,700
nsel,a,loc,z,750
nsel,a,loc,z,800
nsel,a,loc,z,850
nsel,a,loc,z,900
nsel,a,loc,z,950
nsel,a,loc,z,1000
nsel,a,loc,z,1050
nsel,a,loc,z,1100
nsel,u,,,1,120
eintf
type,7    !加钢筋和混凝土弹簧单元
real,9
nsel,s,loc,y,0
nsel,a,loc,y,1311.5
eintf
type,7
real,10
nsel,s,loc,y,186.5
eintf
type,7
real,11
nsel,s,loc,y,236.5
nsel,a,loc,y,286.5
nsel,a,loc,y,1236.5
eintf
type,7
real,12
nsel,s,loc,y,336.5
eintf
type,7
real,13
nsel,s,loc,y,436.5
nsel,a,loc,y,536.5
nsel,a,loc,y,636.5
nsel,a,loc,y,736.5
nsel,a,loc,y,836.5
nsel,a,loc,y,936.5
nsel,a,loc,y,1036.5
nsel,a,loc,y,1136.5
eintf
type,8
real,14
nsel,s,loc,y,0
nsel,a,loc,y,186.5
nsel,a,loc,y,236.5
nsel,a,loc,y,286.5
nsel,a,loc,y,336.5
nsel,a,loc,y,436.5
nsel,a,loc,y,536.5
nsel,a,loc,y,636.5
nsel,a,loc,y,736.5
nsel,a,loc,y,836.5
nsel,a,loc,y,936.5
nsel,a,loc,y,1036.5
nsel,a,loc,y,1136.5
nsel,a,loc,y,1236.5
nsel,a,loc,y,1311.5
eintf
/SOLU   
nsel,s,loc,z,0  !加钢管底部约束
d,all,uz
d,all,ux
d,all,uy
allsel,all  
asel,s,loc,y,0  !加钢管对称面对称约束
da,all,symm
dtran   
allsel,all  
asel,s,loc,z,1100  !耦合钢管加载端z方向约束，使之未刚性面，减少应立集中
nsla,s,all，
cp,1,uz,all
d,all,ux
d,all,uy
allsel,all  
NSUBST,50,50,10
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
LNSRCH,1
TIME,1  
asel,s,loc,z,1100   
sfa,all,1,pres,-60    !加钢管顶面力
allsel,all  
cnvtol,f, ,0.005,2,0.01
/STATUS,SOLU
SOLVE

[ 本帖最后由 xinyuxf 于 2007-7-22 16:02 编辑 ]

andyqiujie

这是计算完得到的图，显示是一个颜色，力和位移都是10e-13级别的，而且钢管柱加载端有侧

移，我加的是z方向的轴向力啊，这么小的应力怎么就有x，y平面的侧移了呢？

kunull

楼主？这个问题解决没有？另外问一下，你combin39单元的F-D曲线如何确定的？