×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

DO loop that modifies loads (APDL) (HELP)

DO loop that modifies loads (APDL) (HELP)

DO loop that modifies loads (APDL) (HELP)

(OP)
Hello,

Im trying to make a transient simulation in wich I introduce a heatflux proportional to the reactionforces in each node.
How can i get the reactin forces of each node manipulate them and input it as a new heatflux in each node? thanks !



/BATCH

!Rectangle
/PREP7
RECTNG,0,0.1,0,0.01,

!Create element
ET,1,PLANE223,11

!Material
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e11
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,
MPDATA,ALPX,1,,0.00015
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,45
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,460

!Meshing
AESIZE,ALL,0.002,
MSHAPE,0,2D
MSHKEY,1
CM,_Y,AREA
ASEL, , , ,1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
AMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
! /REPLOT,RESIZE
FINISH

!Solution settings
/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0

!Initial temperature
TUNIF,25

!Fixing vertical direction
NSEL,s,loc,y,0
NSEL,a,loc,y,0.01
D,ALL,UY,0
ALLSEL

!Fixing horizontal direction
NSEL,s,loc,y,0
NSEL,r,loc,x,0
D,ALL,Ux,0
ALLSEL

!Loop
*DO, i, 1, 60, 1
TIME,i
KBC,0
NSEL,s,loc,y,0
NSEL,a,loc,y,0.01
SF,ALL,HFLUX,100000,
ALLSEL
SOLVE
FINISH
/POST1
! Here i must get the reaction forces and update the heatflux in some way for the next step
FINISH
/SOLU
*ENDDO


RE: DO loop that modifies loads (APDL) (HELP)

I'm assuming when you say reaction force, you mean at the restrained nodes. I think this is what you are looking for:

/post1
set,last
NSEL,s,loc,y,0
NSEL,a,loc,y,0.01
cm,fluxnodes,node
ncount=ndinqr(0,14)
*dim,fluxes,table,ncount,,,node
nnum=ndnext(0)
*do,loop,1,ncount
fluxes(loop,0)=nnum
nsel,s,node,,nnum
fsum
*get,reactforc,fsum,0,item,fy
fluxes(loop,1)=10000+reactforc/1e3
cmsel,s,fluxnodes
nnum=ndnext(nnum)
*enddo

/solu

cmsel,s,fluxnodes
sf,all,hflux,%fluxes%
allsel

Fill in whatever math you want at fluxes(loop,1)=

You might want to take a look at your reaction forces at the midside nodes. They are an order of magnitude smaller then at the corner nodes, and that may give you an unwanted flux value there. Maybe use PLANE222?

Rick Fischer
Principal Engineer
Argonne National Laboratory

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members! Already a Member? Login


Resources


Close Box

Join Eng-Tips® Today!

Join your peers on the Internet's largest technical engineering professional community.
It's easy to join and it's free.

Here's Why Members Love Eng-Tips Forums:

Register now while it's still free!

Already a member? Close this window and log in.

Join Us             Close