opensees妯℃嫙锛氭湁绮樼粨瀹炰綋妗ュⅸ鍒嗘瀽绠椾緥 - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

# Written by wyy 

# Units: N, m, sec 

wipe;

 #wipe memery for new calculation

 

model basic -ndm 2 -ndf 3;

 #2-d problem, 3-degree of freedom per node

 

## Define constant 

set

 PI [

expr

 2*asin(1.0)]   

 #

绛変簬

蟺

 

set

 g 9.81; 

 

## Data input 

set

 hCol 0.36;

 # Height of Column 

 

set

 bCol 0.50;

 # Width of Column  

 

set

 colAspect 4.86;

 # Column aspect ratio:  L/D  

鏌卞瓙鐨勭旱妯�姣�

 

set

 v [

expr

 0.1];

 # Column axial load rate --> P = v*f'c*Ag  

杞村悜鑽疯浇鐜?/p>

 

set

 LCol [

expr

 $hCol*$colAspect];

 # Column length  

鏌遍暱

 

set

 fc  -4.16E7;

 # Concrete compressive strength

 

set

 fpc -5.747E+7  

set

 epsc0  -0.006328 

set

 fpcu -3.302E+7 

set

 epsu -0.0215 

#set ratio 0.3 

#set ft 5.474E+5 

set

 Ec 2.43E10 

set

 colArea [

expr

 $hCol*$bCol];

 # Column cross-section area 

妯�鎴�闈㈢Н

 

set

 P [

expr

 2*$v*$fc*$colArea];

 # Column axial load  

杞村悜鑽疯浇

 

#set P [expr -576000-691484]  

set

 m [

expr

 -$P/$g];

 # Beam nodal mass  

鑺傜偣璐ㄩ噺

 

set

 beamArea $colArea;

 # Beam cross-section area

 

set

 GJ 1.0e15;

 # Torsional stiffness  

鎵�杞�鍒氬害

 

set

 np 3;

 # Number of integration points

 

 

## Define uniaxial materials 

# Core concrete (confined

绾︽潫娣峰嚌鍦?/p>

) 

uniaxialMaterial

 

Concrete01

 1 [

expr

 $fpc] $epsc0 $fpcu $epsu  

#uniaxialMaterial Concrete02 1 [expr $fpc]  $epsc0 $fpcu $epsu $ratio 

$ft $Ets 

# Cover concrete (unconfined

鏃犵害鏉?/p>

) 

uniaxialMaterial

 

Concrete01

 2 [

expr

 $fc] -0.002 0 -0.005  

# STEEL 

## Cover steel Tube   

#set fy 3.67e8;      # Yield stress    

#set E 2.0E11;    # Young's modulus E=200000000kN/m 

#                        tag  fy E0    b 

#uniaxialMaterial Steel01  3  $fy $E 0.01 

#   STEEL  

set

    fy  3.2E8 

set

    fu  4.59E8 

set

    Es  2.0E11 

set

   Esh  6E9             ;

 # 0.03*Es

 

set

   esh  0.016           ;

 # 8*ey

 

set

   eult  0.08           

;# 40*ey

 

set

   lsr  6.0 

set

   beta  1.0 

set

   r   0.8