Iterationsproblem

[Bjoerk]

Ich hann folschendes Iterationsproblem (Skizze Siehe Anlage). Es soll eine Fläche A gefunden werden, die gewisse Kriterien erfüllt.

Vorgabe ist ein Rechteck b * d sowie eine Belastung Nx, My und Mz.

Das Koordinatensystem x,y,z befindet sich im Schwerpunkt des Querschnitts b * d.

In Richtung der x-Achse drückt jemand (Kraft Nx), an den beiden anderen Achsen drehen welche (Momente My, Mz).

Der Querschnitt soll keine Zugspannungen aufnehmen können, d. h. der Querschnitt reißt auf und es stellt sich eine Fläche A ein, deren Flächeninhalt gleich Nx ist (In Wirklichkeit ein räumlicher Spannungspolyeder, hier aber wurscht).

Damit die Fläche A auch die Momente My und Mz im Gleichgewicht halten kann, muß sie sich wie folgt einstellen (S = Schwerpunkt der Fläche A):

Delphi-Quellcode:

repeat
  A := GetA .. // ??
  S:= GetSchwerpunkt(A);
  C:= GetC; // Konstante
until (A * C = NX) and (A * C * S.y = Mz) and (A * C * S.z = My);

Ich hab das mal vor 20 Jahren mit einer doppelten Repeat Schleife programmiert, kann diese procedure (1000 Zeilen) aber beim besten Willen nicht mehr nachvollziehen.

Anybody?

Crosspost