Hallo Ansysgemeinde,
Ich arbeite erst seit zwei Wochen mit Ansys und wollte einen starren Rahmen, welcher auf sechs Luftfedern gelagert ist, mit Unwuchterregung simulieren. Es klappt auch, aber,wenn ich über eine Schleife meine Unwuchtkraft in Abhängigkeit von der Erregerfrequenz simuliere, bekomme ich im Resonanzbereich mir nicht erklärbare Ergebnisse für die Amplitude. Ich glaube, es liegt an der Schleife.
Kann mir jemand vielleicht helfen. Ich stehle mein ganzes Makro rein.
FINISH
/CLEAR,NOSTART
/PREP7
/TITLE Diplomarbeit Shakerrahmen
EMOD=210e+12
! Masse vom Rahmen mit Shaker[kg]
m11 = 988 !(Konstellation 1)
m12 = 786.5 !(Konstellation 2)
m13 = 700 !(Konstellation 3)
!Massenträgheitsmomente [kg*m²] für m11, m12, m13
Ixx1 = 53.637
Iyy1 = 39.675
Izz1 = 53.637
Ixx2 = 44.6975
Iyy2 = 33.0625
Izz2 = 44.6975
Ixx3 = 35.785
Iyy3 = 26.450
Izz3 = 35.785
! die Zusatzmasse[kg]
m21 = 94.7
m22 = 54.7
m23 = 14.7
! alle Abmassen sind in [m]
!Auswahl des Elementtypes Nr. 1 für Rahmen
!########################################################
!########################################################
ET,1,BEAM188
SECTYPE, 1, BEAM, HREC, Balken, 0
SECOFFSET, CENT
SECDATA,0.06,0.1,0.006,0.006,0.006,0.006,0,0,0,0,0,0
MP,EX,1,EMOD
MP,DENS,1,0.5
!######################################################
!######################################################
!Auswahl des Elementtypes Nr. 2 für Shakerzylinder
!########################################################
!########################################################
ET,2,BEAM188
SECTYPE, 2, BEAM, CSOLID, Zylinder, 0
SECOFFSET, CENT
SECDATA,3.09e-3,0,0,0,0,0,0,0,0,0,0,0
MP,EX,2,1e6
MP,DENS,2,0.1
!######################################################
!######################################################
!Auswahl des Elementtypes Nr.3 für Luftfeder
!######################################################
!######################################################
ET,3,COMBIN14
KEYOPT,3,1,0
KEYOPT,3,2,0
KEYOPT,3,3,0
R,3,118395.63,439.997,0, , ,0, !REAL CONSTANT (SET3)
!#########################################################################
!#########################################################################
!Auswahldes Elementtypes Nr. 4 für Masse m1
!#####################################################
!#####################################################
ET,4,MASS21 !Options (K3= 3-D w rot inert)
KEYOPT,4,1,0
KEYOPT,4,2,0
KEYOPT,4,3,0
R,4,m11,m11,m11,Ixx1,Iyy1,Izz1, !REAL CONSTANT (SET4)
!################################################################
!################################################################
!Auswahldes Elementtypes Nr. 5 für Masse m2
!#####################################################
!#####################################################
ET,5,MASS21 !Options (K3= 3-D w/o rot iner)
KEYOPT,5,1,0
KEYOPT,5,2,0
KEYOPT,5,3,2
R,5,m21 !REAL CONSTANT (SET5)
!################################################################
!################################################################
!Auswahl des Elementtypes Nr.6 für Luftfeder
!######################################################
!######################################################
ET,6,COMBIN14
KEYOPT,6,1,0
KEYOPT,6,2,0
KEYOPT,6,3,0 !REAL CONSTANT (SET6)
R,6,118395.63,0.00001, 0, , , 0,
!#########################################################################
!#########################################################################
!Erzeugung der Knoten für Rahmen
!#################################################################
!#################################################################
N,1,-0.52,0,0.36
N,2,0,0,0.36
N,3,0.52,0,0.36
N,4,0.52,0,-0.36
N,5,0,0,-0.36,
N,6,-0.52,0,-0.36
N,13,0,0,0
!###############################################################
!###############################################################
TYPE, 1 !Es wird in Mesh Attributes eingestellt
MAT, 1
ESYS, 0
SECNUM, 1
E,1,2
E,2,3
E,3,4
E,4,5
E,5,6
E,6,1
E,1,13
E,4,13
E,3,13
E,6,13
!Erzeugung der Knoten für Shakerzylinder
!#################################################################
!#################################################################
N,14,0,0.3,0
!#################################################################
!#################################################################
TYPE, 2
MAT, 2
ESYS, 0
SECNUM,2
E,13,14
/SHRINK,0.05
/ESHAPE,0.1
eplot
!Erzeugung der Knoten für Luftfedern
!###############################################################
!###############################################################
N,7,-0.52,-0.1,0.36
N,8,0,-0.1,0.36
N,9,0.52,-0.1,0.36
N,10,0.52,-0.1,-0.36
N,11,0,-0.1,-0.36
N,12,-0.52,-0.1,-0.36
!##############################################################
!##############################################################
TYPE, 3
ESYS, 0
REAL, 3,
E,1,7
E,2,8
E,3,9
E,4,10
E,5,11
E,6,12
!Erzeugung der Knoten für Luftfedern in x- und z-Richtung ( Dämpfung ist 0)
!###############################################################
!###############################################################
N,15,-0.62,0,0.36
N,16,-0.52,0,0.46
N,17,0,0,0.46
N,18,0.1,0,0.36
N,19,0.52,0,0.46
N,20,0.62,0,0.36
N,21,0.62,0,-0.36
N,22,0.52,0,-0.46
N,23,-0.12,0,-0.36
N,24,0,0,-0.46
N,25,-0.52,0,-0.46
N,26,-0.62,0,-0.36
TYPE,6
ESYS,0
REAL,6 ,
E,1,15
E,1,16
E,2,17
E,2,18
E,3,19
E,3,20
E,4,21
E,4,22
E,5,23
E,5,24
E,6,25
E,6,26
!Konstellationsmasse m(1,2,3) auf den Knoten 13 positionieren
!##############################################################
!##############################################################
TYPE,4
ESYS,0
REAL,4 ,
E,13
!Zusatzmasse m(1,2,3)auf den Knoten 14 positionieren
!##############################################################
!##############################################################
TYPE,5
ESYS,0
REAL,5 ,
E,14
FINISH
/SOLU
!Festlager
!############################################################
!############################################################
D,7,UX
D,7,UY
D,7,UZ
D,7,ROTx
D,7,ROTy
D,7,ROTz
D,8,UX
D,8,UY
D,8,UZ
D,8,ROTx
D,8,ROTy
D,8,ROTz
D,9,UX
D,9,UY
D,9,UZ
D,9,ROTx
D,9,ROTy
D,9,ROTz
D,10,UX
D,10,UY
D,10,UZ
D,10,ROTx
D,10,ROTy
D,10,ROTz
D,11,UX
D,11,UY
D,11,UZ
D,11,ROTx
D,11,ROTy
D,11,ROTz
D,12,UX
D,12,UY
D,12,UZ
D,12,ROTx
D,12,ROTy
D,12,ROTz
D,15,UX
D,15,UY
D,15,UZ
D,15,ROTx
D,15,ROTy
D,15,ROTz
D,16,UX
D,16,UY
D,16,UZ
D,16,ROTx
D,16,ROTy
D,16,ROTz
D,17,UX
D,17,UY
D,17,UZ
D,17,ROTx
D,17,ROTy
D,17,ROTz
D,18,UX
D,18,UY
D,18,UZ
D,18,ROTx
D,18,ROTy
D,18,ROTz
D,19,UX
D,19,UY
D,19,UZ
D,19,ROTx
D,19,ROTy
D,19,ROTz
D,20,UX
D,20,UY
D,20,UZ
D,20,ROTx
D,20,ROTy
D,20,ROTz
D,21,UX
D,21,UY
D,21,UZ
D,21,ROTx
D,21,ROTy
D,21,ROTz
D,22,UX
D,22,UY
D,22,UZ
D,22,ROTx
D,22,ROTy
D,22,ROTz
D,23,UX
D,23,UY
D,23,UZ
D,23,ROTx
D,23,ROTy
D,23,ROTz
D,24,UX
D,24,UY
D,24,UZ
D,24,ROTx
D,24,ROTy
D,24,ROTz
D,25,UX
D,25,UY
D,25,UZ
D,25,ROTx
D,25,ROTy
D,25,ROTz
D,26,UX
D,26,UY
D,26,UZ
D,26,ROTx
D,26,ROTy
D,26,ROTz
D,14,UX
!D,14,UY
D,14,UZ
D,14,ROTx
D,14,ROTy
D,14,ROTz
FINISH
!###################################################################
!###################################################################
*dim,par_sol,array,501,7,
! Initialisierung der Kräfte Spalte 2
*ask,r0,'konstante Anfangsamplitude [m]',
*ask,u,'Zusatzmasse mit Anker [kg]',
start_freq = 2.0 ! Startfrequenz der Simulation
final_freq = 30 ! Endfrequenz der Simulation
step_freq = 5 ! Schrittweite Frequenz
PI = ACOS(-1) ! PI_Berechnung
par_index=1
*do,par_freq,start_freq,final_freq,step_freq
par_sol(par_index,2)= -(((2*PI*par_freq)*(2*PI*par_freq))*r0*u)
par_index = par_index + 1
*enddo
par_index=1
*do,par_freq,start_freq,final_freq,step_freq
par_sol(par_index,5)= (((2*PI*par_freq)*(2*PI*par_freq))*r0*u)
par_index = par_index + 1
*enddo
! *do loop Harmonische Analyse
par_index=1
*do,par_freq,start_freq,final_freq,step_freq
/solu
ANTYPE,3
HROPT,FULL
HARFRQ,par_freq,par_freq
NSUBST,1
KBC,1
dmprat=0.001
par_sol(par_index,1)= par_freq
F,13,FY,par_sol(par_index,2)
F,14,FY,par_sol(par_index,5)
solve
finish
/post1
set,1
temp1 = 0
*get,temp1,node,13,u,y
par_sol(par_index,3) = temp1
set,1
temp2 = 0
*get,temp2,node,14,u,y
par_sol(par_index,4) = temp2
set,1
temp3 = 0
*get,temp3,node,1,u,y
par_sol(par_index,6) = temp3
set,1
temp4 = 0
*get,temp4,node,4,u,y
par_sol(par_index,7) = temp4
finish
par_index=par_index + 1
*enddo ! Ende der Analyseschleife
finish
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP