---

Hallo,

ich versuche grad eine Kontaktfläche in Abhängigkeit der entstehenden Verformungen / Bewegungen anzupassen.
Dies bedeutet, dass der Reibkoeffizient in Abhängigkeit der Verformungen variiert werden soll.

Nachfolgend die Erstellung des Kontaktbereichs:

!----- Elementerstellung
et,1,170
et,2,174
keyopt,2,2,1
keyopt,2,10,5
keyopt,2,9,1
keyopt,2,12,1

!---- Real-Constants
r,1,,,0.1,0.1,0,-4
!---- Material
mptemp,,,,,,,,
mptemp,1,0
mpdata,mu,1,,0.5        !Reibkoeffizient 0,5
!---- Kontakt_erstellen
!----- Target
cmsel,s,CONTACT_TARGET
type,1
real,1
mat,1
esurf
alls
!----- Surface
cmsel,s,CONTACT_SURFACE
type,2
real,1
mat,1
esurf
alls

Im nächsten Schritt veruche ich vereinfachend den Reibkoeffizient elementweise zu zuweisen (hier vereinfachend auf den Wert 1.0 zu setzen):

alls
esel,s,type,,1,,
esel,a,type,,2,,
cm,AKT_ELEM,elem
*get,ELEMNUM,elem,0,count,,,            !Ermittlung der Elementanzahl
*get,STARTELEM,elem,0,num,min            !Startelement bestimmen
*get,ENDELEM,elem,0,num,max              !letztes Element bestimmen
ELEMNUM=ENDELEM-STARTELEM
alls
!--- Materialzuweisung
*del,E_MU,,nopr
*dim,E_MU,array,ENDELEM,1
/prep7
*do,i,STARTELEM,ENDELEM
  E_MU(i)=1.0
  mp,mu,i,E_MU(i)
*enddo

Nun funktioniert das leider scheinbar nicht.
Kann mir jemand helfen und sieht den Fehler bzw. hat vielleicht selbst ein ähnliches Vorhaben erfolgreich umgesetzt.

Besten Dank im Voraus für die Unterstützung.

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

---

Du musst noch den einzelnen Elementen die entsprechende Materialnummer zuweisen, bisher hast Du nur ganz viele Materialien definiert.

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

---

Hallo ChristophN,

erst einmal besten Dank für die Antwort.

Ich habe versucht Deinen Ratschlag zu befolgen. Dementsprechend habe ich den Befehl "emodif,i,mat,i" eingefügt. Jetzt ist nur das Problem, dass der Reibkoeffizient scheinbar nicht in den einzelnen Iterationen elementweise verwendet wird, da die entstehenden Bewegungen immer gleich sind.

alls
esel,s,type,,1,,
esel,a,type,,2,,
cm,AKT_ELEM,elem
*get,ELEMNUM,elem,0,count,,,            !Ermittlung der Elementanzahl
*get,STARTELEM,elem,0,num,min            !Startelement bestimmen
*get,ENDELEM,elem,0,num,max              !letztes Element bestimmen
ELEMNUM=ENDELEM-STARTELEM
alls
!--- Materialzuweisung
*del,E_MU_0,,nopr
*dim,E_MU_0,array,ENDELEM,1
/prep7
*do,i,STARTELEM,ENDELEM
  E_MU(i)=1.0
  mp,mu,i,E_MU_0(i)
  emodif,i,mat,i
*enddo

Hier die Schleife für die iterationsweise Anpassung:
T_PREV=IT-1
*do,i,STARTELEM,ENDELEM
esel,s,type,,1,,
esel,a,type,,2,,
cm,AKT_ELEM,elem
*if,SYS_E_MMOT_%IT%(i),gt,0.020,then
  E_MU_%IT%(i)=E_MU_%IT_PREV%(i)-0.1
  mp,mu,i,E_MU_%IT%(i)
*elseif,SYS_E_MMOT_%IT%(i),le,0.020
  E_MU_%IT%(i)=E_MU_%IT_PREV%(i)+0.1
  mp,mu,i,E_MU_%IT%(i)
*endif
  emodif,i,mat,i
*enddo

Das Array "SYS_E_MMOT_%IT%" handelt es sich um ein aus den Bewegungen im Kontaktbereich erstelltes Array.

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