Hallo Leute !
Hab mir hier ein Linienmodell einer Feder erstellt.
Nun möchte ich das gern mit rechteckigen querschnitt als Volumenmodell darstellen, kann mir einer helfen!?
INPUT:
finish
/clear,nostart
/filename,input_Beleg_Schraubenfeder_neu
/title,Schraubenfeder
/eshape,1
/pnum,line,1
!ABFRAGE
!===============================================================
!Parameter
*ask,KRAFT,WIE IST DIE KRAFT,1
*ask,Rad_fed,Wie ist der Radius der Feder?,10
*ask,r_draht,Wie ist der Radius des Drahtes?,1
*ask,steig,Wie ist die Steigung der Feder?,4
*ask,wickl,Wie ist die Wicklungsanzahl?,10
*ask,ele,Wie viele Elemente sollen verwendet werden?,200
ele2=(3/4)*(ele/wickl)
ro=7.81e-9 !Dichte roh
emod=205000 !E-modul
nxy=0.28 !Poisson-Zahl
!ELEMENTS CONSTANTS NODES
!===============================================================
/prep7 !PREPROZESSOR
et,1,pipe16
r,1,r_draht,0.5*r_draht !real constant mit Außendurchmesser,Wandstärke(0,5*d==>Vollkörper)
et,3,mass21 !Elementtyp 3 definieren
Keyopt,3,1,0
keyopt,3,2,0
keyopt,3,3,0 !REAL 3 definieren
R,3,0,0,0,0,0,0
mp,ex,1,emod !Material Properties: Elastizitätsmodul
mp,nuxy,1,nxy !Material Properties: Poisson-Zahl
mp,dens,1,ro !Material Properties: Dichte
! Kontakt Element Definition
*do,p,1,400
et,4,CONTA176,,,0 !3D-Point-Point-Contact
r,4,r_draht/2,r_draht/2,1 !REAL 4
keyopt,4,10,2 !KEYOPT 4
et,5,targe170 !TARGE170
r,4+p,0,r_draht/2,1 !REAL CONSTANT
*enddo
! Knoten Massepunkt
n,1000,0,0,0 !Master NODE
!FEDER
!===============================================================
! FEDER IM STEIGUNGSBEREICH
csys,1 !Zylindrisches Koordinatensystem mit z als rotationsachse
n,1,Rad_fed,0 !Knotendefinition für 1.Knoten
*repeat,ele,1,0,wickl*360/ele,steig*wickl/ele !Knotengenerierung für Feder im Bereich mit Steigung
type,1 !Def: TYP
real,1 !Def: REAL
mat,1 !Def: MAT
e,1,2 !Elementenerzeugung für 1.Element zw. Knoten 1 u. 2
*repeat,ele-1,1,1 !Eementengenerierung für Feder im Bereich mit Steigung
! FEDERENDE 1
n,ele+1,rad_fed,0 !Knotendefinition für 1.Knoten des 1.Federendes
*repeat,ele2,1,0,-270/ele2,0*wickl/ele2 !Knotengenerierung für Feder im Bereich ohne Steigung
type,1 !Def: TYP
real,1 !Def: REAL
mat,1 !Def: MAT
e,ele+1,ele+2 !Elementenerzeugung für 1.Element des 1.Federendes
*repeat,ele2-1,1,1 !Eementengenerierung für Feder im Bereich ohne Steigung
! FEDERENDE 2
n,ele+ele2+1,Rad_fed,0,(steig*(wickl*((ele-1)/ele))) !Knotenerzeugung für nächsten Knoten mit Steigung 0
*repeat,ele2,1,0,270/ele2,0*wickl/ele2 !Knotengenerierung für Feder im Bereich ohne Steigung
type,1 !Def: TYP
real,1 !Def: REAL
mat,1 !Def: MAT
e,ele,ele+ele2+1 !Elementenerzeugung für 1.Element des 2.Federendes
e,ele+ele2+1,ele+ele2+2 !Elementenerzeugung für 2.Element des 2.Federendes
*repeat,ele2-1,1,1 !Knotengenerierung für Feder im Bereich ohne Steigung
!KRAFTEINLEITUNG
!===============================================================
csys,0 !Kartesisches Koordinatensysten
type,3 !MASSE-Punkt Elementtypaktivierung
real,3
e,1000
csys,0
nsel,s,loc,z,0, !Selektierung der NODES für z=0
cerig,1000,all,all !Lastverteilung wird auf alle Selektierten NODES verteilt
f,1000,fz,kraft !Krafteinleitung
!LAGERUNG
!===============================================================
k,5000,0,0,(steig*(wickl*((ele-1)/ele))) ! KEYPOINT x,y,z z==Variabel
*get,p,kp,5000,loc,x,y,z !Koordinatenermittlung für KEYPOINT
local,13,0,0,0,(steig*(wickl*((ele-1)/ele))) !Lokales Koordinatensystem
wpcsys,13 !Working Plane = CSYS13
csys,13 !Selektierung der NODES in z-Richtung
nsel,s,loc,z,-0.0005,10
d,all,all !Allseitige Einspannung
!Kontakt
!===============================================================
type,5 !Aktivierung Targeelement
esel,s,elem,,1,1 !Elementselektierung
aslv,s !AREA-SELECT
nsla,s,1 !NODE-SELECT
esurf,, !Elementgenerierung für Targe170
allsel,all !Deselektierung
type,4 !Aktivierung Kontaktelement
esel,s,elem,,1,1 !Elementselektierung
aslv,s !AREA-SELECT
nsla,s,1 !NODE-SELECT
esurf,, !Elementgenerierung für Targe176
allsel,all !Deselektierung
!ALLSEL
!===============================================================
/REPLOT,RESIZE
ALLSEL,ALL
GPLOT
FINISH
!SOLUTION
!===============================================================
/nerr,20,99999999
!Festlegung der zulässigen Fehleranzahl auf 99999999
/solu
antype,static,new
solcontrol,on
time,1
nsubst,kraft/30,10000,kraft/30
!Anzahl der substeps pro Rechnung
autots,on
! auto time stepping
pred,on
! predictor on
nropt,full,,on
! Newton-Raphson on
neqit,100
!maximal 100 Iterationsschritte pro substep
!cnvtol,f,kraft/50*1000,,0
!Angabe des Abbruchkriteriums
solve
finish
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP