/TITLE, Your name, ME 421, and date here, Electr2a.a
 ! Electric bridge circuit, O'Malley
 ! diamond shape

/PREP7             ! preprocessor section
volt1 = 100
volt5 = 0

ET,1,Link68        ! thermal-electrical line element
R,1,(1/2)          ! length / resistance (area), 2 ohms
R,2,(SQRT(2)/5)    ! diag length / 5 ohms
R,3,(SQRT(2)/3)    ! diag length / 3 ohms
R,4,(2/10)         ! length / 10 ohms
R,5,(SQRT(2)/4)    ! diag length / 4 ohms
MP,KXX             ! Thermal cond not used but is required
MP,RSVX,1,1        ! electrical resistivity

N,1,0,2
N,2,1,2
N,3,0,1
N,4,2,1
N,5,1,0
/PNUM,node,1       ! plot the nodes
APLOT,all          ! plot all areas
/PLOPTS,INFO,OFF   ! Removes entire legend column
/PLOPTS,FRAME,OFF  ! Removes frame

TYPE,1             ! element type (ET) 68
MAT,1              ! material property MP 1
REAL,1             ! 2 ohms
E,1,2
REAL,2             ! 5 ohms
E,2,3
E,2,4
REAL,3             ! 3 ohms
E,4,5
REAL,4             ! 10 ohms
E,3,4
REAL,5             ! 4 ohms
E,3,5
FINISH             ! preprocessor section

/SOLUTION          ! solution section
antype,static
outpr,basic,1
outpr,nload,1
KRF,1              ! print nodal forces
KBC,1              ! step boundary conditions
D,1,VOLT,volt1     ! nodal voltage, top
D,5,VOLT,volt5     ! nodal voltage, bottom
/output,Volts,txt,D:\Ansys\   !  volts, node currents
! /output,Volts,txt,U:\win\Ansys\   !  volts, node currents

solve
FINISH           ! solution section

/POST1           ! post processor
etable,cur,nmisc,5,4
pldisp,1             ! plot displacements
/output,Curr,txt,D:\Ansys\   !  branch currents
! /output,Curr,txt,U:\win\Ansys\   !  branch currents
pretab
! FINISH           ! post processor

!    NODE       VOLT
!      1      100.000
!      2      67.8930
!      3      29.4314
!      4      26.0870
!      5       0.

!    NODE       AMPS
!      1      16.0535
!      5     -16.0535

!   ELEM   CURRENT
!      1  16.054
!      2   7.6923
!      3   8.3612
!      4   8.6957
!      5   0.33445
!      6   7.3579