ZYF

/COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0
/VERIFY,VM4
/PREP7
/TITLE, VM4, DEFLECTION OF A HINGED SUPPORT
C***     STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE 10, PROB. 2
L=15*12                         ! LENGTH OF BAR IN INCHES
*AFUN,DEG                       ! TRIG FUNCTIONS IN DEGREES 
THETA=30                        ! ANGLE TO BE USED TO CALCULATE A AND B
A=2*L*COS(THETA)                ! CALCULATED X LOCATION - NODE 3
B=L*SIN(THETA)                  ! CALCULATED Y LOCATION - NODE 2
ET,1,LINK1
R,1,.5
MP,EX,1,30E6
N,1
N,2,A/2,-B                      ! X LOCATION = A/2; A AND B AS ABOVE
N,3,A                           
E,1,2
E,2,3
D,1,ALL,,,3,2
F,2,FY,-5000
OUTPR,,1
FINISH
/SOLU    
SOLVE
FINISH
/POST1 
MID_NODE = NODE (A/2,-B,0 )
*GET,DISP,NODE,MID_NODE,U,Y 
LEFT_EL = ENEARN (MID_NODE) 
ETABLE,STRS,LS,1  
*GET,STRSS,ELEM,LEFT_EL,ETAB,STRS 

*DIM,LABEL,CHAR,2
*DIM,VALUE,,2,3
LABEL(1) = 'STRS_psi','DEF_in'
*VFILL,VALUE(1,1),DATA,10000,-0.120
*VFILL,VALUE(1,2),DATA,STRSS,DISP
*VFILL,VALUE(1,3),DATA,ABS(STRSS /10000 ) ,ABS( DISP /0.120 )
/OUT,vm4,vrt
/COM
/COM,------------------- VM4 RESULTS COMPARISON ---------------------
/COM,
/COM,         |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,'   ',F10.3,'  ',F10.3,'   ',1F5.3)
/COM,----------------------------------------------------------------
/OUT

FINISH
*LIST,vm4,vrt

finish
/clear
/COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0
/VERIFY,VM3
/PREP7
/TITLE, VM3, THERMALLY LOADED SUPPORT STRUCTURE
C***      STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE 30, PROB. 9
ANTYPE,STATIC             ! STATIC ANALYSIS
ET,1,LINK1
R,1,.1
MP,EX,1,16E6
MP,ALPX,1,92E-7           ! alpx, alpy, alpz 热膨胀系数？
MP,EX,2,30E6
MP,ALPX,2,70E-7
TREF,70                   ! REFERENCE TEMPERATURE
N,1,-10                   ! DEFINE NODES AND ELEMENTS
N,3,10   
FILL
N,4,-10,-20
N,6,10,-20   
FILL
E,1,4   
E,3,6
MAT,2
E,2,5
CP,1,UY,5,4,6             ! 耦合4, 5, 6三个节点的Y向位移
D,1,ALL,,,3               ! BOUNDARY CONDITIONS AND LOADING
F,5,FY,-4000
BFUNIF,TEMP,80            ! UNIFORM TEMPERATURE (TREF+10)
FINISH
/SOLU    
OUTPR,BASIC,1       
OUTPR,NLOAD,1             
NSUBST,1
SOLVE
FINISH
/POST1                                
STEEL_N  = NODE (,,,)                      ! 选择节点
COPPER_N = NODE (10,0,0)      
STEEL_E = ENEARN (STEEL_N)                 ! 选择单元
COPPER_E = ENEARN (COPPER_N)
ETABLE,STRS_ST,LS,1                        ! LINK1单元的ls,1表示SAXL: Axial stress in the element
ETABLE,STRS_CO,LS,1   
*GET,STRSS_ST,ELEM,STEEL_E,ETAB,STRS_ST  
*GET,STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO

*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'STRSS_ST','STRSS_CO'
LABEL(1,2) = ' (psi)  ',' (psi)  '
*VFILL,VALUE(1,1),DATA,19695,10152 
*VFILL,VALUE(1,2),DATA,STRSS_ST,STRSS_CO
*VFILL,VALUE(1,3),DATA,ABS(STRSS_ST/19695 ) ,ABS( STRSS_CO/10152 )
/COM
/OUT,vm3,vrt
/COM,------------------- VM3 RESULTS COMPARISON ---------------------
/COM,
/COM,                 |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,'   ',F10.0,'  ',F10.0,'   ',1F5.3)
/COM,----------------------------------------------------------------
/OUT

FINISH
*LIST,vm3,vrt

• 分析使用2D梁单元（beam3: 两节点，每个节点3个自由度ux, uy, rotz)
  
• 注意均布载荷为线载荷 (N/M)
  
• 单元输出位置是at 9 intermediate locations
  
• 根据beam3的单元说明，是用于etable的输出选项，确定table里各输出项的意义。
  

COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0
/VERIFY,VM2
JPGPRF,500,100,1             ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7    
MP,PRXY,,0.3
/TITLE, VM2, BEAM STRESSES AND DEFLECTIONS
C***        STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE 98, PROB. 4
ANTYPE,STATIC                
ET,1,BEAM3
KEYOPT,1,9,9                  ! OUTPUT AT 9 INTERMEDIATE LOCATIONS
R,1,50.65,7892,30
MP,EX,1,30E6
N,1                           ! DEFINE NODES AND ELEMENTS
N,5,480
FILL
E,1,2
EGEN,4,1,1
D,2,UX,,,,,UY                 ! BOUNDARY CONDITIONS AND LOADING
D,4,UY
SFBEAM,1,1,PRES,(10000/12)      
SFBEAM,4,1,PRES,(1E4/12)
FINISH
/SOLU    
OUTPR,BASIC,1
SOLVE
FINISH
/POST1
SET,1,1
PRNSOL,U,COMP
PRNSOL,ROT,COMP
PLDISP,1
MID_NODE = NODE (240,,, )               
*GET,DISP,NODE,MID_NODE,U,Y         
MID_ELM = ENEARN (MID_NODE)       
ETABLE,STRS,LS,3              
*GET,STRSS,ELEM,MID_ELM,ETAB,STRS

*DIM,LABEL,CHAR,2
*DIM,VALUE,,2,3
LABEL(1) = 'STRS_psi','DEF_in'
*VFILL,VALUE(1,1),DATA,-11400,0.182
*VFILL,VALUE(1,2),DATA,STRSS,DISP
*VFILL,VALUE(1,3),DATA,ABS(STRSS /11400 ) ,ABS( DISP /0.182 )
/OUT,vm2,vrt
/COM
/COM,-------------------VM2 RESULTS COMPARISON ---------------------
/COM,
/COM,         |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,'   ',F10.3,'  ',F10.3,'   ',1F5.3)
/COM,---------------------------------------------------------------
/OUT
FINISH
*LIST,vm2,vrt