pdsqian0319.mac

钢管桁架随机缺陷模态分析

!定义参数
DIA355=0.35426
DIA219=0.21866
THICKNESS16=0.01587
THICKNESS12=0.011483
THICKNESS10=0.010213
KZGD=5.3990
QG4=1.2484
QG3=1.5523
QG2=1.7989
YS=405.227e6
RAD355= (DIA355)/2
INRAD12=RAD355-THICKNESS12
INRAD16=RAD355-THICKNESS16
RAD219= (DIA219)/2
INRAD10=RAD219-THICKNESS10

/PREP7
ET,1,BEAM188
ET,2,LINK10
ET,3,LINK8

SECTYPE,   1, BEAM, CTUBE, , 0  
SECOFFSET, CENT 
SECDATA,INRAD12,RAD355,0,0,0,0,0,0,0,0
SECTYPE,   2, BEAM, CTUBE, , 0  
SECOFFSET, CENT 
SECDATA,INRAD16,RAD355,0,0,0,0,0,0,0,0
SECTYPE,   3, BEAM, CTUBE, , 0  
SECOFFSET, CENT 
SECDATA,INRAD10,RAD219,0,0,0,0,0,0,0,0

R,4,0.000314,0,
R,5,0.03141,,

!弦杆D355*12的材料
MPTEMP,1,0  
MPDATA,DENS,1,,7800  
MPDATA,EX,1,,2.06e11 
MPDATA,PRXY,1,,0.3  
TB,BISO,1,1,2,
TBDATA,,YS,0,,,,
TBPLOT，BIOS，1

!弦杆D355*16的材料
MPTEMP,1,0  
MPDATA,DENS,2,,7800  
MPDATA,EX,2,,2.06e11 
MPDATA,PRXY,2,,0.3  
TB,BISO,2,1,2,
TBDATA,,YS,0,,,,
TBPLOT，BIOS，2

!腹杆D219*10的材料
MPTEMP,1,0  
MPDATA,DENS,3,,7800  
MPDATA,EX,3,,2.06e11 
MPDATA,PRXY,3,,0.3  
TB,BISO,3,1,2,
TBDATA,,YS,0,,,,
TBPLOT，BIOS，3

!拉杆及耳板的材料
MPTEMP,1,0  
MPDATA,DENS,4,,7800 
MPDATA,EX,4,,2.06e11 
MPDATA,PRXY,4,,0.3  

K,1,3,0,0,,,
K,2,6,0,0.958000000000,,,
K,3,12,0,2.57900000000,,,
K,4,18,0,3.82200000000,,,
K,5,24,0,4.70100000000,,,
K,6,30,0,5.22400000000,,,
K,7,36,0,KZGD,,,
K,8,42,0,5.22400000000,,,
K,9,48,0,4.70100000000,,,
K,10,54,0,3.82200000000,,,
K,11,60,0,2.57900000000,,,
K,12,66,0,0.958000000000,,,
K,13,69,0,0,,,
K,14,6,0.443000000000,0.269400000000,,,
K,15,12,1.32800000000,0.808100000000,,,
K,16,18,1.94450000000,QG4,,,
K,17,24,2.40400000000,QG3,,,
K,18,30,2.67150000000,1.73920000000,,,
K,19,36,2.76200000000,QG2,,,
K,20,42,2.67150000000,1.73920000000,,,
K,21,48,2.40400000000, QG3,,,
K,22,54,1.94450000000, QG4,,,
K,23,60,1.32800000000,0.808100000000,,,
K,24,66,0.442700000000,0.269400000000,,,
K,25,3,6,0,,,
K,26,6,5.55730000000,0.269400000000,,,
K,27,12,4.67200000000,0.808100000000,,,
K,28,18,4.05550000000, QG4,,,
K,29,24,3.59600000000, QG3,,,
K,30,30,3.32850000000,1.73920000000,,,
K,31,36,3.23800000000,QG2,,,
K,32,42,3.32850000000,1.73920000000,,,
K,33,48,3.59600000000, QG3,,,
K,34,54,4.05550000000, QG4,,,
K,35,60,4.67200000000,0.808100000000,,,
K,36,66,5.55730000000,0.269400000000,,,
K,37,69,6,0,,,
K,38,6,6,0.958000000000,,,
K,39,12,6,2.57900000000,,,
K,40,18,6,3.82200000000,,,
K,41,24,6,4.70100000000,,,
K,42,30,6,5.22400000000,,,
K,43,36,6,KZGD,,,
K,44,42,6,5.22400000000,,,
K,45,48,6,4.70100000000,,,
K,46,54,6,3.82200000000,,,
K,47,60,6,2.57900000000,,,
K,48,66,6,0.958000000000,,,
K,49,15.50000,0.000000,3.304083,,,
K,50,15.50000,6.000000,3.304083,,,
K,51,27.50000,0.000000,5.006083,,,
K,52,27.50000,6.000000,5.006083,,,
K,53,39.50000,0.000000,5.296500,,,
K,54,39.50000,6.000000,5.296500,,,
K,55,51.50000,0.000000,4.188250,,,
K,56,51.50000,6.000000,4.188250,,,
K,57,61.00000,0.000000,2.308833,,,
K,58,61.00000,6.000000,2.308833,,,
K,59,5.388900,0.3527609,0.2145232,,,
K,60,5.388900,5.647478,0.2145232,,,
K,61,14.00000,1.533500,0.9551000,,,
K,62,14.00000,4.466500,0.9551000,,,
K,63,23.00000,2.327417,1.502933,,,
K,64,23.00000,3.672583,1.502933,,,
K,65,33.00000,2.716750,1.768600,,,
K,66,33.00000,3.283250,1.768600,,,
K,67,45.00000,2.537750,1.646450,,,
K,68,45.00000,3.462250,1.646450,,,
K,69,57.00000,1.636250,1.028600,,,
K,70,57.00000,4.363750,1.028600,,,
K,71,66.61110,0.352522,0.2145232,,,
K,72,66.61110,5.647478,0.2145232,,,
N,3,12,0,2.57900000000,,,
N,4,18,0,3.82200000000,,,
N,5,24,0,4.70100000000,,,
N,6,30,0,5.22400000000,,,
N,7,36,0,KZGD,,,
N,8,42,0,5.22400000000,,,
N,9,48,0,4.70100000000,,,
N,10,54,0,3.82200000000,,,
N,11,60,0,2.57900000000,,,
N,39,12,6,2.57900000000,,,
N,40,18,6,3.82200000000,,,
N,41,24,6,4.70100000000,,,
N,42,30,6,5.22400000000,,,
N,43,36,6, KZGD,,,
N,44,42,6,5.22400000000,,,
N,45,48,6,4.70100000000,,,
N,46,54,6,3.82200000000,,,
N,47,60,6,2.57900000000,,,
N,16,18,1.94450000000,QG4,,,
N,17,24,2.40400000000,QG3,,,
N,18,30,2.67150000000,1.73920000000,,,
N,19,36,2.76200000000,QG2,,,
N,20,42,2.67150000000,1.73920000000,,,
N,21,48,2.40400000000, QG3,,,
N,22,54,1.94450000000, QG4,,,
N,28,18,4.05550000000, QG4,,,
N,29,24,3.59600000000, QG3,,,
N,30,30,3.32850000000,1.73920000000,,,
N,31,36,3.23800000000,QG2,,,
N,32,42,3.32850000000,1.73920000000,,,
N,33,48,3.59600000000, QG3,,,
N,34,54,4.05550000000, QG4,,,

!D355*16
L,1,59
L,59,14
L,1,2
L,25,60
L,60,26
L,12,13
L,25,38
L,24,71
L,71,13
L,36,72
L,72,37
L,48,37

!D355*12前上弦
L,2,3
L,3,49
L,49,4
L,4,5
L,5,51
L,51,6
L,6,7
L,7,53
L,53,8
L,8,9
L,9,55
L,55,10
L,10,11
L,11,57
L,57,12

!D355*12前下弦
L,14,15
L,15,61
L,61,16
L,16,63
L,63,17
L,17,18
L,18,65
L,65,19
L,19,20
L,20,67
L,67,21
L,21,22
L,22,69
L,69,23
L,23,24

!D355*12后下弦
L,26,27
L,27,62
L,62,28
L,28,64
L,64,29
L,29,30
L,30,66
L,66,31
L,31,32
L,32,68
L,68,33
L,33,34
L,34,70
L,70,35
L,35,36

!D355*12后上弦
L,38,39
L,39,50
L,50,40
L,40,41
L,41,52
L,52,42
L,42,43
L,43,54
L,54,44
L,44,45
L,45,56
L,56,46
L,46,47
L,47,58
L,58,48

!腹杆
L,3,15
L,4,16
L,5,17
L,6,18
L,7,19
L,8,20
L,9,21
L,10,22
L,11,23
L,3,39
L,4,40
L,5,41
L,6,42
L,7,43
L,8,44
L,9,45
L,10,46
L,11,47
L,15,27
L,16,28
L,17,29
L,18,30
L,19,31
L,20,32
L,21,33
L,22,34
L,23,35
L,27,39
L,28,40
L,29,41
L,30,42
L,31,43
L,32,44
L,33,45
L,34,46
L,35,47

!两头腹杆
L,1,25
L,13,37

!拉杆
L,3,16
L,4,17
L,5,18
L,6,19
L,8,19
L,9,20
L,10,21
L,11,22
L,39,28
L,40,29
L,41,30
L,42,31
L,44,31
L,45,32
L,46,33
L,47,34

!类似耳板
L,2,14
L,26,38
L,12,24
L,36,48

!划分弦杆D355*12
TYPE,1   
MAT,1
REAL,   
ESYS,0   
SECNUM,1
LMESH,	13
LMESH,	14
LMESH,	15
LMESH,	16
LMESH,	17
LMESH,	18
LMESH,	19
LMESH,	20
LMESH,	21
LMESH,	22
LMESH,	23
LMESH,	24
LMESH,	25
LMESH,	26
LMESH,	27
LMESH,	28
LMESH,	29
LMESH,	30
LMESH,	31
LMESH,	32
LMESH,	33
LMESH,	34
LMESH,	35
LMESH,	36
LMESH,	37
LMESH,	38
LMESH,	39
LMESH,	40
LMESH,	41
LMESH,	42
LMESH,	43
LMESH,	44
LMESH,	45
LMESH,	46
LMESH,	47
LMESH,	48
LMESH,	49
LMESH,	50
LMESH,	51
LMESH,	52
LMESH,	53
LMESH,	54
LMESH,	55
LMESH,	56
LMESH,	57
LMESH,	58
LMESH,	59
LMESH,	60
LMESH,	61
LMESH,	62
LMESH,	63
LMESH,	64
LMESH,	65
LMESH,	66
LMESH,	67
LMESH,	68
LMESH,	69
LMESH,	70
LMESH,	71
LMESH,	72

!划分弦杆D355*16
TYPE,1   
MAT,2
REAL,  
ESYS,0   
SECNUM,2
LMESH,	1
LMESH,	2
LMESH,	3
LMESH,	4
LMESH,	5
LMESH,	6
LMESH,	7
LMESH,	8
LMESH,	9
LMESH,	10
LMESH,	11
LMESH,	12

!划分腹杆D219*10
TYPE,1   
MAT,3
REAL,   
ESYS,0   
SECNUM,3
LMESH,	73
LMESH,	74
LMESH,	75
LMESH,	76
LMESH,	77
LMESH,	78
LMESH,	79
LMESH,	80
LMESH,	81
LMESH,	82
LMESH,	83
LMESH,	84
LMESH,	85
LMESH,	86
LMESH,	87
LMESH,	88
LMESH,	89
LMESH,	90
LMESH,	91
LMESH,	92
LMESH,	93
LMESH,	94
LMESH,	95
LMESH,	96
LMESH,	97
LMESH,	98
LMESH,	99
LMESH,	100
LMESH,	101
LMESH,	102
LMESH,	103
LMESH,	104
LMESH,	105
LMESH,	106
LMESH,	107
LMESH,	108
LMESH,	109
LMESH,	110

!划分拉杆
TYPE,2
MAT,4
REAL,4
ESYS,0
SECNUM,
E,3,16
E,4,17
E,5,18
E,6,19
E,8,19
E,9,20
E,10,21
E,11,22
E,39,28
E,40,29
E,41,30
E,42,31
E,44,31
E,45,32
E,46,33
E,47,34

!划分耳板
TYPE,3
MAT,4
REAL,5
ESYS,0
SECNUM,
ESIZE,0,1,
LMESH,	127
LMESH,	128
LMESH,	129
LMESH,	130
NUMMRG,NODE,,,,LOW
NUMCMP,NODE
FINISH


!静力分析
/SOLU
ANTYPE,STATIC
D,429,UX,0,,,,UY,UZ,ROTX,ROTZ
D,442,UX,0,,,,UY,UZ,ROTX,ROTZ
D,450,UY,0,,,,UZ,ROTX,ROTZ
D,471,UY,0,,,,UZ,ROTX,ROTZ
F,429,FZ,-19660.67	
F,1,FZ,-32528.21
F,2,FZ,-25602.03
F,3,FZ,-25382.27
F,4,FZ,-25253.07
F,5,FZ,-25210.59
F,6,FZ,-25253.07
F,7,FZ,-25382.27
F,8,FZ,-25602.03
F,9,FZ,-32528.21
F,450,FZ,-19660.67
F,442,FZ,-19660.67
F,24,FZ,-32528.21
F,25,FZ,-25602.03
F,26,FZ,-25382.27
F,27,FZ,-25253.07
F,28,FZ,-25210.59
F,29,FZ,-25253.07
F,30,FZ,-25382.27
F,31,FZ,-25602.03
F,32,FZ,-32528.21
F,471,FZ,-19660.67
ACEL,0,0,9.8,
PSTRES,ON
NLGEOM,OFF
NROPT,FULL,,
EQSLV,SPAR,,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
TOFFST,0,
NSUBST,500,0,0
Solve
FINISH





!提取结果数并处存到参数中
/POST1
SET,FIRST
NSORT,U,Z
*GET,DMAX,SORT,,MIN
*GET,D5,NODE,5,U,Z
*GET,D28,NODE,5,U,Z
*GET,D2,NODE,2,U,Z
*GET,D8,NODE,8,U,Z
*GET,D25,NODE,25,U,Z
*GET,D31,NODE,31,U,Z
FINISH