!!!FEA NL2: CONTACT PROBLEM OF A BALL AND A DISK

!!Initial setup

! White canvas
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,100,100,100,0
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,0,0,0,15
! Define constant parameters
! Units: [N],[Nmm],[N/mm^2=MPa],[N/mm]
ELAST = 210E3
NU = 0.3
F = 3000
ELEMNUM=6

!!SETUP: element types, materials, layup

! Axysimmetric plane element
KEYOPTKEYOPT, ITYPE, KNUM, VALUE Sets element key options.,1,3,1
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,2,TARGE169
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,3,CONTA175
! Linear, isotropic, elastic material model
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EX,1,,ELAST
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,1,,NU

!!Modeling

! glue areas
AGLUEAGLUE, NA1, NA2, NA3, NA4, NA5, NA6, NA7, NA8, NA9 Generates new areas by "gluing" areas.,1,2
!
KK, NPT, X, Y, Z Defines a keypoint.,5,0,10,0
KK, NPT, X, Y, Z Defines a keypoint.,8,10,10,0
KK, NPT, X, Y, Z Defines a keypoint.,9,0,0,0
LARCLARC, P1, P2, PC, RAD Defines a circular arc.,8,9,5,10
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,9,5
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,5,8
! divide the arc (L5)
LDIVLDIV, NL1, RATIO, PDIV, NDIV, KEEP Divides a single line into two or more lines.,5,2/3
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,5,10
!
ALAL, L1, L2, L3, L4, L5, L6, L7, L8, L9, L10 Generates an area bounded by previously defined lines.,7,11,12
ALAL, L1, L2, L3, L4, L5, L6, L7, L8, L9, L10 Generates an area bounded by previously defined lines.,12,5,8

!!Meshing

! Set the number of elements along lines
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,11,,,2*ELEMNUM
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,5,,,2*ELEMNUM
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,7,,,2*ELEMNUM
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,8,,,2*ELEMNUM
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,12,,,2*ELEMNUM
! horizontal line
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,1,,,ELEMNUMELEMNUM = 6
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,3,,,ELEMNUMELEMNUM = 6
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,9,,,2*ELEMNUM
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,10,,,2*ELEMNUM
! vertical line
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,4,,,ELEMNUMELEMNUM = 6
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,2,,,ELEMNUMELEMNUM = 6
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,6,,,ELEMNUMELEMNUM = 6
! mesh areas using mapped algorithm
MSHKEYMSHKEY, KEY Specifies whether free meshing or mapped meshing should be used to mesh a model.,1
AMESHAMESH, NA1, NA2, NINC Generates nodes and area elements within areas.,ALL
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
! Add Contact and Target elements of surface lines
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,,,11
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
TYPETYPE, ITYPE Sets the element type attribute pointer.,3
ESURFESURF, XNODE, Tlab, Shape Generates elements overlaid on the free faces of selected nodes.,ALL
!
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,,,3
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
TYPETYPE, ITYPE Sets the element type attribute pointer.,2
ESURFESURF, XNODE, Tlab, Shape Generates elements overlaid on the free faces of selected nodes.,ALL
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
GPLOTGPLOT Controls general plotting.! Click here to view image

!!Apply kinematic and dynamic loads on structure.

! kinematic constraints
DKDK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6 Defines DOF constraints at keypoints.,6,UY,0
! impose identical UY displacement on line 8 through constraint
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,,,8
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
CPCP, NSET, Lab, NODE1, NODE2, NODE3, NODE4, NODE5, NODE6, NODE7, NODE8, NODE9, NODE10, NODE11, NODE12, NODE13, NODE14, NODE15, NODE16, NODE17 Defines (or modifies) a set of coupled degrees of freedom.,1,UY,ALL
!*
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL

!!Define loadcases

! Loadcase 1 displacement
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,,,8
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UY,-0.005
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
FINISHFINISH Exits normally from a processor.
! Solve loadcase 1
/SOL/SOLU Enters the solution processor.
SOLVESOLVE, Action Starts a solution.
FINISHFINISH Exits normally from a processor.
! Loadcase 2
/PREP7/PREP7 Enters the model creation preprocessor.
! Delete displacement constraint
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,,,8
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
DDELEDDELE, NODE, Lab, NEND, NINC, Rkey Deletes degree-of-freedom constraints.,ALL,ALL
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
! Apply force on keypoint
FKFK, KPOI, Lab, VALUE, VALUE2 Defines force loads at keypoints.,5,FY,-F
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
! Set solution options
NLGEOMNLGEOM, Key Includes large-deflection effects in a static or full transient analysis.,1
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,2,10,1
SOLVESOLVE, Action Starts a solution.

!!Postprocessing

! Loadcase 2
! Displacement fields
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,X,0,1.0! Click here to view image
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,Y,0,1.0! Click here to view image
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,SUM,0,1.0! Click here to view image
! Stress fields
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,S,Y,0,1.0! Click here to view image
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,S,EQV,0,1.0! Click here to view image
! Stress distribution along the symmetry axis
PATHPATH, NAME, nPts, nSets, nDiv Defines a path name and establishes parameters for the path.,SYMM,2,30,100
PPATHPPATH, POINT, NODE, X, Y, Z, CS Defines a path by picking or defining nodes, or locations on the currently active working plane, or by entering specific coordinate locations.,1,NODE(0,-5,0)
PPATHPPATH, POINT, NODE, X, Y, Z, CS Defines a path by picking or defining nodes, or locations on the currently active working plane, or by entering specific coordinate locations.,2,NODE(0,10,0)
PDEFPDEF, Lab, Item, Comp, Avglab Interpolates an item onto a path.,SIGY,S,Y,NOAVG
!
/AXLAB/AXLAB, Axis, Lab Labels the X and Y axes on graph displays.,X,'y [mm]'
/AXLAB/AXLAB, Axis, Lab Labels the X and Y axes on graph displays.,Y,'sigY [MPa]'
PLPATHPLPATH, Lab1, Lab2, Lab3, Lab4, Lab5, Lab6 Displays path items on a graph.,SIGY! Click here to view image

! This webpage was generated with APDL highlighter (version: 0.1-beta) on: 09.05.2018 11:16:37, from APDL script: script.txt