!Apdl-Code zum Mappen !Knoten des Netzes ermitteln *get,mpmax,node,,num,max !Knotenanzahl *dim,nadata,array,mpmax,3 !Def. von Array ndata *vget,nadata(1,1),node,0,loc,x !laden der x koordinaten *vget,nadata(1,2),node,0,loc,y !laden der y koordinaten *vget,nadata(1,3),node,0,loc,z !laden der z koordinaten *dim,yielstr,array,mpmax,1 !Def. von Array yielstr !Einlesen der externen Koordinaten defanzahl=8278 *dim,flgr,array,defanzahl,3 !Def. von Array flgr *vread,flgr(1,1),koord,txt,,,,defanzahl,1 !laden der x koordinaten aus Datei (F8.6) *vread,flgr(1,2),koord,txt,,,,defanzahl,1 !laden der y koordinaten aus Datei (9xF8.6) *vread,flgr(1,3),koord,txt,,,,defanzahl,1 !laden der z koordinaten aus Datei (18xF8.6) !Einlesen der externen Fliessgrenzen *dim,yiel,array,mpmax,1 !Def. von Array yiel *vread,yiel(1,1),fliessgrenze,txt,,,,defanzahl,1 !Laden der Fließgrenzwerte aus Datei (27xF8.6) !Mappen der Daten von flgr auf yielstr *moper,yielstr(1,1),nadata(1,1),map,yiel(1,1),flgr(1,1),3 !Mappen der Daten von flgr auf ndata *get,enum,elem,,num,max !Elementgesamtanzahl !DO-SCHLEIFE !von Knoten auf Elemente Schliessen *DO,J,1,enum,1 !DO-Schleife,beginnt bei J=1,elemno-fach,J+1/Durchgang *dim,elemno,array,8,1 !Element-Knoten-Zusammenhang elemno = knotennummer der Knoten !Knotennummer von ELement J !Material J Definieren a=yielstr(1,elemno(1,1)) !Fliessgrenze in jedem Knoten aus mapping b=yielstr(1,elemno(1,2)) c=yielstr(1,elemno(1,3)) d=yielstr(1,elemno(1,4)) e=yielstr(1,elemno(1,5)) f=yielstr(1,elemno(1,6)) g=yielstr(1,elemno(1,7)) h=yielstr(1,elemno(1,8)) m=(a+b+c+d+e+f+g+h)/8 !gemittelte Fliessgrenze von Material J in ELement J MPTEMP,J,20 !Temperatur von Material J MP,EX,J,210000 !E-Modul von Material J MP,PRXY,J,0.3 !Poisson-Zahl von Material J !Material mit Bilinearen isotropen Verfestigung MP,tangmod,J,2954.2708 !Tangentenmodul von Material J MP,yldstrs,J,m !Streckgrenze/Fließgrenze von Material J !Materialzeiger zu Element zuweisen MAT,1,,0,, !Materialdatenzeiger auf Element J *Enddo