Nichtkonstante Randbedingungen "teilweise" an einer Kante definieren - MATLAB

Question

Ich muss Dirichlet- und Neumann-Bedingungen in einem einfachen Spannungsproblem definieren nur auf Teil einer der Kanten einer Platte. Matlab help für nicht konstante Randbedingungen definiert, zeigt an, dass Funktionen wie geschrieben werden müssen:

applyBoundaryCondition(model,'edge',1,'r',@myrfun); 
applyBoundaryCondition(model,'face',2,'g',@mygfun,'q',@myqfun); 
applyBoundaryCondition(model,'edge',[3,4],'u',@myufun,'EquationIndex',[2,3]); 

Er sagt außerdem, dass jede Funktion die folgende Syntax haben muss.

function bcMatrix = myfun(region,state) 

und schließlich, dass " Region" ist eine Struktur, die Felder region.x (x-Koordinate der Punkte) enthält, region.y (y-Koordinate der Punkte), usw. Wenn es Neumann-Bedingungen gibt (was mein Fall ist), dann geben Solver die folgenden Daten in der Regionsstruktur weiter: region.nx - (x-Komponente des Normalenvektors an den Auswertungspunkten), usw. Meine Fragen sind:

1. Fr om wo bekomme ich die Struktur region?
2. Wie kann ich das Argument übergeben, dass Randbedingungen in einem Teil einer der Kanten gelten?

Danke !!

Answer 1

@Oliver,

1) Ich denke, dass Sie nicht die Struktur region machen müssen, aber Funktionen machen, die sie verwenden können. Da die Randbedingung im Allgemeinen vom Standort abhängt, benötigen Sie region.x und region.y.

2) Sie können region.x und region.y verwenden, um Randbedingungen abhängig von der Position zu machen. Dies ist eine Möglichkeit, sie "teilweise" anzuwenden, wenn die Art der Randbedingung gleich ist. Andernfalls müssen Sie die Aufteilung in der Grenze explizit definieren. Dies geschieht wahrscheinlich während der Definition der Geometrie des Problems.