!1.3  超靜定拉壓桿的反力計算
! 本程序來源于邢靜忠等編著的《ANSYS分析實例與工程應用》，機械工業出版社（2004年）
FINISH
/CLEAR,START
! (1) 設置工程的選項, 進入前處理模塊,聲明單元類型,實常數和材料參數
/FILNAME, EX1.3
/PREP7
/TITLE, EX1.3, STATICALLY INDETERMINATE REACTION FORCE ANALYSIS
ANTYPE,STATIC
ET,1,LINK1
R,1,1
MP,EX,1,210E9
! (2) 定義節點
N,1
N,2,,0.4
N,3,,0.7
N,4,,1.0
! (3) 定義單元
E,1,2
EGEN,3,1,1
! (4) 定義位移約束和荷載
D,1,ALL,,,4,3
F,2,FY,-500
F,3,FY,-1000
FINISH
! (5) 進入求解模塊SOLU
/SOLU
OUTPR,BASIC,1
OUTPR,NLOAD,1
SOLVE
FINISH
! (6) 進入一般后處理模塊POST1
/POST1
NSEL,S,LOC,Y,1.0
FSUM
*GET,REAC_1,FSUM,,ITEM,FY
NSEL,S,LOC,Y,0
FSUM
*GET,REAC_2,FSUM,,ITEM,FY
! (7) 申明數組,輸出計算結果并比較
*DIM,LABEL,CHAR,2
*DIM,VALUE,,2,3
LABEL(1) = 'R1, N ','R2, N  '
*VFILL,VALUE(1,1),DATA,900,600
*VFILL,VALUE(1,2),DATA,ABS(REAC_1),ABS(REAC_2)
*VFILL,VALUE(1,3),DATA,ABS(REAC_1/900) ,ABS(REAC_2/600)
/OUT,EX1_3,out
/COM
/COM,------------------- EX1.3 RESULTS COMPARISON ---------------------
/COM,
/COM,         |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,'   ',F10.1,'  ',F10.1,'   ',1F5.3)
/COM,----------------------------------------------------------------
/OUT
FINISH
*LIST,EX1_3,out