VM140 (Stretching, Twisting and Bending of a Long, Solid Shaft) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM140
/PREP7
/TITLE, VM140, STRETCHING, TWISTING AND BENDING OF A LONG SOLID SHAFT
C*** STR. OF MATL., TIMOSHENKO, PART 1, 3RD ED., PAGE 296
ET,1,PLANE83,,,,2,,1
MP,EX,1,30E6
MP,NUXY,1,0 ! ZERO POISSONS RATIO
N,1
N,25,,24
FILL
NGEN,3,25,1,25,1,.25
E,51,53,3,1,52,28,2,26
EGEN,12,2,1
MODE,0,1 ! AXIAL + TORSION MODE
D,1,ALL,,,51,25 ! FIX SUPPORT
D,2,UZ,,,25 ! CENTERLINE CONSTRAINTS AGAINST TORSION
F,75,FY,100 ! APPLY AXIAL FORCE
F,75,FZ,400 ! APPLY TORSION
FINISH
/SOLU
OUTPR,BASIC,LAST ! PRINTOUT SOLUTION
SOLVE
MODE,1,1 ! BENDING MODE
DDELE,2,UZ,25 ! DELETE PREVIOUS UZ CONSTRAINTS
FDEL,75,ALL ! DELETE PREVIOUS FORCES
D,2,UY,,,25 ! CENTERLINE CONSTRAINTS AGAINST BENDING
F,75,FX,-50 ! APPLY VERTICAL FORCE
SOLVE
FINISH
/POST1
SET,1,1,,,,0.0 ! GET LOAD STEP 1 AT 0.0 DEGREES
LCWRITE,1 ! WRITE OUT AS LOAD CASE 1
SET,2,1,,,,0.0 ! GET LOAD STEP 2 AT 0.0 DEGREES
LCOPER,ADD,1 ! ADD LOAD CASE 1 TO LOAD CASE 2
ESEL,S,ELEM,,1 ! SELECT ELEMENT 1
NSLE
PRNSOL,S,COMP
PRNSOL,S,PRIN ! PRINT PRINCIPLE STRESSES
*GET,TOR,NODE,51,S,YZ
*GET,AXBND,NODE,51,S,Y
*GET,COMB,NODE,51,S,1
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'TORSION ','AXIAL+BEND ','COMBINED '
LABEL(1,2) = ' psi',' psi',' psi'
*VFILL,VALUE(1,1),DATA,1018.6,6238.9,6401
*VFILL,VALUE(1,2),DATA,TOR,AXBND,COMB
*VFILL,VALUE(1,3),DATA,ABS(TOR/1018.6),ABS(AXBND/6238.9),ABS(COMB/6401)
/COM
/OUT,vm140,vrt
/COM,------------------- VM140 RESULTS COMPARISON -------------
/COM,
/COM, LOAD STEP 3 | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F5.3)
/COM,-----------------------------------------------------------
*STATUS
/OUT
FINISH
*LIST,vm140,vrt
VM141 (Diametral Compression of a Disk) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM141
/SHOW
/PREP7
smrt,off
/DEVICE,VECTOR,ON
/TITLE, VM141, DIAMETRAL COMPRESSION OF A DISK
! THEORY OF ELASTICITY, TIMOSHENKO AND GOODIER, 2ND ED., PG 107
! PLANE STRESS ELEMENTS (PLANE82 AND PLANE2)
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,PLANE2,,,3,,,,1 ! THICKNESS INPUT, SUPPRESS SOLUTION PRINTOUT
ET,2,PLANE82,,,3,,,,1
ET,3,PLANE82,,,3,,2 ! NODAL STRESS PRINTOUT SELECTED
R,1,.2 ! THICKNESS
MP,EX,1,30E6
MP,NUXY,1,0.3
CSYS,1 ! CYLINDRICAL COORDINATES
K,1,1,90
K,2,.5,90
K,4,1
K,5,1,50
L,1,5
LESIZE,1,,,7,5
L,5,4
LESIZE,2,,,4,2
CSYS,0 ! CARTESIAN COORDINATES
K,3
L,3,4
LESIZE,3,,,5
L,2,3
LESIZE,4,,,4,2
L,2,5
LESIZE,5,,,5
L,1,2
LESIZE,6,,,7,5
A,1,2,5,5
A,2,3,4,5
TYPE,2
MSHK,1 ! MAPPED AREA MESH
MSHA,0,2D ! USING QUADS
AMESH,2 ! QUADRILATERAL MESHING (PLANE82)
EPLOT
TYPE,1
MSHK,0 ! FREE AREA MESH
MSHA,1,2D ! USING TRIS
AMESH,1 ! TRIANGLE MESHING (PLANE2)
EPLOT
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
ESLN
TYPE,3
EMODIF,ALL ! MODIFY ALL SELECTED ELEMENTS
SAVE,VM141,DB
NSEL,S,LOC,Y,0
DSYM,SYMM,Y ! SYMMETRY ALONG X AXIS
NSEL,S,LOC,X,0
DSYM,SYMM,X ! SYMMETRY ALONG Y AXIS
NSEL,ALL
ESEL,ALL
OUTPR,NSOL,NONE ! NODAL DISPL. & REACTION FORCES PRINTOUT CONTROL
OUTPR,ESOL,ALL ! ELEMENTAL PRINTOUT CONTROL
FK,1,FY,-1000 ! APPLY HALF OF FORCE (SYMMETRY)
FINISH
/SOLU
SOLVE
FINISH
/POST1
NSEL,S,LOC,X,0 ! SELECT ONLY THE NODES OF INTEREST
NSEL,A,LOC,X,0.1
NSEL,R,LOC,Y,0
PRNSOL,S,COMP ! PRINT COMPONENT NODAL STRESSES
NSEL,R,LOC,X,0
NSEL,A,LOC,X,0.2
NSEL,R,LOC,Y,0
*GET,SNOD,NODE,,NUM,MIN ! GET STARTING NODE FOR PATH (X=0)
*GET,FNOD,NODE,,NUM,MAX ! GET END NODE FOR PATH (X=0.2)
NSEL,ALL
ESEL,ALL
PATH,STRESS1,2,,48 ! DEFINE PATH WITH NAME = "STRESS1"
PPATH,1,SNOD ! DEFINE PATH POINTS BY NODE
PPATH,2,FNOD
PDEF,SY,S,Y ! INTERPOLATE SY STRESS ON PATH
PRANGE,24 ! PRINT EVERY 24TH POINT
PRPATH,SY ! PRINT SY STRESS ALONG THE PATH
*GET,S1,PATH,0,MIN,SY
NSEL,R,LOC,X,0
NSEL,A,LOC,X,0.1
NSEL,R,LOC,Y,0
*GET,FNOD,NODE,,NUM,MAX
PATH,STRESS2,2,,48 ! DEFINE PATH WITH NAME = "STRESS2"
PPATH,1,SNOD
PPATH,2,FNOD
PDEF,SY,S,Y ! INTERPOLATE SY STRESS ON PATH
PRANGE,24 ! PRINT EVERY 24TH POINT
PRPATH,SY ! PRINT SY STRESS ALONG THE PATH
*GET,S2,PATH,0,MAX,SY
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'P (psi)','P (psi)'
LABEL(1,2) = ' X=0 ',' X=.1'
*VFILL,VALUE(1,1),DATA,-9549,-9298
*VFILL,VALUE(1,2),DATA,S1,S2
*VFILL,VALUE(1,3),DATA,ABS(S1/9549),ABS(S2/9298)
SAVE,TABLE_1
FINISH
/CLEAR, NOSTART
/FILNAM,GEN
/PREP7
smrt,off
! S.E. GENERATION PASS (SHELL ELEMENTS ,SHELL93)
RESUME,VM141,DB
ANTYPE,SUBST ! SUBSTRUCTURE GENERATION PASS
SEOPT,GEN
! NOTE: SINCE PLANE2 AND SHELL93 HAVE DIFFERENT NODE ORDER, DELETE OLD MESH
ACLEAR,1,2 ! DELETE NODES AND ELEMENTS
ETDELE,1,3 ! DELETE PREVIOUS ELEMENT TYPES
ET,1,SHELL93,,,,,,,1 ! THICKNESS INPUT, SUPPRESS SOLUTION PRINTOUT
ET,2,SHELL93,,,,,,,1
ET,3,SHELL93,,,,,2 ! NODAL STRESS PRINTOUT
/OUTPUT,SCRATCH
NUMCMP,NODE ! COMPRESS NODE NUMBER TO ZERO
/OUTPUT
TYPE,2
MSHK,1 ! MAPPED AREA MESH
MSHA,0,2D ! USING QUADS
AMESH,2 ! MESH AREA 2 WITH QUADRILATERALS
EPLOT
TYPE,1
MSHK,0 ! FREE AREA MESH
MSHA,1,2D ! USING TRIS
AMESH,1 ! MESH AREA 1 WITH TRIANGLES
EPLOT
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
ESLN
TYPE,3
EMODIF,ALL
NSEL,S,LOC,Y,0
DSYM,SYMM,Y ! SYMMETRY ALONG X AXIS
NSEL,S,LOC,X,0
DSYM,SYMM,X ! SYMMETRY ALONG Y AXIS
NSEL,ALL
ESEL,ALL
D,ALL,UZ,,,,,ROTX,ROTY ! CONSTRAIN UNNEEDED DOF'S
NSEL,S,LOC,X
NSEL,R,LOC,Y,1
*GET,NDOF,NODE,,NUM,MAX ! GET NODE NUMBER FOR MASTER DOF
M,NDOF,UY ! SELECT MASTER DOF AT LOAD APPLICATION POINT
NSEL,ALL
FINISH
/SOLU
SOLVE
SAVE ! SAVE SUBSTRUCTURE DATA BASE FOR EXPANSION PASS
PARSAV,SCALAR,GEN,PARM
FINISH
/CLEAR, NOSTART
/FILNAM,USE
PARRES,,GEN,PARM
/PREP7
smrt,off
/TITLE, VM141, DIAMETRAL COMPRESSION OF A DISK (S.E. USE PASS)
ET,1,MATRIX50
SE,GEN
F,NDOF,FY,-1000
FINISH
/SOLU
SOLVE
FINISH
/CLEAR, NOSTART
/FILNAM,GEN
RESUME
/SOLU
EXPASS,ON,YES ! EXPANSION PASS WITH ELEMENT SOLUTION
SEEXP,GEN,USE
/TITLE, VM141, DIAMETRAL COMPRESSION OF A DISK (S.E. EXPANSION PASS)
OUTPR,NSOL,NONE ! DISPLACEMENT PRINTOUT CONTROL
OUTPR,ESOL,ALL
EXPSOL,1,1
SOLVE
FINISH
/POST1
NSEL,R,LOC,X,0
NSEL,A,LOC,X,0.2
NSEL,R,LOC,Y,0
*GET,SNOD,NODE,,NUM,MIN ! GET STARTING NODE FOR PATH (X=0)
*GET,FNOD,NODE,,NUM,MAX ! GET END NODE FOR PATH (X=0.2)
NSEL,ALL
ESEL,ALL
PATH,STRESS3,2,,48 ! DEFINE PATH WITH NAME = "STRESS3"
PPATH,1,SNOD
PPATH,2,FNOD
PDEF,SY,S,Y ! INTERPOLATE SY STRESS ON PATH
PRANGE,24 ! PRINT EVERY 24TH POINT
PRPATH,SY ! PRINT SY STRESS ALONG THE PATH
*GET,S1,PATH,0,MIN,SY
NSEL,R,LOC,X,0
NSEL,A,LOC,X,0.1
NSEL,R,LOC,Y,0
*GET,FNOD,NODE,,NUM,MAX
NSEL,ALL
ESEL,ALL
PATH,STRESS4,2,,48 ! DEFINE PATH WITH NAME = "STRESS4"
PPATH,1,SNOD
PPATH,2,FNOD
PDEF,SY,S,Y ! INTERPOLATE SY STRESS ON PATH
PRANGE,24 ! PRINT EVERY 24TH POINT
PRPATH,SY ! PRINT SY STRESS ALONG THE PATH
*GET,S2,PATH,0,LAST,SY
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'P (psi)','P (psi)'
LABEL(1,2) = ' X=0 ',' X=.1'
*VFILL,VALUE(1,1),DATA,-9549,-9298
*VFILL,VALUE(1,2),DATA,S1,S2
*VFILL,VALUE(1,3),DATA,ABS(S1/9549),ABS(S2/9298)
SAVE,TABLE_2
FINISH
FINISH
/CLEAR, NOSTART
/PREP7
smrt,off
/TITLE, VM141, DIAMETRAL COMPRESSION OF A DISK
ET,1,PLANE145,4,5,3 ! DEFINE ELEMENT WITH P-LEVELS & THICKNESS INPUT
R,1,.2 ! THICKNESS
MP,EX,1,30E6 ! DEFINE MATERIAL PROPERTIES
MP,NUXY,1,0.3
PCIRC,,1,0,90 ! CREATE MODEL GEOMETRY
ESIZE,,5
AMESH,ALL ! MESH ALL AREAS
EPLOT
FINISH
/SOLU
NSEL,S,LOC,Y,0 ! APPLY BOUNDARY CONDITIONS
DSYM,SYMM,Y ! SYMMETRY ALONG X AXIS
NSEL,S,LOC,X,0
DSYM,SYMM,X ! SYMMETRY ALONG Y AXIS
NSEL,ALL
ESEL,ALL
F,NODE(0,1,0),FY,-1000 ! APPLY HALF OF FORCE (SYMMETRY)
PCONV,1,S,Y,NODE(0,0,0) ! SET CONVERGENCE CRITERIA FOR SY AT CENTER
PCONV,1,S,Y,NODE(.10,0,0) ! SET CONVERGENCE CRITERIA FOR SY AT X=0.1,Y=0
SOLVE
PRCONV ! PRINT P-CONVERGENCE HISTORY
/POST1
*GET,S1,NODE,12,S,Y
*GET,S2,NODE,22,S,Y
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'P (psi)','P (psi)'
LABEL(1,2) = ' X=0 ',' X=.1'
*VFILL,VALUE(1,1),DATA,-9549,-9298
*VFILL,VALUE(1,2),DATA,S1,S2
*VFILL,VALUE(1,3),DATA,ABS(S1/9549),ABS(S2/9298)
SAVE,TABLE_3
RESUME,TABLE_1
FINISH
/DELETE,SCRATCH
/DELETE,GEN,PARM
/DELETE,GEN,db
/DELETE,GEN,emat
/DELETE,GEN,esav
/DELETE,GEN,redm
/DELETE,GEN,rst
/DELETE,GEN,seld
/DELETE,GEN,sub
/DELETE,GEN,tri
/DELETE,USE,dsub
/DELETE,USE,emat
/DELETE,USE,esav
/DELETE,USE,rst
/DELETE,USE,sord
/DELETE,USE,tri
/DELETE,VM141,DB
/COM
/OUT,vm141,vrt
/COM,------------------- VM141 RESULTS COMPARISON ------------
/COM,
/COM, TARGET | ANSYS | RATIO
/COM, PLANE82 AND PLANE2
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.0,' ',F11.0,' ',1F6.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM, SHELL93
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.0,' ',F11.0,' ',1F6.3)
/NOPR
RESUME,TABLE_3
/GOPR
/COM,
/COM, PLANE145
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.0,' ',F11.0,' ',1F6.3)
/COM,---------------------------------------------------------
/OUT
FINISH
! CLEANUP FILES
/DELETE,TABLE_1
/DELETE,TABLE_2
/DELETE,TABLE_3
*LIST,vm141,vrt
VM142 (Stress Concentration At a Hole in a Plate) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM142
/FILNAM,vm142 ! DEFINE JOBNAME FOR THE COARSE MODEL
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
smrt,off
/TITLE, VM142, STRESS CONCENTRATION AT A HOLE IN A PLATE
C*** ROARK 4TH EDITION, PAGE 384.
/NOPR
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,PLANE2
MP,EX,1,30E6
MP,NUXY,1,0.3
K,1,6 ! KEYPOINTS
K,2,6,6
K,3,,6
K,4,,.5
K,5,.5
K,6
L,1,2 ! LINE SEGMENTS
L,2,3
L,3,4
LESIZE,3,,,4,.25 ! DEFINE NO. OF DIVISONS AND SPACING FOR LINE 3
LARC,4,5,6,0.5
LESIZE,4,,,6 ! DEFINE NO. OF DIVISONS FOR LINE 4
L,5,1
LESIZE,5,,,4,4 ! DEFINE NO. OF DIVISONS AND SPACING FOR LINE 5
AL,1,2,3,4,5 ! AREA DEFINED BY 5 BOUNDING LINES
ESIZE,,4 ! 4 DIVISIONS PER LINE
AMESH,ALL
/AUTO,1 ! AUTO SCALE ON WINDOW 1
/PLOPTS,INFO,0 ! TURN WINDOW DOCUMENTATION OFF
/PLOPTS,WINS,0
/WINDOW,,LTOP ! WINDOW 1 AT LEFT TOP CORNER
EPLOT
LSEL,S,LINE,,3,5,2
DL,ALL,,SYMM
LSEL,S,LINE,,1
NSLL,,1
SF,ALL,PRES,-1000. ! APPLY TENSION ON PLATE
LSEL,ALL
NSEL,ALL
CSYS,1
FINISH
/SOLU
SOLVE
FINISH
SAVE ! SAVE FILE AS VM142.DB
/POST1
SET,1,1
NSORT,S,X,,,3 ! SORT BASED ON SX, RETAIN ONLY THE HIGHEST 3
PRNSOL,S,COMP
/WINDOW,1,OFF ! TURN WINDOW 1 OFF
/NOERASE ! OVERLAY DISPLAYS
/DSCALE,2,1
/WINDOW,2,RTOP ! WINDOW 2 AT RIGHT TOP CORNER
PLNSOL,S,X
*GET,CRSESTR,NODE,18,S,X
*STATUS
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'MX STR '
LABEL(1,2) = 'CRS MODEL'
*VFILL,VALUE(1,1),DATA,3018
*VFILL,VALUE(1,2),DATA,CRSESTR
*VFILL,VALUE(1,3),DATA,ABS(CRSESTR/3018)
SAVE,TABLE_1
FINISH
/CLEAR, NOSTART ! CLEAR THE DATABASE
/FILNAM,SUBMODEL ! DEFINE JOBNAME FOR THE SUBMODEL
/PREP7
smrt,off
/NOPR
/TITLE, VM142, STRESS CONCENTRATION AT A HOLE IN A PLATE
C*** BOUNDARY INTERPOLATION, MODIFIED SECTION
ANTYPE,STATIC
ET,1,PLANE42
MP,EX,1,30E6
MP,NUXY,1,0.3
CSYS,1
K,10,.5,45 ! DEFINE KEYPOINTS FOR SUBMODEL
K,11,.5,90
K,12,1.5,45
K,13,1.5,90
A,10,12,13,11
ESIZE,,8 ! 8 DIVISIONS
MSHK,1 ! MAPPED AREA MESH
MSHA,0,2D ! USING QUADS
AMESH,1
/WINDOW,2,OFF
/NOERASE
/PLOPTS,INFO,0
/PLOPTS,WINS,0
/WINDOW,1,LTOP ! WINDOW 1 AT LEFT TOP CORNER
/USER ! USER DEFINED SCALE
/DIST,1,3.3 ! USER DEFINED VIEW DISTANCE
/FOCUS,1,3,3 ! USER DEFINED FOCUS POINT
EPLOT
LSEL,S,LINE,,1,2
NSLL,,1 ! SELECT NODES OF CUT BOUNDARY
NWRITE ! WRITE GEOMETRY TO SUBMODEL.NODE
LSEL,ALL
NSEL,ALL
FINISH
SAVE ! SAVE SUBMODEL DATA IN FILE SUBMODEL.DB
/POST1
RESUME,vm142,db ! RESUME FROM FILE VM142.DB
FILE,vm142,rst ! DEFINE RESULTS FILE NAME
CBDOF,,,,,,,0,,0 ! ACTIVATE CUT BOUNDARY INTERPOLATION
FINISH
/PREP7
smrt,off
RESUME ! RESUME SUBMODEL FROM FILE SUBMODEL.DB
/NOPR
/INPUT,,cbdo,,:cb1 ! READ IN INTERPOLATED B.C.'S FROM SUBMODEL.CBDO
/GOPR
LSEL,S,LINE,,3 ! APPLY REMAINING BOUNDARY CONDITIONS
DL,ALL,,SYMM
FINISH
/SOLU
SOLVE
FINISH
/POST1
SET,1,1
NSORT,S,X,,,3
PRNSOL,S,COMP
/WINDOW,1,OFF
/AUTO,3
/WINDOW,3,BOT ! WINDOW 3 AT LOWER HALF
/NOERASE
/PLOPTS,MINM,1
/USER,3
/DIST,3,.2
/FOCUS,3,.2,.5
/CONTOUR,3,,AUTO ! AUTO CONTOUR SCALE FOR WINDOW 3
PLNSOL,SX
*GET,SUBSTR,NODE,18,S,X
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'MAX STRS'
LABEL(1,2) = ' SUBMOD'
*VFILL,VALUE(1,1),DATA,3018
*VFILL,VALUE(1,2),DATA,SUBSTR
*VFILL,VALUE(1,3),DATA,ABS(SUBSTR/3018 )
SAVE,TABLE_2
FINISH
/CLEAR,NOSTART
/FILNAM,SUBMODEL
/PREP7
smrt,off
/TITLE, VM142, STRESS CONCENTRATION AT A HOLE IN A PLATE
C*** ROARK 4TH EDITION, PAGE 384.
/NOPR
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,PLANE146
MP,EX,1,30E6
MP,NUXY,1,0.3
K,1,6 ! KEYPOINTS
K,2,6,6
K,3,,6
K,4,,.5
K,5,.5
K,6
L,1,2 ! LINE SEGMENTS
L,2,3
L,3,4
LESIZE,3,,,4,.25 ! DEFINE NO. OF DIVISONS AND SPACING FOR LINE 3
LARC,4,5,6,0.5
LESIZE,4,,,6 ! DEFINE NO. OF DIVISONS FOR LINE 4
L,5,1
LESIZE,5,,,4,4 ! DEFINE NO. OF DIVISONS AND SPACING FOR LINE 5
AL,1,2,3,4,5 ! AREA DEFINED BY 5 BOUNDING LINES
ESIZE,,4 ! 4 DIVISIONS PER LINE
AMESH,ALL
/AUTO,1 ! AUTO SCALE ON WINDOW 1
/PLOPTS,INFO,0 ! TURN WINDOW DOCUMENTATION OFF
/PLOPTS,WINS,0
/WINDOW,,LTOP ! WINDOW 1 AT LEFT TOP CORNER
EPLOT
LSEL,S,LINE,,3,5,2
DL,ALL,,SYMM
LSEL,S,LINE,,1
NSLL,,1
SF,ALL,PRES,-1000. ! APPLY TENSION ON PLATE
LSEL,ALL
NSEL,ALL
CSYS,1
FINISH
/SOLU
SOLVE
FINISH
/POST1
SET,1,1
!PRNSOL,S,COMP
/WINDOW,1,OFF ! TURN WINDOW 1 OFF
/NOERASE ! OVERLAY DISPLAYS
/DSCALE,2,1
/WINDOW,2,RTOP ! WINDOW 2 AT RIGHT TOP CORNER
PLNSOL,S,X
*GET,SUBSTR,NODE,18,S,X
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'MAX STRS'
LABEL(1,2) = 'CRS MOD '
*VFILL,VALUE(1,1),DATA,3018
*VFILL,VALUE(1,2),DATA,SUBSTR
*VFILL,VALUE(1,3),DATA,ABS(SUBSTR/3018 )
SAVE,TABLE_3
FINISH
/DEL,SUBMODEL,cbdo
/DEL,SUBMODEL,db
/DEL,SUBMODEL,emat
/DEL,SUBMODEL,esav
/DEL,SUBMODEL,node
/DEL,SUBMODEL,rst
/DEL,SUBMODEL,tri
RESUME,TABLE_1
/OUT,vm142,vrt
/COM,-------------------(VM142)RESULTS COMPARISON---------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,PLANE2
*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)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,PLANE42
*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)
/NOPR
RESUME,TABLE_3
/GOPR
/COM,
/COM,PLANE146
*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,-----------------------------------------------------------
/COM,
/OUT
FINISH
/DEL,TABLE_1
/DEL,TABLE_2
/DEL,TABLE_3
*LIST,vm142,vrt
VM143 (Fracture Mechanics Stress Intensity for a Crack in a Finite Width Plate) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM143
*CREATE,FRACT,MAC
! MACRO TO CREATE 3D SOLID95 CRACK TIP ELEMENTS FROM 3D SOLID45 ELEMENTS
! MAKE A COMPONENT CONTAINING THE CRACK TIP NODES (CRACKTIP)
! THE CRACK TIP IS BETWEEN NODES K AND O
! SET ELEMENT TYPE TO POINT TO SOLID95
! SET ARG1 TO N (THE TYPE OF THE ELEMENTS AROUND THE CRACK TIP)
!
/NOPR
NSEL,ALL
*GET,N,NODE,,NUM,MAX ! CURRENT MAXIMUM NODE NUMBER
CMSEL,S,CRACKTIP ! SELECT THE TIP NODES
ESLN ! ANY ELEMENTS ATTACHED
*GET,ELMAX,ELEM,,NUM,MAX ! CURRENT MAXIMUM ELEMENT NUMBER
*DO,IEL,1,ELMAX ! LOOP ON MAX ELEMENT
ELMI=IEL
*IF,ELMI,LE,0,EXIT ! NO MORE SELECTED
*GET,ELTYPE,ELEM,ELMI,ATTR,TYPE ! GET ELEMENT TYPE
*IF,ELTYPE,NE,ARG1,CYCLE ! CHECK FOR SELECTED ELEMENT
N3 = NELEM(ELMI,3) ! GET NODE 3 (K)
*IF,NSEL(N3),LE,0,CYCLE ! IT MUST BE SELECTED
N7 = NELEM(ELMI,7) ! GET NODE 7 (L)
*IF,NSEL(N7),LE,0,CYCLE ! IT MUST ALSO BE SELECTED
N1 = NELEM(ELMI,1) ! GET NODE 1 (I)
N2 = NELEM(ELMI,2) ! GET NODE 2 (J)
N5 = NELEM(ELMI,5) ! GET NODE 5 (M)
N6 = NELEM(ELMI,6) ! GET NODE 6 (N)
X3 = 0.75*NX(N3) ! WEIGHTED POSITION OF N3
Y3 = 0.75*NY(N3)
Z3 = 0.75*NZ(N3)
X = 0.25*NX(N2) + X3 ! QUARTER POINT LOCATION ( NODE (R) )
Y = 0.25*NY(N2) + Y3
Z = 0.25*NZ(N2) + Z3
N = N + 1 ! NEXT NODE
N10 = N
N,N10,X,Y,Z ! MIDSIDE NODE LOCATION
X = 0.25*NX(N1) + X3
Y = 0.25*NY(N1) + Y3
Z = 0.25*NZ(N1) + Z3
N = N + 1
N12= N
N,N12,X,Y,Z
X7 = 0.75*NX(N7)
Y7 = 0.75*NY(N7)
Z7 = 0.75*NZ(N7)
X = 0.25*NX(N6) + X7
Y = 0.25*NY(N6) + Y7
Z = 0.25*NZ(N6) + Z7
N = N + 1
N14 = N
N,N14,X,Y,Z
X = 0.25*NX(N5) + X7
Y = 0.25*NY(N5) + Y7
Z = 0.25*NZ(N5) + Z7
N = N + 1
N16 = N
N,N16,X,Y,Z
N4=N3
N8=N7
NSEL,ALL
TYPE,3
EN,ELMI,N1,N2,N3,N4,N5,N6,N7,N8 ! REDEFINE THE ELEMENT
EMORE,0,N10,0,N12,0,N14,0,N16
EMORE,
*ENDDO
CMSEL,U,CRACKTIP ! UNSELECT THE TIP NODES
NUMMRG,NODE ! MERGE MIDSIDE NODES
NSEL,ALL ! SELECT ALL ELEMENTS
ESEL,ALL ! SELECT ALL ELEMENTS
/GOPR
*END
/PREP7
SMRT,OFF
/TITLE, VM143, FRACTURE MECHANICS STRESS INTENSITY - CRACK IN A FINITE WIDTH PLATE
C*** BROWN AND SRAWLEY, ASTM SPECIAL TECHNICAL PUBLICATION NO. 410.
/COM, ****** CRACK IN 3-DIMENSIONS USING SOLID45 AND SOLID95
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,SOLID45
ET,2,SOLID45 ! ELEMENTS AROUND THE CRACK TIP
ET,3,SOLID95 ! CRACK TIP ELEMENTS CREATED USING MACRO FRACT
MP,EX,1,3E7
MP,NUXY,1,.3
CSYS,1 ! CYLINDRICAL COORDINATE SYSTEM
N,1
NGEN,9,20,1
N,11,.8
N,171,.8,180
FILL,11,171,7,31,20
CSYS,0 ! CARTESIAN COORDINATE SYSTEM
FILL,1,11,9,2,1,9,20,3
N,15,4
N,75,4,5
FILL,15,75,2,35,20
N,155,-1,5
FILL,75,155,3,95,20
N,172,-1
FILL,155,172,5,177,-1,,,.15
FILL,11,15,3,,,7,20,3
NGEN,2,200,1,177,,,,.25
E,2,22,1,1,202,222,201,201
EGEN,8,20,-1
E,2,3,23,22,202,203,223,222
EGEN,8,20,-1
EGEN,9,1,-8
EGEN,5,1,73,78
E,171,151,173,172,371,351,373,372
E,151,131,174,173,351,331,374,373
E,131,132,175,174,331,332,375,374
EGEN,3,1,-1
E,134,135,155,177,334,335,355,377
TYPE,2
EMODIF,1 ! MODIFY ELEMENTS 1 TO 8 FROM TYPE,1 TO TYPE,2
*REPEAT,8,1
NUMMRG,NODE ! MERGE COINCIDENT NODES
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
CM,CRACKTIP,NODE
/NERR,0 ! TEMPORARILY NO WARNINGS OR ERRORS PRINTOUT
! (IN ORDER TO AVOID WARNING MESSAGES DUE TO
! MIDSIDE NODES LOCATION)
FRACT,2 ! CONVERSION MACRO, TYPE 2 IS SOLID45
! ELEMENTS AROUND THE CRACK TIP
/NERR,DEFA ! TURN ON THE WARNINGS OR ERRORS PRINTOUT
/OUTPUT
OUTPR,,ALL
OUTPR,VENG,ALL ! STORE STRAIN ENERGY FOR J-INTEGRAL EVALUATION
NSEL,S,LOC,X,-1
DSYM,SYMM,X ! SYMMETRIC B.C.'S AT X = -1
NSEL,S,LOC,X,0,4
NSEL,R,LOC,Y,0
DSYM,SYMM,Y ! SYMMETRIC B.C.'S AT Y = 0 EXCEPT CRACK NODES
NSEL,ALL
D,ALL,UZ ! Z CONSTRAINTS FOR PLANE STRAIN PROBLEM
NSEL,S,LOC,Y,5
SF,ALL,PRES,-.5641895
NSEL,ALL
ESEL,ALL
FINISH
/OUTPUT,SCRATCH
/SOLU
SOLVE
FINISH
/OUTPUT
/POST1
ETABLE,SENE,SENE ! RETRIEVE STRAIN ENERGY PER ELEMENT
ETABLE,VOLU,VOLU ! RETRIEVE VOLUME PER ELEMENT
C*** IN POST1 DETERMINE KI (STRESS INTENSITY FACTOR) USING KCALC !**
PATH,KI1,3,,48 ! DEFINE PATH WITH NAME = "KI1"
PPATH,1,1 ! DEFINE PATH POINTS BY NODE
PPATH,2,406
PPATH,3,162
KCALC,,,1 ! COMPUTE KI FOR A HALF-MODEL WITH SYMM. B.C.
*GET,KI1,KCALC,,K,1 ! GET KI AS PARAMETER KI1
!********************************************************************************
!************************** J-INTEGRAL USER FILE *****************************
!*******************************************************************************
! ****NOTE:- IN GENERAL USAGE, THE USER FILE WOULD BE AVAILABLE IN THE
! LOCAL DIRECTORY RATHER THAN BEING CREATED IN THE INPUT
!*******************************************************************************
*CREATE,JIN1
STINFC ! DATA BLOCK NAME
SEXP,W,SENE,VOLU,1,-1 ! CALCULATE STRAIN ENERGY DENSITY
PATH,JINT,4,50,48 ! DEFINE PATH WITH NAME = "JINT"
PPATH,1,ARG1 ! DEFINE PATH POINTS BY NODE
PPATH,2,ARG2
PPATH,3,ARG3
PPATH,4,ARG4
PDEF,W,ETAB,W ! PUT STRAIN ENERGY DENSITY ON THE PATH
PCALC,INTG,J,W,YG ! INTEGRATE ENERGY W.R.T. GLOBAL Y
*GET,JA,PATH,,LAST,J ! GET FINAL VALUE OF INTEGRAL FOR 1ST TERM OF J
PDEF,CLEAR ! CLEAR OLD PATH VARIABLES
PVECT,NORM,NX,NY,NZ ! DEFINE THE PATH UNIT NORMAL VECTOR
PDEF,INTR,SX,SX ! PUT STRESS SX ON THE PATH
PDEF,INTR,SY,SY ! PUT STRESS SY ON THE PATH
PDEF,INTR,SXY,SXY ! PUT STRESS SXY ON THE PATH
PCALC,MULT,TX,SX,NX ! CALCULATE TRACTION TX
PCALC,MULT,C1,SXY,NY ! TX = SX*NX + SXY*NY
PCALC,ADD,TX,TX,C1
PCALC,MULT,TY,SXY,NX ! CALCULATE TRACTION TY
PCALC,MULT,C1,SY,NY ! TY = SXY*NX + SY*NY
PCALC,ADD,TY,TY,C1
*GET,DX,PATH,,LAST,S ! DEFINE PATH SHIFT AS 1% OF PATH LENGTH
DX=DX/100
PCALC,ADD,XG,XG,,,,-DX/2 ! SHIFT PATH FROM X TO X-DX/2 (GLOBAL X DIR.)
PDEF,INTR,UX1,UX ! DEFINE UX AT X-DX
PDEF,INTR,UY1,UY ! DEFINE UY AT X-DX
PCALC,ADD,XG,XG,,,,DX ! SHIFT PATH FROM X-DX/2 TO X+DX/2
PDEF,INTR,UX2,UX ! DEFINE UX AT X+DX
PDEF,INTR,UY2,UY ! DEFINE UY AT X+DX
PCALC,ADD,XG,XG,,,,-DX/2 ! SHIFT PATH BACK TO ORIGINAL POSITION
C=(1/DX)
PCALC,ADD,C1,UX2,UX1,C,-C ! CALCULATE DERIVATIVE DUX/DX
PCALC,ADD,C2,UY2,UY1,C,-C ! CALCULATE DERIVATIVE DUY/DX
PCALC,MULT,C1,TX,C1 ! DEFINE INTEGRAND
PCALC,MULT,C2,TY,C2 ! = TX*DUX/DX + TY*DUY/DX
PCALC,ADD,C1,C1,C2
PCALC,INTG,J,C1,S ! FORM SECOND INTEGRAL (W.R.T. PATH LENGTH S)
*GET,JB,PATH,,LAST,J ! GET FINAL VALUE OF INTEGRAL FOR 2ND TERM OF J
JINT=2*(JA-JB) ! ADD BOTH TERMS AND DOUBLE FOR HALF MODELS
PDEF,CLEAR ! CLEAR PATH VARIABLES
*END
C**************** IN POST1 DETERMINE KI FROM J-INTEGRAL !***********************
CON1=30E6/(1-(0.3*0.3)) ! J-TO-KI CONVERSION FACTOR
*ULIB,JIN1 ! ASSIGN LOCAL FILE JIN1 AS USER FILE
*USE,STINFC,4,45,125,164 ! USE DATA BLOCK STINFC AND GIVE PATH NODES
KI2=SQRT(CON1*JINT) ! CALCULATE KI FROM J
*STATUS,KI1 ! VIEW RESULTS
*STATUS,KI2
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'BY DISP ','BY J-'
LABEL(1,2) = 'EXTRP ','INT'
*VFILL,VALUE(1,1),DATA,1.0249,1.0249
*VFILL,VALUE(1,2),DATA,KI1,KI2
*VFILL,VALUE(1,3),DATA,ABS(KI1/1.0249),ABS(KI2/1.0249)
SAVE,TABLE_1
FINISH
/CLEAR, NOSTART ! CLEAR DATABASE FOR 2ND SOLUTION
/PREP7
SMRT,OFF
/TITLE, VM143, FRACTURE MECHANIC STRESS INTENSITY - CRACK IN A FINITE WIDTH PLATE
/COM, ****** CRACK IN 2-DIMENSIONS USING 2-D PLANE82 ******
ET,1,PLANE82,,,2 ! PLANE82 (PLANE STRAIN)
MP,EX,1,30E6
MP,NUXY,1,0.3
K,1 ! DEFINE KEYPOINTS AND LINE SEGMENTS
K,2,4
K,3,4,5
K,4,-1,5
K,5,-1
L,1,2
L,2,3
LESIZE,2,,,4
L,3,4
LESIZE,3,,,4
L,4,5,
LESIZE,4,,,6,.2
L,5,1
ESIZE,,5
KSCON,1,.15,1,8 ! DEFINE CRACK TIP ELEMENT SIZE
AL,1,2,3,4,5
DL,1,1,SYMM ! APPLY SOLID MODEL BOUNDARY CONDITIONS
DL,4,1,SYMM
SFL,3,PRES,-.5641895
AMESH,1
OUTPR,ALL
FINISH
/COM
/OUTPUT,SCRATCH
/SOLU
SOLVE
FINISH
/OUTPUT
/POST1
ETABLE,SENE,SENE ! RETRIEVE STRAIN ENERGY PER ELEMENT
ETABLE,VOLU,VOLU ! RETRIEVE VOLUME PER ELEMENT
C*** IN POST1 DETERMINE KI (STRESS INTENSITY FACTOR) USING KCALC !**
NSEL,S,LOC,Y,0 ! SELECT NODES FOR LPATH COMMAND
NSEL,R,LOC,X,0
*GET,NOD1,NODE,,NUM,MIN
NSEL,A,LOC,Y
NSEL,R,LOC,X,-.005,-.145
*GET,NOD2,NODE,,NUM,MIN
NSEL,A,LOC,Y
NSEL,R,LOC,X,-.145,-.155
*GET,NOD3,NODE,,NUM,MIN
NSEL,ALL
PATH,KI2,3,,48 ! DEFINE PATH WITH NAME = "KI2"
PPATH,1,NOD1 ! DEFINE PATH POINTS BY NODE
PPATH,2,NOD2
PPATH,3,NOD3
KCALC,,,1 ! COMPUTE KI FOR A HALF-MODEL WITH SYMM. B.C.
*GET,KI1,KCALC,,K,1 ! GET KI AS A PARAMETER KI1
C**************** IN POST1 DETERMINE KI FROM J-INTEGRAL !***********************
CSYS,1
NSEL,S,LOC,X,.5,.8 ! SELECT NODES FOR LPATH COMMAND IN STINFC
NSEL,R,LOC,Y,-1,1
*GET,NOD4,NODE,,NUM,MAX
NSEL,S,LOC,X,.5,.8
NSEL,R,LOC,Y,35,55
*GET,NOD5,NODE,,NUM,MAX
NSEL,S,LOC,X,.5,.8
NSEL,R,LOC,Y,120,145
*GET,NOD6,NODE,,NUM,MAX
NSEL,S,LOC,X,.5,.8
NSEL,R,LOC,Y,179,181
*GET,NOD7,NODE,,NUM,MIN
NSEL,ALL
CSYS,0
*USE,STINFC,NOD4,NOD5,NOD6,NOD7 ! USE DATA BLOCK STINFC AND GIVE PATH NODES
CON1=30E6/(1-(0.3*0.3)) ! J-TO-KI CONVERSION FACTOR
KI2=SQRT(CON1*JINT) ! CALCULATE KI FROM J
*STATUS,KI1 ! VIEW RESULTS
*STATUS,KI2
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'BY DISP ','BY J-'
LABEL(1,2) = 'EXTRP ','INT'
*VFILL,VALUE(1,1),DATA,1.0249,1.0249
*VFILL,VALUE(1,2),DATA,KI1,KI2
*VFILL,VALUE(1,3),DATA,ABS(KI1/1.0249),ABS(KI2/1.0249)
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm143,vrt
/COM,------------------- VM143 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,USING SOLID95 AND SOLID45 (3-D ANALYSIS)
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.4,' ',F10.4,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,USING PLANE82 (2-D ANALYSIS)
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.4,' ',F10.4,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm143,vrt
/DELETE,FRACT,MAC
VM144 (Bending of a Composite Beam) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM144
/PREP7
/TITLE, VM144, BENDING OF A COMPOSITE BEAM
C*** FORMULAS FOR STRESS AND STRAIN, ROARK, 5TH ED.
C*** USING LAYERED SHELL ELEMENT (SHELL99)
ANTYPE,STATIC
ET,1,SHELL99 ! 8 NODE LAYERED SHELL ELEMENT
R,1,2 ! 2 LAYERS PER ELEMENT
RMORE
RMORE,1,0,.2,2,0,.1 ! LAYER 1: 0.2 THK; LAYER 2: 0.1 THK
MP,EX,1,1.2E6 ! MATERIAL 1 PROPERTIES
MP,NUXY,1,0
MP,ALPX,1,18E-5
MP,ALPY,1,0.0
MP,EX,2,0.4E6 ! MATERIAL 2 PROPERTIES
MP,NUXY,2,0
MP,ALPX,2,6E-5
MP,ALPY,2,0
N,1
N,9,8
FILL
NGEN,3,10,1,9,,,.25
E,1,3,23,21,2,13,22,11
EGEN,4,2,-1 ! 4 ELEMENTS ALONG BEAM LENGTH
CP,1,ROTY,9,19,29 ! COUPLE FREE END NODES FOR ROTATION
D,1,ALL,,,21,10 ! FIXED END
F,19,MY,10 ! APPLY BENDING MOMENT AT FREE EDGE
BFUNIF,TEMP,100 ! ELEVATED TEMPERATURE LOAD
FINISH
/SOLU
SOLVE
FINISH
/POST1
NSEL,S,LOC,X,8 ! SELECT FREE EDGE
PRNSOL,U,Z ! PRINT DISPLACEMENTS
*GET,U1,NODE,9,U,Z
NSEL,S,LOC,Y ! SELECT NODES ALONG LENGTH
SHELL,TOP
PRNSOL,S,COMP ! PRINT TOP STRESSES
*GET,ST1,NODE,1,S,X
SHELL,BOT
PRNSOL,S,COMP ! PRINT BOTTOM STRESSES
*GET,SB1,NODE,1,S,X
ALLSEL
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DISP ','PRS TP ','PRS BTM '
LABEL(1,2) = 'in','psi','psi'
*VFILL,VALUE(1,1),DATA,.832,2258,1731
*VFILL,VALUE(1,2),DATA,U1,ST1,SB1
*VFILL,VALUE(1,3),DATA,ABS(U1/.832),ABS(ST1/2258),ABS(SB1/1731)
SAVE,TABLE_1
FINISH
/PREP7
C*** USING SHELL99 ELEMENTS WITH NODE OFFSET OPTION
ET,1,SHELL99 ! 8 NODE LAYERED SHELL ELEMENT
KEYOPT,1,11,2 ! WITH NODAL PLANE ON TOP FACE
R,1,1 ! BOTTOM LAYER
RMORE
RMORE,1,0,.2
ET,2,SHELL99 ! 8 NODE LAYERED SHELL ELEMENT
KEYOPT,2,11,1 ! WITH NODAL PLANE ON BOTTOM FACE
R,2,1 ! TOP LAYER
RMORE
RMORE,2,0,.1
EGEN,2,0,1,4,1,,1,1 ! GENERATE OVERLAPPING ELEMENTS
FINISH
/SOLU
SOLVE
FINISH
/POST1
NSEL,S,LOC,X,8 ! SELECT FREE EDGE
PRNSOL,U,Z ! PRINT DISPLACEMENTS
*GET,U2,NODE,9,U,Z
NSEL,S,LOC,Y ! SELECT NODES ALONG LENGTH
SHELL,TOP
ESEL,S,TYPE,,2 ! ELEMENTS REPRESENTING TOP LAYER
PRNSOL,S,COMP ! PRINT TOP STRESSES
*GET,ST2,NODE,1,S,X
SHELL,BOT
ESEL,S,TYPE,,1 ! ELEMENTS REPRESENTING BOTTOM LAYER
PRNSOL,S,COMP ! PRINT BOTTOM STRESSES
*GET,SB2,NODE,1,S,X
*VFILL,VALUE(1,1),DATA,.832,2258,1731
*VFILL,VALUE(1,2),DATA,U2,ST2,SB2
*VFILL,VALUE(1,3),DATA,ABS(U2/.832),ABS(ST2/2258),ABS(SB2/1731)
SAVE,TABLE_2
FINISH
/CLEAR,NOSTART
/PREP7
/TITLE, VM144, BENDING OF A COMPOSITE BEAM
C*** FORMULAS FOR STRESS AND STRAIN, ROARK, 5TH ED.
C*** USING LAYERED SOLID ELEMENTS (SOLID46)
ANTYPE,STATIC
ET,1,SOLID46,,,,,2,4 ! LAYERED SOLID ELEMENTS
R,1,2 ! 2 LAYERS
RMORE
RMORE,1,0,2,2,0,1 ! RELATIVE LAYER THICKNESSES
MP,EX,1,1.2E6 ! MATERIAL 1 PROPERTIES
MP,NUXY,1,0
MP,ALPX,1,1.8E-4
MP,ALPY,1,0.0
MP,ALPZ,1,0.0
MP,EX,2,0.4E6 ! MATERIAL 2 PROPERTIES
MP,NUXY,2,0
MP,ALPX,2,0.6E-4
MP,ALPY,2,0.0
MP,ALPZ,2,0.0
N,1
N,9,8
FILL
NGEN,2,10,1,9,1,,.5
NGEN,2,20,1,19,1,,,.3
E,1,2,12,11,21,22,32,31
EGEN,8,1,-1 ! 8 ELEMENTS ALONG LENGTH
D,1,ALL,,,31,10 ! FIXED END
F,9,FX,-(50/3),,19,10 ! APPLY NODAL FORCES TO GENERATE MOMENT
F,29,FX,(50/3),,39,10
BFUNIF,TEMP,100 ! ELEVATED TEMPERATURE LOAD
FINISH
/SOLU
SOLVE
FINISH
/POST1
NSEL,S,LOC,X,8
PRNSOL,U,Z ! PRINT FREE END DISPLACEMENTS
*GET,U3,NODE,9,U,Z
NSEL,S,LOC,Z,0.3
PRNSOL,S,COMP ! PRINT STRESSES ALONG TOP SURFACE
*GET,ST3,NODE,21,S,X
NSEL,S,LOC,Z
PRNSOL,S,COMP ! PRINT STRESSES ALONG BOTTOM SURFACE
*GET,SB3,NODE,1,S,X
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DISP ','PRS TP ','PRS BTM '
LABEL(1,2) = 'in','psi','psi'
*VFILL,VALUE(1,1),DATA,.832,2258,1731
*VFILL,VALUE(1,2),DATA,U3,ST3,SB3
*VFILL,VALUE(1,3),DATA,ABS(U3/.832),ABS(ST3/2258),ABS(SB3/1731)
SAVE, TABLE_3
FINISH
/CLEAR,NOSTART
/PREP7
/TITLE, VM144, BENDING OF A COMPOSITE BEAM
C*** FORMULAS FOR STRESS AND STRAIN, ROARK, 5TH ED.
C*** USING LAYERED SOLID ELEMENTS (SOLID191)
ANTYPE,STATIC
ET,1,SOLID191,,,,,2, ! LAYERED SOLID ELEMENTS
R,1,2 ! 2 LAYERS
RMORE
RMORE,1,0,2,2,0,1 ! RELATIVE LAYER THICKNESSES
MP,EX,1,1.2E6 ! MATERIAL 1 PROPERTIES
MP,NUXY,1,0
MP,ALPX,1,1.8E-4
MP,ALPY,1,0.0
MP,ALPZ,1,0.0
MP,EX,2,0.4E6 ! MATERIAL 2 PROPERTIES
MP,NUXY,2,0
MP,ALPX,2,0.6E-4
MP,ALPY,2,0.0
MP,ALPZ,2,0.0
N,1
N,9,8
FILL
NGEN,2,10,1,9,1,,.5
NGEN,2,20,1,19,1,,,.3
E,1,2,12,11,21,22,32,31
EGEN,8,1,-1 ! 8 ELEMENTS ALONG LENGTH
emid
nsel,s,loc,x
D,all,ALL, ! FIXED END
nsel,all
nlist,all
sfe,8,3,pres,,4000/3,4000/3,-4000/3,-4000/3 ! tapered pressure to apply moment on face
BFUNIF,TEMP,100 ! ELEVATED TEMPERATURE LOAD
FINISH
/SOLU
outpr,,1
SOLVE
FINISH
/POST1
NSEL,S,LOC,X,8
PRNSOL,U,Z ! PRINT FREE END DISPLACEMENTS
*GET,U3,NODE,9,U,Z
NSEL,S,LOC,Z,0.3
PRNSOL,S,COMP ! PRINT STRESSES ALONG TOP SURFACE
*GET,ST3,NODE,21,S,X
NSEL,S,LOC,Z
PRNSOL,S,COMP ! PRINT STRESSES ALONG BOTTOM SURFACE
*GET,SB3,NODE,1,S,X
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DISP ','PRS TP ','PRS BTM '
LABEL(1,2) = 'in','psi','psi'
*VFILL,VALUE(1,1),DATA,.832,2258,1731
*VFILL,VALUE(1,2),DATA,U3,ST3,SB3
*VFILL,VALUE(1,3),DATA,ABS(U3/.832),ABS(ST3/2258),ABS(SB3/1731)
SAVE, TABLE_4
VM176 (Frequency Response of Electrical Input Admittance for a Piezoelectric Transducer) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM176
/PREP7
SMRT,OFF
/TITLE, VM176, FREQUENCY RESPONSE OF ELECTRICAL INPUT ADMITTANCE FOR A
/COM PIEZOELECTRIC TRANSDUCER
/COM KAGAWA AND YAMABUCHI, FINITE ELEMENT SIMULATION OF A COMPOSITE
/COM PIEZOELECTRIC ULTRASONIC TRANSDUCER, IEEE TRANS. SONICS AND
/COM ULTRASONICS, VOL. SU-26, NO.2, MARCH 1979
ET,1,SOLID5,0 ! 3-D COUPLED-FIELD SOLID
MP,DENS,3,7730 ! NEPEC DENSITY
MP,EX,2,7.03E10 ! ALUMINUM MODULUS OF ELASTICITY
MP,NUXY,2,.345 ! ALUMINUM POISSON RATIO
MP,DENS,2,2690 ! ALUMINUM DENSITY
MP,EX,4,10E9 ! ADHESIVE MODULUS OF ELASTICITY
MP,DENS,4,1700 ! ADHESIVE DENSITY
MP,NUXY,4,.38 ! ADHESIVE POISSON RATIO
/COM ** DEFINE GEOMETRIC PARAMETERS **
R=27.5E-3 ! DISK RADIUS
HA=15.275E-3 ! ALUMINUM 1/2 HEIGHT LOCATION
HN=5E-3 ! NEPEC 1/2 HEIGHT
HB=5.275E-3 ! ADHESIVE MATERIAL HEIGHT
RDIV=5 ! NO. ELEMENTS ALONG RADIUS
HADV=3 ! NO. ELEMENTS ALONG ALUMINUM HEIGHT
HNDV=2 ! NO. ELEMENTS ALONG NEPEC HEIGHT
HBDV=1 ! NO. ELEMENTS ALONG ADHESIVE HEIGHT
ZRO=1E-5 ! DEFINE ZERO FOR KEYPOINT LOCATION
CSYS,1
K,1,ZRO,-5 ! DEFINE KEYPOINTS FOR MESH, WEDGE ELEMENTS
K,2,R,-5 ! ARE NOT ALLOWED BY MESH MODULE SO KEYPOINTS
K,3,R,-5,HN ! ARE DEFINED NEAR ZERO AND LATER MERGED
K,4,ZRO,-5,HN
K,5,R,-5,HB
K,6,ZRO,-5,HB
K,7,R,-5,HA
K,8,ZRO,-5,HA
KGEN,2,1,8,1,,10
L,2,3 ! DEFINE LINE SEGMENTS
LESIZE,1,,,HNDV
L,3,5
LESIZE,2,,,HBDV
L,5,7
LESIZE,3,,,HADV
L,2,10
LESIZE,4,,,1
V,11,3,4,12,10,2,1,9 ! CREATE NEPEC VOLUME
VATT,3 ! ASSIGN MATERIAL ATTRIBUTES
V,13,5,6,14,11,3,4,12 ! CREATE ADHESIVE VOLUME
VSEL,S,VOLU,,2
VATT,4 ! ASSIGN MATERIAL ATTRIBUTES
V,15,7,8,16,13,5,6,14 ! CREATE ALUMINUM VOLUME
VSEL,S,VOLU,,3
VATT,2 ! ASSIGN MATERIAL ATTRIBUTES
VSEL,ALL
MSHK,1 ! MAPPED VOLUME MESH
MSHA,0,3D ! USING HEX
ESIZE,,RDIV
SHPP,OFF ! TURN OFF SHAPE CHECKING TO ALLOW
! FOR WEDGE SHAPE MESH OF SOLIDS
VMESH,ALL ! MESH ALL VOLUMES
NUMMRG,NODE ! MERGE NODES TO CREATE WEDGE ELMENTS
! AT AXIS
NSEL,S,LOC,Y,-5
DSYM,SYMM,Y,1 ! SYMMETRY B.C. AT THETA=-5 DEG.
NSEL,S,LOC,Y,5
DSYM,SYMM,Y,1 ! SYMMETRY B.C. AT THETA=5 DEG.
NSEL,S,LOC,X,0,.001
DSYM,SYMM,X,1 ! SYMMETRY B.C. AT X=0
NSEL,S,LOC,Z,0
DSYM,SYMM,Z,1 ! SYMMETRY B.C. AT Z=0
NSEL,S,LOC,Z,HN ! SELECT NODES AT TOP ELECTRODE
CP,1,VOLT,ALL ! COUPLE VOLT DOF ON ELECTRODE
*GET,N1,NODE,,NUM,MIN ! GET NODE ON ELECTRODE
D,N1,VOLT,-0.5 ! APPLY -0.5 VOLT TO TOP ELECTRODE
NSEL,S,LOC,Z,0 ! SELECT NODES AT SYMMETRY PLANE
D,ALL,VOLT,0.0 ! SET VOLT TO ZERO AT SYMMETRY PLANE
NSEL,ALL
FINISH
/SOLU
EQSLV,SPARSE ! USING SPARSE MATRIX SOLVER
ANTYPE,HARMIC ! PERFORM HARMONIC ANALYSIS
OUTRES,ALL,ALL ! STORE EVERY SUBSTEP
HARFRQ,5000,35000 ! SOLVE FOR FREQ=20KHZ AND 35KHZ
NSUBST,2
KBC,1 ! STEP BOUNDARY CONDITIONS
EQSLV,ICCG ! ICCG SOLVER
SOLVE
HARFRQ,39000,45000 ! SOLVE FOR FREQ=42KHZ AND 45KHZ
SOLVE
HARFRQ,46000,54000 ! SOLVE FOR FREQ=50KHZ AND 54KHZ
SOLVE
FINISH
/POST26
RFORCE,2,N1,AMPS ! STORE CHARGES ON ELECTRODE
PI2=(3.14159*2.)
PROD,3,2,1,,MHOS,,,PI2 ! CALCULATE ADMITTANCE (10 DEG. SLICE)
PROD,4,3,,,MMHO,,,36000 ! CALCULATE TOTAL ADMITTANCE (MMHOS)
PRVAR,4 ! PRINT ELECTRICAL ADMITTANCE VS. FREQUENCY
*GET,F1,VARI,4,RTIME,20000
*GET,F2,VARI,4,RTIME,35000
*GET,F3,VARI,4,RTIME,42000
*GET,F4,VARI,4,RTIME,45000
*GET,F5,VARI,4,RTIME,50000
*GET,F6,VARI,4,RTIME,54000
*STATUS
*DIM,LABEL,CHAR,6,2
*DIM,VALUE,,6,3
LABEL(1,1) = 'Y MMHOS ','Y MMHOS','Y MMHOS ','Y MMHOS ','Y MMHOS ','Y MMHOS '
LABEL(1,2) = '@20 kHz','@35 kHz','@42 kHz','@45 kHz','@50 kHz','@54 kHz'
*VFILL,VALUE(1,1),DATA,.41,.9,2,0,.39,.65
*VFILL,VALUE(1,2),DATA,F1,F2,F3,F4,F5,F6
*VFILL,VALUE(1,3),DATA,ABS(F1/.41),ABS(F2/0.9),ABS(F3/2),0,ABS(F5/.39),ABS(F6/0.65)
/COM
/OUT,vm176,vrt
/COM,------------------- VM176 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.2,' ',F10.2,' ',1F5.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm176,vrt
VM177 (Natural Frequency of a Submerged Ring) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM177
/PREP7
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM177, NATURAL FREQUENCY OF SUBMERGED RING
/COM, REF. "FINITE ELEMENT SOLUTION OF FLUID STRUCTURE
/COM, INTERACTION PROBLEMS" SCHROEDER & MARCUS
/COM, SHOCK & VIBRATION SYMPOSIUM, SAN DIEGO, 1975
/COM, USING FLUID30 AND SHELL63 AND FULL HARMONIC ANALYSIS (ANTYPE=3)
ET,1,FLUID30 ! FLUID ELEMENTS INTERFACING WITH STRUCTURE
ET,2,SHELL63 ! SHELL ELEMENTS TO MODEL STEEL RING
ET,3,FLUID30,,1 ! NON-INTERFACING FLUID ELEMENTS
R,1,1
R,2,0.25
MP,DENS,1,9.6333E-5
MP,SONC,1,57480.0 ! SPEED OF SOUND IN WATER
MP,EX,2,30.E6
MP,DENS,2,7.4167E-4
MP,NUXY,2,0.3
CSYS,1
N,1,10.0
N,7,30.0
FILL
NGEN,9,10,1,10,1,0,(90/8)
NGEN,2,100,1,99,1,0,0,1 ! DEFINE UPPER PLANE OF NODES
E,1,101,111,11,2,102,112,12
EGEN,6,1,-1
EGEN,8,10,-6
TYPE,2
REAL,2
MAT,2
E,1,101,111,11
EGEN,8,10,-1
NSEL,S,LOC,X,10.0
ESLN
ESEL,INVE
TYPE,3
REAL,1
MAT,1
EMODIF,ALL ! CHANGE ELEMENT TYPE TO TYPE 3
ESEL,ALL
NSEL,S,LOC,X,10
SF,ALL,FSI ! COUPLE STRUCTURAL MOTION & FLUID PRESS.
NSEL,ALL
D,7,PRES,0.0,,87,10 ! SET PRESSURE AT OUTER RADIUS TO ZERO
D,107,PRES,0.0,,187,10
NSEL,S,LOC,X,10.0
NSEL,R,LOC,Y,90.0
DSYM,SYMM,X
NSEL,S,LOC,X,10.0
NSEL,R,LOC,Y,0.0
DSYM,SYMM,Y
NSEL,ALL
D,1,UZ,0.0,,7,6
D,81,UZ,0.0,,87,6
F,1,FX,1.0 ! EXCITE THE EVEN MODES OF VIBRATION
F,81,FY,0,1.0
FINISH
/SOLU
ANTYPE,HARMIC ! FULL HARMONIC ANALYSIS
HARFRQ,36.59,36.64 ! SELECT FREQUENCY RANGE
NSUBST,5
KBC,1
OUTRES,,1
EQSLV,ICCG ! ICCG SOLVER
/OUT,SOLVE,LOG ! DIVERT PLATFORM DEPENDENT SOLVER OUTPUT
SOLVE
/OUT
FINISH
/POST26
NSOL,2,1,U,X
NSOL,3,41,U,X
NSOL,4,81,U,Y
PRVAR,2,3,4 ! PRINT DISPLACEMENTS OF RING VS. FREQUENCY
FINISH
/POST1
SET,1,2
PLDISP,1
*GET,F1,ACTIVE,0,SET,FREQ
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'F, '
LABEL(1,2) = ' HZ'
*VFILL,VALUE(1,1),DATA,35.62
*VFILL,VALUE(1,2),DATA,F1
*VFILL,VALUE(1,3),DATA,ABS(F1/35.62 )
SAVE,TABLE_1
SET,1,2,,1
PLDISP,1
FINISH
/CLEAR, NOSTART ! CLEAR DATABASE BEFORE STARTING PART 2
/PREP7
/COM, USING FLUID29 AND BEAM3 AND UNSYMMETRIC MATRIX MODAL ANALYSIS (ANTYPE=2)
ET,1,FLUID29 ! FLUID ELEMENTS INTERFACING WITH STRUCTURE
ET,2,BEAM3 ! BEAM ELEMENTS TO MODEL STEEL RING
ET,3,FLUID29,,1 ! NON-INTERFACING FLUID ELEMENTS
MP,EX,1,1
MP,DENS,1,9.6333E-5
MP,SONC,1,57480.0
MP,EX,2,30.E6
MP,DENS,2,7.4167E-4
MP,NUXY,2,0.3
R,1,1
R,2,.25,(.25**3)/12,.25
CSYS,1
N,1,10.0
N,7,30.0
FILL
NGEN,9,10,1,10,1,0,(90/8)
E,1,2,12,11 ! DEFINE FLUID ELEMENTS INTERFACE WITH STRUCTURE
EGEN,6,1,-1
EGEN,8,10,-6
TYPE,2
MAT,2
REAL,2
E,1,11 ! DEFINE BEAM ELEMENTS TO MODEL STEEL RING
EGEN,8,10,-1
ESEL,,TYPE,,2
NSLE
NSEL,R,LOC,Y,90
DSYM,SYMM,X ! APPLY SYMMETRY BOUNDARY CONDITION ON RING
NSLE
NSEL,R,LOC,Y,0
DSYM,SYMM,Y
ESEL,ALL
NSEL,S,LOC,X,10
ESLN
ESEL,INVE
TYPE,3
MAT,1
REAL,1
EMODIF,ALL ! DEFINE UNCOUPLED FLUID ELEMENTS
NSEL,ALL
ESEL,S,TYPE,,1 ! DEFINE FLUID-STRUCTURE INTERFACE
NSEL,R,LOC,X,10
SF,ALL,FSI ! COUPLE STRUCTURAL MOTION & FLUID PRESS.
ESEL,ALL
NSEL,S,LOC,X,30
D,ALL,PRES,0.0 ! SET PRESSURE AT OUTER RADIUS TO ZERO
NSEL,ALL
FINISH
/SOLU
ANTYPE,MODAL ! MODAL ANALYSIS
MODOPT,UNSYM,4,-100 ! SELECT UNSYMMETRIC MATRIX MODE EXTRACTION
MXPAND,,36.50,36.60 ! EXPAND THE MODES OVER A SELECTED FREQUENCY
SOLVE ! RANGE
FINISH
/POST1
SET,1,1
PLDISP,1
*GET,F2,ACTIVE,0,SET,FREQ
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'F, '
LABEL(1,2) = ' Hz'
*VFILL,VALUE(1,1),DATA,35.62
*VFILL,VALUE(1,2),DATA,F2
*VFILL,VALUE(1,3),DATA,ABS(F2/35.62 )
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm177,vrt
/COM,------------------- VM177 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS USING FLUID30 & SHELL63
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS USING FLUID29 & BEAM3
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/DELETE,TABLE_1
/DELETE,TABLE_2
FINISH
*LIST,vm177,vrt
VM178 (Plane Poiseuille Flow) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM178
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM178, PLANE POISEULEE FLOW
! FOUNDATIONS OF FLUID MECHANICS, S.W. YUAN, SECTION 8.36
/PREP7
smrt,off
ET,1,FLUID141 ! FLUID ELEMENT
FLDATA,NOMI,DENS,1 ! DENSITY
FLDATA,NOMI,VISC,1 ! VISCOSITY
K,1,0 ! GEOMETRY
K,2,10.
K,3,10.,2.
K,4, 0.,2.
L,1,2
LESIZE,1,,,8
L,2,3
LESIZE,2,,,6
L,3,4
LESIZE,3,,,8
L,4,1
LESIZE,4,,,6
A,1,2,3,4
AMESH,1
LSEL,S,LINE,,4 ! SELECT LINE SEGMENT AT INLET
NSLL,S,1 ! SELECT NODES ON INLET LINE SEGMENT
D,ALL,VY ! NO VELOCITY IN Y DIRECTION
D,ALL,PRES,0.1 ! INLET PRESSURE = 0.1 PSI
LSEL,S,LINE,,2 ! SELECT LINE SEGMENT AT OUTLET
NSLL,S,1
D,ALL,VY ! NO VELOCITY IN Y DIRECTION
D,ALL,PRES,0.0 ! OUTLET PRESSURE = 0.0 PSI
LSEL,S,LINE,,1,3,2
NSLL,S,1
D,ALL,VX,,,,,VY ! NO SLIP BC AT CHANNEL EDGES
LSEL,S,LINE,,ALL
NSEL,S,NODE,,ALL
NSEL,S,LOC,X,10. ! GET NODE NUMBERS FOR PATH ENDPOINTS
NSEL,R,LOC,Y, 0.
*GET,NOD1,NDMX
NSEL,S,LOC,X,10.
NSEL,R,LOC,Y, 2.
*GET,NOD2,NDMX
NSEL,ALL
FINISH
/SOLU
/OUT,SCRATCH ! REDIRECT CONVERGENCE OUTPUT
FLDATA,ITER,EXEC,50 ! # OF ITERATIONS
SOLVE
/OUT
FINISH
/POST1
SET,1,1
/EDGE,1,1
PLVECT,V ! DISPLAY VELOCITY VECTORS
/EDGE,1,0
PATH,VXX,2,,48 ! DEFINE PATH WITH NAME = "VXX"
PPATH,1,NOD1 ! DEFINE PATH POINTS BY NODE
PPATH,2,NOD2
PDEF,VX,V,X
PRPATH,VX
/AXLAB,Y,VELOCITY
PLPATH,VX ! DISPLAY VELOCITY PROFILE
LSEL,S,LINE,,ALL
NSEL,S,NODE,,ALL
NSEL,S,LOC,X,10
NSEL,R,LOC,Y,1
*GET,NOD2,NDMX
*GET,V1,NODE,NOD2,V,X
LSEL,S,LINE,,ALL
NSEL,S,NODE,,ALL
NSEL,S,LOC,X,10
NSEL,R,LOC,Y,(1+(1/3))
*GET,NOD3,NDMX
*GET,V2,NODE,NOD3,V,X
LSEL,S,LINE,,ALL
NSEL,S,NODE,,ALL
NSEL,S,LOC,X,10
NSEL,R,LOC,Y,(1+(2/3))
*GET,NOD4,NDMX
*GET,V3,NODE,NOD4,V,X
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = '@ Y=','@ Y=','@ Y='
LABEL(1,2) = '1 in ','1.33 in','1.67 in'
*VFILL,VALUE(1,1),DATA,.005,.0044,.0028
*VFILL,VALUE(1,2),DATA,V1,V2,V3
*VFILL,VALUE(1,3),DATA,ABS(V1/.005),ABS(V2/.0044),ABS(V3/.0028)
FINISH
/DELETE,SCRATCH
/DELETE,vm178,pfl
/DELETE,vm178,rsw
/OUT,vm178,vrt
/COM,------------------- VM178 RESULTS COMPARISON --------------
/COM,
/COM,
/COM,VX (in/sec) | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.4,' ',F10.4,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm178,vrt
VM179 (Dynamic Double Rotation of a Jointed Beam) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM179
/PREP7
/TITLE, VM179, DYNAMIC DOUBLE ROTATION OF A JOINTED BEAM
C*** REFERENCE -- ANY BASIC MECHANICS TEXT
ANTYPE,TRANS ! FULL TRANSIENT DYNAMIC ANALYSIS
NLGEOM,ON ! LARGE DEFLECTION
SSTIF,OFF ! TURN OFF STRESS STIFFENING
NROPT,FULL,,OFF ! USE FULL NEWTON-RAPHSON WITHOUT ADAPTIVE DESCENT
ET,1,COMBIN7,,1 ! JOINT LOCK-UP WHEN STOP ENGAGED
ET,2,BEAM4
R,1,1E6,1E6,1E6,1E7,0,.05 ! TF = .05
RMORE,0,.5,,,,0.08727 ! 5 DEG. UPPER LIMIT
R,2,1,1,1,1,1
MP,EX,1,70E9 ! ALUMINUM MATERIAL PROPERTIES
MP,DENS,1,1E-6
MP,NUXY,1,.35
BETAD,3.1831E-6 ! 1% DAMPING AT 1000 HZ
N,1
N,2,1
N,3,1
N,4,2
N,5,1,-1 ! NODE 5 IS USED FOR ORIENTING JOINT
TYPE,1
REAL,1
E,2,3,5,5,1 ! JOINT ELEMENT
TYPE,2
REAL,2
E,1,2 ! BEAMS
E,3,4
D,1,UX,0,,,,UY,UZ,ROTX,ROTY ! ALLOW ONLY Z-ROTATION AT PINNED END
KBC,1
FINISH
/OUTPUT,SCRATCH
/SOLU
PRED,ON
NSUBST,5
F,1,MZ,.7854 ! BEGIN ROTZ WITH NON-ZERO ACCELERATION
TIME,.05
SOLVE
NSUBST,60
TIME,1 ! DELTA THETA-Z = 45 DEG. ( AT 1 SEC )
SOLVE
NSUBST,5
F,1,MZ,-.7854 ! SLOW ROTZ W/ "INSTANTANEOUS" REVERSAL OF APPLIED MOMENT
TIME,1.05
SOLVE
NSUBST,60
TIME,2 ! DELTA THETA-Z = 90 DEG. ( AT 2 SEC )
SOLVE
FINISH
/OUTPUT
/POST1
SET,2
PRNSOL,DOF
*GET,DX1,NODE,4,U,X
*GET,DY1,NODE,4,U,Y
*GET,RZ1,NODE,4,ROT,Z
ESEL,S,TYPE,,1
ETABLE,ROT_STAT,NMISC,1
ETABLE,ROTATE,NMISC,9
ESEL,ALL
PRETAB,ROTATE,ROT_STAT
SET,4
PRNSOL,DOF
*GET,DX2,NODE,4,U,X
*GET,DY2,NODE,4,U,Y
*GET,RZ2,NODE,4,ROT,Z
ESEL,S,TYPE,,1
ETABLE,REFL
ESEL,ALL
PRETAB,ROTATE,ROT_STAT
FINISH
PARSAV
/COPY,,rdb,,rest,rdb ! SAVE THE FILES NEEDED FOR RESTART (NOT NEEDED
/COPY,,ldhi,,rest,ldhi ! FOR STRAIGHT-THRU RUN)
/COPY,,r001,,rest,r001
/COPY,,rst,,rest,rst ! NEEDED FOR CONTINUITY OF THE RESULTS FILE
/CLEAR,NOSTART ! CLEAR THE DATABASE (TO SIMULATE RESTART)
/COM,
/COM, ----- RESTART ANALYSIS -----
/COPY,rest,rdb,,file,rdb ! COPY THE FILES NEEDED FOR RESTART (NOT NEEDED
/COPY,rest,ldhi,,file,ldhi ! FOR STRAIGHT-THRU RUN)
/COPY,rest,r001,,file,r001
/COPY,rest,rst,,file,rst
/OUTPUT,SCRATCH
/SOLU
ANTYPE,,REST ! USE RESTART ANALYSIS TO DEFINE MORE LOADSTEPS
NSUBST,5
F,1,MZ,0 ! REMOVE M1, ALLOW ROTZ MOTION TO STABILIZE
TIME,2.05
SOLVE
NSUBST,60
F,3,MX,.5 ! DELTA THETA-X = 5 DEG., THEN STOP ( AT 3 SEC )
TIME,3
SOLVE
FINISH
/OUTPUT
PARRES ! RESTORE PARAMETERS FROM INITIAL RUN
/POST1
SET,LAST
PRNSOL,DOF
*GET,DX3,NODE,4,U,X
*GET,DY3,NODE,4,U,Y
*GET,DZ3,NODE,4,U,Z
*GET,RX3,NODE,4,ROT,X
*GET,RZ3,NODE,4,ROT,Z
ESEL,S,TYPE,,1
ETABLE,ROT_STAT,NMISC,1
ETABLE,ROTATE,NMISC,9
ESEL,ALL
PRETAB,ROTATE,ROT_STAT
*STATUS
*DIM,LABEL,CHAR,11,2
*DIM,VALUE,,11,3
LABEL(1,1) = 'DX T1 ','DY T1 ','AZ T1 ','DX T2 ','DY T2 ','AZ T2 ','DX T3 ','DY T3 ','DZ T3 '
LABEL(10,1) = 'AX T3 ','AZ T3 '
LABEL(1,2) = 'in','in','rad','in','in','rad','in','in','in'
LABEL(10,2) = 'RAD','RAD'
*VFILL,VALUE(1,1),DATA,-.5858,1.4142,.7854,-2,2,1.5708,-2,1.9962,.08716
*VFILL,VALUE(10,1),DATA,.08727,1.5708
*VFILL,VALUE(1,2),DATA,DX1,DY1,RZ1,DX2,DY2,RZ2,DX3,DY3,DZ3
*VFILL,VALUE(10,2),DATA,RX3,RZ3
V1 = ABS(DX1/.5858)
V2 = (DY1/1.4142)
V3 = (RZ1/.7854)
V4 = ABS(DX2/2)
V5 = (DY2/2)
V6 = (RZ2/1.5708)
V7 = ABS(DX3/2)
V8 = (DY3/1.9962)
*VFILL,VALUE(1,3),DATA,V1,V2,V3,V4,V5,V6,V7,V8
*VFILL,VALUE(9,3),DATA,(DZ3/.08716),(RX3/.08727),(RZ3/1.5708)
/COM
/OUT,vm179,vrt
/COM,------------------- VM179 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.5,' ',F10.5,' ',1F5.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm179,vrt
VM180 (Bending of a Curved Beam) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM180
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
smrt,off
/TITLE, VM180, BENDING OF A CURVED BEAM
C*** THEORY OF ELASTICITY, TIMOSHENKO & GOODIER, 3RD ED., P. 78
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,PLANE2
ET,2,BEAM3
R,2,1,1,1 ! ARBITRARY BEAM ELEMENT REAL CONSTANTS
MP,EX,1,30E6
MP,NUXY,1,0
CSYS,1
K,1,3.5
K,2,3.5,90
KGEN,2,1,2,1,1.0
L,2,4 ! CREATE STIFFENING BEAM ELEMENTS
LESIZE,1,,,4
TYPE,2
REAL,2
LMESH,1
L,1,2 ! CREATE CURVED BEAM
LESIZE,2,,,25
L,3,4
LESIZE,3,,,25
ESIZE,,4
A,1,2,4,3
AATT,1,1,1
AMESH,1
NROTAT,ALL
DK,1,ALL,,,1 ! DEFINE KEYPOINT CONSTRAINTS AND FORCES
DK,3,UY,,,1 ! IN GLOBAL CYLINDRICAL COORDINATE SYSTEM
FK,2,FY,100
FK,4,FY,-100
FINISH
/SOLU
SOLVE
FINISH
/POST1
RSYS,1
NSEL,R,LOC,X,3.5 ! SELECT INNER RADIUS NODES
NSEL,U,LOC,Y,90
PRNSOL,S,COMP
NSORT,S,Y
*GET,SI,SORT,,MAX
NSEL,ALL
NSEL,R,LOC,X,4.5 ! SELECT OUTER RADIUS NODES
NSEL,U,LOC,Y,90
PRNSOL,S,COMP
NSORT,S,Y
*GET,SO,SORT,,MIN
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'INR STR ','OTR STR '
LABEL(1,2) = 'psi ','psi '
*VFILL,VALUE(1,1),DATA,655,-555
*VFILL,VALUE(1,2),DATA,SI,SO
*VFILL,VALUE(1,3),DATA,ABS(SI/655),ABS(SO/555)
/COM
/OUT,vm180,vrt
/COM,------------------- VM180 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.1,' ',F10.1,' ',1F5.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm180,vrt
VM181 (Natural Frequency of a Flat Circular Plate with a Clamped Edge) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM181
/PREP7
MP,PRXY,,0.3
SMRT,OFF
/TITLE, VM181, NATURAL FREQUENCY OF A FLAT CIRCULAR PLATE WITH A CLAMPED EDGE
C*** FORMULAS FOR NATURAL FREQUENCY AND MODE SHAPE, BLEVINS, PAGE 241
ET,1,PLANE2,,,1 ! AXISYMMETRIC ELEMENTS
MP,EX,1,3E7
MP,DENS,1,0.00073
MP,PRXY,1,0.3
K,1
K,2,17
KGEN,2,1,2,1,,.5
L,1,2
L,3,4
LESIZE,ALL,,,10 ! TEN DIVISIONS ALONG LENGTH
ESIZE,,1
A,1,2,4,3
AMESH,1
NSEL,R,LOC,X,0
D,ALL,UX
NSEL,ALL
DK,2,ALL,,,1
DK,4,UX,,,1
FINISH
/SOLU
ANTYPE,MODAL ! MODAL ANALYSIS
MODOPT,SUBSP,9
EQSLVE,PCG ! POWERSOLVER
SOLVE
*GET,F1,MODE,1,FREQ
*GET,F2,MODE,2,FREQ
*GET,F3,MODE,3,FREQ
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'F (0,0) ','F (0,1) ','F (0,2) '
LABEL(1,2) = 'Hz','Hz','Hz'
*VFILL,VALUE(1,1),DATA,172.64,671.79,1505.7
*VFILL,VALUE(1,2),DATA,F1,F2,F3
*VFILL,VALUE(1,3),DATA,ABS(F1/172.64),ABS(F2/671.79),ABS(F3/1505.7)
/COM
/OUT,vm181,vrt
/COM,------------------- VM181 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.2,' ',F15.2,' ',1F10.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm181,vrt
VM182 (Transient Response of a Spring-mass System) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM182
/PREP7
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM182; TRANSIENT RESPONSE OF A SPRING-MASS SYSTEM
C*** R. K. VIERCK, "VIBRATION ANALYSIS", 2ND EDITION, SEC.5-8
ANTYPE,MODAL ! MODE-FREQUENCY ANALYSIS
MODOPT,REDUC,2,,,2 ! PRINT TWO REDUCED MODE SHAPES
ET,1,COMBIN40,,,2 ! UY DOF
R,1,6,,2 ! K1=6 N/M M1=2 KG
R,2,16,,2 ! K2=16 N/M M2=2 KG
N,1
N,2,0,1
N,3,0,2
REAL,1
E,1,2
REAL,2
E,2,3
M,1,UY,2 ! UY MASTERS AT NODES 1 & 2
D,3,ALL
OUTPR,,ALL
FINISH
/SOLU
SOLVE
FINISH
/SOLU
ANTYPE,TRANS ! TRANSIENT DYNAMIC ANALYSIS
TRNOPT,MSUP,2 ! MODE SUPERPOSITION, BOTH MODES
DELTIM,0.01 ! INTEGRATION TIME STEP = .01
OUTPR,,NONE
OUTRES,,1
KBC,1 ! STEP BOUNDARY CONDITIONS
F,1,FY,0 ! FORCE = 0 AT TIME = 0 (INIT. CONDITIONS)
SOLVE
TIME,1.8
F,1,FY,50 ! FORCE = 50N FROM TIME = 0 TO 1.8 SEC
SOLVE
TIME,2.4
F,1,FY,0 ! FORCE = 0 FROM TIME = 1.8 TO 2.4 SEC
SOLVE
FINISH
/POST26
FILE,,rdsp ! REDUCED DISPLACEMENTS FILE
NSOL,2,1,U,Y,UY1
NSOL,3,2,U,Y,UY2
/GRID,1
/AXLAB,Y,DISP
PLVAR,2,3 ! DISPLAY DISPLACEMENT RESPONSE VS. TIME
PRVAR,2,3
*GET,Y1,VARI,2,RTIME,1.3
*GET,Y2,VARI,3,RTIME,1.3
*GET,Y3,VARI,2,RTIME,2.4
*GET,Y4,VARI,3,RTIME,2.4
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'DISP 1 ','DISP 2 '
LABEL(1,2) = 'm','m'
*VFILL,VALUE(1,1),DATA,14.48,3.99
*VFILL,VALUE(1,2),DATA,Y1,Y2
*VFILL,VALUE(1,3),DATA,ABS(Y1/14.48),ABS(Y2/3.99)
SAVE,TABLE_1
*VFILL,VALUE(1,1),DATA,18.32,6.14
*VFILL,VALUE(1,2),DATA,Y3,Y4
*VFILL,VALUE(1,3),DATA,ABS(Y3/18.32),ABS(Y4/6.14),
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm182,vrt
/COM,------------------- VM182 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS AT T=1.3 S
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS AT T=2.4 S
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/DELETE,TABLE_1
/DELETE,TABLE_2
FINISH
*LIST,vm182,vrt
VM183 (Harmonic Response of a Spring-mass System) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM183
/PREP7
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM183, HARMONIC RESPONSE OF A SPRING-MASS SYSTEM
C*** R.K. VIERCK, "VIBRATION ANALYSIS", 2ND EDITION, SECTION 4-2
ANTYPE,MODAL ! MODE - FREQUENCY ANALYSIS
MODOPT,REDUC,2,,,2 ! PRINT TWO REDUCED MODE SHAPES
ET,1,COMBIN40,,,2 ! UY DOF
R,1,6,,2 ! K1=6 N/M M1=2 KG
R,2,16,,2 ! K2=16 N/M M2=2 KG
N,1
N,2,0,1
N,3,0,2
REAL,1
E,1,2
REAL,2
E,2,3
M,1,UY,2 ! UY MASTERS AT NODES 1 & 2
OUTPR,,ALL
D,3,ALL
FINISH
/SOLU
SOLVE
FINISH
/SOLU
ANTYPE,HARMIC ! HARMONIC ANALYSIS
HROPT,MSUP,2 ! MODE SUPERPOSITION USING TWO MODES
HARFRQ,0.1,1.0 ! RANGE OF FREQUENCIES FROM 0.1 TO 1.0 HZ
F,1,FY,50
KBC,1 ! STEP BOUNDARY CONDITIONS
NSUBST,50
OUTPR,,NONE
OUTRES,,1
SOLVE
FINISH
/POST26
FILE,,rfrq ! REDUCED FREQUENCIES FILE
NSOL,2,1,U,Y,UY1
NSOL,3,2,U,Y,UY2
/GRID,1
/AXLAB,Y,DISP
PLVAR,2,3 ! DISPLAY DISPLACEMENT RESPONSE VS. FREQUENCY
PRVAR,2,3
*GET,Y1,VARI,2,RTIME,.226
*GET,Y2,VARI,3,RTIME,.226
*GET,Y3,VARI,2,RTIME,.910
*GET,Y4,VARI,3,RTIME,.910
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'DISP 1 ','DISP 2 '
LABEL(1,2) = 'm','m'
*VFILL,VALUE(1,1),DATA,-1371.7,-458.08
*VFILL,VALUE(1,2),DATA,Y1,Y2
*VFILL,VALUE(1,3),DATA,ABS(Y1/1371.7),ABS(Y2/458.08)
SAVE,TABLE_1
*VFILL,VALUE(1,1),DATA,-.8539,.1181
*VFILL,VALUE(1,2),DATA,Y3,Y4
*VFILL,VALUE(1,3),DATA,ABS(Y3/.8539),ABS(Y4/.1181),
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm183,vrt
/COM,------------------- VM183 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS AT .226 Hz
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS AT .910 Hz
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.4,' ',F10.4,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm183,vrt
/DELETE,TABLE_1
/DELETE,TABLE_2
FINISH
VM184 (Straight Cantilever Beam) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM184
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
smrt,off
/TITLE, VM184, STRAIGHT CANTILEVER BEAM
C*** ANY BASIC MECHANICS OF MATERIALS TEXT
C*** USING SOLID5 HEXAHEDRONS
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,SOLID5,2 ! MULTI-FIELD SOLID5
MP,EX,1,1E7
MP,NUXY,1,0.3
K,1
K,2,6
KGEN,2,1,2,1,,.2
KGEN,2,1,4,1,,,.1
L,1,2
LESIZE,ALL,,,10
ESIZE,,1
V,1,2,4,3,5,6,8,7
VMESH,1
NSEL,S,LOC,X,0
D,ALL,ALL ! CONSTRAIN LEFT END
NSEL,ALL
FK,2,FX,0.25 ! APPLY AXIAL FORCES
*REPEAT,4,2
NOORDER
FINISH
/SOLU
SOLVE
FKDELE,ALL,FX ! DELETE AXIAL FORCES
FK,2,FY,0.25 ! APPLY IN-PLANE LOADS
*REPEAT,4,2
SOLVE
FKDELE,ALL,FY ! DELETE IN-PLANE LOADS
FK,2,FZ,0.25 ! APPLY OUT-OF-PLANE LOADS
*REPEAT,4,2
SOLVE
FINISH
/POST1 ! PRINT END DISPLACEMENTS AS RATIO OF ANSYS:TARGET
CSYS,0
*CREATE,MAC ! DEFINE MACRO TO CALCULATE ANSYS:TARGET RATIOS
SET,ARG1,1
LCDEF,ARG1,ARG1
NSEL,S,LOC,X,6 ! SELECT NODE AT END OF BEAM
PRNSOL,U,COMP ! PRINT DISPLACEMENTS
LCFACT,ARG1,ARG2 ! APPLY SCALE FACTOR "ARG2" TO LOAD CASE 1
LCASE,ARG1
PRNSOL,U,COMP ! PRINT DISPLACEMENTS
*END
/COM *** USE MACRO TO PROCESS ALL 3 LOADCASES ***
*USE,MAC,1,(1/3E-5)
SET,1,1
LCSEL,S,0,0
N1=NODE(6,0,0)
*GET,UX1,NODE,N1,U,X
*USE,MAC,2,(1/.108)
SET,2,1
LCSEL,S,0,0
*GET,UY1,NODE,N1,U,Y
*USE,MAC,3,(1/.432)
SET,3,1
LCSEL,S,0,0
*GET,UZ1,NODE,N1,U,Z
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z '
LABEL(1,2) = 'in','in','in'
*VFILL,VALUE(1,1),DATA,3E-5,.108,.432
*VFILL,VALUE(1,2),DATA,UX1,UY1,UZ1
*VFILL,VALUE(1,3),DATA,ABS(UX1/(3E-5)),ABS(UY1/.108),ABS(UZ1/.432)
SAVE,TABLE_1
FINISH
/CLEAR, NOSTART
/PREP7
smrt,off
MOPT,VMESH,MAIN
MOPT,AMESH,ALTE
/TITLE, VM184, STRAIGHT CANTILEVER BEAM
C*** USING SOLID92 TETRAHEDRONS
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,SOLID92 ! STRUCTURAL SOLID92
MP,EX,1,1E7
MP,NUXY,1,0.3
K,1
K,2,,.2
KGEN,2,1,2,1,,,.1
KGEN,2,1,4,1,(2/3)
ESIZE,,1
V,1,2,4,3,5,6,8,7
VMESH,1
VGEN,9,1,1,1,(2/3),,,4 ! GENERATE 9 VOLUMES TO COMPLETE BEAM
NSEL,S,LOC,X,0
D,ALL,ALL ! CONSTRAIN LEFT END
SAVE
/COM *** GET NODE NUMBERS FOR LOAD APPLICATION ***
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.1
NSEL,R,LOC,Z,.05
*GET,MIDD,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.2
NSEL,R,LOC,Z,.05
*GET,TOPP,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.0
NSEL,R,LOC,Z,.05
*GET,BOTT,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.1
NSEL,R,LOC,Z,.10
*GET,LFT,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.1
NSEL,R,LOC,Z,.00
*GET,RGHT,NDMX
NSEL,ALL
/COM *** APPLY LOADS TO PARAMETRIC NODE NUMBERS ***
F,RGHT,FX,(1/6)
F,LFT,FX,(1/6)
F,TOPP,FX,(1/6)
F,BOTT,FX,(1/6)
F,MIDD,FX,(1/3)
FINISH
/SOLU
SOLVE
FDELE,ALL ! REMOVE ALL FORCES
F,RGHT,FY,(1/6)
F,LFT,FY,(1/6)
F,TOPP,FY,(1/6)
F,BOTT,FY,(1/6)
F,MIDD,FY,(1/3)
SOLVE
FDELE,ALL ! REMOVE ALL FORCES
F,RGHT,FZ,(1/6)
F,LFT,FZ,(1/6)
F,TOPP,FZ,(1/6)
F,BOTT,FZ,(1/6)
F,MIDD,FZ,(1/3)
SOLVE
FINISH
/POST1 ! PRINT END DISPLACEMENTS AS RATIO OF ANSYS:TARGET
/COM *** USE MACRO TO PROCESS ALL 3 LOADCASES ***
*USE,MAC,1,(1/3E-5)
SET,1,1
LCSEL,S,0,0
N1=NODE(6,0,0)
*GET,UX2,NODE,N1,U,X
*USE,MAC,2,(1/.108)
SET,2,1
LCSEL,S,0,0
*GET,UY2,NODE,N1,U,Y
*USE,MAC,3,(1/.432)
SET,3,1
LCSEL,S,0,0
*GET,UZ2,NODE,N1,U,Z
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z '
LABEL(1,2) = 'in','in','in'
*VFILL,VALUE(1,1),DATA,3E-5,.108,.432
*VFILL,VALUE(1,2),DATA,UX2,UY2,UZ2
*VFILL,VALUE(1,3),DATA,ABS(UX2/(3E-5)),ABS(UY2/.108),ABS(UZ2/.432)
SAVE,TABLE_2
FINISH
/CLEAR, NOSTART
/PREP7
smrt,off
MOPT,VMESH,MAIN
MOPT,AMESH,ALTE
/COM *** REPEAT USING SOLID98 TETRAHEDRONS ***
RESUME ! RESTORE PREP7 DATABASE
ET,1,SOLID98,2 ! MULTI-FIELD SOLID98
/COM *** GET NODE NUMBERS FOR LOAD APPLICATION ***
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.1
NSEL,R,LOC,Z,.05
*GET,MIDD,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.2
NSEL,R,LOC,Z,.05
*GET,TOPP,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.0
NSEL,R,LOC,Z,.05
*GET,BOTT,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.1
NSEL,R,LOC,Z,.10
*GET,LFT,NDMX
NSEL,S,LOC,X,6
NSEL,R,LOC,Y,.1
NSEL,R,LOC,Z,.0
*GET,RGHT,NDMX
NSEL,ALL
/COM *** APPLY LOADS TO PARAMETRIC NODE NUMBERS ***
F,RGHT,FX,(1/6)
F,LFT,FX,(1/6)
F,TOPP,FX,(1/6)
F,BOTT,FX,(1/6)
F,MIDD,FX,(1/3)
FINISH
/SOLU
SOLVE
FDELE,ALL ! REMOVE ALL FORCES
F,RGHT,FY,(1/6)
F,LFT,FY,(1/6)
F,TOPP,FY,(1/6)
F,BOTT,FY,(1/6)
F,MIDD,FY,(1/3)
SOLVE
FDELE,ALL ! REMOVE ALL FORCES
F,RGHT,FZ,(1/6)
F,LFT,FZ,(1/6)
F,TOPP,FZ,(1/6)
F,BOTT,FZ,(1/6)
F,MIDD,FZ,(1/3)
SOLVE
FINISH
/POST1 ! PRINT END DISPLACEMENTS AS RATIO OF ANSYS:THEORY
/COM *** USE MACRO TO PROCESS ALL 3 LOADCASES ***
*USE,MAC,1,(1/3E-5)
SET,1,1
LCSEL,S,0,0
N1=NODE(6,0,0)
*GET,UX3,NODE,N1,U,X
*USE,MAC,2,(1/.108)
SET,2,1
LCSEL,S,0,0
*GET,UY3,NODE,N1,U,Y
*USE,MAC,3,(1/.432)
SET,3,1
LCSEL,S,0,0
*GET,UZ3,NODE,N1,U,Z
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z '
LABEL(1,2) = 'in','in','in'
*VFILL,VALUE(1,1),DATA,3E-5,.108,.432
*VFILL,VALUE(1,2),DATA,UX3,UY3,UZ3
*VFILL,VALUE(1,3),DATA,ABS(UX3/(3E-5)),ABS(UY3/.108),ABS(UZ3/.432)
SAVE,TABLE_3
FINISH
/CLEAR,NOSTART
/PREP7
smrt,off
MOPT,VMESH,ALTE
MOPT,AMESH,ALTE
ANTYPE,STATIC ! STATIC ANALYSIS
ET,1,SOLID147
MP,EX,1,1E7
MP,NUXY,1,0.3
K,1
K,2,6
KGEN,2,1,2,1,,.2
KGEN,2,1,4,1,,,.1
L,1,2
LESIZE,ALL,,,10
ESIZE,,1
V,1,2,4,3,5,6,8,7
VMESH,1
NSEL,S,LOC,X,0
D,ALL,ALL ! CONSTRAIN LEFT END
NSEL,ALL
FK,2,FX,0.25 ! APPLY AXIAL FORCES
*REPEAT,4,2
NOORDER
FINISH
/SOLU
SOLVE
FKDELE,ALL,FX ! DELETE AXIAL FORCES
FK,2,FY,0.25 ! APPLY IN-PLANE LOADS
*REPEAT,4,2
SOLVE
FKDELE,ALL,FY ! DELETE IN-PLANE LOADS
FK,2,FZ,0.25 ! APPLY OUT-OF-PLANE LOADS
*REPEAT,4,2
SOLVE
FINISH
/POST1 ! PRINT END DISPLACEMENTS AS RATIO OF ANSYS:TARGET
CSYS,0
*CREATE,MAC ! DEFINE MACRO TO CALCULATE ANSYS:TARGET RATIOS
SET,ARG1,1
LCDEF,ARG1,ARG1
NSEL,S,LOC,X,6 ! SELECT NODE AT END OF BEAM
PRNSOL,U,COMP ! PRINT DISPLACEMENTS
LCFACT,ARG1,ARG2 ! APPLY SCALE FACTOR "ARG2" TO LOAD CASE 1
LCASE,ARG1
PRNSOL,U,COMP ! PRINT DISPLACEMENTS
*END
/COM *** USE MACRO TO PROCESS ALL 3 LOADCASES ***
*USE,MAC,1,(1/3E-5)
SET,1,1
LCSEL,S,0,0
N1=NODE(6,0,0)
*GET,UX1,NODE,N1,U,X
*USE,MAC,2,(1/.108)
SET,2,1
LCSEL,S,0,0
*GET,UY1,NODE,N1,U,Y
*USE,MAC,3,(1/.432)
SET,3,1
LCSEL,S,0,0
*GET,UZ1,NODE,N1,U,Z
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z '
LABEL(1,2) = 'in','in','in'
*VFILL,VALUE(1,1),DATA,3E-5,.108,.432
*VFILL,VALUE(1,2),DATA,UX1,UY1,UZ1
*VFILL,VALUE(1,3),DATA,ABS(UX1/(3E-5)),ABS(UY1/.108),ABS(UZ1/.432)
SAVE,TABLE_4
RESUME,TABLE_1
/COM *** CLIPPED AND CAPPED DISPLAY OF STRESS CONTOURS ***
NSEL,ALL
/VIEW,1,2,1,1
EPLOT ! ELEMENT PLOT
/TYPE,1,CAP ! DISPLAY TYPE CAP
/DIST,1,.2
/FOCUS,1,.3,.15,.09 ! SET FOCUS FOR SECTION LOCATION
PLNSOL,S,X ! STRESS CONTOUR PLOT
/COM
/OUT,vm184,vrt
/COM,------------------- VM184 RESULTS COMPARISON ------------
/COM,
/COM, TARGET | ANSYS | RATIO
/COM, SOLID5
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.6,' ',F11.6,' ',1F6.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM, SOLID92
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.6,' ',F11.6,' ',1F6.3)
/NOPR
RESUME,TABLE_3
/GOPR
/COM,
/COM, SOLID98
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.6,' ',F11.6,' ',1F6.3)
/NOPR
RESUME,TABLE_4
/GOPR
/COM
/COM, SOLID147
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F11.6,' ',F11.6,' ',1F6.3)
/COM,---------------------------------------------------------
/OUT
FINISH
/DELETE,MAC
/DELETE,TABLE_1
/DELETE,TABLE_2
/DELETE,TABLE_3
/DELETE,TABLE_4
*LIST,vm184,vrt
VM185 (AC Analysis of a Slot Embedded Conductor) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM185
/PREP7
SMRT,OFF
/TITLE, VM185, AC ANALYSIS OF A SLOT EMBEDDED CONDUCTOR
C*** KONRAD, A., "INTEGRODIFFERENTIAL FINITE ELEMENT FORMULATION
C*** OF TWO-DIMENSIONAL STEADY-STATE SKIN EFFECT PROBLEMS",
C*** IEEE TRANS. MAGNETICS, VOL. MAG-18, NO. 1, JAN. 1982
C*** PP. 284-292.
C***
ET,1,PLANE13 ! PLANE13, AZ DOF, (FOR AIR)
ET,2,PLANE13,6 ! PLANE13, AZ VOLT DOF, (FOR CONDUCTOR)
EMUNIT,MKS ! DEFINE SYSTEM UNITS
MP,MURX,1,1 ! RELATIVE PERMEABILITY
MP,MURX,2,1
RES=1.724E-8 ! DEFINE RESISTIVITY OF CONDUCTOR
MP,RSVZ,2,RES ! CONDUCTOR RESISTIVITY
A=6.45E-3 ! DEFINE GEOMETRY IN TERMS OF PARAMETERS
B=8.55E-3
C=8.45E-3
D=18.85E-3
E=8.95E-3
F=(D-E)/2
PTXY,0,0,D,0,D,C,D-F,C ! CREATE POLYGON AREA OF CONDUCTOR
PTXY,D-F,B+C,F,B+C,F,C,0,C
POLY
RECTNG,F,F+E,B+C,B+C+A ! CREATE AIR AREA
AGLUE,1,2 ! GLUE AREAS TOGETHER
ASEL,S,AREA,,3
AATT,1,,1 ! SET ATTRIBUTES FOR AIR
ASEL,S,AREA,,1
AATT,2,,2 ! SET ATTRIBUTES FOR CONDUCTOR
ASEL,ALL
ESIZE,D/15 ! SET ELEMENT EDGE LENGTH
MSHK,0 ! FREE AREA MESH
MSHA,1,2D ! USING TRIS
AMESH,ALL ! MESH AREAS
ESEL,S,MAT,,2 ! SELECT ALL NODES IN CONDUCTOR
NSLE,S
CP,1,VOLT,ALL ! COUPLE ALL NODES IN VOLT
I=1.0 ! DEFINE TOTAL CURRENT
ASUM ! CALCULATE AREA ATTRIBUTES
*GET,AREA,AREA,1,AREA ! GET AREA OF CONDUCTOR
*GET,N1,NODE,,NUM,MIN ! SELECT A NODE IN THE CONDUCTOR
F,N1,AMPS,I ! APPLY 1 AMP TOTAL CURRENT
ESEL,ALL
NSEL,S,LOC,Y,.02345 ! SELECT NODES AT TOP PLANE
D,ALL,AZ,0 ! SET FLUX PARALLEL B.C.
NSEL,ALL
FINISH
/SOLU
ANTYPE,HARMIC ! HARMONIC ANALYSIS
HARFRQ,45 ! SET OPERATING FREQUENCY
SOLVE
FINISH
/POST1
SET,1,1 ! RETRIEVE REAL SOLUTION
ETABLE,JT,NMISC,7 ! STORE TOTAL CURRENT DENSITY
ETABLE,JS,SMISC,1 ! STORE SOURCE CURRENT DENSITY
ETABLE,JE,NMISC,6 ! STORE EDDY CURRENT DENSITY
/PNUM,MAT,1
/EDGE,1,1
/NUM,1
/SHOW,JPEG
/GFILE,500
JPEG,QUAL,100
/TRIAD,OFF
/PLOPTS,LOGO,0
/PLOPTS,INFO,2
/PLOPTS,WP,0
/RGB,INDEX,100,100,100,0
/RGB,INDEX,80,80,80,13
/RGB,INDEX,60,60,60,14
/RGB,INDEX,0,0,0,15
PLNSOL,A,Z ! DISPLAY FLUX LINES
/NUM,0
ESEL,MAT,2 ! SELECT COPPER ONLY
PLETAB,JT,1 ! DISPLAY TOTAL CURRENT DENSITY
PLETAB,JE,1 ! DISPLAY EDDY CURRENT DENSITY
*GET,JSR,ELEM,1,ETAB,JS ! GET REAL COMPONENT OF JS
ACRE=JSR*RES/I ! CALCULATE AC RESISTANCE/LENGTH
SET,1,1,,1 ! READ IN IMAGINARY DATA
ETABLE,REFL ! REFILL ELEMENT TABLE WITH IMAG. DATA
*GET,JSI,ELEM,1,ETAB,JS ! GET IMAGINARY COMPONENT OF JS
ACRA=JSI*RES/I ! CALCULATE AC REACTANCE/LENGTH
DCRE=RES/AREA ! CALCULATE DC RESISTANCE/LENGTH
RAT=ACRE/DCRE ! AC/DC LOSS RATIO
/OUTPUT
/GOPR
*STATUS ! SHOW PARAMETER STATUS
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'JS ','JS ','LOSS '
LABEL(1,2) = '(RE)','(IM)','RATIO'
*VFILL,VALUE(1,1),DATA,10183,27328,2.33
*VFILL,VALUE(1,2),DATA,JSR,JSI,RAT
*VFILL,VALUE(1,3),DATA,ABS(JSR/10183),ABS(JSI/27328),ABS(RAT/2.33)
/COM
/OUT,vm185,vrt
/COM,------------------- VM185 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,' ',F11.2,' ',F11.2,' ',1F6.3)
/COM,-----------------------------------------------------------
/COM,
/COM,-----------------------------------------------------------
/COM,NOTE: THERE ARE VERIFIED RESULTS IN VM185 NOT CONTAINED IN
/COM,THIS TABLE
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm185,vrt
VM186 (Transient Analysis of a Slot Embedded Conductor) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM186
/PREP7
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM186, TRANSIENT ANALYSIS OF A SLOT EMBEDDED CONDUCTOR
C*** KONRAD, IEEE TRANS., MAGNETICS, VOL. MAG-18, NO. 1, JAN. 1982
/NOPR
ANTYPE,TRANS ! NONLINEAR TRANSIENT DYNAMIC ANALYSIS
ET,1,PLANE13 ! PLANE13, AZ DOF, (FOR AIR)
ET,2,PLANE13,6 ! PLANE13, AZ VOLT DOF, (FOR CONDUCTOR)
EMUNIT,MUZRO,1 ! SET MUZERO=1
MP,MURX,1,1 ! RELATIVE PERMEABILITY
MP,MURX,2,1 ! RELATIVE PERMEABILITY (CONDUCTOR)
MP,RSVX,2,1 ! RESISTIVITY (CONDUCTOR)
N,1
N,8,,7
FILL
NGEN,2,8,1,8,1,1
MAT,2
TYPE,2
E,1,2,10,9
EGEN,4,1,-1
MAT,1
TYPE,1
E,5,6,14,13
EGEN,3,1,-1
CP,1,AZ,1,9 ! COUPLE AZ TO ENSURE 1-D SOLUTION
*REPEAT,5,1,,1,1
ESEL,,MAT,,2
NSLE
CP,6,VOLT,ALL ! COUPLE VOLT IN CONDUCTOR
ESEL,ALL
NSEL,S,LOC,Y,7
D,ALL,AZ,0 ! FLUX-PARALLEL B.C.
NSEL,ALL
FINISH
/SOLU
EQSLV,JCG,1E-9 ! USE THE JACOBI CONJUGATE GRADIENT SOLVER
T=1E-8 ! INITIALIZE TIME PARAMETER
C=0 ! INITIALIZE COUNTER PARAMETER
N=80 ! NUMBER OF TIME INCREMENTS PER TURN
PI=2*ASIN(1) ! VALUE OF PI
CON=2*PI/N ! SET TIME INCREMENT
NEQIT,1 ! 1 ITERATION PER TIME STEP
*CREATE,LOAD ! CREATE MACRO TO SET UP LOAD STEPS
TIME,T
I=4*SIN(T) ! CALCULATE CURRENT
F,1,AMPS,I ! APPLY CURRENT TO A NODE IN CONDUCTOR
T=T+CON ! INCREMENT TIME
C=C+1 ! INCREMENT COUNTER
OUTRES,ALL,1
*IF,C,EQ,((N*.75)+1),THEN ! SET FOR PRINTOUT AT DESIRED TIME POINTS
OUTPR,,1
*ELSEIF,C,EQ,(N+1),THEN
OUTPR,,1
*ELSE
OUTPR,,0
*ENDIF
SOLVE
*END
*DO,I,1,81 ! REPEAT MACRO EXECUTION
*USE,LOAD ! EXECUTE MACRO
*ENDDO
FINISH
/POST26
NUMVAR,12 ! INCREASE STORAGE ARRAY SIZE
ESOL,2,1,,NMISC,6,JE ! STORE JE
*REPEAT,4,1,1
ESOL,6,1,,VOLUME ! STORE VOLUME
*REPEAT,4,1,1
PROD,2,2,6 ! CALCULATE IE=JE*VOLUME
*REPEAT,4,1,1,1
ADD,2,2,3,4,IE ! SUMM IE OVER ALL CONDUCTOR ELEMENTS
ADD,10,2,5,,IE ! IE TOTAL
ESOL,2,1,,SMISC,1,JS ! STORE JS
*REPEAT,4,1,1
PROD,2,2,6 ! CALCULATE IS=JS*VOLUME
*REPEAT,4,1,1,1
ADD,2,2,3,4,IS ! SUM IS OVER ALL CONDUCTOR ELEMENTS
ADD,11,2,5,,IS ! IS TOTAL
ESOL,2,1,,NMISC,7,JT ! STORE JT
*REPEAT,4,1,1
PROD,2,2,6 ! CALCULATE IT=JT*VOLUME
*REPEAT,4,1,1,1
ADD,2,2,3,4,IT ! SUM IT OVER ALL CONDUCTOR ELEMENTS
ADD,12,2,5,,IT ! IT TOTAL
/AXLAB,Y,CURRENT
/GROPT,AXNSC,2.0
PRVAR,10,11,12 ! PRINT EDDY, SOURCE, AND TOTAL CURRENT
PLVAR,10,11,12 ! DISPLAY EDDY, SOURCE, AND TOTAL CURRENT
FINISH
/POST1
SET,61,1,,,4.7124
*GET,A1,NODE,1,A,Z
*GET,A2,NODE,4,A,Z
*GET,A3,NODE,7,A,Z
SET,81,1,,,6.2832
*GET,A4,NODE,1,A,Z
*GET,A5,NODE,4,A,Z
*GET,A6,NODE,7,A,Z
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'NODE ','NODE ','NODE '
LABEL(1,2) = '1','4','7'
*VFILL,VALUE(1,1),DATA,-15.18,-14.68,-4
*VFILL,VALUE(1,2),DATA,A1,A2,A3
*VFILL,VALUE(1,3),DATA,ABS(A1/15.18),ABS(A2/14.68),ABS(A3/4)
SAVE,TABLE_1
*VFILL,VALUE(1,1),DATA,-3.26,-.92,0
*VFILL,VALUE(1,2),DATA,A4,A5,A6
*VFILL,VALUE(1,3),DATA,ABS(A4/3.26),ABS(A5/.92),
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm186,vrt
/COM,------------------- VM186 RESULTS COMPARISON --------------
/COM,
/COM, VECTOR POTENTIAL | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS AT T=(3*PI/2)
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS AT T=(2*PI)
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/COM,-----------------------------------------------------------
/COM,
/OUT
FINISH
*LIST,vm186,vrt
VM188 (Force Calculation on a Current Carrying Conductor) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM188
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
smrt,off
/TITLE, VM188, FORCE CALCULATION ON A CURRENT CARRYING CONDUCTOR
/COM, REF: MOON, FRANCIS C., MAGNETO-SOLID MECHANICS, PG. 418, 1984
ET,1,PLANE53 ! 8-NODE QUADRILATERAL MAGNETICS ELEMENT
ET,2,INFIN9 ! 2-D INFINITE BOUNDARY ELEMENT
EMUNIT,MKS ! MKS UNITS
MP,MURX,1,1 ! MATERIAL 1 RELATIVE PERMEABILITY=1.0
MP,MURX,2,1 ! MATERIAL 2 RELATIVE PERMEABILITY=1.0
D=.01 ! DEFINE GEOMETRY IN TERMS OF PARAMETERS
A=.012
T=.002
OB=.04 ! OUTER BOUNDARY SIZE
X1=D/2-T/2
X2=D/2+T/2
GP=.0002 ! GAP FOR THIN AIR LAYER NEXT TO CONDUCTOR
RECTNG,0,OB,0,OB ! DEFINE BOOLEAN AREAS
RECTNG,0,.012,0,.012
RECTNG,X1,X2,0,A/2
RECTNG,X1-GP,X2+GP,0,A/2+GP
AOVLAP,ALL ! OVERLAP AREAS
ASEL,S,AREA,,3
AATT,2 ! ASSIGN MATERIAL ATTRIBUTE TO CONDUCTOR
ASEL,ALL
KSEL,S,LOC,X,0,.012 ! SELECT KEYPOINTS FOR KESIZE SPEC.
KSEL,R,LOC,Y,0,.012
KESIZE,ALL,A/8 ! ASSIGN ELEMENT SIZE AT KEYPOINTS
KSEL,INVE
KESIZE,ALL,OB/5
KSEL,ALL
LSEL,S,LOC,X,OB ! SELECT FAR-FIELD BOUNDARY LINES
LSEL,A,LOC,Y,OB
TYPE,2
LMESH,ALL ! MESH WITH INFINITE LINE ELMENTS
LSEL,ALL
MSHK,0 ! FREE AREA MESH
MSHA,0,2D ! USING QUADS
TYPE,1
AMESH,ALL ! MESH AREAS WITH PLANE13
FINISH
/SOLU
ANTYPE,STATIC ! STATIC MAGNETICS ANALYSIS
ESEL,S,MAT,,2 ! SELECT CONDUCTOR ELEMENTS
BFE,ALL,JS,,,,1E6 ! APPLY CURRENT DENSITY TO CONDUCTOR
NSLE,S ! SELECT NODES IN CONDUCTOR REGION
BF,ALL,MVDI,1 ! APPLY VIRTUAL WORK DISPLACEMENT = 1
NSEL,INVE ! SELECT ALL OTHER NODES
BF,ALL,MVDI,0 ! APPLY VIRTUAL WORK DISPLACEMENT = 0
NSEL,ALL
ESEL,ALL
SOLVE
FINISH
/POST1
ETABLE,FMAGX,FMAG,X ! STORE J*B FORCE INFORMATION
ETABLE,FVWX,NMISC,3 ! STORE VIRTUAL WORK FORCE
SSUM ! SUM TABLE ENTRIES
*GET,FXL,SSUM,,ITEM,FMAGX ! GET J*B FORCE AS PARAMETER
FXL=FXL*2 ! TOTAL LORENTZ FORCE
*GET,FXVW,SSUM,,ITEM,FVWX ! GET VIRTUAL WORK FORCE AS PARAMETER
FXVW=FXVW*2 ! TOTAL VIRTUAL WORK FORCE
PATH,MAXWELL,4,,48 ! DEFINE PATH WITH NAME = "MAXWELL"
PPATH,1,,.012,0,0 ! DEFINE PATH POINTS BY LOCATION
PPATH,2,,.012,.012,0
PPATH,3,,0,.012,0
PPATH,4,,0,0,0
FOR2D ! COMMAND MACRO FOR MAXWELL STRESS FORCE CALC
FXM=FX*2 ! TOTAL MAXWELL FORCE (SYMMETRY)
*STATUS ! SHOW RESULTS
/PBC,PATH,1 ! ACTIVATE PATH B.C. FOR DISPLAY
PLF2D ! DISPLAY FLUX LINES
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'F (LRNZ) ','F (MAXW) ','F (VW) '
LABEL(1,2) = 'N/m','N/m','N/m'
*VFILL,VALUE(1,1),DATA,-9.684E-3,-9.684E-3,-9.684E-3
*VFILL,VALUE(1,2),DATA,FXL,FXM,FXVW
*VFILL,VALUE(1,3),DATA,ABS(FXL/(9.684E-3)),ABS(FXM/(9.684E-3)),ABS(FXVW/(9.684E-3))
/COM
/OUT,vm188,vrt
/COM,------------------- VM188 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.6,' ',F15.6,' ',1F10.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/nopr
/DELETE,for2d,out
/DELETE,scratch
/DELETE,vm188,gsav
FINISH
*LIST,vm188,vrt
VM189 (Hollow Sphere in a Uniform Magnetic Field) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM189
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
SMRT,OFF
/TITLE, VM189, HOLLOW SPHERE IN A UNIFORM MAGNETIC FIELD
/COM, COMPEL, VOL 7, NOS. 1&2, 89-101, 1988
/COM, TEAM WORKSHOP PROBLEM NO. 6
ET,1,SOLID97,1 ! MAGNETIC SOLID
ET,2,PLANE13 ! 2-D MAGNETIC OPTION
ET,3,SOLID96 ! MAGNETIC SCALAR SOLID
ET,4,INTER115 ! INTERFACE ELEMENT
EMUNIT,MKS ! MKS UNITS
MP,MURX,1,1 ! RELATIVE PERMEABILITY, AIR
MP,MURX,2,1 ! RELATIVE PERMEABILITY, ALUMINUM
MP,RSVX,2,2E-9 ! RESISTIVITY, ALUMINUM
RECTNG,0,MOD,0,MOD ! CREATE 2-D GEOMETRY
PCIRC,CIR,COR,,90
AOVLAP,ALL
CSYS,1
LSEL,S,LOC,Y,45 ! SET LINE SEGMENT DIVISIONS
LESIZE,ALL,,,12
LSEL,S,LOC,X,MOD,1
LESIZE,ALL,,,6
LSEL,S,LOC,X,(CIR+COR)/2
LESIZE,ALL,,,6
LSEL,S,LOC,X,CIR/2
LESIZE,ALL,,,12
KSEL,S,LOC,X,COR
KESIZE,ALL,.01
LSEL,S,LOC,X,COR+(MOD-COR)/2
LSEL,ALL
KSEL,ALL
CSYS,0
MSHK,1 ! MAPPED AREA MESH
MSHA,0,2D ! USING QUADS
ET,2,PLANE13 ! USE PLANE13 TO MESH AREA
TYPE,2
AMAP,4,4,2,5,6 ! MAP MESH AREA 4 BY CORNERS
AMESH,ALL ! MAP MESH REMAINING AREAS
ALLSEL,ALL
TYPE,1 ! SWITCH TO 3-D SOLID96 ELEMENT
ESIZE,,2 ! SPECIFIY CIRCUMFERENTIAL DIVISIONS
MAT,2
VROTAT,2,,,,,,4,6,20 ! CREATE 20 DEG. SLICE OF SHELL
TYPE,3 ! SCALAR POTENTIAL ELEMENTS
MAT,1
VROTAT,3,4,,,,,4,6,20 ! CREATE 20 DEG. SLICE OF AIR
MODMSH,DETACH ! DETACH SOLID MODEL
NUMMRG,NODE ! MERGE NODES
ESEL,S,TYPE,,2 ! SELECT 2-D ELEMENTS
EDELE,ALL ! DELETE 2-D ELEMENTS
ESEL,S,TYPE,,1
NSLE,S
ESEL,INVE
TYPE,4
ESURF ! MESH INTERFACE ELEMENTS
ALLSEL,ALL
LOCAL,11,1,,,,,-90 ! CREATE LOCAL CYLINDRICAL C.S.
NROTAT,ALL ! ROTATE NODES
ESEL,S,TYPE,,4 ! SELECT INTERFACE ELEMENTS
NSLE,S ! SELECT NODES AT INTERFACE
CSYS,2
NROTAT,ALL ! ROTATE INTERFACE NODES
FINISH
/SOLU
ANTYPE,HARMIC ! FULL HARMONIC ANALYSIS
HARFRQ,50 ! 50 HERTZ FREQUENCY
D,ALL,AX,0 ! SET A*N=0 AT INTERFACE
CSYS,11
ESEL,S,TYPE,,1
NSLE,S
NSEL,R,LOC,Y,0 ! SELECT NODES AT ONE PLANE
DSYM,ASYM,Y,11 ! APPLY FLUX-PARALLEL CONDITIONS
NSLE,S
NSEL,R,LOC,Y,-20 ! SELECT NODES AT OTHER PLANE
DSYM,ASYM,Y,11 ! APPLY FLUX-PARALLEL CONDITIONS
NSLE,S
NSEL,R,LOC,Z,0
NSEL,R,LOC,Y,-10
D,ALL,AY,0 ! APPLY FLUX-NORMAL CONSTRAINT
ESEL,S,TYPE,,1
NSLE,S
NSEL,R,LOC,X,0
D,ALL,AX,0,,,,AY,AZ ! FLUX-PARALLEL AT SYMMETRY AXIS
NSEL,ALL
MU=16*ATAN(1)*1E-7 ! DEFINE FREE-SPACE PERMEABILITY
CSYS,0
NSEL,S,LOC,Y,MOD
D,ALL,MAG,-MOD/MU ! APPLY POTENTIAL TO GIVE 1 TESLA FIELD
NSEL,S,LOC,Y,0
D,ALL,MAG,0
CSYS,11
ESEL,S,MAT,,2 ! SELECT SHELL ELEMENTS
NSLE,S ! SELECT NODES ATTACHED TO SHELL
NSEL,R,LOC,Y,0 ! SELECT NODES AT ONE PLANE
D,ALL,VOLT,0 ! SET TIME-INTEGRATED POTENTIAL TO 0
NSLE,S
NSEL,R,LOC,Y,-20 ! SELECT NODES AT OTHER PLANE
D,ALL,VOLT,0 ! SET TIME-INTEGRATED POTENTIAL TO 0
NSEL,ALL
ESEL,ALL
EQSLV,JCG,1e-6 ! JCG SOLVER
SOLVE ! SOLVE HARMONIC ANALYSIS
FINISH
/POST1
SET,1,1
/VIEW,,-1,.5,.5
/VUP,1,Y
/EDGE,,1
ESEL,S,MAT,,2 ! SELECT SHELL ELEMENTS
NSLE,S
PLVECT,JT,,,JT,VECT,ELEM,ON ! PLOT REAL EDDY CURRENTS
SET,1,1,,1
PLVECT,JT,,,JT,VECT,ELEM,ON ! PLOT IMAGINARY EDDY CURRENTS
POWERH ! CALCULATE POWER LOSS
PAVG=PAVG*8*4.5 ! CALCULATE FOR TOTAL SHELL
NSEL,ALL
ESEL,ALL
NORIG=NODE(0,0,0) ! GET NODE NUMBER AT ORIGIN
SET,1
*GET,BYR,NODE,NORIG,B,Y ! GET BY (REAL)
SET,1,1,,1
*GET,BYI,NODE,NORIG,B,Y ! GET BY (IMAGINARY)
BYM=SQRT(BYR**2+BYI**2) ! CALCULATE BY (MAGNITUDE)
*STATUS
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'BY (0,0)','POW LS '
LABEL(1,2) = ' T',' W'
*VFILL,VALUE(1,1),DATA,.0524,10062
*VFILL,VALUE(1,2),DATA,BYM,PAVG
*VFILL,VALUE(1,3),DATA,ABS(BYM/.0524),ABS(10062/PAVG)
/COM
/OUT,vm189,vrt
/COM,------------------- VM189 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.4,' ',F15.4,' ',1F10.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm189,vrt
VM190 (Ferromagnetic Inductor) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM190
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
smrt,off
/TITLE, VM190, FERROMAGNETIC INDUCTOR
/COM, CHAPMAN, "ELECTRIC MACHINERY FUNDAMENTALS", MCGRAW-HILL,
/COM, 1985, EXAMPLE 1-1, PG. 14
ET,1,SOLID98,10 ! 10-NODE TETRAHEDRAL, MAG OPTION
ET,2,INFIN47 ! INFINITE ELEMENT
ET,3,SOURC36 ! CURRENT ELEMENT
EMUNIT,MKS ! MKS UNITS
MP,MURX,1,1 ! RELATIVE PERMEABILITY OF AIR
MP,MURX,2,2500 ! RELATIVE PERMEABILITY OF IRON
R,1,1,200,.02,.25 ! COIL DIMENSIONS AND CURRENT
LOCAL,11,0,-.325 ! SHIFT ORIGIN TO CENTER OF MODEL
WPCSYS,,11 ! WORKING PLANE FOR SOLID MODELLING
N,1,.125,0,0 ! CREATE NODES TO LOCATE COIL
N,2,.235,0,0
N,3,.235,.235,0
TYPE,3
E,2,3,1 ! DEFINE COIL
BLOCK,.05,.20,0,.05,0,.45 ! CREATE SOLID MODEL OF IRON
BLOCK,.20,.50,0,.05,.30,.45
BLOCK,.50,.60,0,.05,0,.45
VGLUE,ALL
BLOCK,0,.65,0,.10,0,.50 ! CREATE SOLID MODEL OF AIR
VOVLAP,ALL ! OVERLAP AIR AND IRON
ASEL,S,AREA,,8,10
ASEL,A,AREA,,4
MSHK,0 ! FREE MESH
MSHA,1,3D ! USING TETS
MSHA,1,2D ! USING TRIS
ESIZE,.10
TYPE,2
AMESH,ALL ! MESH EXTERIOR BOUNDARY WITH INFIN47
TYPE,1
MAT,2
VMESH,4,6 ! MESH IRON
MAT,1
VMESH,2 ! MESH AIR
NSEL,S,LOC,Z,0
D,ALL,MAG,0 ! SET FLUX-NORMAL SYMMETRY CONDITION
NSEL,ALL
FINISH
/SOLU
ANTYPE,STATIC ! STATIC MAGNETIC FIELD ANALYSIS
MAGSOLV,4
FINISH
/POST1
PATH,IRON,7,,48 ! DEFINE PATH WITH NAME = "IRON"
PPATH,1,,-.2,0,0 ! DEFINE PATHS POINTS BY
PPATH,2,,-.2,0,.20
PPATH,3,,-.2,0,.375
PPATH,4,,.025,0,.375
PPATH,5,,.225,0,.375
PPATH,6,,.225,0,.20
PPATH,7,,.225,0,0
*CREATE,MAC ! CREATE MACRO FOR MMF CALCULATION
PDEF,HX,H,X ! INTERPOLATE H FIELD TO PATH
PDEF,HY,H,Y
PDEF,HZ,H,Z
PVECT,TANG,TX,TY,TZ ! INTERPOLATE UNIT TANGENTS
PDOT,D,HX,HY,HZ,TX,TY,TZ ! PERFORM DOT PRODUCT
PCALC,INTG,MMF,D,S ! INTEGRATE OVER PATH
*GET,MMF,PATH,,LAST,MMF ! GET MMF
MMF=MMF*2 ! MULTIPLY BY 2 FOR SYMMETRY
*STATUS,MMF
*END
ESEL,S,MAT,,2 ! SELECT IRON ELEMENTS
*USE,MAC ! USE MACRO TO CALCULATE MMF
/VIEW,,1,-3,1
/VUP,1,Z
/TRIAD,OFF
/PBC,PATH,1 ! SHOW PATH ON DISPLAY
/AUTO
WPSTYL,,,,,,,OFF
NSLE,S ! SELECT NODES ATTACHED TO IRON
/COM *** THE FOLLOWING ANNOTATION COMMANDS ARE ***
/COM *** TYPICALLY GENERATED INTERACTIVELY ***
/ANUM,1,12,-.28056,.71310 ! ANNOTATION NUMBER, TYPE, AND HOT SPOT
/LSYMBOL,-.282,.511,269,4,1.0 ! ANNOTATION SYMBOL DEFINITION - ARROW
/LINE,-.282,.511,-.279,.915 ! ANNOTATION LINE DEFINITION
/ANUM,2,4,-.21690,.91150
/LINE,-.282,.911,-.151,.911
/ANUM,3,1.14734,.93021
/TLABEL,-.133,.930,CONTOUR PATH FOR ! ANNOTATION LOCATION AND TEXT
/ANUM,4,1,.91097E-01,.87406
/TLABEL,-.136,.874,LINE INTEGRAL
/ANNOT,ON
/TITLE, VM190: MAGNETIC FLUX DENSITY
PLNSOL,B,SUM ! DISPLAY FLUX DENSITY IN IRON
FINISH
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'MMF DRP '
LABEL(1,2) = 'A-t'
*VFILL,VALUE(1,1),DATA,200
*VFILL,VALUE(1,2),DATA,MMF
*VFILL,VALUE(1,3),DATA,ABS(MMF/200),
/COM
/OUT,vm190,vrt
/COM,------------------- VM190 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.2,' ',F10.2,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/DELETE,MAC
/DELETE,scratch
/DELETE,magsolv,out
*LIST,vm190,vrt
VM191 (Hertz Contact Between Two Cylinders) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM191
/TITLE, VM191, HERTZ CONTACT BETWEEN TWO CYLINDERS
/COM "FINITE ELEMENT ANALYSIS OF HERTZ CONTACT PROBLEM"
/COM N. CHANDRASEKARAN, W.E. HAISLER, R.E. GOFORTH,
/COM FINITE ELEMENTS IN ANALYSIS AND DESIGN 3, 1987, PP 39-56.
/PREP7
SMRT,OFF
/COM 2-D ANALYSIS USING PLANE42 AND CONTAC48
ANTYPE,STATIC
ET,1,PLANE42 ! 2-D SOLID ELEMENTS
MP,EX,1,30000 ! SMALLER CYLINDER PROPERTIES
MP,NUXY,1,0.25
MP,EX,2,29120 ! LARGER CYLINDER PROPERTIES
MP,NUXY,2,0.30
ET,2,CONTAC48 ! 2-D CONTACT ELEMENTS
R,2,3E4,,,,0.001 ! KN=30000, TOLS=0.001
CSYS,1
K,1 ! CREATE BIGGER CYLINDER
K,2,13
K,3,13,82
K,4,13,90
K,5,11,90
L,1,5
L,2,3
LESIZE,ALL,,,7
L,3,4 ! TARGET SURFACE (LINE 3)
LOCAL,11,1,,13
L,3,5
CSYS,1
A,1,2,3,5
A,5,3,4,4
MAT,2
MSHK,1 ! MAPPED AREA MESH
MSHA,0,2D ! USING QUADS
ESIZE,,4
AMESH,1,2
LOCAL,12,1,,23-1E-5,,-90 ! INTRODUCE SLIGHT INTERFERENCE
K,11 ! CREATE SMALLER CYLINDER
K,12,10
K,13,10,8
K,14,10,90
K,15,8
L,11,15
L,13,14
LESIZE,7,,,6
LESIZE,8,,,6
L,12,13 ! CONTACTOR SURFACE (LINE 9)
CSYS,11
L,13,15
CSYS,12
MAT,1
A,12,13,15,15
A,15,13,14,11
ESIZE,,6
AMESH,3,4
LSEL,S,LINE,,9 ! SELECT CONTACT NODES ON SMALLER CYLINDER
NSLL,,1
CM,CYL1,NODE ! CONTACT NODES COMPONENT
LSEL,S,LINE,,3 ! SELECT TARGET NODES ON BIGGER CYLINDER
NSLL,,1
CM,CYL2,NODE ! TARGET NODES COMPONENT
NSEL,ALL
TYPE,2
REAL,2
GCGEN,CYL1,CYL2 ! GENERATE ASYMMETRIC CONTACT ELEMENTS
CSYS,0
NSEL,S,LOC,Y,23 ! SELECT TOP EDGE OF MODEL
CP,1,UY,ALL ! COUPLE NODES ON TOP EDGE
*GET,NC,NODE,,NUM,MIN ! GET LOWEST NODE NUMBER (MASTER)
NSEL,S,LOC,X ! SYMMETRY CONSTRAINTS
D,ALL,UX
NSEL,S,LOC,Y
D,ALL,UY
NSEL,ALL
FINISH
/SOLU
D,NC,UY,-0.005 ! APPLY SMALL DISPLACEMENT TO ENGAGE CONTACT
SOLVE ! SOLVE FIRST LOAD STEP
DDELE,NC,UY ! DELETE IMPOSED DISPLACEMENT
F,NC,FY,-1600 ! APPLY HALF LOAD ON (SYMMETRY) MODEL
SOLVE ! SOLVE SECOND LOAD STEP
FINISH
/POST1
NSEL,,LOC,Y,23 ! SELECT TOP EDGE OF SMALLER CYLINDER
*GET,D,NODE,NC,U,Y ! GET APPROACH DISTANCE (D)
ESEL,S,TYPE,,2 ! SELECT CONTACT ELEMENTS
ETABLE,NSTAT,NMISC,1 ! STORE CONTACT STATUS
ESEL,R,ETAB,NSTAT,2,2 ! SELECT ELEMENTS WITH CONTACT (STAT=2)
CMSEL,S,CYL1 ! SELECT CONTACT COMPONENT NODES
NSLE,R ! RESELECT NODES WITH CONTACT
NSORT,LOC,X,1 ! SORT CONTACT NODES BY ASCENDING X LOCATION
*GET,B,SORT,,MAX ! GET SEMI-CONTACT LENGTH (B)
*STATUS
*GET,DEF,NODE,1,U,X
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'AP DIS ','S-CON LEN '
LABEL(1,2) = ' mm',' mm'
*VFILL,VALUE(1,1),DATA,-.4181,1.2
*VFILL,VALUE(1,2),DATA,D,B
*VFILL,VALUE(1,3),DATA,ABS(D/.4181),ABS(B/1.2)
SAVE,TABLE_1
FINISH
/CLEAR, NOSTART
/PREP7 $SMRT,OFF
ANTYPE,STATIC
ET,1,SOLID45 ! 3-D SOLID ELEMENTS
MP,EX,1,30000 ! SMALLER CYLINDER PROPERTIES
MP,NUXY,1,0.25
MP,EX,2,29120 ! LARGER CYLINDER PROPERTIES
MP,NUXY,2,0.30
ET,2,CONTAC49 ! 3-D CONTACT ELEMENTS
R,2,3E4,,,,0.001 ! KN=30000, TOLS=0.001
CSYS,1
K,1 ! CREATE LOWER BIGGER CYLINDER
K,2,13
K,3,13,82
K,4,13,90
K,5,11,90
KGEN,2,1,5,1,,,1,100 ! UNIT THICKNESS SLICE
L,1,5
L,2,3
L,101,105
L,102,103
LESIZE,ALL,,,7
L,1,101
*REPEAT,5,1,1
LESIZE,5,,,1
*REPEAT,5,1
LOCAL,11,1,,13
L,3,5
L,103,105
CSYS,1
MAT,2
MSHK,1 ! MAPPED VOLUME MESH
MSHA,0,3D ! USING HEX
ESIZE,,4
V,1,2,3,5,101,102,103,105
V,5,3,4,4,105,103,104,104
VMESH,ALL
LOCAL,12,1,,23-1E-5,,-90 ! INTRODUCE SLIGHT INTERFERENCE
K,11 ! CREATE UPPER SMALLER CYLINDER
K,12,10
K,13,10,8
K,14,10,90
K,15,8
KGEN,2,11,15,1,,,1,100
L,11,15
L,13,14
LESIZE,18,,,6
LESIZE,19,,,6
L,11,111
*REPEAT,5,1,1
LESIZE,20,,,1
*REPEAT,5,1
CSYS,11
L,13,15
L,113,115
CSYS,12
MAT,1
ESIZE,,6
V,12,13,15,15,112,113,115,115
V,15,13,14,11,115,113,114,111
VMESH,3,4
ASEL,S,AREA,,12
NSLA,,1
CM,CYL1,NODE ! CONTACT NODES COMPONENT
ASEL,S,AREA,,8 ! SELECT TARGET NODES ON BIGGER CYLINDER
NSLA,,1
CM,CYL2,NODE ! TARGET NODES COMPONENT
NSEL,ALL
TYPE,2
REAL,2
GCGEN,CYL1,CYL2 ! GENERATE ASYMMETRIC CONTACT ELEMENTS
CSYS,0
NSEL,S,LOC,Y,23 ! SELECT TOP EDGE OF MODEL
CP,1,UY,ALL ! COUPLE NODES ON TOP EDGE
*GET,NC,NODE,,NUM,MIN ! GET LOWEST NODE NUMBER (MASTER)
NSEL,S,LOC,X ! SYMMETRY CONSTRAINTS
D,ALL,UX
NSEL,S,LOC,Y
D,ALL,UY
NSEL,S,LOC,Z
D,ALL,UZ
NSEL,ALL
FINISH
/SOLU
D,NC,UY,-0.001 ! APPLY SMALL DISPLACEMENT TO ENGAGE CONTACT
SOLVE ! SOLVE FIRST LOAD STEP
DDELE,NC,UY ! DELETE IMPOSED DISPLACEMENT
F,NC,FY,-1600 ! APPLY HALF LOAD ON (SYMMETRY) MODEL
SOLVE ! SOLVE SECOND LOAD STEP
FINISH
/POST1
NSEL,,LOC,Y,23 ! SELECT TOP EDGE OF SMALLER CYLINDER
*GET,D,NODE,NC,U,Y ! GET APPROACH DISTANCE (D)
ESEL,S,TYPE,,2 ! SELECT CONTACT ELEMENTS
ETABLE,NSTAT,NMISC,1 ! STORE CONTACT STATUS
ESEL,R,ETAB,NSTAT,2,2 ! SELECT ELEMENTS WITH CONTACT (STAT=2)
CMSEL,S,CYL1 ! SELECT CONTACT COMPONENT NODES
NSLE,R ! RESELECT NODES WITH CONTACT
NSORT,LOC,X,1 ! SORT CONTACT NODES BY ASCENDING X LOCATION
*GET,B,SORT,,MAX ! GET SEMI-CONTACT LENGTH (B)
*STATUS
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'AP DIS ','S-CON LEN '
LABEL(1,2) = ' mm',' mm'
*VFILL,VALUE(1,1),DATA,-.4181,1.2
*VFILL,VALUE(1,2),DATA,D,B
*VFILL,VALUE(1,3),DATA,ABS(D/.4181),ABS(B/1.2)
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm191,vrt
/COM,------------------- VM191 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS USING 2-D ANALYSIS:
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.4,' ',F10.4,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS USING 3-D ANALYSIS:
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.4,' ',F10.4,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm191,vrt
VM192 (Cooling of a Billet by Radiation) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM192
/PREP7
SMRT,OFF
/TITLE, VM192, COOLING OF A BILLET BY RADIATION
/COM THERMAL RADIATION HEAT TRANSFER, SIEGEL AND HOWELL, 2ND EDITION,
/COM PG. 229, PROBLEM NO. 21.
ET,1,SOLID70 ! 3-D THERMAL SOLID ELEMENT
ET,2,SURF152,,,,1,1 ! 3-D THERMAL SURFACE EFFECT ELEMENTS
KEYOPT,2,9,1 ! RADIATION OPTION
R,2,1,0.1712E-8 ! FORM FACTOR = 1, STEFAN-BOLTZMANN CONSTANT
MP,KXX,1,10000 ! ARBITRARY CONDUCTIVITY
MP,C,1,0.11
MP,DENS,1,487.5
MP,EMIS,2,1 ! BLACK BODY EMISSIVITY
BLOCK,,2,,2,,4
ESIZE,,1
VMESH,1 ! MESH WITH A SINGLE SOLID70 ELEMENT
TYPE,2
REAL,2
MAT,2
N,100,5,5,5 ! EXTRA "SPACE" NODE FOR RADIATION
ESURF,100 ! GENERATE SURF152 ELEMENTS
FINISH
/SOLU
SOLCONTROL,0
ANTYPE,TRANS ! TRANSIENT ANALYSIS
D,100,TEMP,530 ! SPECIFY SURROUNDING ABSOLUTE TEMPERATURE
TUNIF,2000 ! INITIAL BILLET ABSOLUTE TEMPERATURE
AUTOTS,ON
KBC,1 ! STEP SURROUNDING TEMPERATURE IN FIRST TIME STEP
DELTIM,0.005 ! INITIAL (MINIMUM) INTEGRATION TIME STEP
OUTRES,,ALL
OUTPR,NSOL,LAST
TIME,3.7 ! TRANSIENT TIME SPAN
SOLVE
FINISH
VM193 (Adaptive Analysis of 2-D Heat Transfer with Convection) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM193
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
SMRT,OFF
/TITLE, VM193, TWO DIMENSIONAL HEAT TRANSFER WITH CONVECTION
C*** "THE STANDARD NAFEMS BENCHMARKS", TEST NO. T4,
C*** NAFEMS, REV 3, OCTOBER 1990.
ANTYPE,STATIC
ET,1,PLANE55
MP,KXX,1,52.0
K,1
K,2,.6
K,3,.6,1.0
K,4,,1.0
K,5,.6,.2
L,1,2
L,2,5
L,5,3
L,3,4
L,4,1
AL,ALL
DK,1,TEMP,100,,1
DK,2,TEMP,100,,1
SFL,2,CONV,750.0,,0.0
SFL,3,CONV,750.0,,0.0
SFL,4,CONV,750.0,,0.0
FINISH
ADAPT,10,,5,0.2,1 ! FINAL PERCENT ERROR NEAR 5% WITHIN 10 LOOPS
/POST1
PLNSOL,TEMP ! DISPLAY TEMP CONTOURS IN FINAL MESH
*GET,TEPC,PRERR,,TEPC
KSEL,,,,5
NSLK
*GET,N1,NODE,,NUM,MAX
*GET,TEMP1,NODE,N1,TEMP
*STATUS
*DIM,VALUE,,1,3
*DIM,LABEL,CHAR,1,2
LABEL(1,1) = 'TEMP '
LABEL(1,2) = 'DEG C'
*VFILL,VALUE(1,1),DATA,18.3
*VFILL,VALUE(1,2),DATA,TEMP1
*VFILL,VALUE(1,3),DATA,ABS(TEMP1/18.3)
/COM
/OUT,vm193,vrt
/COM,------------------- VM193 RESULTS COMPARISON --------------
/COM,
/COM, LOAD STP 4 | TARGET | ANSYS | RATIO
/COM,
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F8.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm193,vrt
VM194 (Element Birth/Death in a Fixed Bar with Thermal Loading) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM194
/PREP7
/TITLE, VM194, ELEMENT BIRTH/DEATH IN A FIXED BAR WITH THERMAL LOADING
/COM ANY STANDARD MECHANICS OF MATERIALS TEXT
ET,1,LINK1 ! 2-D SPAR ELEMENTS
MP,EX,1,30E6 ! BAR MATERIAL PROPERTIES
MP,ALPX,1,.00005
MP,EX,2,30E6
MP,ALPX,2,.00005 ! MATERIAL PROPERTIES FOR RE-BORN ELEMENT
MP,REFT,2,100 ! REFERENCE TEMPERATURE FOR ELEMENT BIRTH
R,1,1.0
N,1
N,4,10
FILL
E,1,2
EGEN,3,1,-1 ! GENERATE THREE ELEMENTS
FINISH
/SOLU
ANTYPE,STATIC
D,1,ALL,,,4,3 ! FIX BOTH ENDS OF THE BAR
TREF,0 ! ZERO REFERENCE TEMPERATURE
TUNIF,100 ! UNIFORM TEMPERATURE THERMAL LOAD
NROPT,FULL
OUTPR,BASIC,ALL
SOLVE
EKILL,2 ! KILL CENTER ELEMENT
SOLVE
EALIVE,2 ! RESURRECT CENTER ELEMENT
MPCHG,2,2 ! AND CHANGE TO MATERIAL 2 FOR STRAIN-FREE BIRTH
SOLVE
TUNIF,0 ! REMOVE THERMAL LOADING
SOLVE
/POST1
ESEL,S,ELEM,,1,2
ETABLE,FO,SMISC,1
ESORT,FO
*GET,F,SORT,,MAX
ETABLE,STR,LEPTH,1
ESORT,STR
*GET,S1,SORT,,MAX
*GET,S2,SORT,,MIN
*STATUS
*DIM,VALUE,,3,3
*DIM,LABEL,CHAR,3,2
LABEL(1,1) = 'PRESS ','MAX STRS ','MAX STRS '
LABEL(1,2) = 'psi','EL(1)','EL(2)'
*VFILL,VALUE(1,1),DATA,150000,0,-.005
*VFILL,VALUE(1,2),DATA,F,S1,S2
*VFILL,VALUE(1,3),DATA,ABS(F/150000),0,ABS(S2/.005)
/COM
/OUT,vm194,vrt
/COM,------------------- VM194 RESULTS COMPARISON --------------
/COM,
/COM, LOAD STP 4 | TARGET | ANSYS | RATIO
/COM,
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.3,' ',F10.3,' ',1F8.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm194,vrt
VM196 (Counter-Balanced Loads on a Block) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM196
/PREP7
MP,PRXY,,0.3
SMRT,OFF
/TITLE, VM196, COUNTER-BALANCED LOADS ON A BLOCK
C*** ANY BASIC MECHANICS BOOK
ET,1,SOLID45
MP,EX,1,70E9 ! ALUMINUM
MP,DENS,1,2712
MP,PRXY,1,0.3
WPOFFS,,,300 ! AXIS OF ROTATION 300 M
BLOCK,-1,1,-1,1,0,3 ! 3 M HIGH BY 2 M SQUARE
ESIZE,1
VMESH,ALL
DK,1,ALL ! CONSTRAIN 6 DOF SUCH THAT NO ROTATIONS OCCUR
DK,4,UX
DK,6,UY
DK,7,UZ
FINISH
/SOLU
ANTYPE,STATIC
IRLF,1 ! INERTIA RELIEF CALCULATIONS
FK,5,FX,-2000
FK,5,FY,3000
OUTPR,RSOL,1 ! PRINT REACTION SOLUTION
SOLVE
IRLIST ! LIST INERTIA RELIEF LOADS AND ACCELERATIONS
*GET,MX,ELEM,,MMOR,X
*GET,MY,ELEM,,MMOR,Y
*GET,MZ,ELEM,,MMOR,Z
*DIM,LABEL,CHAR,3,2
*DIM,VALUE,,3,3
LABEL(1,1) = 'MOM X ','MOM Y','MOM Z'
LABEL(1,2) = 'N-m','N-m','N-m'
*VFILL,VALUE(1,1),DATA,-909000,-606000,-5000
*VFILL,VALUE(1,2),DATA,MX,MY,MZ
*VFILL,VALUE(1,3),DATA,ABS(MX/909000),ABS(MY/606000),ABS(MZ/5000)
/COM
/OUT,vm196,vrt
/COM,------------------- VM196 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,-----------------------------------------------------------
/COM,
/COM,-----------------------------------------------------------
/COM,NOTE: THERE ARE VERIFIED RESULTS IN VM196 NOT CONTAINED IN
/COM,THIS TABLE
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST,vm196,vrt
VM197 (GES Write/Read for Thick-Walled Cylinder with Spherical End Caps) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM197
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
/TITLE, VM197, IGES WRITE/READ FOR THICK-WALLED CYLINDER WITH SPHERICAL END CAPS
/COM, ANY BASIC GEOMETRY TEXT
CYLIND,RADIUS,RADIUS+THICK,0,LENGTH,90,270 ! CREATE HOLLOW CYLINDER
WPROTA,,90 ! ROTATE WORKING PLANE FOR CREATION OF END CAPS
SPHERE,RADIUS,RADIUS+THICK,180,270 ! CREATE END CAP AT ONE END OF CYLINDER
WPOFFS,,LENGTH ! MOVE WORKING PLANE TO OTHER END OF CYL
SPHERE,RADIUS,RADIUS+THICK,90,180 ! CREATE END CAP AT OTHER END OF CYLINDER
/VIEW,1,1,1,1 ! SET VIEW DIRECTION FOR DISPLAY
/TYPE,1,4 ! SET TYPE OF DISPLAY TO PRECISE HIDDEN
/TRIAD,OFF ! TURN OFF COORDINATE SYSTEM TRIAD ON DISPLAYS
WPSTYL,,,,,,,OFF ! TURN OFF WORKING PLANE TRIAD ON DISPLAYS
VPLOT
IGESOUT,VM197,IGS ! WRITE IGES FILE NAMED VM197.IGS
FINISH
PARSAV,ALL ! SAVE PARAMETERS (TO BE AVAILABLE AFTER CLEAR)
/CLEAR, NOSTART
/AUX15
IOPTN,IGES,NODEFE ! SELECT NO DEFEATURING IGES IMPORT
IOPTN,MERGE,YES
IOPTN,SOLID,NO
IGESIN,VM197,IGS ! READ IN IGES FILE NAMED VM197.IGS
FINISH
/PREP7
VSUM ! PERFORM VSUM FOR SUBSEQUENT *GET OF TOT VOLUME
PARRES ! RESUME PARAMETERS SAVED PRIOR TO CLEAR
*GET,VOLUME,VOLU,,VOLU ! VOLUME = TOTAL VOLUME OF GEOMETRY THAT WAS READ
! IN FROM IGES FILE (FROM VSUM)
PI=(4.0)*ATAN(1.0) ! HAVE ANSYS CALCULATE VALUE FOR PI
CAPVOL=2/3*PI*(((RADIUS+THICK)**3)-(RADIUS**3)) ! CALC VOL OF TWO END CAPS
CYLVOL=.5*LENGTH*PI*(((RADIUS+THICK)**2)-(RADIUS**2)) ! CALC VOL OF THE CYLINDER
CALCVOL=CAPVOL+CYLVOL ! TOTAL CALC'D VOL = CYL VOL + END CAPS VOL
NORMVOL=VOLUME/CALCVOL ! NORMALIZE ANSYS'S VOLUME BY CALCULATED VOLUME
CAPVOL= ! DELETE PARAMETERS USED IN CALC. OF VOLUME
CYLVOL=
PI=
RADIUS= ! DELETE PARAMETERS USED IN MODEL GENERATION
THICK=
LENGTH=
*DIM,VALUE,,1,3
*DIM,LABEL,CHAR,1,2
LABEL(1,1) = 'VOLUME '
LABEL(1,2) = 'NO UNTS'
*VFILL,VALUE(1,1),DATA,79063
*VFILL,VALUE(1,2),DATA,VOLUME
*VFILL,VALUE(1,3),DATA,ABS(VOLUME/79063)
/COM
/OUT,vm197,vrt
/COM,------------------- VM197 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.0,' ',F10.0,' ',1F8.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/NOPR
/DELETE,VM197,IGS
FINISH
*LIST,vm197,vrt
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM198
/SHOW,,,1
EL='VISCO106'
*DIM,A,CHAR,2 ! DEFINE AND DIMENSION THE CHARACTER PARAMETER A
A(1)='UX','UY' ! SET CHARACTER STRINGS AS VALUES OF CHARACTER
! PARAMETER A
! ANALYZE THE PROBLEM WITH ELEMENT VISCO106
/PREP7
/TITLE, VM198, LARGE STRAIN IN-PLANE TORSION TEST (%EL%)
! "SOME COMPUTATIONAL ASPECTS ....",NAGTEGAAL ET AL.
ET,1,VISCO106 ! 2-D 4-NODE VISCOPLASTIC ELEMENT
MP,EX,1,7200
MP,NUXY,1,0.33
TB,BISO,1 ! BILINEAR ISOTROPIC HARDENING
TBDATA,1,10 ! YIELD STRESS
TBDATA,2,40 ! TANGENT MODULUS
SAVE ! SAVE DATABASE FOR LATER USE
CSYS,1
N,1,10
N,11,20
FILL
NGEN,2,100,1,11,1,,3
E,1,2,102,101
EGEN,10,1,1,1
! APPLY BOUNDARY CONDITIONS
NSEL,S,LOC,X,20
D,ALL,ALL ! CLAMP OUTER SURFACE
NSEL,ALL
! ROTATE APPROPRIATE NODES AND APPLY COUPLING
LOCAL,11,0,,,,3
NROTAT,101,111
CSYS,0
CP,1,A(1),2,102 ! USE A(1) VALUE FOR DOF LABEL
CP,2,A(2),2,102 ! USE A(2) VALUE FOR DOF LABEL
CPSGEN,9,1,1,2,1
FINISH
/OUTPUT,SCRATCH ! DIVERT VOLUMINOUS SOLUTION OUTPUT
/SOLU
! CREATE SOLUTION OPTIONS AND LOADING MACRO FOR MULTIPLE USE
*CREATE,SOLD,MAC
SOLCONTROL,0
ANTYPE,STATIC
NLGEOM,ON ! LARGE STRAIN OPTION ACTIVATED
NEQIT,100 ! 100 EQUILIBRIUM ITERATIONS ALLOWED
CUTC,PLSLIMIT,0.5 ! RESET MAXIMUM PLASTIC INCREMENTAL STRAIN
D,ARG1,ARG2,,,ARG3,ARG4,ARG5
SOLVE
OUTRES,ESOL,1 ! STORE RESULTS FOR EVERY SUBSTEP
NSTP = 10 ! NO. OF LOAD STEPS USED
T1 = 60/NSTP ! ROTATION PER LOAD STEP
T2 = 3.1415927/180 ! PARAMETER FOR FURTHER CALCULATIONS
T33 = T1*T2 ! DEGREES TO RADIANS CONVERSION
*DO,I,1,NSTP ! USE DO LOOP FOR LOADING
T3 = (I*T33) ! CURRENT ANGLE
T4 = (10*SIN(T3)) ! UY DISPLACEMENT
T5 = (10*COS(T3))
T5 = (T5-10) ! UX DISPLACEMENT
D,ARG1,ARG2,T5,,ARG3,ARG4
D,ARG1,ARG5,T4,,ARG3,ARG4
SOLVE
*ENDDO
*END
SOLD,1,A(1),101,100,A(2) ! USE A(1) AND A(2) AS ARG2 AND ARG5 VALUES FOR
FINISH ! DOF LABELS
/OUTPUT
! CREATE POST PROCESSING MACRO FOR MULTIPLE USE
*CREATE,POSP,MAC
!COM PLOT THE FINAL DISPLACED GEOMETRY USING POST1
/POST1
SET,NSTP+1
/DSCALE,1,1
/DIST,1,13
PLDISP,1
/DIST,1
FINISH
! OBTAIN MAXIMUM SHEAR STRESS USING POST26
/POST26
ESOL,2,1,1,S,1
ESOL,3,1,1,S,3
FILLDATA,5,1,,,0,6
ADD,4,2,3,,SHEAR,,,-1/2,1/2 ! COMPUTE MAX. SHEAR USING
! PRINCIPAL STRESSES
PRVAR,4,5
/GRID,1
/XRANGE,0,60
/YRANGE,-60,0
/AXLAB,X,ROTATION (DEGREES)
/AXLAB,Y,SHEAR STRESS (PSI)
XVAR,5
PLVAR,4
*GET,P1,VARI,4,RTIME,11
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'PRS MAX '
LABEL(1,2) = 'psi'
*VFILL,VALUE(1,1),DATA,-48
*VFILL,VALUE(1,2),DATA,P1
*VFILL,VALUE(1,3),DATA,ABS(P1/48)
FINISH
*END
POSP
SAVE,TABLE_1
/CLEAR, NOSTART
! ANALYZE THE SAME PROBLEM WITH ELEMENT VISCO108
/PREP7
RESUME
EL='VISCO108'
/STATUS,TITLE ! DISPLAY TITLE WITH NEW VALUE OF PARAMETER EL
ET,1,VISCO108 ! 2-D, 8-NODE VISCOPLASTIC ELEMENT
N,1,10
N,21,20
FILL
CSYS,1
NGEN,3,30,1,21,1,,1.5
E,1,3,63,61,2,33,62,31
EGEN,10,2,1,1,1
NSLE
NSEL,INVE
NDELE,ALL
! APPLY BOUNDARY CONDITIONS
NSEL,S,LOC,X,20
D,ALL,ALL
! ROTATE APPROPRIATE NODES AND APPLY COUPLING
NSEL,S,NODE,,61,81
LOCAL,11,0,,,,3
NROTAT,ALL
NSEL,S,NODE,,31,51,2
LOCAL,12,0,,,,1.5
NROTAT,ALL
NSEL,ALL
CSYS,0
CP,1,A(1),3,33,63
CP,2,A(2),3,33,63
CPSGEN,9,2,1,2,1
CP,21,A(1),2,62
CP,22,A(2),2,62
CPSGEN,10,2,21,22,1
FINISH
/OUTPUT,SCRATCH ! DIVERT VOLUMINOUS SOLUTION OUTPUT
/SOLU
SOLD,1,A(1),61,30,A(2)
FINISH
/OUTPUT
! POSTPROCESS SOLUTION RESULTS USING POSP MACRO
POSP
SAVE,TABLE_2
/CLEAR, NOSTART
! ANALYZE THE SAME PROBLEM WITH ELEMENT VISCO107
/PREP7
RESUME
EL='VISCO107'
/STATUS,TITLE ! DISPLAY TITLE WITH NEW VALUE OF PARAMETER EL
ET,1,VISCO107
CSYS,1
N,1,10
N,11,20
FILL
NGEN,2,100,1,11,1,,3
NGEN,2,200,1,111,1,,,1
E,1,2,102,101,201,202,302,301
EGEN,10,1,1,1
! APPLY BOUNDARY CONDITIONS
D,ALL,UZ
NSEL,S,LOC,X,20
D,ALL,ALL ! CLAMP OUTER SURFACE
NSEL,ALL
! ROTATE APPROPRIATE NODES AND APPLY COUPLING
LOCAL,11,0,,,,3
NROTAT,101,111
NROTAT,301,311
CSYS,0
CP,1,A(1),2,102,202,302
CP,2,A(2),2,102,202,302
CPSGEN,9,1,1,2,1
FINISH
/OUTPUT,SCRATCH ! DIVERT VOLUMINOUS SOLUTION OUTPUT
/SOLU
SOLD,1,A(1),301,100,A(2)
FINISH
/OUTPUT
POSP ! POSTPROCESS RESULTS USING POSP MACRO
SAVE,TABLE_3
RESUME,TABLE_1
/COM
/OUT,vm198,vrt
/COM,------------------- VM198 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS USING VISCO106
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F5.3)
/NOPR
RESUME,TABLE_3
/GOPR
/COM,
/COM,RESULTS USING VISCO107
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS USING VISCO108
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/NOPR
/DELETE,POSP,MAC
/DELETE,SOLD,MAC
/DELETE,SCRATCH
/DELETE,TABLE_1
/DELETE,TABLE_2
/DELETE,TABLE_3
FINISH
*LIST,vm198,vrt
VM199 (Viscoplastic Analysis of a Body Undergoing Shear Deformation) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM199
/PREP7
/TITLE, VM199, VISCOPLASTIC ANALYSIS OF A BODY UNDERGOING SHEAR DEFORMATION
! " AN IMPLICIT STRESS UPDATE ....", LWO ET AL.
! VISCO106 ELEMENT
ET,1,VISCO106 ! 2-D, 4-NODE VISCOPLASTIC ELEMENT
MP,EX,1,60.6E9
MP,NUXY,1,0.4999
TB,ANAND,1 ! VISCOPLASTIC MODEL BY ANAND
TBDATA,1,29.7E6
TBDATA,2,21.08999E3
TBDATA,3,1.91E7
TBDATA,4,7.0
TBDATA,5,0.23348
TBDATA,6,1115.6E6
TBDATA,7,18.92E6
TBDATA,8,0.07049
TBDATA,9,1.3
SAVE ! SAVE DATABASE FOR LATER USE
N,1
N,2,1E-2
N,3,1E-2,1E-2
N,4,,1E-2
E,1,2,3,4
D,1,ALL,,,2
FINISH
/OUTPUT,SCRATCH ! DIVERT VOLUMINOUS SOLUTION OUTPUT
/SOLU
SOLCONTROL,0
ANTYPE,STATIC
NLGEOM,ON ! LARGE DEFORMATION ACTIVATED
OUTRES,RSOL,ALL ! STORE REACTION RESULTS FOR ALL SUBSTEPS
BFUNIF,TEMP,673 ! UNIFORM TEMPERATURE OF 673 K
D,3,ALL,0.0,,4
TIME,0.000001 ! NEAR ZERO TIME FOR FIRST LOAD STEP
SOLVE
NSUBST,20
D,3,UX,0.2E-2,,4
TIME,20
SOLVE
FINISH
/OUTPUT
/POST26
RFORCE,2,3,F,X
RFORCE,3,4,F,X
ADD,4,2,3,,LOAD,,,(1/100),(1/100)
! THE FX FORCE IS DIVIDED BY 100 BECAUSE THE DEFAULT OUT-
! OF-PLANE THICKNESS IS 1 METER WHILE PROBLEM CONSIDERED
! HAS OUT-OF-PLANE THICKNESS OF 0.01 METER WHICH IS 100
! TIMES LESS THAN 1 METER.
PRVAR,4
*GET,F1,VARI,4,RTIME,20
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'Fx '
LABEL(1,2) = 'N'
*VFILL,VALUE(1,1),DATA,845
*VFILL,VALUE(1,2),DATA,F1
*VFILL,VALUE(1,3),DATA,ABS(F1/845)
SAVE,TABLE1
FINISH
/CLEAR, NOSTART
! VISCO108 ELEMENT TYPE
/PREP7
RESUME
ET,1,VISCO108 ! 2-D, 8-NODE VISCOPLASTIC ELEMENT
N,1
N,2,1E-2
N,3,1E-2,1E-2
N,4,,1E-2
N,5,0.5E-2
N,6,1E-2,0.5E-2
N,7,0.5E-2,1E-2
N,8,,0.5E-2
E,1,2,3,4,5,6,7,8
D,1,ALL
D,2,ALL
D,5,ALL
FINISH
/OUTPUT,SCRATCH ! DIVERT VOLUMINOUS SOLUTION OUTPUT
/SOLU
SOLCONTROL,0
NLGEOM,ON ! LARGE DEFORMATION ACTIVATED
OUTRES,RSOL,ALL ! STORE REACTION RESULTS FOR ALL SUBSTEPS
BFUNIF,TEMP,673 ! UNIFORM TEMPERATURE OF 673 K
D,3,ALL,0.0,,8
TIME,0.000001 ! NEAR ZERO TIME FOR FIRST LOAD STEP
SOLVE
NSUBST,20
D,3,UX,0.2E-2,,4
D,7,UX,0.2E-2
D,6,UX,0.1E-2,,8,2
TIME,20
SOLVE
FINISH
/OUTPUT
/POST26
RFORCE,2,3,F,X
RFORCE,3,4,F,X
RFORCE,4,7,F,X
ADD,5,2,3,4,LOAD,,,(1/100),(1/100),(1/100)
PRVAR,5
*GET,F1,VARI,5,RTIME,20
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'Fx '
LABEL(1,2) = 'N'
*VFILL,VALUE(1,1),DATA,845
*VFILL,VALUE(1,2),DATA,F1
*VFILL,VALUE(1,3),DATA,ABS(F1/845)
SAVE,TABLE2
FINISH
/CLEAR, NOSTART
! VISCO107 ELEMENT
/PREP7
RESUME
ET,1,VISCO107 !3-D, 8-NODE VISCOPLASTIC ELEMENT
N,1
N,2,1E-2
N,3,1E-2,1E-2
N,4,,1E-2
NGEN,2,4,1,4,1,,,1E-2
E,1,2,3,4,5,6,7,8
D,1,ALL,,,2
D,5,ALL,,,6
FINISH
/OUTPUT,SCRATCH ! DIVERT VOLUMINOUS SOLUTION OUTPUT
/SOLU
SOLCONTROL,0
NLGEOM,ON ! LARGE DEFORMATION ACTIVATED
OUTRES,RSOL,ALL ! STORE REACTION RESULTS FOR ALL SUBSTEPS
BFUNIF,TEMP,673 ! UNIFORM TEMPERATURE OF 673 K
D,3,ALL,0.0,,8
TIME,0.000001 ! NEAR ZERO TIME FOR FIRST LOAD STEP
SOLVE
NSUBST,20
D,3,UX,0.2E-2,,4
D,7,UX,0.2E-2,,8
TIME,20
SOLVE
FINISH
/OUTPUT
/POST26
RFORCE,2,3,F,X
RFORCE,3,4,F,X
RFORCE,4,7,F,X
RFORCE,5,8,F,X
ADD,6,2,3,4
ADD,7,6,5,,LOAD
PRVAR,7
*GET,F1,VARI,7,RTIME,20
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'Fx '
LABEL(1,2) = 'N'
*VFILL,VALUE(1,1),DATA,845
*VFILL,VALUE(1,2),DATA,F1
*VFILL,VALUE(1,3),DATA,ABS(F1/845)
/COM
/OUT,vm199,vrt
/COM,------------------- VM199 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,VISCO107 RESULTS:
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F8.3)
/NOPR
RESUME,TABLE1
/GOPR
/COM,
/COM,VISCO106 RESULTS:
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F8.3)
/NOPR
RESUME,TABLE2
/GOPR
/COM,
/COM,VISCO108 RESULTS:
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F8.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm199,vrt
VM200 (Viscoelastic Sandwich Seal Analysis) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM200
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
/TITLE, VM200, VISCOELASTIC SANDWICH SEAL ANALYSIS
/COM, ----- 2-D ANALYSIS -----
/COM, "FE CALCULATIONS OF RESIDUAL STRESSES ....",SOULES ET AL.
ET,1,VISCO88,,,1 ! AXISYMMETRIC 2-D VISCOELASTIC ELEMENT
*CREATE,MAC1 ! MACRO FOR MATERIAL PROPERTIES
/COM, MATERIAL ONE IS G-11 GLASS AND MATERIAL TWO IS ALUMINA
/COM, NOTE THAT ALUMINA IS AN ELASTIC MATERIAL THEREFORE IT
/COM, DOES NOT HAVE VISCOELASTICITY AND STRUCTURAL RELAXATION
/COM, MATERIAL PROPERTIES OF ALUMINA
TB,EVISC,2
TBDATA,31,52.6E-7 ! THERMAL EXPANSION COEFFICENTS FOR ALUMINA
TBDATA,32,.119E-7 ! LOCATIONS 31-33
TBDATA,33,-1.0E-11
TBDATA,46,1.435E5 ! SHEAR MODULUS AT TIME = 0
TBDATA,47,1.435E5 ! SHEAR MODULUS AT TIME = INFINITY
TBDATA,48,3.11E5 ! BULK MODULUS AT TIME = 0
TBDATA,49,3.11E5 ! BULK MODULUS AT TIME = INFINITY
/COM, MATERIAL PROPERTIES OF G-11 GLASS
TB,EVISC,1
TBDATA,1,6.45D4 ! H/R
TBDATA,2,0.53 ! VALUE OF X
TBDATA,3,6 ! NO. OF MAXWELL ELEMENTS FOR STRUCTURAL RELAXATION
TBDATA,6,0.108 ! COEFFICIENTS OF THE MAXWELL ELEMENTS FOR
TBDATA,7,0.443 ! VOLUME DECAY (STRUCTURAL RELAXATION)
TBDATA,8,0.166 ! LOCATIONS 6-11
TBDATA,9,0.161
TBDATA,10,0.046
TBDATA,11,0.076
TBDATA,16,3.00 ! RELAXATION TIMES FOR VOLUME DECAY FUNCTION
TBDATA,17,0.671 ! LOCATIONS 16-21
TBDATA,18,0.247
TBDATA,19,0.091
TBDATA,20,0.033
TBDATA,21,0.008
TBDATA,26,3.43E-5 ! THERMAL EXPANSION COEFF. IN LIQUID STATE
TBDATA,31,64.7E-7 ! THERMAL EXPANSION COEFF. IN GLASS (SOLID) STATE
TBDATA,32,.02E-7 ! LOCATIONS 31-32
TBDATA,36,618 ! FICTIVE TEMPERATURE LOCATIONS 36-41
TBDATA,37,618
TBDATA,38,618
TBDATA,39,618
TBDATA,40,618
TBDATA,41,618
TBDATA,46,2.79D4 ! SHEAR MODULUS AT TIME = 0
TBDATA,47,0.0D0 ! SHEAR MODULUS AT TIME = INFINITY
TBDATA,48,6.05D4 ! BULK MODULUS AT TIME = 0
TBDATA,49,6.05D4 ! BULK MODULUS AT TIME = INFINITY
TBDATA,50,3 ! THREE MAXWELL ELEMENTS FOR SHEAR RELAXATION
TBDATA,51,0.422 ! COEFF. FOR SHEAR RELAXATION LOCATIONS 51-53
TBDATA,52,0.423
TBDATA,53,0.155
TBDATA,61,0.0689 ! RELAXATION TIMES FOR SHEAR RELAXATION
TBDATA,62,0.0065 ! LOCATIONS 61-63
TBDATA,63,0.0001
TBDATA,71,0 ! NO BULK MODULUS RELAXATION
*END
*USE,MAC1 ! EXECUTE MACRO FOR MATERIAL PROPERTIES
:COM, CREATE FINITE ELEMENT MODEL
N,1,
N,3,,0.00025
FILL
N,5,0,(0.00025+0.00325)
FILL
NGEN,3,10,1,5,1,.001
MAT,2
E,1,21,23,3,11,22,13,2
MAT,1
E,3,23,25,5,13,24,15,4
/COM, APPLY BOUNDARY CONDITIONS AND COUPLING
NSEL,S,LOC,Y
DSYM,SYMM,Y
NSEL,S,LOC,X
DSYM,SYMM,X
NSEL,ALL
D,1,ALL
CP,1,UX,21,22,23,24,25
CP,2,UY,2,22
CP,3,UY,3,13,23
CPSGEN,2,2,2,3,1
FINISH
/COM SINCE THE SOLUTION OUTPUT IS VOLUMINOUS IT IS DIVERTED TO A
/COM SCRATCH FILE
/OUTPUT,SCRATCH
*CREATE,MAC2 ! CREATE MACRO FOR ANALYSIS TYPE AND LOADING
/SOLU
SOLCONTROL,0
ANTYPE,STATIC
/COM, TEMPERATURE SET UP
TREF,618
TOFFST,273
TUNIF,618
TIME,1E-5
CNVTOL,F,,,,.00001 ! VERY SMALL MINIMUM ENFORCED
! FOR CONVERGENCE
SOLVE
OUTRES,ESOL,1 ! STORE RESULTS FOR EVERY SUBSTEP
NSUBST,200
TUNIF,460 ! COOLING
TIME,3160
SOLVE
TIME,(14400+3160) ! ISOTHERMAL HOLD
SOLVE
TUNIF,18 ! FURTHER COOLING
TIME,(14400+12000)
SOLVE
*END
FINISH
*USE,MAC2 ! EXECUTE ANALYSIS AND LOADING MACRO
/OUTPUT
/POST26
ESOL,2,2,,BFE,TEMP
ESOL,3,2,3,S,X,STRESS
*CREATE,MAC3 ! MACRO FOR PROCESSING RESULTS
XVAR,2
/GRID,1
/AXLAB,X,TEMPERATURE
/AXLAB,Y,IN-PLANE STRESS (MPA)
PLVAR,3
*GET,MXSX,VARI,3,EXTREM,VMAX ! MAXIMUM IN-PLANE STRESS
NSTORE,20 ! STORE EVERY 20TH TIME POINT RESULTS
PRVAR,2,3
*END
*USE,MAC3 ! EXECUTE POSTPROCESSING MACRO
*SET,P1,(MXSX)
*GET,T1,VARI,3,EXTREM,TMAX
*GET,TE,VARI,2,RTIME,T1
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'PRES MX ','TEMP '
LABEL(1,2) = 'MPa','DEG C'
*VFILL,VALUE(1,1),DATA,12.5,460
*VFILL,VALUE(1,2),DATA,P1,TE
*VFILL,VALUE(1,3),DATA,ABS(P1/12.5),ABS(TE/460)
SAVE,TABLE_1
FINISH
/CLEAR,NOSTART ! CLEAR THE DATABASE
/PREP7
/TITLE, VM200, VISCOELASTIC SANDWICH SEAL ANALYSIS
/COM, ----- 3-D ANALYSIS -----
ET,1,89 ! 3-D VISCOELASTIC ELEMENT
*USE,MAC1 ! EXECUTE MACRO FOR MATERIAL PROPERTIES
/COM, CREATE FINITE ELEMENT MODEL
N,1,
N,3,0.00025
FILL
N,5,(0.00025+0.00325)
FILL
NGEN,3,10,1,5,1,,.001
NGEN,3,100,1,25,1,,,0.001
MAT,2
E,1,3,23,21,201,203,223,221
EMORE,2,13,22,11,202,213,222,211
EMORE,101,103,123,121
EGEN,2,2,1,1,1,-1
NSLE,S
NSEL,INVE
NDELE,ALL
NSLE,S
/COM, APPLY BOUNDARY CONDITIONS AND COUPLING
NSEL,S,LOC,Y
DSYM,SYMM,Y
NSEL,S,LOC,X
DSYM,SYMM,X
NSEL,S,LOC,Z
DSYM,SYMM,Z
NSEL,S,LOC,Y,0.002
CP,1,UY,ALL
NSEL,S,LOC,Z,0.002
CP,2,UZ,ALL
NSEL,S,LOC,X,0.00025
CP,3,UX,ALL
NSEL,S,LOC,X,0.0035
CP,4,UX,ALL
NSEL,ALL
FINISH
/COM, SINCE THE SOLUTION OUTPUT IS VOLUMINOUS IT IS DIVERTED TO A
/COM, SCRATCH FILE
/OUTPUT,SCRATCH
*USE,MAC2 ! EXECUTE ANALYSIS AND LOADING MACRO
/OUTPUT
/POST26
ESOL,2,2,,BFE,TEMP
ESOL,3,2,3,S,Y,STRESS
*USE,MAC3 ! EXECUTE POSTPROCESSING MACRO
*SET,P2,(MXSX)
*GET,T2,VARI,3,EXTREM,TMAX
*GET,TE2,VARI,2,RTIME,T2
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'PRES MX ','TEMP '
LABEL(1,2) = 'MPa','DEG C'
*VFILL,VALUE(1,1),DATA,12.5,460
*VFILL,VALUE(1,2),DATA,P2,TE2
*VFILL,VALUE(1,3),DATA,ABS(P2/12.5),ABS(TE2/460)
SAVE,TABLE_2
RESUME,TABLE_1
/COM
/OUT,vm200,vrt
/COM,------------------- VM200 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS USING VISCO88
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS USING VISCO89
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.1,' ',F10.1,' ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
/NOPR
/DELETE,TABLE_1
/DELETE,TABLE_2
/DELETE,MAC1
/DELETE,MAC2
/DELETE,MAC3
/DELETE,SCRATCH
FINISH
*LIST,vm200,vrt
VM201 (Rubber Cylinder Pressed Between Two Plates) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM201
/config,nproc,4
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
R = 0.2 ! RADIUS OF CYLINDER
/PREP7
smrt,off
/TITLE, VM201, RUBBER CYLINDER PRESSED BETWEEN TWO PLATES
/COM REF: T. SUSSMAN, K.J. BATHE, "A FE FORMULATION FOR NONLINEAR ..."
/COM COMPUTERS & STRUCTURES, VOL. 26, NOS. 1/2, 1987
ET,1,HYPER56 ! 2-D 4-NODE MIXED U-P HYPERELASTIC SOLID
ET,2,CONTAC26 ! 2-D POINT-TO-GROUND CONTACT ELEMENT
R,2,2000 ! SET SURFACE STIFFNESS
MP,EX,1,2.82 ! YOUNG'S MODULUS [MPA]
MP,NUXY,1,0.49967 ! POISSON'S RATIO
TB,MOONEY,1 ! MOONEY-RIVLIN CONSTANTS [MPA]
TBDATA,1,0.293
TBDATA,2,0.177
CSYS,1 ! SWITCH TO CYLINDRICAL C.S.
K,1 ! DEFINE KEYPOINTS
K,2,R,-90
K,3,R
K,4,(0.5*R),-90
K,5,(0.6*R),-45
K,6,(0.5*R)
K,7,R,-45
L,2,7
L,7,3
CSYS,0 ! SWITCH TO CARTESIAN C.S.
A,2,7,5,4
A,7,3,6,5
A,4,5,6,1
ESIZE,,4 ! SET ELEMENT DIVISION SIZE
AMESH,ALL ! MESH ALL AREAS
SAVE ! SAVE MODEL FOR SECOND ANALYSIS
*CREATE,CONTAC2D,MAC ! MACRO TO GENERATE CONTACT ELEMENTS
N,1001,(-2*R),-R ! TARGET SURFACE NODES
N,1002,(2*R) ,-R
CSYS,1 ! SWITCH TO CYLINDRICAL C.S.
TYPE,2
REAL,2
*DO,J,1,ARG1
NODC = NODE(R,90*((J-1)/ARG1-1),0)
E,NODC,1001,1002
*ENDDO
CSYS,0 ! SWITCH TO CARTESIAN C.S.
NSEL,S,LOC,X ! SELECT LEFT EDGE
D,ALL,UX ! CONSTRAIN LEFT EDGE IN UX
NSEL,S,LOC,Y ! SELECT TOP EDGE
CP,1,UY,ALL ! COUPLE TOP EDGE IN UY
*GET,NCEN,NODE,,NUM,MIN ! GET MINIMUM NODE NUMBER FROM SELECTED SET
NSEL,ALL
*END
CONTAC2D,8 ! USE MACRO CONTAC2D
FINISH
*CREATE,SOLVIT,MAC ! MACRO TO SOLVE MODEL
/SOLU
SOLCONTROL,0
ANTYPE,STATIC
CNVTOL,F,,,,-1
NLGEOM,ON ! INCLUDE LARGE DEFORMATION EFFECTS
NSUBST,6 ! SPECIFY NUMBER OF SUBSTEPS IN LOAD STEP
OUTRES,,1 ! WRITE SOLUTION FOR EVERY SUBSTEP
D,NCEN,UY,-0.1 ! APPLY DISPLACEMENT UY = -0.1 TO COUPLED NODES
SOLVE
FINISH
*END
SOLVIT ! USE MACRO SOLVE2D
*CREATE,PLOTS,MAC ! MACRO FOR POST-PROCESSING
/POST1
/DSCALE,1,1
PLDISP,1 ! PLOT DISPLACED SHAPE
FINISH
/POST26
/AXLAB,Y,FORCE
/AXLAB,X,DISPLACEMENT
NSOL,2,NCEN,U,Y
RFORCE,3,NCEN,F,Y
PROD,2,2,,,,,,-2
PROD,3,3,,,,,,-2
XVAR,2
PLVAR,3 ! PLOT DISPLACEMENT VS FORCE
PRVAR,2,3 ! PRINT DISPLACEMENT, FORCE
*GET,F1,VARI,3,RTIME,.5
*GET,F2,VARI,3,RTIME,1
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'F @ ','F @ '
LABEL(1,2) = '.1','.2'
*VFILL,VALUE(1,1),DATA,.25,1.4
*VFILL,VALUE(1,2),DATA,F1,F2
*VFILL,VALUE(1,3),DATA,ABS(F1/.25),ABS(F2/1.4)
FINISH
*END
PLOTS ! USE MACRO PLOTS
SAVE,TABLE_1
RESUME
/PREP7
smrt,off
ESEL,S,TYPE,,1
ET,3,HYPER74 ! 2-D 8-NODE MIXED U-P HYPERELASTIC SOLID
EMODIF,ALL,TYPE,3 ! SWITCH ELEMENTS TO HIGHER ORDER
EMID ! ADD MIDSIDE NODES
ESEL,ALL
CONTAC2D,16 ! USE MACRO TO GENERATE CONTACT ELEMENTS
FINISH
SOLVIT ! USE MACRO TO OBTAIN SOLUTION
PLOTS ! USE MACRO TO POSTPROCESS
SAVE,TABLE_2
RESUME
/PREP7
smrt,off
ET,5,HYPER58 ! 3-D 8-NODE MIXED U-P HYPERELASTIC SOLID
ET,6,CONTAC49 ! 3-D POINT-TO-SURFACE CONTACT ELEMENT
ET,7,MESH200,6 ! 2-D 4-NODED QUAD
R,6,2000,,,,0.1
R,7,0.05
ALLSEL
TYPE,5
ESIZE,,1
VEXT,ALL,,,,,1
N,1001,,-R ! CREATE TARGET PLANE OF NODES
N,1002,2*R,-R
N,1003,2*R,-R,8*R
N,1004,,-R,8*R
TYPE,7
REAL,7
E,1002,1001,1004,1003
NSEL,S,NODE,,1001,1004
CM,TARGET,NODE ! DEFINE TARGET COMPONENT NODES
D,ALL,ALL,0
CSYS,1 ! SWITCH TO CYLINDRICAL C.S.
ESEL,S,TYPE,,5
NSLE
NSEL,R,LOC,X,R
CM,CONTAC,NODE ! DEFINE CONTACT COMPONENT NODES
ESEL,S,TYPE,,5,7
TYPE,6
REAL,6
GCGEN,CONTAC,TARGET ! GENERATE 3-D CONTACT ELEMENTS
CSYS,0 ! SWITCH TO CARTESIAN C.S.
NSEL,ALL
D,ALL,UZ ! CONSTRAIN ALL NODES IN Z (PLANE STRAIN)
NSEL,S,LOC,X
D,ALL,UX
NSEL,S,LOC,Y
CP,1,UY,ALL ! COUPLE TOP NODES IN Y
*GET,NCEN,NODE,,NUM,MIN
ESEL,S,TYPE,,5,6
NSLE
FINISH
SOLVIT ! USE MACRO TO OBTAIN SOLUTION
PLOTS ! USE MACRO TO POSTPROCESS
SAVE,TABLE_3
RESUME,TABLE_1
/COM
/OUT,vm201,vrt
/COM,------------------- VM201 RESULTS COMPARISON --------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,RESULTS USING HYPER56:
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/NOPR
RESUME,TABLE_2
/GOPR
/COM,
/COM,RESULTS USING HYPER74:
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/NOPR
RESUME,TABLE_3
/GOPR
/COM,
/COM,RESULTS USING HYPER78:
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.2,' ',F10.2,' ',1F5.3)
/COM,-----------------------------------------------------------
/COM,
/OUT
FINISH
/NOPR
/DELETE,CONTAC2D,MAC
/DELETE,PLOTS,MAC
/DELETE,SOLVIT,MAC
/DELETE,TABLE_1
/DELETE,TABLE_2
/DELETE,TABLE_3
*LIST,vm201,vrt
VM202 (Transverse Vibrations of a Shear Beam) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM202
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
/TITLE, VM202, TRANSVERSE VIBRATIONS OF A SHEAR BEAM
/COM, REF: R. BLEVINS, FORMULAS FOR NATURAL FREQUENCY AND MODE SHAPE, 1979.
ANTYPE,MODAL
ET,1,SHELL28
R,1,0.1 ! DEFINE REAL CONSTANT SET 1; THICKNESS OF SHELL
MP,EX, 1,200E9
MP,NUXY,1,0.27
MP,DENS,1,7860.0
N,1 ! DEFINE NODES
N,2,10
N,3,10,10
N,4, ,10
NGEN,9,4,1,4,1,,,3.75
E,1,2,6,5 ! DEFINE ELEMENTS
EGEN,3,1,-1
E,4,1,5,8
EGEN,8,4,-4
D,ALL,UZ ! DEFINE BOUNDARY CONDITIONS
D,1,UX,,,4,1,UY
NSEL,S,LOC,Z,3.75 ! SELECT NODES AT COMMON Z LOCATION
CP,1,UX,ALL ! COUPLE NODES IN X AND Y DOF
CP,2,UY,ALL
NSEL,S,LOC,Z,7.5
CP,3,UX,ALL
CP,4,UY,ALL
NSEL,S,LOC,Z,11.25
CP,5,UX,ALL
CP,6,UY,ALL
NSEL,S,LOC,Z,15
CP,7,UX,ALL
CP,8,UY,ALL
NSEL,S,LOC,Z,18.75
CP,9,UX,ALL
CP,10,UY,ALL
NSEL,S,LOC,Z,22.5
CP,11,UX,ALL
CP,12,UY,ALL
NSEL,S,LOC,Z,26.25
CP,13,UX,ALL
CP,14,UY,ALL
NSEL,S,LOC,Z,30
CP,15,UX,ALL
CP,16,UY,ALL
NSEL,ALL
FINISH
/SOLU
MXPAND,4 ! EXPAND FIRST TWO REPEATED MODES
MODOPT,REDUC
TOTAL,8 ! SET AUTOMATIC MASTER DOF GENERATION TO 8
SOLVE
*GET,F1,MODE,1,FREQ
*GET,F2,MODE,3,FREQ
FINISH
/POST1
SET,1,1
/VUP,,Z
/VIEW,,,1
PLDISP,1 ! DISPLAY DISPLACEMENTS (FIRST MODE)
SET,1,3
PLDISP,1 ! DISPLAY DISPLACEMENTS (SECOND MODE)
*STATUS
*DIM,VALUE,,2,3
*DIM,LABEL,CHAR,2,2
LABEL(1,1) = 'F1 ','F2 '
LABEL(1,2) = 'Hz ','Hz '
*VFILL,VALUE(1,1),DATA,17.375,52.176
*VFILL,VALUE(1,2),DATA,F1,F2
*VFILL,VALUE(1,3),DATA,ABS(F1/17.375),ABS(F2/52.176)
/COM
/OUT,vm202,vrt
/COM,------------------- VM202 RESULTS COMPARISON -------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,' ',F10.3,' ',F10.3,' ',1F8.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm202,vrt
N1=NODE(S,0,0) ! GET A NODE ON THE COIL
N2=NODE(0,0,0) ! GET NODE AT ORIGIN
ESEL,S,MAT,,2 ! GET COIL ELEMENTS
NSLE,S
CP,1,CURR,ALL ! COUPLE CURR DOF IN COIL
*GET,ELM,ELEM,,NUM,MIN ! GET AN ELEMENT NUMBER IN THE COIL REGION
NSEL,ALL
ESEL,ALL
CSYS,1
NSEL,S,LOC,X,12*S
SF,ALL,INF
NSEL,S,LOC,X,0
D,ALL,AZ,0
NSEL,ALL
FINISH
/POST1
SET,1,1 ! READ IN REAL SOLUTION
PRESOL,ELEM ! PRINT CIRCUIT SOLUTION PER ELEMENT
SET,1,1,,1 ! READ IN IMAGINARY SOLUTION
PRESOL,ELEM ! PRINT CIRCUIT SOLUTION PER ELEMENT
FINISH
/POST26
NSOL,2,4,VOLT
PRVAR,2
*GET,REAL,VARI,2,RTIME,15915
*GET,IMAG,VARI,2,ITIME,15915
FINISH
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'REAL VOL','IMAG VOL'
LABEL(1,2) = 'TAGE, V ','TAGE, V '
*VFILL,VALUE(1,1),DATA,16.44,-1.41
*VFILL,VALUE(1,2),DATA,REAL,IMAG
*VFILL,VALUE(1,3),DATA,ABS(REAL/16.44),ABS(IMAG/1.41)
/COM
/OUT,vm208,vrt
/COM,------------------- VM208 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.5,' ',F10.5,' ',1F5.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm208,vrt
VM209 (Multiple Species Flow Entering a Circular Pipe) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM209
/PREP7
/TITLE,VM209, ENTRANCE FLOW OF SPECIES THROUGH A CIRCULAR PIPE
/NOPR
/COM
/COM -----------------------------------------------------------------------
/COM
/COM TEST THE SPECIES TRANSPORT EQUATION.
/COM
/COM FLOW OF TWO SPECIES THROUGH A CIRCULAR PIPE IS CONSIDERED.
/COM THE FLOW IS HYDRODYNAMICALLY DEVELOPED. THIS IS ACHIEVED BY SPECIFYING
/COM THE PRESURE AT THE INLET AND OUTLET. UNIFORM BUT DISSIMILIAR
/COM MASS FRACTIONS ARE SPECIFIED AT THE INLET AND WALL. SYMMETRY BOUNDARY
/COM CONDITIONS ARE APPLIED AT THE PIPE CENTERLINE. BOTH SPECIES HAVE THE
/COM SAME DENSITY, VISCOSITY, MASS DIFFUSION COEFFICIENT. RESULTING MASS
/COM FRACTION DISTRUBUTIONS ARE ANALOGOUS TO THAT OF TEMPERATURE
/COM IN GREATZ PROBLEM ( HYDRODYNAMICALLY FULLY DEVELOPED, THERMALLY
/COM DEVELOPING FLOW IN A PIPE WITH UNIFORM WALL TEMPERATURE.) THE RESULTS
/COM ARE COMPARED WITH ANALYTICAL SOLUTION ( REF: CONVECTIVE HEAT AND MASS
/COM TRANSFER, THIRD ED., KAYS, W.M., CRAWFORD,M.E,MCGRAW-HILL,1993.)
/COM
/COM ------------------------------------------------------------------------
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
ET,1,FLUID141,,,2 ! 2D AXISYMMETRIC XR SYSTEM
MSHK,1
MSHA,0,2D
PI = ACOS( -1.0 )
L = 0.1 ! PIPE LENGHT (m)
R = 0.0025 ! PIPE RADIUS (m)
NRD = 12 ! NO. OF DIVISIONS IN THE RADIAL DIRECTION
RS = 4 ! SPACING RATIO IN THE RADIAL DIRECTION
NXD = 100 ! NO. OF DIVISIONS IN THE AXIAL DIRECTION
XS = 1 ! SPACING RATIO IN THE AXIAL DIRECTION
PIN = 1.28 ! INLET PRESSURE
YBIN = 1.0 ! MASS FRACTION OF SPECIES B AT THE INLET
YBWALL = 0.0 ! MASS FRACTION OF SPECIES B AT THE WALL
RHO = 1.0 ! DENSITY (kg/m**3)
MU = 1.0E-5 ! VISCOSITY (Pa*s)
DIFF = 1.43E-5 ! DIFFUSION COEFFICIENT (m**2/s)
!
! DEFINE GEOMETRY AND GENERATE MESH
!
RECT,,L,,R
LESIZE,2,,,NRD,-RS
LESIZE,4,,,NRD,-RS
LESIZE,1,,,NXD,XS
LESIZE,3,,,NXD,1/XS
AMESH,1
KEYOPT,1,1,2 ! SET THE NUMBER OF SPECIES TO TW0
!
! APPLY THE BOUNDARY CONDITIONS
!
LSEL,S,,,4 ! INLET BC'S
NSLL,S,1
D,ALL,PRES,PIN ! PRESSURE
D,ALL,SP01,1.0-YBIN ! MASS FRACTION OF A
D,ALL,SP02,YBIN ! MASS FRACTION OF B
ALLS
LSEL,S,,,2 ! OUTLET BC'S
NSLL,S,1
D,ALL,PRES,0.0 ! PRESSURE
ALLS
LSEL,S,,,3 ! WALL BC'S
NSLL,S,1
D,ALL,VX,0.0 ! AXIAL VELOCITY IS SET TO ZERO
D,ALL,VY,0.0 ! RADIAL VELOCITY IS SET TO ZERO
D,ALL,SP01,1.0-YBWALL ! MASS FRACTION OF A
D,ALL,SP02,YBWALL ! MASS FRACTION OF B
ALLS
LSEL,S,,,1 ! SYMMETRY BC'S
NSLL,S,1
D,ALL,VY,0.0 ! RADIAL VELOCITY IS SET TO ZERO
ALLS
/SOLU
!
! SOLUTION OPTIONS
!
FLDA,ITER,EXEC,700 ! NO. OF GLOBAL ITERATIONS
FLDATA,ITER,CHEC,50 ! CHECKPOINT FREQUENCY
FLDA,PROT,VISC,CONSTANT
FLDA,PROT,DENS,CONSTANT
FLDA,NOMI,DENS,RHO ! FLUID DENSITY
FLDA,NOMI,VISC,MU ! FLUID VISCOSITY
!
FLDA,SOLU,SPEC,T ! ACTIVATE SPECIES TRANSPORT EQUATIONS
MSDATA,1 ! DESIGNATE ALGEBRAIC SPECIES
MSSPEC,1,SP01 ! ESTABLISH SPECIES 1
MSSPEC,2,SP02 ! ESTABLISH SPECIES 2
MSMETH,2,1 ! SPECIFY THE METHOD OF SOLUTION
MSNOMF,1,(1.0-YBWALL) ! INITIAL CONDITION FOR SPECIES 1
MSNOMF,2,YBWALL ! INITIAL CONDITION FOR SPECIES 2
MSRELAX,2,1.0 ! NO RELAXATION FOR SPECIES
FLDATA,PROT,VISC,CMIX ! RE-SPECIFY VISCOSITY AS OF MIXTURE
FLDATA,PROT,DENS,CMIX ! RE-SPECIFY DENSITY AS OF MIXTURE
FLDATA,VARY,VISC,T ! ALLOW MIXTURE VISCOSITY TO VARY
FLDATA,VARY,DENS,T ! ALLOW MIXTURE DENSITY TO VARY
MSPROP,1,VISC,CONSTANT,MU ! VISCOSITY OF SP01
MSPROP,2,VISC,CONSTANT,MU ! VISCOSITY OF SP02
MSPROP,1,DENS,CONSTANT,RHO ! DENSITY OF SP01
MSPROP,2,DENS,CONSTANT,RHO ! DENSITY OF SP02
MSPROP,1,MDIF,CONSTANT,DIFF ! MASS DIFFUSION COEFFICIENT OF SP01
MSPROP,2,MDIF,CONSTANT,DIFF ! MASS DIFFUSION COEFFICIENT OF SP02
!
SOLVE
FINISH
SAVE
!
! POST PROCESSING
!
/POST1
SET,LAST
/TRIAD,OFF ! DO NOT SHOW TRIAD
/RATIO,1,1,10 ! SCALE THE R_COORDINATE FOR CLEARITY
EPLOT ! PLOT THE MESH
/EDGE,1,1 ! SHOW THE EDGES ONLY
/CONTOUR,1,21 ! SET THE NUMBER OF CONTOUR LEVELS
/TITLE,CONTOURS OF AXIAL VELOCITY
PLNSOL,VX ! PLOT CONTOURS OF AXIAL VELOCITY
/TITLE,CONTOURS OF MASS FRACTION OF SPECIES A
PLNSOL,SP01 ! PLOT CONTOURS OF MASS FRACTION (A)
/TITLE,CONTOURS OF MASS FRACTION OF SPECIES B
PLNSOL,SP02 ! PLOT CONTOURS OF MASS FRACTION (B)
VAVE = 0.125*PIN*R**2/(L*MU) ! CALC. THE AVERAGE VELOCITY FROM INPUT
DATA
RE = 2*RHO*R*VAVE/MU ! CALC. THE REYNOLDS NUMBER
SC = (MU/RHO)/DIFF ! CALC. THE SCHMIDT NUMBER
*DIM,XCORD,,10 ! DECLARE AN ARRAY TO STORE AXIAL COORDINATES
*DIM,AVEYB,,10,2 ! DECLARE AN ARRAY TO STORE AVERAGE MASS
*DIM,RAT,,10
! FRACTIONS OF SPECIES B AT VARIOUS AXIAL
! LOCATIONS ( BOTH CALCULATED AND ANALYTICAL)
/NOPR
*DIM,pth,char,10
pth(1)='pth1'
pth(2)='pth2'
pth(3)='pth3'
pth(4)='pth4'
pth(5)='pth5'
pth(6)='pth6'
pth(7)='pth7'
pth(8)='pth8'
pth(9)='pth9'
pth(10)='pth10'
/GOPR
*DO,I,1,10
XCORD(I) = L*I/10 ! DETERMINE THE AXIAL LOCATION
N1 = NODE(XCORD(I),0,0)
N2 = NODE(XCORD(I),R,0)
PATH,pth(i),2,,48 ! DEFINE THE PATH
PPATH,1,N1 ! DEFINE THE PATH POINT
PPATH,2,N2 ! DEFINE THE PATH POINT
PDEF,SP02,SP02
PCALC,INTG,YBC,SP02,S ! INTEGRATE ALONG THE PATH TO CALC. TOTAL
! MASS FRACTION
*GET,YBC,PATH,,LAST,YBC
AVEYB(I,1) = YBC/R ! CALCULATE AVERAGE MASS FRACTION OF B
ARG = XCORD(I)/(R*RE*SC)
TERM0 = 0.1024204*EXP(-7.313*ARG)
TERM1 = 0.0121946*EXP(-44.61*ARG)
TERM2 = 0.0040650*EXP(-113.9*ARG)
TERM3 = 0.0019284*EXP(-215.2*ARG)
TERM4 = 0.0010987*EXP(-348.6*ARG)
TERM = TERM0 + TERM1 + TERM2 + TERM3 + TERM4
! CALCULATE THE AVERAGE MASS FRACTION
! ANALYTICALLY
AVEYB(I,2) = YBWALL - 8.0*(YBWALL-YBIN)*TERM
RAT(I) = AVEYB(I,1)/AVEYB(I,2)
*ENDDO
*vlen,10
/COM
/OUT,vm209,vrt
/COM,------------------- VM209 RESULTS COMPARISON ---------------
/COM,
/COM,Axial location(m)| TARGET | ANSYS | RATIO
/COM,
*VWRITE,XCORD(1),AVEYB(1,2),AVEYB(1,1),RAT(1)
(4(5x,f8.4))
/COM,
/COM,-------------------------------------------------------------
/OUT
FINISH
*LIST,vm209,vrt
VM210 (Pyramid Validation of Tetrahedron to Hexahedron Interface) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM210
/TITLE,VM210, BENDING OF HEX-TO-TET INTERFACE, FORMATION OF PYRAMIDS
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
SMRT,OFF
ET,1,95 !ELEMENT TYPE SOLID95
ET,2,95 !ELEMENT TYPE SOLID95
MP,EX,1,30E6 !ELASTIC MODULUS
MP,NUXY,1,0.3 !POISSON RATIO
P = 200 !FORCE
W = 31.071 !WIDTH
H = 33.917 !HEIGHT
L = 37.264 !LENGTH
I = 1/12*(W)*(H**3) !MOMENT OF INERTIA
SURF = P/(W*H) !SURFACE FORCE
BLOCK,-W/2, 0 , 0 , H/2, 0 , L
BLOCK, 0 , W/2, 0 , H/2, 0 , L !CREATE BLOCKS
BLOCK, 0 , W/2,-H/2, 0 , 0 , L !FORMING MIDPLANES
BLOCK,-W/2, 0 ,-H/2, 0 , 0 , L
WPOFF,,,L/2 !OFFSET WORKPLANE
WPLANE,-1,0,0,18.632,.87653,.402592,18.896,-.437455,.894952,18.72
WPLANE,-1,0,0,18.632,.890222,.338091,18.937,-.405661,.893273,18.826
VSBW,ALL,,
VGLUE,ALL
NUMCMP,ALL
/VIEW, 1, 0.9227 , 0.3132 ,-0.2246
/ANG, 1, 4.473
/PNUM,LINE,1
/NUM,-1
LPLO
/PNUM,LINE,0
/NUM,0
LESIZE,ALL,,,2
MOPT,PYRA,ON !MESH SET TO TRANSITIONAL PYRAMID
MSHK,1
MSHA,0
MSHM,0
TYPE,1
VSEL,S,VOLU,,1,2,1
VSEL,A,VOLU,,5,7,2
VMESH,ALL !VMESH ELEMENT TYPE 1
MSHK,0
MSHA,1,3D
MSHM,0
TYPE,2
VSEL,S,VOLU,,3,4,1
VSEL,A,VOLU,,6,8,2
VMESH,ALL !VMESH ELEMENT TYPE 2
VSEL,ALL
CSYS,4 ! USE WP AS LOCAL COOR. SYS
NSEL,S,LOC,Z,0 ! SELECT ALL NODES ON WP
CSYS,0 ! REVERT COOR. SYS TO CART.
ESLN ! ALL ELEMENTS ATTACHED TO WP NODES
NSLE ! AND ALL NODES ATT. TO ELEMS
/SHRINK,0.5
/VIEW,1, 0.51440 , -0.35450 , -0.78090
/ANG,1 ,1.41
EPLO ! ADJUST VIEW AND CAPTURE INTERFACE
/SHRINK,0
!
! THIS IS THE SECTION CONTAINING THE LOADING. THE MOMENT
! USED IN THE VERIFICATION EQUATION IS THE SUM OF ALL
! OF THE MOMENTS ON THE AREA LOCATED AT Z=0 OR Z=L. NOTE
! THAT THE LOADS ON EITHER AREA FORM A COUPLE.
!
SFGRAD,PRES,,Y,H/2,2*SURF/H !APPLY GRADIENT SURFACE LOADS
ASEL,S,LOC,Z,0
ASEL,A,LOC,Z,L
SFA,ALL,1,PRES,SURF,
ALLS
SFTRAN
! LOADS END
NSEL,S,LOC,X,W/2 !DEFINE DOF CONSTRAINTS
NSEL,R,LOC,Y,0 !SIMILAR TO PATCH TEST
NSEL,R,LOC,Z,0 !CONSTRAINTS
D,ALL, UX ,
D,ALL, UY ,
D,ALL, UZ ,
NSEL,S,LOC,X,-W/2
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,0
D,ALL, UY ,
D,ALL, UZ ,
NSEL,S,LOC,X,W/2
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,L
D,ALL, UY ,
! !CONSTRAINTS END
ALLSEL
! !CALCULATION : TOP AVERAGE OF STRESS
!
NSEL,S,LOC,Y,H/2 !SELECT NODES IN TOP AREA
*GET,NUMNOD1,NODE,,COUNT !OBTAIN NUMBER OF TOP SURFACE NODES
ZSTR1=0
TOTAL=0
COUNT=0
*DO,J,1,NUMNOD1,1
NDE1=NDNEXT(NDE1)
NSEL,,,,NDE1
ESLN
*GET,NUMELM,ELEM,0,COUNT !GET NUMBER OF ELEMS CONNECTED TO NDE1
ELNUM=0
TRIP=0
*DO,K,1,NUMELM,1 !LOOP ON ELEMS CONNECTED TO NDE1
ELNUM=ELNEXT(ELNUM)
*DO,L,1,8,1 !VOLUME : SOLID95 CORNER NODES
POS=NELEM(ELNUM,L) !CHECK POS 1-8 ON SOLID95 FOR
! NODE NUMBER
*IF,POS,EQ,NDE1,THEN
TRIP=1 !SET TRIP IF OUR CURRENTLY SELECTED
!NODE IS CORNER OF ELEMENT ELNUM
*ENDIF
*IF,TRIP,EQ,1,EXIT
*ENDDO
*IF,TRIP,EQ,1,EXIT
*ENDDO
*IF,TRIP,NE,1,THEN
NSEL,S,LOC,Y,H/2 !IF ENTRY IS GAINED TO HERE
ESEL,ALL !THEN IT MUST BE THAT NDE1 IS A
*CYCLE !MIDNODE : NO STRESS CALC WANTED!
*ENDIF
ALLSEL
*GET,ZSTR1,NODE,NDE1,S,Z !IF ENTRY GAINED TO THIS POINT OF
TOTAL=TOTAL+ZSTR1 !THE LOOP, THEN GRAB Z-STRESS AT
COUNT=COUNT+1 !NODE NDE1 AND ADD IT TO TOTAL.
NSEL,S,LOC,Y,H/2 !RESET SELECTED SETS TO WHAT IS
ESEL,ALL !NEEDED.
*ENDDO
MST1=(TOTAL/COUNT) !AVERAGE OF Z-STRESS ON TOP SURFACE
RESULTS(1)=MST1
NSEL,S,LOC,Y,0 !SELECT NODES ALONG Z AXIS
NSEL,R,LOC,X,0 ! ( THE NURTRAL AXIS )
*GET,NUMNOD2,NODE,,COUNT !OBTAIN NUMBER OF TOP SURFACE NODES
VM211 (Rubber Cylinder Pressed Between Two Plates) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
!! CPU 500
/VERIFY,VM211
/CONFIG,NPROC,4
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
R = 0.2 ! RADIUS OF CYLINDER
/PREP7
SMRT,OFF
/TITLE, VM211, RUBBER CYLINDER PRESSED BETWEEN TWO PLATES
/COM REF: T. SUSSMAN, K.J. BATHE, "A FE FORMULATION FOR NONLINEAR ..."
/COM COMPUTERS & STRUCTURES, VOL. 26, NOS. 1/2, 1987
ET,1,HYPER56 ! 2-D 4-NODE MIXED U-P HYPERELASTIC SOLID
ET,2,TARGE169 ! 2-D TARGET ELEMENT
ET,3,CONTA171 ! 2-D CONTACT ELEMENT
R,2,,,2 ! SET STIFFNESS
MP,EX,1,2.82 ! YOUNG'S MODULUS [MPA]
MP,NUXY,1,0.49967 ! POISSON'S RATIO
TB,MOONEY,1 ! MOONEY-RIVLIN CONSTANTS [MPA]
TBDATA,1,0.293
TBDATA,2,0.177
CSYS,1 ! SWITCH TO CYLINDRICAL C.S.
K,1 ! DEFINE KEYPOINTS
K,2,R,-90
K,3,R
K,4,(0.5*R),-90
K,5,(0.6*R),-45
K,6,(0.5*R)
K,7,R,-45
L,2,7
L,7,3
CSYS,0 ! SWITCH TO CARTESIAN C.S.
A,2,7,5,4
A,7,3,6,5
A,4,5,6,1
TSHAPE,LINE
K,1001,-2*R,-R
K,1002,2*R,-R
L,1002,1001
SAVE,temp,db ! SAVE MODEL FOR SECOND ANALYSIS
TYPE,1
AMESH,ALL ! MESH ALL AREAS
REAL,2
TYPE,2
LMESH,10
LSEL,S,LINE,,1,2,1
TYPE,3
LMESH,ALL
FINI
*CREATE,SOLV2D,MAC ! MACRO TO SOLVE MODEL
/SOLU
SOLCONTROL,0
ANTYPE,STATIC
NSEL,S,LOC,X
D,ALL,UX
NSEL,S,LOC,Y
CP,1,UY,ALL
*GET,NCEN,NODE,,NUM,MIN
ALLSEL
CNVTOL,F,,,,-1
NLGEOM,ON ! INCLUDE LARGE DEFORMATION EFFECTS
NSUBST,6 ! SPECIFY NUMBER OF SUBSTEPS IN LOAD STEP
OUTRES,,1 ! WRITE SOLUTION FOR EVERY SUBSTEP
D,NCEN,UY,-0.1 ! APPLY DISPLACEMENT UY = -0.1 TO COUPLED NODES
NROPT,FULL,,OFF
SOLVE
FINISH
*END
VM212 (Modal Analysis of a Rectangular Cavity) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM212
!! MEM 128
/PREP7
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM212, MODAL ANALYSIS OF A RECTANGULAR CAVITY
/COM,********************************************************************************
/COM,
/COM, INTRO TO GUIDED WAVES AND MICROWAVE CIRCUITS
/COM, ROBERT S. ELLIOT, PAGE 264
/COM,
/COM,********************************************************************************
ET,1,HF120,2 ! HF SOLID BRICK ELEMENT, 2ND ORDER
MP,MURX,1,1. ! RELATIVE PERMEABILITY
MP,PERX,1,1. ! RELATIVE PERMITTIVITY
CH=0.3 ! CAVITY HEIGTH
CW=0.4 ! CAVITY WIDTH
CL=1.0 ! CAVITY LENGTH
FRQSTR=1.0E6 ! SHIFT POINT FOR EIGENVALUE EXTRACTION
FRQEND=100E9 ! UPPER BOUND FOR FREQUNECY EXTRACTION
MODES=1 ! NUMBER OF MODES TO EXTRACT
BLOCK,0,CL,0,CW,0,CH ! CREATE CAVITY
LSEL,S,LOC,X,0
LSEL,A,LOC,X,CL
LESIZE,ALL,,,5 ! SET LINE DIVISIONS
LSEL,INVE
LESIZE,ALL,,,10
VMESH,1 ! MESH VOLUME
DA,ALL,AX,0 ! SET ELECTRIC WALL CONDITION (TANGENTIAL E=0)
ASEL,ALL
FINISH
/SOLUTION
ANTYPE,MODAL ! MODAL ANALYSIS
MODOPT,LANB,MODES,FRQSTR,FRQEND,,ON ! BLOCK LANCZOS SOLVER
MXPAND,,,,YES ! EXPAND MODE
SOLVE
FINISH
/POST1
SET,LAST
*GET,MODFRQ,ACTIVE,,SET,FREQ
MODFRQ=MODFRQ/1E9 ! GHZ.
/VIEW,,.75,.5,.6
/VUP,1,Z
/DEVICE,VECTOR,1
EPLOT
PLVECT,H,,,,VECT,NODE,ON ! DISPLAY H FIELD
PLVECT,EF,,,,VECT,NODE,ON ! DISPLAY E FIELD
*DIM,LABEL,CHAR,1
*DIM,VALUE,,1,3
LABEL(1,1)='FREQUENCY'
*VFILL,VALUE(1,1),DATA,.40389
*VFILL,VALUE(1,2),DATA,MODFRQ
*VFILL,VALUE(1,3),DATA,ABS(MODFRQ/.40389)
/OUT,vm212,vrt
/COM,------------------- VM212 RESULTS COMPARISON ----------------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1,1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,' ',F10.7,' ',F10.7,' ',1F5.3)
/COM,--------------------------------------------------------------------
/OUT
FINISH
*LIST,vm212,vrt
VM213 (Harmonic Response Analysis of a Coaxial Cable) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM213
/PREP7
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE,VM213, HARMONIC RESPONSE ANALYSIS OF A COAXIAL CABLE
/COM,********************************************************************************
/COM,
/COM, INTRO TO GUIDED WAVES AND MICROWAVE CIRCUITS
/COM, ROBERT S. ELLIOT
/COM, SECTION 3.3, THE COAXIAL CABLE, PAGE 48
/COM,
/COM,********************************************************************************
ET,1,HF120,2 !HIGH FREQUENCY SOLID BRICK, 2ND ORDER
MP,MURX,1,1. !RELATIVE PERMEABILITY
MP,PERX,1,1. !RELATIVE PERMITTIVITY
FREQ=8E8 !FREQUENCY (HZ.)
CA=0.025 !INNER RADIUS
CB=0.075 !OUTER RADIUS
CL=0.375 !CABLE LENGTH
V0=1.0 !VOLTAGE DIFFERENCE BETWEEN CONDUCTORS
NL=15 !NUMBER OF ELEMENTS ALONG LENGTH
NR=8 !NUMBER OF ELEMENTS ALONG RADIUS
NTH=1 !NUMBER OF ELEMENT THROUGH CIRCUMFERENCE
ANG=5 !CIRCUMFERENTIAL MODEL ANGLE
CYLIND,CA,CB,0,CL,0,5 !CREATE WEDGE MODEL
CSYS,1
LSEL,S,LOC,Z,CL/2
LESIZE,ALL,,,NL,-10 !NL ELEMENTS ALONG LENGTH OF COAX, CLUSTER AT PORTS
LSEL,S,LOC,X,(CA+CB)/2
LESIZE,ALL,,,NR !NR ELEMENTS ALONG THE RADIUS
LSEL,S,LOC,Y,ANG/2
LESIZE,ALL,,,NTH !NTH ELEMENTS AROUND THE CIRCUMFERENCE
VMESH,1 !MESH THE VOLUME
ASEL,S,LOC,X,CA
ASEL,A,LOC,X,CB
DA,ALL,AX,0 !SET ELECTRIC WALL BC (TANGENTIAL E = 0)
LOCAL,11,1
CSYS,0
ASEL,S,LOC,Z,0 !SELECT AREA AT PORT 1 LOCATION
SFA,ALL,,PORT,1 !DEFINE AS PORT 1
PORTOPT,1,COAX,11,CA,CB,V0 !SPECIFY PORT OPTIONS
ASEL,S,LOC,Z,CL !SELECT AREA AT PORT 2 LOCATION
SFA,ALL,,PORT,2 !DEFINE PORT 2
PORT,2,COAX,11, CA,CB, !SPECIFY PORT OPTIONS (MATCHING PORT)
ASEL,ALL
CSYS,1
NBI=NODE(CA,0,CL) !RETRIEVE NODE AT INNER RADIUS
NBO=NODE(CB,0,CL) !RETRIEVE NODE AT OUTER RADIUS
NBA=NODE(CB,ANG,CL) !RETRIVE NODE AT OUTER RADIUS, ANGLE "ANG"
FINISH
/SOLUTION
ANTYPE,HARMIC !FULL HARMONIC ANALYSIS
HARFRQ,FREQ
SOLVE
FINISH
/POST1
SPARM,1,2 !CALCULATE S-PARAMETERS
SET,1,1,
/VIEW,,1,.5,.8
EPLOT
PLVECT,H,,,H,VECT,NODE,ON !DISPLAY H FIELD
PLVECT,EF,,,EF,VECT,NODE,ON !DISPLAY E FIELD
RSYS,1 !RETRIEVE RESULTS IN CYLINDRICAL C.S.
*GET,EFX,NODE,NBO,EF,X !RETRIEVE EF(X) AT NODE NBO
*GET,HY,NODE,NBO,H,Y !RETRIEVE H(Y) AT NODE NBO
RSYS,0
ESEL,S,SFE,PORT,1 !SELECT ELEMENTS AT PORT1
ETABLE,PS1,NMISC,5 !STORE INCIDENT POWER
SSUM !SUM FOR ALL SELECTED ELEMENTS
*GET,PINC,SSUM,,ITEM,PS1 !RETRIEVE INCIDENT POWER
PINC=PINC*360/ANG !CALCULATE POWER FOR FULL 360 DEGREES
ESEL,ALL
!
PATH,VLTG,2 !CREATE PATH FOR VOLTAGE CALC.
PPATH,1,NBI !DEFINE PATH POINTS VIA NODES
PPATH,2,NBO
PATH,CURR,2 !CREATE CURRENT PATH FOR CURRENT CALC.
PPATH,1,NBO !DEFINE PATH POINTS VIA NODES
PPATH,2,NBA
IMPD,'VLTG','CURR',1,360/ANG !CALCULATE IMPEDANCE
EX=12.14
H=.0322
S12=1.0
S11=0.0
P=.01517/2
IMP=65.87
*DIM,LABEL,CHAR,6,2
*DIM,VALUE,,6,3
LABEL(1,1) = 'EX ','HY ','S11 ','S12 ','P ','IMPD'
LABEL(1,2) = '= ','= ','= ','= ','= ','= '
*VFILL,VALUE(1,1),DATA,12.14,.0322,0,1.0,ABS(.01517/2),65.87
*VFILL,VALUE(1,2),DATA,EFX,HY,SII,SIJ,PINC,ZRE
*VFILL,VALUE(1,3),DATA,(EFX/12.14),(HY/.0322),0,(SIJ/1.0),(PINC/(.01517/2)),(ZRE/65.87)
/COM
/OUT,vm213,vrt
/COM,------------------- VM213 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.2,' ',F10.2,' ',1F5.3)
/COM,
/COM,-------------------------------------------------------------
/OUT
FINISH
*LIST,vm213,vrt
VM214 (Harmonic Response of a Rectangular Waveguide) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM214
/PREP7 $ SMRT,OFF
/TITLE,VM214, HARMONIC ANALYSIS OF A TE10-LOADED RECTANGULAR WAVEGUIDE
/COM,********************************************************************************
/COM,
/COM, INTRO TO GUIDED WAVES AND MICROWAVE CIRCUITS
/COM, ROBERT S. ELLIOT
/COM, SECTION 4.3, THE RECTANGULAR WAVEGUIDE, PAGE 82
/COM,
/COM,********************************************************************************
/PREP7 $ SMRT,OFF
ET,1,HF119 ! HIGH FREQUENCY SOLID TETRAHEDRAL
MP,MURX,1,1. ! REALTIVE PERMEABILITY
MP,PERX,1,1. ! RELATIVE PERMITTIVITY
CH=0.01 ! WIDTH
CW=0.03 ! HEIGHT
CL=0.048 ! LENGTH EQUAL TO ONE WAVELENGTH
FREQ=8E9 ! FREQUENCY (HZ.)
LAMDA=3E8/FREQ ! CALCULATE WAVELENGTH
BLOCK,-CW/2,CW/2,0,CH,0,CL ! CREATE WAVEGUIDE
ESIZE,.003 ! SET ELEMENT DIVISIONS
VMESH,1 ! MESH VOLUME
ASEL,S,LOC,X,-CW/2
ASEL,A,LOC,X,CW/2
ASEL,A,LOC,Y,0
ASEL,A,LOC,Y,CH
DA,ALL,AX,0 ! SET ELECTRIC WALL CONDITION (TANGENTIAL E=0)
LOCAL,11
PORT,1,TE10,11,CW,CH,1.0 ! INPUT PORT, 1 VOLT
PORT,2,TE10,11,CW,CH, ! OUTPUT PORT, MATCHED
ASEL,S,LOC,Z,0 ! SELECT AREAS TO ASSIGN INPUT PORT
SFA,ALL,,PORT,1 ! ASSIGN INPUT PORT TO PORT 1
ASEL,S,LOC,Z,CL ! SELECT AREAS TO ASSIGN OUTPUT PORT
SFA,ALL,,PORT,2 ! ASSIGN OUTPUT PORT TO PORT 2
ASEL,ALL
FINISH
/SOLUTION
ANTYPE,HARMIC ! FULL HARMONIC ANALYSIS
EQSLV,ICCG ! SELECT ICCG SOLVER
HARFRQ,FREQ ! SET FREQUENCY
SOLVE
FINISH
/POST1
SPARM,1,2
S11C = SII
S12C = SIJ
SET,1,1
ETABLE,PS1,NMISC,5 ! STORE INCIDENT POWER
/OUT,SCRATCH
SSUM ! SUM FOR ALL SELECTED ELEMENTS
/OUT
*GET,PINC,SSUM,,ITEM,PS1 ! RETRIEVE INCIDENT POWER
N1=NODE(0,0,0) ! INPUT PORT AT CENTERLINE
*GET,EYC,NODE,N1,EF,Y
*GET,HXC,NODE,N1,H,X
S11T=0.0
S12T=1.0
PIN=.1554E-6
EY=1.0
HX=-.2072E-2
*DIM,LABEL,CHAR,5,2
*DIM,VALUE,,5,3
LABEL(1,1) = 'S11 ','S12 ','POWER ','EY ','HX '
LABEL(1,2) = '= ','= ','= ','= ','= '
*VFILL,VALUE(1,1),DATA,S11T,S12T,PIN,EY,HX
*VFILL,VALUE(1,2),DATA,S11C,S12C,PINC,EYC,HXC
*VFILL,VALUE(1,3),DATA,0,(S12C/S12T),(PINC/PIN),(EYC/EY),(HXC/HX)
/COM
/OUT,vm214,vrt
/COM,------------------- VM214 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.8,' ',F10.8,' ',1F5.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,vm214,vrt
VM215 (Thermal-Electric Hemispherical Shell with Hol) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM215
/PREP7
/TITLE,VM215, CONDUCTING SPHERE
ET,1,157
R,1,0.2 ! THICKNESS = 0.2
MP,RSVX,,7 ! DEFINE ELECTRICAL RESISTIVITIES PROPERTY
MP,KXX,,3 ! DEFINE THERMAL CONDUCTIVITIES PROPERTY
CSYS,2 ! SPHERICAL COORDINATE SYSTEM
N,1,10 ! 10 UNIT RADIUS SPHERE
N,21,10,,80 ! 10 DEGREE HOLE AT THE TOP
FILL,,, ,,, ,,0.1 ! SHIFT ELEMENTS TOWARD HOLE
NGEN,2,30,1,21,1,,3 ! ANALYSE A 3 DEGREE SECTOR
E,1,2,32,31 ! DEFINE ELEMENT
EGEN,20,1,-1 ! USE 20 ELEMENTS
CP,1,VOLT,1,31
CP,2,TEMP,1,31
CP,3,VOLT,21,51
CP,4,TEMP,21,51
FINISH
/SOLU
OUTPR,,1
D,21,ALL ! SET ALL VOLTAGES AND TEMPERATURES AT THE HOLE TO ZERO
D,1,VOLT,100 ! SET VOLTAGE AT THE EQUATOR TO 100
SOLVE
FINISH
/POST1
PRRSOL,AMPS ! PRINTS THE CONSTRAINED NODE REACTION AT CURRENT FLOW
PRRSOL,HEAT ! PRINTS THE CONSTRAINED NODE REACTION AT HEAT FLOW
NSEL,S,NODE,,21,51,30
FSUM
*GET,I21,FSUM,0,ITEM,AMPS
*GET,H21,FSUM,0,ITEM,HEAT
*DIM,LABEL,CHAR,2,2
*DIM,VALUE,,2,3
LABEL(1,1) = 'I AMPS ','Q WATT'
LABEL(1,2) = 'NODE 21','NODE 21'
*VFILL,VALUE(1,1),DATA,0.0614,6.14058
*VFILL,VALUE(1,2),DATA,I21,H21
*VFILL,VALUE(1,3),DATA,(I21/0.0614),(H21/6.14058)
/COM
/OUT,vm215,vrt
/COM,------------------- VM215 RESULTS COMPARISON --------------------------
/COM,
/COM, SHELL157 | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A9,A8,' ',F7.5,' ',F7.5,' ',1F5.3)
/COM,------------------------------------------------------------------------
/OUT
FINISH
*LIST,vm215,vrt
VM216 (Lateral Buckling of a Right-Angle Frame) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM216
/TITLE,VM216, LATERAL BUCKLING OF RIGHT-ANGLE FRAME
! J.C. SIMO, L. VU-QUOC, "THREE-DIMENSIONAL FINITE-STRAIN ROD
! MODEL",PART II CMAME, VOL 58, 1986, PP 79-116
/PREP7
N,1 , ,240 !DEFINE NODES
N,21,240,240
N,41,240,0
FILL,1,21
FILL,21,41
N,101,120,360
N,201,360,120
SECNUM,1
ET,1,188 !USE BEAM188 ELEMENT TYPE
*DO,I,1,20,1
E,I,I+1,101
*ENDDO
*DO,I,21,40,1
E,I,I+1,201
*ENDDO
SECTYPE,1,BEAM,ASEC
SECDATA,18.0,1350.0,0.0,0.54,0.0,2.16
MP,EX,1,71240
MP,NUXY,1,0.31
D,1,ALL
FINISH
/SOLU
NLGEOM,ON
NSUBST,2
F,41,FZ,1.0E-3 !APPLY PERTURBATION FORCE
/OUT,SCRATCH, !SUPPRESS SOLUTION DATA
SOLVE
/OUT
OUTRES,ALL,ALL
ARCLEN,ON
ARCTRM,U,60,41,UZ
F,41,FX,1.485 !APPLY END FORCE
NSUBST,10
/OUT,SCRATCH, !SUPPRESS SOLUTION DATA
SOLVE
/OUT
FINISH
/POST26
/AXLAB,X,TIP DISPLACEMENT
/AXLAB,Y,END FORCE
NSOL,2,41,U,Z,DISP
RFORCE,3,1,F,X,FORCE
PROD,4,3, , ,FORCE , , ,-1.0,1,1,
XVAR,2
/COM, THE LOAD DEFLECTION CURVE SHOWN IN VM216.GRPH
/COM, SHOW A CRITICAL LOAD OF APPROX. 1.09
PRVAR,2,4
PLVAR,4
FINISH
/POST1
SET,2,7,1
NSEL,S,LOC,X,0
FSUM
*GET,CP1,FSUM,,ITEM,FX
*DIM,LABEL1,CHAR,1
*DIM,VALUE1,,1,3
LABEL1(1) = 'FX_CRLD'
*VFILL,VALUE1(1,1),DATA,1.09
*VFILL,VALUE1(1,2),DATA,ABS(CP1)
*VFILL,VALUE1(1,3),DATA,ABS(CP1 / 1.09)
/OUT,vm216,vrt
/COM,
/COM,------------------- VM216 RESULTS COMPARISON ---------------------
/COM,
/COM,BEAM188 | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL1(1),VALUE1(1,1),VALUE1(1,2),VALUE1(1,3)
(1X,A15,' ',F10.3,' ',F10.3,' ',1F5.3)
/OUT
FINISH
VM217 (Portal Frame Under Symmetric Loading) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM217
/TITLE,VM217, PORTAL FRAME UNDER SYMMETRIC LOADING
! N. J. HOFF, THE ANALYSIS OF STRUCTURES, PG 115
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
W=-500
A=400
EX=30E6
IO=20300
MROT=((W*A*A*A)/(EX*IO))*(1/27) !CALCULATE MAX ROT TARGET VALUE
BNDM=(W*A*A)*(19/54) !CALCULATE MAX BEND MOMENT
/GRA,POWER
/GST,ON
/TRIAD,OFF
/ESHAPE,1
/PREP7
ET,1,BEAM188
ET,2,BEAM189
SECTYPE,1,BEAM,I
SECDATA,16.655,16.655,36.74,1.68,1.68,.945
SECPLOT,1
C = 1.49535
SECTYPE,2,BEAM,I
SECDATA,C*16.655,C*16.655,C*36.74,C*1.68,C*1.68,C*.945
SECPLOT,2
MP,EX,1,30E6
MP,NUXY,1,0.3
MP,EX,2,30E6
MP,NUXY,2,0.3
A = 400
COLUMDIV = 4
SPANDIV = 16
K,1
K,2,,A
K,3,2*A
K,4,2*A,A
L,2,1
L,3,4
L,4,2
LSEL,,,,1
LATT,,,,,3
LSEL,,,,2
LATT,,,,,1
LSEL,,,,3
LATT,,,,,1
ALLSEL
LESIZE,1,,,COLUMDIV
LESIZE,2,,,COLUMDIV
LESIZE,3,,,SPANDIV
TYPE,1
SECNUM,1
REAL,1
LMESH,1,2
TYPE,2
SECNUM,2
REAL,2
LMESH,3
ALLSEL
DK,1,ALL
DK,3,ALL
LSEL,,,,3
ESLL
SFBEAM,ALL,1,PRESS,-500,-500
ALLSEL
/VIEW,1,1,1,1
/ANG,1
EPLOT
FINISH
/SOLUTION
SOLVE
FINISH
/POST1
ETABLE,SMIS2,SMISC,2
ETABLE,SMIS15,SMISC,15
ETABLE,SMIS5,SMISC,5
ETABLE,SMIS18,SMISC,18
PRNSOL,DOF
PRETAB,SMIS2,SMIS15,SMIS5,SMIS18
/VIEW,1,,,1
PLLS,SMIS2,SMIS15
PLDISP,0
*GET,A,ROTZ,48 !GET ANSYS VALUE FOR MAX ROT
*GET,B,SMIS15,16 !GET ANSYS VALUE FOR MAX BEN MOMENT
/FORMAT,10,E,16
AUTOTS,ON
NSUBST, 400, 1200,25
NLGEOM,ON
NROPT,FULL, ,OFF
OUTRES, ALL, ALL,
SF,ALL,PRES,50.0
NEQITR,20
/AUTO,1
/VIEW,1,,,1
/ANG,1
/ESHAPE,1
EPLOT
FINISH
/SOLUTION
/OUT,SCRATCH
SOLVE
/OUT
FINISH
/POST1
/NOPR !SUPPRESS GRAPHING DATA
/VIEW,1,,-1
/ANG,1
/USER
/FOCUS,1,4,,8,0 !SET UP CENTER OF GRAPHICS SCREEN FOR DISPLACEMENT PLOT
/DIST,,12 !SET DISTANCE TO ZOOM OUT
/TRIAD,OFF
SET,FIRST !SET DISPLACEMENT DATA FOR FIRST SUBSTEP
PLDISP,0 !PLOT DISPLACEMENT DATA
/NOERASE !SET DISPLAY TO OVERLAY PLOTS
SET,,10
PLDISP,0
SET,,20
PLDISP,0
SET,,25
PLDISP,0
SET,LAST
PLDISP,1 !PLOT FINAL DISPLACEMENT WITH ORIGINAL POSITION
/ERASE
/TRIAD,ON
/GOPR
/ESHAPE,0
FINISH
/POST1
/NOPR
SET,NEAR,,,,0.08 !SELECT UZ VALUE FOR NODE 1 AT T=0.08
*GET,VR1,NODE,1,U,Z
SET,NEAR,,,,0.48 !UZ VALUE FOR NODE 1 AT T=0.48
*GET,VR2,NODE,1,U,Z
SET,NEAR,,,,0.76 !UZ VALUE FOR NODE 1 AT T=0.76
*GET,VR3,NODE,1,U,Z
PRES1 = (0.08*50) !SOLVE FOR PRES IN RESULTS TABLE
PRES2 = (0.48*50)
PRES3 = (0.76*50)
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM219
/TITLE,VM219, PRESSURE DRIVEN NON-NEWTONIAN FLOW
! FLUID DYNAMICS BY HUGHES & BRIGHTON- SCHAUM SERIES
R = 1.0 ! PIPE RADIUS
DELTAP= 0.2 ! PRESSURE DROP ACROSS PIPE
L = 0.2 ! LENGTH OF PIPE
NDX = 040 ! NUMBER OF X DIVISIONS
NDY = 002 ! NUMBER OF Y DIVISIONS
RHO = 1.0 ! FLUID DENSITY
MU0 = 1.29684 ! POWER LAW COEFFICIENTS
K = 1.0
D0 = 1E-6
N = 0.8
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
DPDX = DELTAP/L ! PRESSURE GRADIENT
POWER= (N+1)/N
FACTOR = (1/POWER)*( (DPDX/(2*MU0*K))**(1/N) )
/PREP7 $SMRT,OFF
ET,1,FLUID141,,,1 ! 2D YR SYSTEM
ESHAPE,2 ! QUAD ELEMENTS
RECTNG,,R,,L
LESIZE,3,,,NDX,-10
LESIZE,1,,,NDX,-10
LESIZE,4,,,NDY
LESIZE,2,,,NDY
AMESH,ALL
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,vm222
/TITLE,vm222, Warping Torsion Bar
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
!
/PREP7
!
!BUILD MODEL USING BEAM188 ELEMENTS
!
K,,
K,,1000
K,,500,100
L,1,2
ET,1,BEAM188
KEYOPT,1,1,1
SECT,1,BEAM,I
SECD,40,40,80,2,2,2
SECPLOT,1
!
MP,EX, 1, 217396.3331684 !SET MATERIAL PROPS
MP,NUXY,1, 0.335579823862
NEL=50
LATT,1,,,,3 !SET DEFAULT ATTRIBUTES
LESIZE,1,,,NEL
LMESH,ALL !MESH MODEL
DK,1,ALL !PIN MODEL AT KEYPOINT 1
SAVE,PREPDATA,DB !SAVE DB FOR FUTURE PARTS OF TEST
DK,2,ALL !PIN MODEL AT KEYPOINT 2
!
MLOAD1=(1000/NEL)/2
MLOAD2=MLOAD1*2
/ESHAPE,1 !TURN ESHAPE ON TO SHOW 3-D IMAGE OF MODEL
FINISH
/SOLU
OUTRES,ALL,ALL
F,ALL,MX,MLOAD2
FDEL,1,ALL
FDEL,2,ALL
F,1,MX,MLOAD1
F,2,MX,MLOAD1
SOLVE
FINISH
/POST1
PRNSOL,DOF
PRRSOL
NSORT,ROT,X
*GET,ROTX1,SORT, ,MAX !GET VERIFIED RESULT FOR PART 1
PRES,SMISC,27
ETABLE,ROTX,ROT,X
/ESHAPE,1
PLDISP !PLOT DISPLACED SHAPE
PLETAB,ROTX,NOAV
ETABLE,BIMOMENT,SMISC,27
PLETAB,BIMOMENT,NOAV
PARSAV !SAVE PARAMETERS FOR NEXT PART OF TEST
FINISH
/CLEAR,NOSTART
PARRES,CHANGE !RESTORE SAVED PARAMETERS
/PREP7
!
! REBUILD AND REMESH MODEL USING BEAM189 ELEMENTS
!
K,,
K,,1000
K,,500,100
L,1,2
ET,1,BEAM189
KEYOPT,1,1,1
SECT,1,BEAM,I
SECD,40,40,80,2,2,2
MP,EX, 1, 217396.3331684
MP,NUXY,1, 0.335579823862
NEL=50
LATT,1,,,,3
LESIZE,1,,,NEL
LMESH,ALL
!
! APPLY BOUNDRY CONDITIONS AND SET UP NEW PARAMETERS
!
DK,1,ALL
DK,2,ALL
MLOAD=(1000/NEL)/6.0
MLOAD1=MLOAD
MLOAD2=MLOAD*2
MLOAD4=MLOAD*4
FINISH
/SOLU
OUTRES,ALL,ALL
!
! SET UP LOOP TO APPLY LOADS
!
*DO,I,1,NEL,1
J=NELEM(I,1)
K=NELEM(I,2)
L=NELEM(I,3)
F,J,MX,MLOAD2
F,K,MX,MLOAD2
F,L,MX,MLOAD4
*ENDDO
F,1,MX,MLOAD1
F,2,MX,MLOAD1
FLIST
SOLVE
FINISH
/POST1
PRNSOL,DOF
NSORT,ROT,X
*GET,ROTX2,SORT, ,MAX !GET VERIFIED RESULTS FOR PART 2
PRRSOL
PRES,SMISC,4
PRES,SMISC,17
PRES,SMISC,27
T1=ROTX1
T2=ROTX2
!
!SET UP AND FILL VM RATIO TABLE
!
RAT_1 = ABS(T1/0.3292617E-3)
RAT_2 = ABS(T2/0.3292617E-3)
*DIM,LABEL,CHAR,2
*DIM,VALUE,,2,2
*DIM,RESULTS,,2,1
LABEL(1) = 'BEAM188','BEAM189'
*VFILL,VALUE(1,1),DATA,0.3292617E-3,0.3292617E-3
*VFILL,RESULTS(1,1),DATA,T1,T2
*VFILL,VALUE(1,2),DATA,RAT_1,RAT_2
/OUT,vm222,vrt
/COM
/COM,------------------- VM222 RESULTS COMPARISON ---------------------
/COM, MX TWIST
/COM, IN X-DIR | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),RESULTS(1,1),VALUE(1,2)
(1X,A10,' ',E10.4,' ',E10.4,' ',1F5.3)
/COM,-------------------------------------------------------------------
/out
*LIST,vm222,vrt
FINISH
VM223 (Frequency Response of a Pre-stressed Beam) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM223
! SET UP JPEG GRAPHIC SETTINGS TO CREATE IMAGES FOR MANUAL !
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/PREP7
MP,PRXY,,0.3
/TITLE, VM223, FREQ. RESPONSE OF A PRESTRESSED BEAM
/COM REFERENCE: "FORMULAS FOR NATURAL FREQUENCY AND MODE SHAPE"
/COM PG. 144, EQN 8-20, R.D. BLEVINS, VAN NOSTRAND REINHOLD CO.
/COM 1979
/NOPR
GAPI=1 ! INITIAL GAP (MICRONS)
GAP=GAPI-P*L/E/B/H ! APPROXIMATE DEFLECTED GAP (IGNORING CAP STIFFNESS) (MICROMETER)
C3=PER0*PLATEA ! TRANDUCER REAL CONTANT
C3P=C3/(GAP**2) ! DERIVITIVE OF C3
C3PP=2*C3/(GAP**3) ! SECOND DERIVITIVE OF C3
VLT=SQRT(2*P/C3P) ! APPLIED VOLTAGE TO PLATE
KUU=C3PP*VLT**2/2 ! GAP STIFFNESS
KBEAM=E*PLATEA/L ! BEAM STIFFNESS (NOTE: GAP STIFFNESS ASSUMED << BEAM STIFFNESS)
UX2=P*L/B/H/E ! DESIRED DEFLECTION
*DIM,FREQ,,5 ! ARRAY PARAMETER FOR BEAM FREQUENCY
*DIM,PFREQ,,5 ! ARRAY PARMETER FOR BEAM PRE-STRESSED FREQUENCY
PI=4*ATAN(1)
*DO,J,1,5
LAMDA=J*PI
LAMDAP2=LAMDA**2*SQRT((1 + P*L**2/(E*I*LAMDA**2)))
LAMDAP=SQRT(LAMDAP2)
!
FREQ(J) = LAMDA**2/(2*PI*L**2)*SQRT(E*I/M)
PFREQ(J) = LAMDAP**2/(2*PI*L**2)*SQRT(E*I/M)
!
*ENDDO
!
ET,1,3
R,1,B*H,I,H ! BEAM PROPERTIES
MP,EX,1,E
MP,DENS,1,DENS
!
ET,2,126
C3=PER0*PLATEA
R,2,0,0,1,0
rmore,c3
!
! BUILD MODEL
!
N,1,-10
N,2,0
N,22,L
FILL
TYPE,2
REAL,2
E,1,2
TYPE,1
REAL,1
E,2,3
*REPEAT,20,1,1
!
!APPLY BOUNDRY CONDITIONS AND LOADS
!
NSEL,S,LOC,X,-10
NSEL,A,LOC,X,L
D,ALL,UX,0,,,,UY
NSEL,S,LOC,X,0
D,ALL,UY,0
D,ALL,VOLT,VLT
IC,ALL,VOLT,VLT
NSEL,S,LOC,X,-10
D,ALL,VOLT,0
NSEL,ALL
FINI
!
! ENTER SOLUTION TO SOLVE STATIC PART OF ANALYSIS
!
/SOLU
ANTYP,STATIC
PSTRES,ON
SOLVE
FINI
!
! NOW SOLVING FOR MODAL SOLUTION
!
/SOLU
ANTYP,MODAL
MODOPT,UNSYM,3 ! EXTRACT 3 MODES
MXPAND
PSTRES,ON
SOLVE
FINISH
!
! ENTER GENERAL POSTPROCESSOR TO DISPLAY RESULTS
!
/POST1
SET,1,1
*GET,FP1,ACTIVE,,SET,FREQ
SET,1,2
*GET,FP2,ACTIVE,,SET,FREQ
SET,1,3
*GET,FP3,ACTIVE,,SET,FREQ
! PLOT DISPLACED SHAPE AT THE THREE FREQUENCIES
SET,,,1,,FP1
/TITLE,VM223 RESULT PLOT FREQ1
PLDISP,1
SET,,,1,,FP2
/TITLE,VM223 RESULT PLOT FREQ2
PLDISP,1
SET,,,1,,FP3
/TITLE,VM223 RESULT PLOT FREQ3
PLDISP,1
FINISH
!
! SET UP ANALYTICAL SOLUTION
!
PFREQ1=PFREQ(1)
PFREQ2=PFREQ(2)
PFREQ3=PFREQ(3)
!
! SET UP VM RATIO TABLE
!
RAT_1 = ABS(FP1/PFREQ1)
RAT_2 = ABS(FP2/PFREQ2)
RAT_3 = ABS(FP3/PFREQ3)
*DIM,LABEL,CHAR,3
*DIM,VALUE,,3,2
*DIM,RESULTS,,3,1
LABEL(1) = 'FREQ1','FREQ2','FREQ3'
*VFILL,VALUE(1,1),DATA,PFREQ1,PFREQ2,PFREQ3
*VFILL,RESULTS(1,1),DATA,FP1,FP2,FP3
*VFILL,VALUE(1,2),DATA,RAT_1,RAT_2,RAT_3
/OUT,vm223,vrt
/COM
/COM,------------------- VM223 RESULTS COMPARISON ---------------------
/COM,
/COM, PRE-STRESS | TARGET | ANSYS | RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),RESULTS(1,1),VALUE(1,2)
(1X,A10,' ',E10.4,' ',E10.4,' ',1F5.3)
/COM,-------------------------------------------------------------------
/OUT
FINISH
*LIST,vm223,vrt
VM224 (Implicit Creep under Biaxial Load) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,vm224
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE,VM224, Implicit Creep Under Biaxial Load
/COM, NAFEMS Fundamental Tests of Creep Behavior, Becker and Hyde
/NOPR
/COM,
/COM, 2D CREEP TESTS WITH BIAXIAL CONSTANT LOAD,
/COM, REFERENCE: TEST 10(A) FROM NAFEMS R0027.
/COM,
/COM, EXPECTED RESULTS:
/COM, TIME | EPCRX | EPCRY
/COM, ----------------------------
/COM, 0.0 | 0.0 | -0.0
/COM, 0.1 | 0.0427 | -0.0427
/COM, 1.0 | 0.135 | -0.135
/COM, 5.0 | 0.3019 | -0.3019
/COM, 10.0 | 0.4269 | -0.4269
/COM, 50.0 | 0.9546 | -0.9546
/COM, 100.0 | 1.35 | -1.35
/COM, 500.0 | 3.019 | -3.019
/COM, 1000.0 | 4.2691 | -4.2691
/COM,
/PREP7
C*** PARAMETRIC INPUT FOR CREEP CONSTANTS (PRIMARY CREEP)!***
*SET,C1,1.5625E-14 !ASSIGN VALUE
*SET,C2,5.0 !ASSIGN VALUE
*SET,C3,-0.5 !ASSIGN VALUE
*SET,C4,0 !ASSIGN VALUE
C*** TIME PARAMETER
*SET,HOUR,1000 !ASSIGN VALUE
C*** ELASTIC CONSTANT
MP,EX,1,200E3 !DEFINE YOUNG'S MODULUS
MP,NUXY,1,0.3 !DEFINE POISON'S RATIO
TUNIF,HOT !ASSIGN TEMP TO NODES
TOFF,OFFS !SPECIFY TEMP RELATIVE TO ABSOLUTE VALUES
TB,CREEP,1,,,6 !ACTIVATE DATA TABLE
TBDATA,1,C1,C2,C3,C4 !DEFINE DATA FOR TABLE
SAVE !SAVE
/PREP7
N,1,0,0,0
N,2,100,0
N,3,100,100
N,4,0,100
N,5,50,0
N,6,100,50
N,7,50,100
N,8,0,50
N,9,50,50
ET,1,PLANE182
KEYOPT,1,1,1
KEYOPT,1,3,0
E,1,5,9,8
E,5,2,6,9
E,9,6,3,7
E,8,9,7,4
NALL
NSEL,S,LOC,X
D,ALL,UX
NSEL,S,LOC,Y
D,ALL,UY
NALL
/SOLU
NSEL,S,LOC,X,100
SF,ALL,PRES,-200
NALL
NSEL,S,LOC,Y,100
SF,ALL,PRES,200
NALL
RATE, OFF
DELT,1.0E-8,1.0E-9,1.0E-8
TIME, 1.0E-8
/OUT,SCRATCH
OUTRES,ESOL,ALL
SOLVE
/OUT
RATE, ON, ON
DELT,1E-5,1E-5,100
TIME,1000
/OUT,SCRATCH
SOLVE
/OUT
FINISH
/POST26
ESOL,2,1,,EPCR,X
ESOL,3,1,,EPCR,Y
PRVAR,2,3
PLVAR,2,3
*GET,RES1X,VARI,2,RTIME,1000
*GET,RES1Y,VARI,3,RTIME,1000
FINISH
PARSAV,ALL
RESUME
PARRES,CHANGE
/PREP7
N,1,0,0,0
N,2,100,0
N,3,100,100
N,4,0,100
N,5,50,0
N,6,100,50
N,7,50,100
N,8,0,50
ET,1,PLANE183
KEYOPT,1,3,0
E,1,2,3,4,5,6,7,8
NALL
NSEL,S,LOC,X
D,ALL,UX
NSEL,S,LOC,Y
D,ALL,UY
NALL
/SOLU
NSEL,S,LOC,X,100
SF,ALL,PRES,-200
NALL
NSEL,S,LOC,Y,100
SF,ALL,PRES,200
NALL
RATE, OFF
DELT,1.0E-8,1.0E-9,1.0E-8
TIME, 1.0E-8
/OUT,SCRATCH
SOLVE
/OUT
RATE, ON, ON
DELT,1E-5,1E-5,100
TIME,1000
/OUT,SCRATCH
OUTRES,ESOL,ALL
SOLVE
/OUT
FINISH
/POST26
ESOL,2,1,,EPCR,X
ESOL,3,1,,EPCR,Y
PRVAR,2,3
PLVAR,2,3
*GET,RES2X,VARI,2,RTIME,1000
*GET,RES2Y,VARI,3,RTIME,1000
VM225 (Rectangular Cross-Section Bar with Preload) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,VM225
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE, VM225, BEAM WITH PRETENSION LOAD
/COM,
/PREP7 !ENTER PREPROCESSOR
MP,PRXY,,0.3
BLOCK,,12,,1,,.5 !CREATE BLOCK
/VIEW,1,1,2,3 !CHANGE VIEW
ET,1,SOLID45 !SET ELEMENT TYPE
MP,EX,1,30E6, !DEFINE YOUNG'S MODULUS
MP,PRXY,1,0.3 !DEFINE POISSON'S RATIO
MP,DENS,1,.283, !DEFINE DENSITY
DA,2,SYMM !DEFINE SYMMETRY BC ON AREA
DA,3,SYMM !DEFINE SYMMETRY BC ON AREA
DA,5,UX !DEFINE DOF ON AREA
DA,6,UX !DEFINE DOF ON AREA
VMESH,ALL !MESH VOLUME
PSMESH, , ,1000,ALL,,0,X,6,,,,EEE,NNN !APPLY LOAD VIA PRE-TENSION ELEMENT
SLOAD,1,9,LOCK,F,125,1,2
EPLOT
FINI !EXIT PREP7
/SOLU !ENTER SOLVER
SOLVE !SOLVE
FINISH !EXIT SOLVER
/POST1 !ENTER POST PROCESSOR
NSORT,S,INT,1,0,, !SORT STRESS RESULTS
*GET,MAXNFEA2,SORT,,IMAX !GET NODE VALUE
*GET,SIGFEA2,NODE,MAXNFEA2,S,INT !GET MAXIMUM VON MISSES STRESS
NSORT,U,X,SUM,1,0,,, !SORT DEFLECTION RESULTS
*GET,MAXNFEA2,SORT,,IMAX !GET NODE VALUE
*GET,UXFEA2,NODE,MAXNFEA2,U,X,SUM !GET MAXIMUM UX VALUE
*STAT,UXFEA2 !LIST PARAMETER VALUE
*STAT,SIGFEA2 !LIST PARAMETER VALUE
/COM,**********************************
/COM,*** CLASSICAL ANALYSIS RESULTS ***
/COM,**********************************
SIGCA = 250
UXCA = 0.0001
VM226 (Fourier Series Analysis of a Diode Rectified Circuit) Input Listing
/COM,ANSYS MEDIA REL. 5.7.1 (4/9/01) REF. VERIF. MANUAL: REL. 5.7.1
/VERIFY,vm226
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TITLE,VM226, FOURIER SERIES ANALYSIS OF A DIODE RECTIFIED CIRCUIT
/COM, REF: SEDRA/SMITH "MICROELECTRONIC CIRCUITS 4TH ED." SEC. 3.7
/COM,
/COM FIRST PART: NO CAPACITANCE
/COM
PI = 4*ATAN(1)
R1=2500 ! RESISTOR VALUE
OMEGA=2*PI*60
IFINAL=3 ! NUMBER OF MODES ( 3 -> A0,A1,B1)
IFIN=IFINAL+20 ! WE MUST COMPUTE MORE COEFF THAN WE NEED
U=135
! FOR THE SECOND PART
*DIM,CAPA,ARRAY,3
CAPA(1)=1E-6,10E-6,1E-3 ! CAPACITOR VALUES
EPS=1E-09 ! ERROR CRITERIA FOR TO1
/PREP7
R,1,,U,OMEGA/2/PI, !SET UP SINUSOIDAL VOLTAGE SOURCE
N,1,-0.85,0.4,0
N,2,-0.85,0.25,0
RMOD,1,15,0,1
ET,1,CIRCU124,4,1
TYPE,1
REAL,1
MAT,1
!
N,3,-0.85,0.325,0
E,1,2,3 !CREATE IND. SINUSOIDAL VOLT SOURCE
R,2,R1, !SET UP 2500 OHM RESISTOR
N,4,-0.75,0.4,0
N,5,-0.75,0.25,0
RMOD,2,15,0,2
ET,2,CIRCU124,0,0
TYPE,2
REAL,2
MAT,1
E,4,5 !CREATE 2500 OHM RESISTOR
!
! THE FOLLOWING COMMANDS ARE USED TO SET UP THE IDEAL DIODE
!
ET,3,125,
R,3
TYPE,3
REAL,3
E,1,4
!
! APPLY GROUND TO CIRCUIT
!
D,2,VOLT,0
D,5,VOLT,0
SAVE
ALLS
EPLOT
FINISH
!
! SOLVE NON-LINEAR CIRCUIT WITH T = 0 TO 0.025
! USING A TIME STEP OF 0.001 FOR EACH ITERATION
!
/SOLU
ANTYPE,TRANS
OUTRES,ALL,ALL,
TIME,0.025
AUTOTS,-1
DELTIM,0.0001, , ,1
CNVTOL,VOLT,,0.0001,2,1.0E-6 !CONVERGANCE CRITERIA
/OUT,SCRATCH
SOLVE
/OUT
FINISH
/POST26
NSOL,2,4,VOLT,,
/COLOR,CURVE,BLUE,1
/TITLE,VM226, LOAD VOLTAGE WAVEFORM WITH NO CAPACITANCE
/AXLAB,Y,OUTPUT POTENTIAL (VOLT)
PLVAR,2,
/TITLE,VM226, FOURIER SERIES ANALYSIS OF A DIODE RECTIFIED CIRCUIT
!
! SET UP TABLE ARRAYS TO DISPLAY OUTPUT RESULTS
!
*DIM,VOLTG,TABLE,251
*DIM,TARGET,ARRAY,251
*DO,INC,1,251,1
T = INC*0.0001
!TIME(INC) = T
*GET,V,VARI,2,RTIME,T
VOLTG(INC,1) = V
VOLTG(INC,0) = T
ANAL = U*SIN(OMEGA*T)
*IF,ANAL,LT,0,THEN !SET TARGET TO ZERO IF ANALYTICAL SOLUTION
TARGET(INC) = 0 ! IS NEGATIVE
*ELSE
TARGET(INC) = ANAL
*ENDIF
*ENDDO
FINISH
*DIM,COEFF,,IFIN
*DIM,MODE,TABLE,IFIN
*DIM,ISYM,TABLE,IFIN
*DIM,THETA,TABLE,121
*DIM,CURVEI,TABLE,121 ! CURVE INPUT TO PROGRAM
*DIM,CURVEO,TABLE,121
*VFILL,THETA(1),RAMP,0,3 ! THETA VALUES INCREMENT 3 DEGREES
*DO,INC,1,121,1
T=(INC-1)*3/360*2*PI/OMEGA
CURVEI(INC)=VOLTG(T,1)
*ENDDO
! CALCULATE FOURIER COEFFICIENT
MODE(1)=0
ISYM(1)=1
ISTART=2
/COM
/COM *** *DO *********************************************************
/COM
*DO,I,ISTART,IFIN,2
MODE(I)=I/2 ! FILL EVEN INDICIES OF {MODE}
ISYM(I)=1
*ENDDO
/COM
/COM *** *ENDDO WAS LAST COMMAND USED *************************************
/COM
! ! FILL ODD INDICIES OF {MODE}
ISTART=3
*DO,I,ISTART,IFIN,2
MODE(I)=(I/2)-.5
ISYM(I)=-1
*ENDDO
*MFOURI,FIT, COEFF(1),MODE(1),ISYM(1),THETA(1),CURVEI(1)
!
! CURVE WHICH WILL BE DEVELOPED FROM GENERATED COEFFICIENTS
!
*MFOURI,EVAL,COEFF(1),MODE(1),ISYM(1),THETA(1),CURVEO(1)
! PLOT CURVE
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TRIAD,OFF
/PLOPTS,LOGO,0
/PLOPTS,INFO,2
/PLOPTS,WP,0
/COLOR,CURVE,CBLU
/XRANGE,0,370
/YRANGE,0,140
/TSPEC,15
/TLAB,1,0.75,CAPACITANCE = 0
/TSPEC,4
/TLAB,1,0.7,BLUE->ANSYS
/TSPEC,1
/TLAB,1,0.65,RED->FOURIER
*VPLOT,THETA(1),CURVEI(1) ! PLOT INPUT CURVE VERSUS THETA
/USER
/NOERASE
/COM OVERLAY THE OUTPUT CURVE ON THE INPUT CURVE
/COLOR,CURVE,RED
/COLOR,AXLAB,BLAC
/AXLAB,X,ANGLE IN DEGREE
/AXLAB,Y,VOLT
*VPLOT,THETA(1),CURVEO(1) ! PLOT OUTPUT CURVE VERSUS THETA
/ERASE
! ANALYTICAL FOURIER COEFFICIENT
RESUL(1,2)=2*U/PI/2 ! FIRST FOURIER COEFFICIENT = A0/2
ISTART=4
*DO,I,ISTART,IFINAL,4
RESUL(I,2)=-2*U/(PI*((I/2)**2-1))
*ENDDO
*DO,I,2,IFINAL,4
RESUL(I,2)=0
*ENDDO
RESUL(3,2)=U/2
ISTART=5
*DO,I,ISTART,IFINAL,2
RESUL(I,2)=0
*ENDDO
*DO,I,1,IFINAL
RESUL(I,1)=COEFF(I)
RESUL(I,3)=RESUL(I,2)/RESUL(I,1)
*ENDDO
*DO,CAP,1,3 ! START DO LOOP ON CAPACITANCE
PARSAV,ALL
/CLEAR,NOSTART
/COM SECOND PART: CAPACITANCE VALUE
/COM
/PREP7
RESUME
PARRES,CHANGE
C1=CAPA(CAP)
TAU=TAUARR(1,CAP)
N,6,-0.65,0.4,0
ET,4,CIRCU124,2
R,4,C1
TYPE,4
REAL,4
E,4,6
!
! APPLY GROUND TO CIRCUIT
!
D,6,VOLT,0
ALLS
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
EPLOT
FINISH
!
! SOLVE NON-LINEAR CIRCUIT WITH T = 0 TO 0.025
! USING A TIME STEP OF 0.0001 FOR EACH ITERATION
!
/SOLU
ANTYPE,TRANS
OUTRES,ALL,ALL,
TIME,0.025
AUTOTS,-1
DELTIM,0.0001, , ,1
CNVTOL,VOLT,,0.0001,2,1.0E-6 !CONVERGANCE CRITERIA
/OUT,SCRATCH
SOLVE
/OUT
FINISH
/POST26
NSOL,2,4,VOLT,,
/COLOR,CURVE,BLUE,1
*IF,CAP,EQ,1,THEN
/TITLE,VM226, VLOAD WITH CAPACITANCE OF 1E-6F
*ENDIF
*IF,CAP,EQ,2,THEN
/TITLE,VM226, VLOAD WITH CAPACITANCE OF 10E-6F
*ENDIF
*IF,CAP,EQ,3,THEN
/TITLE,VM226, VLOAD WITH CAPACITANCE OF 1E-3F
*ENDIF
/AXLAB,Y,OUTPUT POTENTIAL (VOLT)
PLVAR,2,
/TITLE,VM226, FOURIER SERIES ANALYSIS OF A DIODE RECTIFIED CIRCUIT
!
! DETERMINE T0
!
T0=1/OMEGA*ATAN(1/(OMEGA*TAU))
!
! DETERMINE T0' : MACRO TO DO A BISECTION BETWEEN THE TWO CURVES
!
TINIT=2*PI/OMEGA/4
TFINAL=2*PI/OMEGA
V1=U*COS(OMEGA*TINIT)-U*COS(OMEGA*T0)*EXP(-1*(TINIT-T0)/TAU)
V2=U*COS(OMEGA*TFINAL)-U*COS(OMEGA*T0)*EXP(-1*(TFINAL-T0)/TAU)
V0=U*COS(OMEGA*T0)
*DO,I,1,10000
T3=(TINIT+TFINAL)/2
V3=U*COS(OMEGA*T3)-V0*EXP(-1*(T3-T0)/TAU)
ERROR=ABS(V1-V2)
*IF,ERROR,LT,EPS,THEN
*EXIT
*ENDIF
*IF,V3,LT,0,THEN
TINIT=T3
V1=V3
*ELSE
TFINAL=T3
V2=V3
*ENDIF
*ENDDO
T01=T3
!
! SET UP TABLE ARRAYS TO DISPLAY OUTPUT RESULTS
!
*DIM,VOLTG,TABLE,250
*DIM,TARGET,ARRAY,250
*DO,INC,1,250,1
T = INC*0.0001
!TIME(INC) = T
*GET,V,VARI,2,RTIME,T
VOLTG(INC,1) = V
VOLTG(INC,0) = T
ANAL = U*SIN(OMEGA*T)
*IF,T,LE,T0+PI/2/OMEGA,THEN !SET TARGET TO ZERO IF ANALYTICAL SOLUTION
TARGET(INC) = ANAL ! IS NEGATIVE
*ELSE
*IF,T,LE,T01+PI/2/OMEGA,THEN
TARGET(INC) = U*COS(OMEGA*T0)*EXP(-(T-T0-PI/2/OMEGA)/TAU)
*ELSE
*IF,T,LE,T0+PI/2/OMEGA+2*PI/OMEGA,THEN
TARGET(INC) = ANAL
*ELSE
TARGET(INC) = U*COS(OMEGA*T0)*EXP(-(T-T0-PI/2/OMEGA-2*PI/OMEGA)/TAU)
*ENDIF
*ENDIF
*ENDIF
*ENDDO
FINISH
*DIM,COEFF,,IFIN
*DIM,MODE,TABLE,IFIN
*DIM,ISYM,TABLE,IFIN
*DIM,THETA,TABLE,121
*DIM,CURVEI,TABLE,121 ! CURVE INPUT TO PROGRAM
*VFILL,THETA(1),RAMP,0,3 ! THETA VALUES INCREMENT 3 DEGREES
*DO,INC,1,121,1
T=(INC-1)*3/360*2*PI/OMEGA+PI/2/OMEGA
CURVEI(INC)=VOLTG(T,1)
*ENDDO
! CALCULATE FOURIER COEFFICIENT
MODE(1)=0
ISYM(1)=1
ISTART=2
/COM
*DO,I,ISTART,IFIN,2
MODE(I)=I/2 ! FILL EVEN INDICIES OF {MODE}
ISYM(I)=1
*ENDDO
/COM
! ! FILL ODD INDICIES OF {MODE}
ISTART=3
*DO,I,ISTART,IFIN,2
MODE(I)=(I/2)-.5
ISYM(I)=-1
*ENDDO
*MFOURI,FIT, COEFF(1),MODE(1),ISYM(1),THETA(1),CURVEI(1)
*IF,CAP,EQ,1,THEN
!
! CURVE WHICH WILL BE DEVELOPED FROM GENERATED COEFFICIENTS
!
*MFOURI,EVAL,COEFF(1),MODE(1),ISYM(1),THETA(1),CURVEO(1)
! PLOT CURVE
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/TRIAD,OFF
/PLOPTS,LOGO,0
/PLOPTS,INFO,2
/PLOPTS,WP,0
/COLOR,CURVE,CBLU
/XRANGE,0,370
/YRANGE,0,140
/TSPEC,15
/TLAB,-0.25,0.75,CAPACITANCE = 1E-06 FARAD
/TSPEC,4
/TLAB,0,0.7,BLUE->ANSYS
/TSPEC,1
/TLAB,0,0.65,RED->FOURIER
*VPLOT,THETA(1),CURVEI(1) ! PLOT INPUT CURVE VERSUS THETA
/USER
/NOERASE
/COM OVERLAY THE OUTPUT CURVE ON THE INPUT CURVE
/COLOR,CURVE,RED
/COLOR,AXLAB,BLAC
/AXLAB,X,ANGLE IN DEGREE
/AXLAB,Y,VOLT
*VPLOT,THETA(1),CURVEO(1) ! PLOT OUTPUT CURVE VERSUS THETA
/ERASE
*ENDIF
!
! ANALYTICAL FOURIER COEFFICIENT
!
*DIM,ANALY,ARRAY,IFINAL
LSEL,S,,,1,4
LSEL,A,,,13,16
SFL,ALL,RDSF,1, ,1, !SET ALL FACING SURFACES TO EMISSIVITY 1
LSEL,S,,,9,12
DL,ALL, ,TEMP,TOUT,1 !APPLY UNIFORM TEMPERATURE TO EXTERIOR
LSEL,S,,,5,8
DL,ALL, ,TEMP,TIN,1 !APPLY UNIFORM TEMPERATURE TO INTERIOR
ASEL,S, , , 1
ESLA,S
CM,INSIDE,ELEM
ASEL,S, , , 2
ESLA,S
CM,OUTSIDE,ELEM
FINI
/AUX12
ALLSEL
!*
STFCONST=0.119E-10
STEF,STFCONST !SET STEFAN-BOLTZMAN CONSTANT FOR MODEL
TOFFST,0.E+00 !SET TEMPERATURE OFFSET
RADOPT,0.1,0.1,0.E+00,1000,0.1,0.1 !SET RADIOSITY OPTIONS
SPCTEMP,1,0.E+00 !SET TEMPERATURE FOR RADIATION TO SPACE
!(NO RADIATION WILL APPEAR IN THIS MODEL)
!*
VFCALC !CALCULATE RADIOSITY VIEW FACTORS
VFQUERY,INSIDE,INSIDE !EXTRACT VIEW FACTOR FROM INTERIOR CYLINDER ELEMENTS DUE TO INTERIOR CYLINDER
*GET,VFAVG1,RAD,,VFAVG !DETERMINE VIEW FACTOR FOR ENTIRE SURFACE
VFQUERY,OUTSIDE,INSIDE !EXTRACT VIEW FACTOR FROM EXTERIOR CYLINDER ELEMENTS DUE TO INTERIOR CYLINDER
*GET,VFAVG2,RAD,,VFAVG !DETERMINE VIEW FACTOR FOR ENTIRE SURFACE
VFQUERY,OUTSIDE,OUTSIDE !EXTRACT VIEW FACTOR FROM EXTERIOR CYLINDER ELEMENTS DUE TO EXTERIOR CYLINDER
*GET,VFAVG3,RAD,,VFAVG !DETERMINE VIEW FACTOR FOR ENTIRE SURFACE
FINISH
ALLSEL
/SOLU
SOLV
FINI
/POST1
LSEL,S, , ,1,4
LSEL,A, , ,13,16
NSLL,S,1