How to simulate contact stress using ANSYS student 19.1 APDL

  • 48 Views
  • Last Post 2 weeks ago
Nuno koh posted this 2 weeks ago

I'm beginner of ANSYS APDL. I'm interested in contact stress.

 

 

I have been trying  for it with the text  "ENGINEERING ANALYSIS with ANSYS SOFTWARE (Author: Y. Nakasone and S. Yoshimoto)". 

The Contact stress problem is described in 120page ~138page in the text.

I have followed  the text,  But I can't get the solution and  I have below error in solution phase. 

 

I have added the capture image and  code source for the contact stress simulation procedure.

 

Thanks. 

Attached Files

Order By: Standard | Newest | Votes
Nuno koh posted this 2 weeks ago

I'm beginner of ANSYS APDL. I'm interested in contact stress.

 

I have been trying  for it with the text  "ENGINEERING ANALYSIS with ANSYS SOFTWARE (Author: Y. Nakasone and S. Yoshimoto)". 

The Contact stress problem is described in 120page ~138page in the text.

I have followed  the text,  But I can't get the solution and  I have below error in solution phase. 

// Error //

"There are no contact elements in the database. Please send the data leading to this operation to your technical support provider, ..." 

 

What is wrong ?  Is there a problem related to making contact element due to ANSYS student version?

I have added the capture image and  code source for the contact stress simulation procedure.

 

Thanks. 

Nuno koh posted this 2 weeks ago

[input file]

 

/PREP7 

BLC4,0,0,500,-500  

CYL4,0,500,500,-90,500,0

SAVE

SAVE

FLST,2,2,5,ORDE,2  

FITEM,2,1  

FITEM,2,-2 

AGLUE,P51X 

!* 

MPTEMP,,,,,,,, 

MPTEMP,1,0 

MPDATA,EX,1,,2.1E5 

MPDATA,PRXY,1,,0.3 

!* 

ET,1,PLANE183  

!* 

KEYOPT,1,1,0

KEYOPT,1,3,2

KEYOPT,1,6,0

!* 

ET,2,TARGE169  

!* 

ET,3,CONTA175 

!* 

FLST,5,2,4,ORDE,2  

FITEM,5,2  

FITEM,5,-3 

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

!* 

LESIZE,_Y1, , ,60,1/10, , , ,1 

!* 

FLST,5,2,4,ORDE,2  

FITEM,5,1  

FITEM,5,4  

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

!* 

LESIZE,_Y1, , ,60,10, , , ,1

!* 

FLST,5,2,4,ORDE,2  

FITEM,5,8  

FITEM,5,-9 

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

!* 

LESIZE,_Y1, , ,40,10, , , ,1

!* 

FLST,5,1,4,ORDE,1  

FITEM,5,6  

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

!* 

LESIZE,_Y1, , ,40,1/10, , , ,1 

!* 

MSHKEY,0

FLST,5,2,5,ORDE,2  

FITEM,5,1  

FITEM,5,3  

CM,_Y,AREA 

ASEL, , , ,P51X

CM,_Y1,AREA

CHKMSH,'AREA'  

CMSEL,S,_Y 

!* 

AMESH,_Y1  

!* 

CMDELE,_Y  

CMDELE,_Y1 

CMDELE,_Y2

 

!!!!!!!!!!  Cylinder btm surface line selection/!!!!!!!!!!

LSEL,S, , ,       8

!!!!!!!!!!  /Cylinder btm surface line selection!!!!!!!!!!

!!!!!!!!!!  Cylinder btm surface node selection/!!!!!!!!!!

NSLL,R,1

!!!!!!!!!!  / Cylinder btm surface node selection!!!!!!!!!!

NPLOT 

 

 !!!!!!!!!!  Cylinder Btm surface contact point  near area selection/!!!!!!!!!!

FLST,5,41,1,ORDE,3 

FITEM,5,242

FITEM,5,11043  

FITEM,5,-11082 

NSEL,R, , ,P51X

!!!!!!!!!!  /Cylinder Btm surface contact point  near area selection!!!!!!!!!!

 

!!!!!!!!!!  Contact element define of cylinder btm surface/ !!!!!!!!!!!!!

TYPE,   3  

MAT,       1

REAL,  

ESYS,       0  

SECNUM,

TSHAP,LINE

!!!!!!!!!!  /Contact element define of cylinder btm surface !!!!!!!!!!!!!

 

!!!!!!!!!!  Contact element attach on cylinder btm surface / !!!!!!!!!!!!!

FLST,5,41,1,ORDE,3 

FITEM,5,242

FITEM,5,11043  

FITEM,5,-11082 

CM,_Y,NODE 

NSEL, , , ,P51X

CM,_Y1,NODE

CMSEL,S,_Y 

!* 

CMSEL,,_Y1 

ESURF, ,BOT,

CMSEL,,_Y  

CMDELE,_Y  

CMDELE,_Y1

!!!!!!!!!!  /Contact element attach on cylinder btm surface !!!!!!!!!!!!!

 

ALLSEL,ALL 

!!!!!!!!!!  Flat plate Top surface line selection/!!!!!!!!!!

LSEL,S, , ,       3

!!!!!!!!!!  /Flat plate Top surface line selection!!!!!!!!!!

 

!!!!!!!!!!  Flat plate Top surface node selection/!!!!!!!!!!

NSLL,R,1

!!!!!!!!!!  /Flat plate Top surface node selection!!!!!!!!!!

NPLOT

 

!!!!!!!!!!  Flat plate Top surface contact point  near area selection/ !!!!!!!!!!

FLST,5,61,1,ORDE,3 

FITEM,5,242

FITEM,5,302

FITEM,5,-361

NSEL,R, , ,P51X

!!!!!!!!!!  /Flat plate Top surface contact point  near area selection!!!!!!!!!!

 

 !!!!!!!!!!  Target element define on Flat plate Top surface/ !!!!!!!!!!!!!

TYPE,   2  

MAT,       1

REAL,  

ESYS,       0  

SECNUM,

TSHAP,LINE

!!!!!!!!!!  /Target element define on Flat plate Top surface !!!!!!!!!!!!!

 

!!!!!!!!!!  Target element attach on Flat plate Top surface/ !!!!!!!!!!!!!

FLST,5,61,1,ORDE,3 

FITEM,5,242

FITEM,5,302

FITEM,5,-361

CM,_Y,NODE 

NSEL, , , ,P51X

CM,_Y1,NODE

CMSEL,S,_Y 

!* 

CMSEL,,_Y1 

ESURF, ,TOP,

CMSEL,,_Y  

CMDELE,_Y  

CMDELE,_Y1 

!!!!!!!!!!  /Target element attach on Flat plate Top surface !!!!!!!!!!!!!

 

!!!!!!!!!!/boundary condition!!!!!!!!!!!

/SOL

FLST,2,2,4,ORDE,2  

FITEM,2,4  

FITEM,2,9  

!* 

/GO

DL,P51X, ,UX,  

FLST,2,1,4,ORDE,1  

FITEM,2,1  

!* 

/GO

DL,P51X, ,UY,  

FLST,2,1,4,ORDE,1  

FITEM,2,6  

/GO

!* 

SFL,P51X,PRES,1,

!!!!!!!!!!/boundary condition!!!!!!!!!!!

 

 

/SOL

/STATUS,SOLU

SOLVE

Nuno koh posted this 2 weeks ago

 

jpasquerell posted this 2 weeks ago

See ANSYS Verification Manual sample model VM211 which is a contact analysis with similar geometry for an example set of commands.

jpasquerell posted this 2 weeks ago

I looked over your input and do see two items that need addressed.  1) There are two lines with real, that should have a value after the comma.  Try changing both to real,1.  2) The other issue is that  everything needs selected prior to the solve command so add an allsel command just before the solve.

  • Liked by
  • SandeepMedikonda
  • Nuno koh
Nuno koh posted this 2 weeks ago

 Thanks jpasquerell !!

Thanks to your comment, I could solve my 1'st problem. (Real, 1  and ALLSEL ALL )

But, I have still encountered another 2 problems.

Problem 1. : timestep problem

Problem 2. : penetration problem

 

[Solution from Text]  No penetration

 

What is wrong ?

 

 

 

[Input file] 

 

/PREP7 

BLC4,0,0,500,-500  

CYL4,0,500,500,-90,500,0

FLST,2,2,5,ORDE,2  

FITEM,2,1  

FITEM,2,-2 

AGLUE,P51X 

MPTEMP,,,,,,,, 

MPTEMP,1,0 

MPDATA,EX,1,,2.1E5 

MPDATA,PRXY,1,,0.3 

ET,1,PLANE82 

KEYOPT,1,1,0

KEYOPT,1,3,2

KEYOPT,1,6,0

ET,2,TARGE169  

ET,3,CONTA172                !!!!!!!!!!!   contact element type !!!!!!!!!!!!

KEYOPT,3,1,0

KEYOPT,3,2,0

KEYOPT,3,3,0

KEYOPT,3,4,0

KEYOPT,3,5,1

KEYOPT,3,6,0

KEYOPT,3,7,0

KEYOPT,3,8,0

KEYOPT,3,9,0

KEYOPT,3,10,0  

KEYOPT,3,11,0  

KEYOPT,3,12,0  

KEYOPT,3,14,0  

KEYOPT,3,15,0  

KEYOPT,3,18,0

R,1,,,4, ,,

FLST,5,2,4,ORDE,2  

FITEM,5,2  

FITEM,5,-3 

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

LESIZE,_Y1, , ,60,1/10, , , ,1 

FLST,5,2,4,ORDE,2  

FITEM,5,1  

FITEM,5,4  

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

LESIZE,_Y1, , ,60,10, , , ,1

FLST,5,2,4,ORDE,2  

FITEM,5,8  

FITEM,5,-9 

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

LESIZE,_Y1, , ,40,10, , , ,1

FLST,5,1,4,ORDE,1  

FITEM,5,6  

CM,_Y,LINE 

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y  

LESIZE,_Y1, , ,40,1/10, , , ,1 

MSHKEY,0

FLST,5,2,5,ORDE,2  

FITEM,5,1  

FITEM,5,3  

CM,_Y,AREA 

ASEL, , , ,P51X

CM,_Y1,AREA

CHKMSH,'AREA'  

CMSEL,S,_Y 

AMESH,_Y1  

CMDELE,_Y  

CMDELE,_Y1 

CMDELE,_Y2

 

!!!!!!!!!!  Cylinder btm surface line selection/!!!!!!!!!!

LSEL,S, , ,       8

!!!!!!!!!!  /Cylinder btm surface line selection!!!!!!!!!!

!!!!!!!!!!  Cylinder btm surface node selection/!!!!!!!!!!

NSLL,R,1

!!!!!!!!!!  / Cylinder btm surface node selection!!!!!!!!!!

NPLOT 

 

!!!!!!!!!!  Cylinder Btm surface contact point  near area selection/!!!!!!!!!!

FLST,5,41,1,ORDE,3 

FITEM,5,242

FITEM,5,11043  

FITEM,5,-11082 

NSEL,R, , ,P51X

!!!!!!!!!!  /Cylinder Btm surface contact point  near area selection!!!!!!!!!!

 

!!!!!!!!!!  Contact element define of cylinder btm surface/ !!!!!!!!!!!!!

TYPE,   3  

MAT,       1

REAL,   1

ESYS,       0  

SECNUM,

TSHAP,LINE

!!!!!!!!!!  /Contact element define of cylinder btm surface !!!!!!!!!!!!!

 

!!!!!!!!!!  Contact element attach on cylinder btm surface / !!!!!!!!!!!!!

FLST,5,41,1,ORDE,3 

FITEM,5,242

FITEM,5,11043  

FITEM,5,-11082 

CM,_Y,NODE 

NSEL, , , ,P51X

CM,_Y1,NODE

CMSEL,S,_Y 

CMSEL,,_Y1 

ESURF, ,BOT,

CMSEL,,_Y  

CMDELE,_Y  

CMDELE,_Y1

!!!!!!!!!!  /Contact element attach on cylinder btm surface !!!!!!!!!!!!!

 

ALLSEL,ALL 

!!!!!!!!!!  Flat plate Top surface line selection/!!!!!!!!!!

LSEL,S, , ,       3

!!!!!!!!!!  /Flat plate Top surface line selection!!!!!!!!!!

 

!!!!!!!!!!  Flat plate Top surface node selection/!!!!!!!!!!

NSLL,R,1

!!!!!!!!!!  /Flat plate Top surface node selection!!!!!!!!!!

NPLOT

 

!!!!!!!!!!  Flat plate Top surface contact point  near area selection/ !!!!!!!!!!

FLST,5,61,1,ORDE,3 

FITEM,5,242

FITEM,5,302

FITEM,5,-361

NSEL,R, , ,P51X

!!!!!!!!!!  /Flat plate Top surface contact point  near area selection!!!!!!!!!!

 

 !!!!!!!!!!  Target element define on Flat plate Top surface/ !!!!!!!!!!!!!

TYPE,   2  

MAT,       1

REAL,   1

ESYS,       0  

SECNUM,

TSHAP,LINE

!!!!!!!!!!  /Target element define on Flat plate Top surface !!!!!!!!!!!!!

 

!!!!!!!!!!  Target element attach on Flat plate Top surface/ !!!!!!!!!!!!!

FLST,5,61,1,ORDE,3 

FITEM,5,242

FITEM,5,302

FITEM,5,-361

CM,_Y,NODE 

NSEL, , , ,P51X

CM,_Y1,NODE

CMSEL,S,_Y 

CMSEL,,_Y1 

ESURF, ,TOP,

CMSEL,,_Y  

CMDELE,_Y  

CMDELE,_Y1 

!!!!!!!!!!  /Target element attach on Flat plate Top surface !!!!!!!!!!!!!

ALLSEL,ALL

 

!!!!!!!!!!/boundary condition!!!!!!!!!!!

/SOL

FLST,2,2,4,ORDE,2  

FITEM,2,4

FITEM,2,6 

/GO

DL,P51X, ,UX,  

FLST,2,1,4,ORDE,1  

FITEM,2,1  

/GO

DL,P51X, ,UY,  

FLST,2,1,4,ORDE,1  

FITEM,2,6  

/GO

SFL,P51X,PRES,1,

!!!!!!!!!!/boundary condition!!!!!!!!!!!

ALLSEL,ALL

 

 

/SOL

/STATUS,SOLU

SOLVE

jj77 posted this 2 weeks ago

As someone mentioned, there is a similar example in the verification manual (VM191). Have a look there.

 

Also as you will see there, first impose a small enforced displacement to active the contacts, and then remove it as you apply the pressure. Otherwise if you apply a pressure it can be difficult to catch the contact. Finally you can do it in workbench which is easier, and then export the .dat file for reading into the APDL. Attached is the apdl code from WB. 

Attached Files

Close