Liebe Forenmitglieder,
vielleicht kann mir jemand bei einem neuen Problem helfen. Wie in meinem letzten Beitrag "Probleme beim meshen", arbeite ich immer noch an der Berechnung der Eigenfrequenzen einer mit Luft gefüllten akustischen Kavität.
Das Problem besteht aus einem inneren, quaderförmigen Volumen, welches mit FLUID220-Elementen gemesht wird. Dieses Volumen wird umschlossen mit einer "Schachtel" aus SOLID95-Elementen, welche die biegeweichen Wände repräsentieren, die die Luft umschließen. Diese äußere Schachtel (SOLID95) ist an den Kanten unverschieblich gelagert. Durchgeführt wird eine Modalanalyse.
Mein Problem ist nun, dass verschiedene Rechendurchgänge verschiedene Ergebnisse für die Eigenfrequenzen erzeugen; diese weichen in zwei nacheinander folgenden Rechnungen bis um Faktor 5 voneinander ab. Ich lasse mir die Ergebnisse meiner Rechnung immer über "General Postprocessing" -> "List Results" -> "Detailed Summary" ausgeben.
Könnte eventuell jemand meine Code laufen lassen und mir mitteilen, ob er den Fehler findet? Über jede Hilfe bin ich sehr dankbar!
Beste Grüße,
Moritz Becker
Der Code lautet:
!!! Präambel =======================================================
FINISH
/CLEAR
!/TITLE, Eigenmoden FSI-Volumen
!!! Präprozessor ===================================================
/PREP7
! Definition der Materialparameter:
Komp = 1420000 ! Kompressionsmodul bei Athmosphärendruck; Einheit: Pa
rho0 = 12
C = SQRT(Komp/rho0)
Ref = 0.00002 ! Referenzschalldruck
E_x = 210000000000
E_y = 210000000000
E_z = 210000000000
nu_x = 0.3
nu_y = 0.3
nu_z = 0.3
rho = 7850
! Definition der Elementtypen:
ET,1,FLUID220,,1 ! ohne FSI
ET,2,FLUID220,,0 ! FSI-Schicht
ET,3,SOLID95
! Definition der Materialien:
!Material 1:
MP,DENS,1,rho0
MP,SONC,1,C
MP,MU,1,0 ! keine Absorption
!Material 2:
MP,DENS,2,rho0
MP,SONC,2,C
MP,MU,2,0 ! keine Absorption
!Material 3:
MP,EX,3,E_x
MP,EY,3,E_y
MP,EZ,3,E_z
MP,PRXY,3,nu_x
MP,PRXZ,3,nu_z
MP,PRYZ,3,nu_y
MP,DENS,3,rho
! Definition der Geometrie:
!Definition der Längen:
L_x = 8
L_y = 3
L_z = 4
!Definition der Elementlängen Fluid und Solid
L_Ele_Fluid = 0.5
L_Ele_Solid = 0.5
! Definition der Keypoints:
X = L_x/2 ! zur Vereinfachung
Y = L_y/2
Z = L_z/2
! Keypoints Innerer Fluidkörper:
K,1,X,-Y,Z
K,2,-X,-Y,Z
K,3,-X,-Y,-Z
K,4,X,-Y,-Z
K,5,X,Y,Z
K,6,-X,Y,Z
K,7,-X,Y,-Z
K,8,X,Y,-Z
! Keypoints äußerer Solidkörper rechts:
K,9,X+L_Ele_Solid,-Y,Z
K,10,X,-Y,Z
K,11,X,-Y,-Z
K,12,X+L_Ele_Solid,-Y,-Z
K,13,X+L_Ele_Solid,Y,Z
K,14,X,Y,Z
K,15,X,Y,-Z
K,16,X+L_Ele_Solid,Y,-Z
! Keypoints äußerer Solidkörper links:
K,17,-X,-Y,Z
K,18,-(X+L_Ele_Solid),-Y,Z
K,19,-(X+L_Ele_Solid),-Y,-Z
K,20,-X,-Y,-Z
K,21,-X,Y,Z
K,22,-(X+L_Ele_Solid),Y,Z
K,23,-(X+L_Ele_Solid),Y,-Z
K,24,-X,Y,-Z
! Keypoints äußerer Soliddkörper vorne:
K,25,(X+L_Ele_Solid),-Y,(Z+L_Ele_Solid)
K,26,-(X+L_Ele_Solid),-Y,(Z+L_Ele_Solid)
K,27,-(X+L_Ele_Solid),-Y,Z
K,28,(X+L_Ele_Solid),-Y,Z
K,29,(X+L_Ele_Solid),Y,(Z+L_Ele_Solid)
K,30,-(X+L_Ele_Solid),Y,(Z+L_Ele_Solid)
K,31,-(X+L_Ele_Solid),Y,Z
K,32,(X+L_Ele_Solid),Y,Z
! Keypoints äußerer Soliddkörper hinten:
K,33,(X+L_Ele_Solid),-Y,-(Z+L_Ele_Solid)
K,34,-(X+L_Ele_Solid),-Y,-(Z+L_Ele_Solid)
K,35,-(X+L_Ele_Solid),-Y,-Z
K,36,(X+L_Ele_Solid),-Y,-Z
K,37,(X+L_Ele_Solid),Y,-(Z+L_Ele_Solid)
K,38,-(X+L_Ele_Solid),Y,-(Z+L_Ele_Solid)
K,39,-(X+L_Ele_Solid),Y,-Z
K,40,(X+L_Ele_Solid),Y,-Z
! Keypoints äußerer Soliddkörper unten:
K,41,(X+L_Ele_Solid),-(Y+L_Ele_Solid),(Z+L_Ele_Solid)
K,42,-(X+L_Ele_Solid),-(Y+L_Ele_Solid),(Z+L_Ele_Solid)
K,43,-(X+L_Ele_Solid),-(Y+L_Ele_Solid),-(Z+L_Ele_Solid)
K,44,(X+L_Ele_Solid),-(Y+L_Ele_Solid),-(Z+L_Ele_Solid)
K,45,(X+L_Ele_Solid),-(Y),(Z+L_Ele_Solid)
K,46,-(X+L_Ele_Solid),-(Y),(Z+L_Ele_Solid)
K,47,-(X+L_Ele_Solid),-(Y),-(Z+L_Ele_Solid)
K,48,(X+L_Ele_Solid),-(Y),-(Z+L_Ele_Solid)
! Keypoints äußerer Solidkörper oben:
K,49,(X+L_Ele_Solid),(Y),(Z+L_Ele_Solid)
K,50,-(X+L_Ele_Solid),(Y),(Z+L_Ele_Solid)
K,51,-(X+L_Ele_Solid),(Y),-(Z+L_Ele_Solid)
K,52,(X+L_Ele_Solid),(Y),-(Z+L_Ele_Solid)
K,53,(X+L_Ele_Solid),(Y+L_Ele_Solid),(Z+L_Ele_Solid)
K,54,-(X+L_Ele_Solid),(Y+L_Ele_Solid),(Z+L_Ele_Solid)
K,55,-(X+L_Ele_Solid),(Y+L_Ele_Solid),-(Z+L_Ele_Solid)
K,56,(X+L_Ele_Solid),(Y+L_Ele_Solid),-(Z+L_Ele_Solid)
!Erstellen der Volumina
V,1,2,3,4,5,6,7,8 ! Solid außen rechts ! Volumen 1
V,9,10,11,12,13,14,15,16 ! Solid außen rechts ! Volumen 2
V,17,18,19,20,21,22,23,24 ! Solid außen links 3
V,25,26,27,28,29,30,31,32 ! Solid außen vorne 4
V,33,34,35,36,37,38,39,40 ! Solid außen hinten 5
V,41,42,43,44,45,46,47,48 ! Solid außen unten 7
V,49,50,51,52,53,54,55,56 ! Solid außen oben 6
! Meshen von Volumen 2 bis 7; SOLID95
TYPE,3
MAT,3
VSEL,S,VOLU,,2,7
ESIZE,L_Ele_Solid
VMESH,ALL
VSEL,ALL
ALLSEL
! Koppeln der Schnittstellen der Solidvolumnina:
ESEL,S,TYPE,,3,3 ! Nur Elemente Typ3
CPINTF,UX,0.0001
CPINTF,UY,0.0001
CPINTF,UZ,0.0001
ESEL,ALL
! Meshen von Volumen 1; FLUID220
TYPE,1
MAT,1
VSEL,S,VOLU,,1
ESIZE,L_Ele_Fluid
VMESH,1,1
VSEL,ALL
! Ändern der Keyoptions von FLUID220 an der Berandung zum Fluid
ESEL,S,TYPE,,1
NSEL,S,EXT
ESLN,ALL
EMODIF,ALL,TYPE,2
ALLSEL
! Definiere FSI-Flag an der Grenze Fluid-Solid:
ESEL,S,TYPE,,2 ! select fluid elements
NSEL,S,EXT ! nodes on the exterior surface
SF,ALL,FSI
ESEL,ALL
NSEL,ALL
!SFELIST,ALL,FSI
! Aufbringen der Lagerungsbedingungen:
NSEL,S,LOC,X,-(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Z,-(L_Z/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,-(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Z,(L_Z/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Z,(-L_Z/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Z,(L_Z/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,-(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Y,-(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,-(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Y,(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Y,-(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,X,(L_X/2 + L_Ele_Solid)
NSEL,R,LOC,Y,(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,Z,-(L_Z/2 + L_Ele_Solid)
NSEL,R,LOC,Y,(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,Z,-(L_Z/2 + L_Ele_Solid)
NSEL,R,LOC,Y,-(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,Z,(L_Z/2 + L_Ele_Solid)
NSEL,R,LOC,Y,-(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
NSEL,S,LOC,Z,(L_Z/2 + L_Ele_Solid)
NSEL,R,LOC,Y,(L_Y/2 + L_Ele_Solid)
D,ALL,ALL,0
NSEL,ALL
FINISH ! /PREP7
!LÖSUNG ============================================================================
/SOLU
ANTYPE,2
MODOPT,UNSYM,40 !Subspace, 20 modes, unsymmertische Matrix
EQSLV,FRONT
MXPAND,40 !Expand 20 modes
SOLVE
FINISH ! /SOLU
(Ich habe diese Frage in keinem anderen Forum gestellt).
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP
Foren auf CAD.de (alle Foren)
| 
