Liebe Nutzerinnen und Nutzer,
im Rahmen meiner Masterarbeit beschäftige ich mich mit der effizienten Modellierung von Schallfeldern im Inneren von akustischen Kavitäten. Diese Kavität hat die Gestalt eines hohlen Quaders, der mit Luft gefüllt ist und soll in Ansys modelliert werden; leider bin ich darin ein absoluter Anfänger.
Das Problem soll folgendermaßen modelliert werden:
Das Innere der Kavität wird mit FLUID220-Elementen modelliert und gemesht. Dieses Innere Volumen soll anschließend mit 6 Wänden aus SOLID95-Elementen umschlossen werden; die Grenze zwischen Festkörper und Fluid wird dann als FSI-Grenze modelliert, da sich Wand und Fluid gegenseitig beeinflussen sollen.
Leider scheitere ich schon an der einfachsten Modellierung, denn ich habe folgendes Problem: Obwohl ich nur das Innere Volumen mit Fluid-Elementen meshen möchte, werden auch die Wandelemente, die mit Solids modelliert werden sollen, ebenfalls mit Fluid-Elementen gemesht. Man sieht das deutlich, wenn man das Skript in Ansys laufen lässt und sich die Geometrie und die Elementtypen anzeigen lässt; die äußere Schale besteht aus Elementtyp 1 (FLUID220), obwohl es ja aus Solids bestehen sollte. Kann mir jemand sagen, was ich da falsch gemacht habe?
Mein Code lautet:
/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 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 = 1
! 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 1; FLUID220
VSEL,S,VOLU,,1
ESIZE,L_Ele_Fluid
ET,1,FLUID220,,1
VMESH,1,1
VSEL,ALL
! Ändern der Keyoptions von FLUID220
NSEL,S,LOC,X,(L_x/2)
NSEL,A,LOC,X,-(L_x/2)
NSEL,A,LOC,Z,(L_z/2)
NSEL,A,LOC,Z,-(L_z/2)
NSEL,A,LOC,Y,-(L_y/2)
NSEL,A,LOC,Y,(L_y/2)
ESLN
ET,2,FLUID220,,0
EMODIF,ALL,TYPE,2
ESEL,ALL
NSEL,ALL
ALLSEL
! Meshen von Volumen 2 bis 7; SOLID95
VSEL,S,VOLU,,2,7
ESIZE,L_Ele_Solid
ET,3,SOLID95
VMESH,ALL
VSEL,ALL
! ! Definiere FSi-Flag
! ESEL,S,TYPE,,2 ! select fluid elements
! NSEL,S,EXT ! nodes on the exterior surface
! SF,ALL,FSIN
! ESEL,ALL
! NSEL,ALL
FINISH ! /PREP7
Über jede Hilfe bin ich sehr dankbar! Ich habe diese Frage in keinem anderen Forum gestellt.
Moritz Becker
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP