Displacement-Limit erhöhen (FEM / Genormte Berechnungen/ANSYS)

Hallo,
ich versuche (immernoch) für meine Diplomarbeit den Eindringvorgang eines konischen Eindrinkörpers in Stahl zu simulieren. Dazu sollen insgesant 15000 N auf den Eindringkörper lasten. In meiner version komme ich nach ein paar loadsteps nur bis auf 250 N und erhalte dann die Nachricht die Displacements eines Knotens sind größer als das Limit. Bei der Simulation erwarte und will ich ja gerade große Displacements (Es soll sich eine Aufwerfung um den Eindringpunkt bilden.
Kann ich dieses Limit irgendwie erhöhen oder die Fehlermeldung umgehen oder stimmt sonst was nicht an der Simulation. Wäre nett wenn mal jemand kurz drüber gucken würde. Mein Abgabetermin ist schon ganz schön nah und ich bin mit meinen Ideen am Ende.

Gruß,

Mitch

!Simulation der Indentation einer planen Oberfläche durch einen konischen Indenter
!mit abgeflachter Spitze

!Überprüfung der Berechnung durch gleichzeitige Simulation eines Zugversuchs

!Alle Werte sind in [N] und [m] angegeben

reib1=0.1 !Reibungskoeffizient für Material1 (S235JR)

hoz=20 !Kantenlänge Probe [20mm]
hoh=2.76165 !Höhe des konischen Indenters [2,76mm]

fkp=0.0005 !Elementgröße um Eindringpunkt, fein
fdge=3 !Elementgröße in der Probe, grob

!/CONFIG, NRES, 10000000

/PREP7 !Preprocessor

!---Elementtypen---

ET,1, plane182,,,2 ! Plane strain (Z strain = 0.0)
REAL, 1 ! Sets the element real constant set attribute pointer.
KEYOPT,1,1,2 ! Enhanced strain formulation, suggestion by ansys

ET,2,CONTA171 !,,,1,,,,,,1,,,2 !Contactelement

KEYOPT, 2, 2,1 !penalty method, requires fkn fks ftoln slto

KEYOPT, 2, 3,0 !no superelements used (or not ???)

KEYOPT, 2, 4,0 !contact detection on gauss points

KEYOPT, 2, 5,3 !close gap/reduce penetration with auto CNOF

KEYOPT, 2, 7,3 ! Change in contact predictions made to achieve the !minimum
!time/load increment whenever a change
!incontact status occurs

KEYOPT, 2, 9,2 ! Include both initial geometrical penetration or gap and offset,
!but with ramped effects

KEYOPT, 2, 10,1 !automatic contact stiffness update

KEYOPT, 2, 12,1 !rough

REAL, 3

R, 3, , , 0.8, ,

ET,3,TARGE169 !Targetelement
REAL, 3

!---Materialkennwert Stahl---

MP, EX,1,210000 !E-Modul
MP,PRXY,1,0.3 !Querkontraktionszahl
MP, MU,1,reib1, !Reibungskoeffizient, vorläufig

!Spannungs-DehnungsDiagramm, vorläufig

TB,MISO,1,1,5 ! Activate a data table
TBTEMP,20 ! Temperature = 20
TBPT,DEFI,1.12E-3, 235 ! Strain, stress at temperature = 20
TBPT,DEFI,5E-2, 260
TBPT,DEFI,10E-2, 310
TBPT,DEFI,15E-2, 350
TBPT,DEFI,17.5E-2, 360
TBPLOT,MISO,1

!---Geometrie----

!---Indenter---
K, 1, 0, 0,
K, 2, 0, hoh,
K, 3, 2.5172989, hoh,
K, 4, 0.2, 0,

LSTR, 1, 4
LSTR, 4, 3
LSTR, 3, 2
LSTR, 2, 1

LFILLT, 1, 2, 0.0001E3,

AL, 1, 2, 3, 4, 5

!---Probe---
rectng,0,hoz,0,-hoz

!----Vernetzung------Probe---

esize, fdge !element size
kesize,10,fkp !element size around first kontakt-node
type,1
real,1
mat,1
amesh,2

!----Randbedingung-----

!----Randbedingung-----Probe
nsel,s,loc,y,-hoz
d,all,uy,0
nsel,s,loc,x,0
d,all,ux,0 !gilt auch für Pilotnode=pilot1

!---Kontakte----

!--------------------TARGETELEMENT
type,3
real,3
TSHAP, LINE
lmesh, 1,2

type,3
real,3
TSHAP, ARC
lmesh, 5

*GET,NMAX,NODE,,NUM,max
Pilot1= NMAX+1
n, Pilot1,, 0.25
TSHAP,PILOT
e,Pilot1
d,pilot1,ROTZ,0
d,pilot1,ux,0

!-------------------DOF-CHECK----
!/PNUM, NODE, 1 !Node numbers on node and element plots.
!nsel,s,node,,all !SELECTS ALL NODES
!DOFSEL, s, disp !SELECTS DOFS
!dlist,all !LIST ALL DOFS

!-------------------CONTACTELEMENT
lsel,s,,,8
nsll,s,1
type,2
real,3
esurf

CNCHECK, detail, 3,

/SOLU !Solution
! Load Step 1:
NROPT,UNSYM
SOLCON,ON,ON
CUTCONTROL, PLSLIMIT, 0.15,
EQSLV,FRONT !Frontal direct equation solver

nsel,s,node,,pilot1
d,pilot1,ux,0
F, pilot1, fy, -50
NSUBST, 40,40,40 ! Load step options

KBC, 0
OUTRES,all, LAST
OUTPR,all, LAST
LSWRITE ! Writes load step file: Jobname.S01

CSYS
ALLSEL, ,
SOLVE

/SOLU !Solution
! Load Step 2:
NROPT,UNSYM
SOLCON,ON
CUTCONTROL, PLSLIMIT, 0.15,
EQSLV,FRONT !Frontal direct equation solver
nsel,s,node,,pilot1
d,pilot1,ux,0
F, pilot1, fy, -100
NSUBST, 40,40,40 ! Load step options
KBC, 0
OUTRES,all, LAST
OUTPR,all, LAST
LSWRITE ! Writes load step file: Jobname.S02
CSYS
ALLSEL, ,
SOLVE

/SOLU !Solution

! Load Step 3:
NROPT,UNSYM
SOLCON,ON,ON
CUTCONTROL, PLSLIMIT, 0.15,
EQSLV,FRONT !Frontal direct equation solver
nsel,s,node,,pilot1
d,pilot1,ux,0
F, pilot1, fy, -150
NSUBST, 300,300,300 ! Load step options
KBC, 0 !ramped loads
OUTRES,all, LAST
OUTPR,all, LAST
LSWRITE
CSYS
ALLSEL, ,
SOLVE

/SOLU !Solution

! Load Step 4:
NROPT,UNSYM
SOLCON,ON,ON
CUTCONTROL, PLSLIMIT, 0.15,
EQSLV,FRONT !Frontal direct equation solver
nsel,s,node,,pilot1
d,pilot1,ux,0
F, pilot1, fy, -200
NSUBST, 300,300,300 ! Load step options
KBC, 0 !ramped loads
OUTRES,all, LAST
OUTPR,all, LAST
LSWRITE
CSYS
ALLSEL, ,
SOLVE

/SOLU !Solution

! Load Step 5:
NROPT,UNSYM
SOLCON,ON,ON
CUTCONTROL, PLSLIMIT, 0.15,
EQSLV,FRONT !Frontal direct equation solver
nsel,s,node,,pilot1
d,pilot1,ux,0
F, pilot1, fy, -250
NSUBST, 700,700,700 ! Load step options
KBC, 0 !ramped loads
OUTRES,all, LAST
OUTPR,all, LAST
LSWRITE
CSYS
ALLSEL, ,
SOLVE

/SOLU !Solution

! Load Step 6:
NROPT,UNSYM
SOLCON,ON,ON
CUTCONTROL, PLSLIMIT, 0.15,
EQSLV,FRONT !Frontal direct equation solver
nsel,s,node,,pilot1
d,pilot1,ux,0
F, pilot1, fy, -500
NSUBST, 10000, 10000, 5000 ! Load step options
KBC, 0 !ramped loads
OUTRES,all, LAST
OUTPR,all, LAST
LSWRITE
CSYS
ALLSEL, ,
SOLVE

/POST1 !Postprocessor
/EDGE,,1
/DSCA,,1
PLNSOL,S,EQV !plnsol,s equivalent stress (wtf),nodal sol v mise

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP