@theerapat-s28/ts-structural-eng-tools 2.3.0 → 3.1.0

package/dist/fea/frame/elements/shell/shell-axes.js

@@ -0,0 +1,102 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.isValidShellQuad = exports.shellAxes = void 0;
+const frame_constant_1 = require("../../../../core/constants/frame.constant");
+const quad_geometry_1 = require("./quad-geometry");
+const sub = (a, b) => [a[0] - b[0], a[1] - b[1], a[2] - b[2]];
+const dot = (a, b) => a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+const cross = (a, b) => [
+    a[1] * b[2] - a[2] * b[1],
+    a[2] * b[0] - a[0] * b[2],
+    a[0] * b[1] - a[1] * b[0],
+];
+const unit = (a) => {
+    const n = Math.hypot(a[0], a[1], a[2]);
+    return [a[0] / n, a[1] / n, a[2] / n];
+};
+/**
+ * Local axes of a four-node shell element, ETABS-style, and its corners on the mean plane.
+ *
+ * - Local 3 is the unit normal `n ∥ d₁ × d₂` of the diagonals `d₁ = x₃ − x₁`, `d₂ = x₄ − x₂`: the
+ *   right-hand rule around the nodes, as for null areas (for four nodes Newell's vector area is
+ *   exactly `½ d₁ × d₂`).
+ * - Non-horizontal element: local 2 is +Z projected onto the plane (upward on a wall) and
+ *   local 1 = 2 × 3 is horizontal.
+ * - Horizontal element (the normal's horizontal part below the vertical-member tolerance):
+ *   local 1 = +X and local 2 = 3 × 1, so local 2 = +Y when local 3 = +Z.
+ * - `angle` then rotates local 1 and 2 about local 3 by the right-hand rule.
+ *
+ * The mean plane passes through the centroid with normal n; a warped element's corners lie
+ * alternately h above and below it.
+ *
+ * @param points - The four corners in global axes (mm), in order around the element
+ * @param angle - Rotation of local 1 and 2 about local 3 in rad (default 0)
+ */
+const shellAxes = (points, angle = 0) => {
+    const e3 = unit(cross(sub(points[2], points[0]), sub(points[3], points[1])));
+    let e1;
+    let e2;
+    if (Math.hypot(e3[0], e3[1]) < frame_constant_1.VERTICAL_MEMBER_TOLERANCE) {
+        e1 = [1, 0, 0];
+        e2 = cross(e3, e1);
+    }
+    else {
+        e2 = unit([-e3[2] * e3[0], -e3[2] * e3[1], 1 - e3[2] * e3[2]]); // Z − (Z·n)·n
+        e1 = cross(e2, e3);
+    }
+    if (angle !== 0) {
+        const c = Math.cos(angle);
+        const s = Math.sin(angle);
+        const r1 = [c * e1[0] + s * e2[0], c * e1[1] + s * e2[1], c * e1[2] + s * e2[2]];
+        const r2 = [-s * e1[0] + c * e2[0], -s * e1[1] + c * e2[1], -s * e1[2] + c * e2[2]];
+        e1 = r1;
+        e2 = r2;
+    }
+    const origin = [0, 1, 2].map((k) => (points[0][k] + points[1][k] + points[2][k] + points[3][k]) / 4);
+    const relative = points.map((p) => sub(p, origin));
+    return {
+        lambda: [e1, e2, e3],
+        origin,
+        corners: relative.map((r) => [dot(r, e1), dot(r, e2)]),
+        offsets: relative.map((r) => dot(r, e3)),
+    };
+};
+exports.shellAxes = shellAxes;
+/**
+ * Whether four corners make a valid shell element: no repeated corner, a non-degenerate
+ * convex quadrilateral in order around it (det J > 0 at every corner of the bilinear map), and
+ * no corner farther than `SHELL_WARPING_LIMIT · √A` from the mean plane.
+ */
+const isValidShellQuad = (points) => {
+    const size = Math.max(...points.map((p, i) => {
+        const q = points[(i + 1) % 4];
+        return Math.hypot(q[0] - p[0], q[1] - p[1], q[2] - p[2]);
+    }));
+    const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * size;
+    for (let i = 0; i < 4; i++) {
+        for (let j = i + 1; j < 4; j++) {
+            const [p, q] = [points[i], points[j]];
+            if (!(Math.hypot(q[0] - p[0], q[1] - p[1], q[2] - p[2]) > tolerance))
+                return false;
+        }
+    }
+    const d1 = [0, 1, 2].map((k) => points[2][k] - points[0][k]);
+    const d2 = [0, 1, 2].map((k) => points[3][k] - points[1][k]);
+    const area = Math.hypot(d1[1] * d2[2] - d1[2] * d2[1], d1[2] * d2[0] - d1[0] * d2[2], d1[0] * d2[1] - d1[1] * d2[0]) / 2;
+    if (!(area > tolerance * size))
+        return false;
+    const axes = (0, exports.shellAxes)(points);
+    for (const [xi, eta] of [
+        [-1, -1],
+        [1, -1],
+        [1, 1],
+        [-1, 1],
+    ]) {
+        if (!((0, quad_geometry_1.quadJacobian)(axes.corners, xi, eta).detJ > (frame_constant_1.AREA_GEOMETRY_TOLERANCE * area) / 4)) {
+            return false;
+        }
+    }
+    return Math.max(...axes.offsets.map(Math.abs)) <= frame_constant_1.SHELL_WARPING_LIMIT * Math.sqrt(area);
+};
+exports.isValidShellQuad = isValidShellQuad;
+//# sourceMappingURL=shell-axes.js.map

package/dist/fea/frame/elements/shell/shell-axes.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"shell-axes.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/shell/shell-axes.ts"],"names":[],"mappings":";;;AAAA,uEAI4C;AAG5C,mDAA0D;AAK1D,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAChF,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1E,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC;IACxC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;CAC1B,CAAC;AACF,MAAM,IAAI,GAAG,CAAC,CAAO,EAAQ,EAAE;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACxC,CAAC,CAAC;AAiBF;;;;;;;;;;;;;;;;;GAiBG;AACI,MAAM,SAAS,GAAG,CAAC,MAAuB,EAAE,KAAK,GAAG,CAAC,EAAa,EAAE;IACzE,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7E,IAAI,EAAQ,CAAC;IACb,IAAI,EAAQ,CAAC;IACb,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,0CAAyB,EAAE,CAAC;QACzD,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACf,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACrB,CAAC;SAAM,CAAC;QACN,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,cAAc;QAC9E,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACrB,CAAC;IACD,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAS,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACvF,MAAM,EAAE,GAAS,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1F,EAAE,GAAG,EAAE,CAAC;QACR,EAAE,GAAG,EAAE,CAAC;IACV,CAAC;IAED,MAAM,MAAM,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAChC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAC/D,CAAC;IACV,MAAM,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;IACnD,OAAO;QACL,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QACpB,MAAM;QACN,OAAO,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAQ,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC5D,OAAO,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;KACzC,CAAC;AACJ,CAAC,CAAC;AA9BW,QAAA,SAAS,aA8BpB;AAEF;;;;GAIG;AACI,MAAM,gBAAgB,GAAG,CAAC,MAAuB,EAAW,EAAE;IACnE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CACnB,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACrB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC,CAAC,CACH,CAAC;IACF,MAAM,SAAS,GAAG,wCAAuB,GAAG,IAAI,CAAC;IACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;gBAAE,OAAO,KAAK,CAAC;QACrF,CAAC;IACH,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,IAAI,GACR,IAAI,CAAC,KAAK,CACR,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAC9B,GAAG,CAAC,CAAC;IACR,IAAI,CAAC,CAAC,IAAI,GAAG,SAAS,GAAG,IAAI,CAAC;QAAE,OAAO,KAAK,CAAC;IAC7C,MAAM,IAAI,GAAG,IAAA,iBAAS,EAAC,MAAM,CAAC,CAAC;IAC/B,KAAK,MAAM,CAAC,EAAE,EAAE,GAAG,CAAC,IAAI;QACtB,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACR,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACP,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;KACR,EAAE,CAAC;QACF,IAAI,CAAC,CAAC,IAAA,4BAAY,EAAC,IAAI,CAAC,OAAO,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,wCAAuB,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACvF,OAAO,KAAK,CAAC;QACf,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,oCAAmB,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAC1F,CAAC,CAAC;AAnCW,QAAA,gBAAgB,oBAmC3B"}

package/dist/fea/frame/index.d.ts

@@ -6,5 +6,5 @@ export { beam3dLocalStiffness, beam3dLocalStiffnessRightHanded, } from "./elemen
 export { memberAxes } from "./transformation/coordinate-transform";
 export { memberDemandEnvelope } from "./post/demand-envelope";
 export type { MemberDemandEnvelope, DemandEnvelopeOptions } from "./post/demand-envelope";
-export type { FrameId, FrameNode, DofRestraint, Restraint, DofRelease, MemberEndRelease, Member, NodalLoad, LoadDirection, UniformMemberLoad, PointMemberLoad, TrapezoidalMemberLoad, TorqueMemberLoad, MemberLoad, AreaDistribution, NullArea, AreaLoadDirection, AreaLoad, LoadCaseType, LoadCase, CombinationTerm, LoadCombination, FrameModel, } from "./types/model.type";
-export type { NodalDisplacement, NodalReaction, MemberEndForce, MemberEndForces, ReleaseDeformation, MemberReleaseDeformations, AutoRestrainedDof, MemberStation, StationComponent, ExtremeValue, MemberDiagram, GoverningValue, GoverningExtremeValue, Bound, FrameCaseResult, FrameEnvelope, FrameResult, } from "./types/results.type";
+export type { FrameId, FrameNode, DofRestraint, Restraint, LineRestraint, DofRelease, MemberEndRelease, Member, NodalLoad, LoadDirection, UniformMemberLoad, PointMemberLoad, TrapezoidalMemberLoad, TorqueMemberLoad, MemberLoad, AreaDistribution, NullArea, ShellType, ShellSection, ShellModifiers, ShellArea, ShellMesh, AreaLoadDirection, AreaLoad, LoadCaseType, LoadCase, CombinationTerm, LoadCombination, FrameModel, } from "./types/model.type";
+export type { NodalDisplacement, NodalReaction, MemberEndForce, MemberEndForces, ReleaseDeformation, MemberReleaseDeformations, AutoRestrainedDof, MemberStation, StationComponent, ExtremeValue, MemberDiagram, ShellCornerForces, ShellForceComponent, ShellElementForces, FrameMesh, GoverningValue, GoverningExtremeValue, Bound, FrameCaseResult, FrameEnvelope, FrameResult, } from "./types/results.type";

package/dist/fea/frame/mesh/expand-model.d.ts

@@ -0,0 +1,29 @@
+import type { Warnings } from "../../../core/types/output-message.type";
+import type { FrameId, FrameModel, Member } from "../types/model.type";
+import type { FrameMesh } from "../types/results.type";
+import { type ShellElement } from "./mesh-shells";
+/** A member with the generated mesh nodes lying inside it, where the solver splits it. */
+export interface SplitMember extends Member {
+    /** Generated node ids strictly between its ends, in order from the start. */
+    through?: FrameId[];
+}
+/** The model the solver analyses: shells meshed, line restraints spread to their nodes. */
+export interface ExpandedModel {
+    /**
+     * Nodes include the generated ones; `shells` are the elements; members list the generated
+     * nodes they pass through; restraints include the lines.
+     */
+    model: FrameModel & {
+        shells: ShellElement[];
+        members: SplitMember[];
+    };
+    /** The shell mesh for the result; undefined without shells. */
+    mesh: FrameMesh | undefined;
+    warnings: Warnings;
+}
+/**
+ * Meshes a valid model's shell areas, spreads its line restraints and marks the members that
+ * generated nodes split (see `meshShells`).
+ * Without shells and line restraints the model is returned as it is.
+ */
+export declare const expandModel: (model: FrameModel) => ExpandedModel;

package/dist/fea/frame/mesh/expand-model.js

@@ -0,0 +1,133 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.expandModel = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+const shell_axes_1 = require("../elements/shell/shell-axes");
+const mesh_shells_1 = require("./mesh-shells");
+const DOF_KEYS = ["ux", "uy", "uz", "rx", "ry", "rz"];
+/**
+ * Restraints of every node on each line restraint's segment (ends included), merged with the
+ * point restraints: a DOF fixed by either is fixed.
+ */
+const spreadLineRestraints = (model, nodes) => {
+    const lines = model.lineRestraints ?? [];
+    if (lines.length === 0)
+        return model.restraints;
+    const position = new Map(nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    const merged = new Map(model.restraints.map((r) => [r.nodeId, { ...r }]));
+    for (const line of lines) {
+        const [P, Q] = line.nodeIds.map((id) => position.get(id));
+        for (const node of nodes) {
+            const on = line.nodeIds.includes(node.id) ||
+                (0, mesh_shells_1.pointOnSegment)([node.x, node.y, node.z], P, Q, true) !== null;
+            if (!on)
+                continue;
+            const restraint = merged.get(node.id) ?? { nodeId: node.id };
+            for (const key of DOF_KEYS)
+                if (line[key] === true)
+                    restraint[key] = true;
+            merged.set(node.id, restraint);
+        }
+    }
+    return [...merged.values()];
+};
+/**
+ * Warnings for member nodes lying inside a shell area but not on its mesh: such a node is not
+ * connected to the shell.
+ */
+const disconnectedNodes = (model, meshNodes) => {
+    const warnings = [];
+    const shells = model.shells ?? [];
+    const position = new Map(model.nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    const memberNodes = new Set(model.members.flatMap((m) => [m.startNodeId, m.endNodeId]));
+    for (const shell of shells) {
+        const points = shell.nodeIds.map((id) => position.get(id));
+        const axes = (0, shell_axes_1.shellAxes)(points);
+        const size = Math.max(...axes.corners.map(([x, y]) => Math.hypot(x, y)));
+        const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * size;
+        for (const id of memberNodes) {
+            if (meshNodes.has(id))
+                continue;
+            const p = position.get(id);
+            const r = [p[0] - axes.origin[0], p[1] - axes.origin[1], p[2] - axes.origin[2]];
+            const [e1, e2, e3] = axes.lambda;
+            const z = r[0] * e3[0] + r[1] * e3[1] + r[2] * e3[2];
+            if (Math.abs(z) > Math.max(tolerance, ...axes.offsets.map(Math.abs)))
+                continue;
+            const x = r[0] * e1[0] + r[1] * e1[1] + r[2] * e1[2];
+            const y = r[0] * e2[0] + r[1] * e2[1] + r[2] * e2[2];
+            const inside = axes.corners.every(([ax, ay], k) => {
+                const [bx, by] = axes.corners[(k + 1) % 4];
+                return (bx - ax) * (y - ay) - (by - ay) * (x - ax) > tolerance * size;
+            });
+            if (inside) {
+                warnings.push({
+                    reference: "Shell area meshing",
+                    message: `Node "${String(id)}" lies inside shell area "${String(shell.id)}" but not on its mesh, so its members are not connected to the shell; divide the area through the node.`,
+                });
+            }
+        }
+    }
+    return warnings;
+};
+/**
+ * The generated nodes lying strictly inside each member (a beam along a meshed edge), so the
+ * member can be split there and share them with the shells.
+ */
+const splitMembers = (model, generated) => {
+    if (generated.length === 0)
+        return model.members;
+    const position = new Map(model.nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    return model.members.map((member) => {
+        const [P, Q] = [position.get(member.startNodeId), position.get(member.endNodeId)];
+        const lo = [0, 1, 2].map((k) => Math.min(P[k], Q[k]));
+        const hi = [0, 1, 2].map((k) => Math.max(P[k], Q[k]));
+        const slack = frame_constant_1.AREA_GEOMETRY_TOLERANCE * Math.hypot(Q[0] - P[0], Q[1] - P[1], Q[2] - P[2]);
+        const through = [];
+        for (const node of generated) {
+            const x = [node.x, node.y, node.z];
+            if (x.some((v, k) => v < lo[k] - slack || v > hi[k] + slack))
+                continue;
+            const t = (0, mesh_shells_1.pointOnSegment)(x, P, Q);
+            if (t !== null)
+                through.push({ id: node.id, t });
+        }
+        return through.length === 0
+            ? member
+            : { ...member, through: through.sort((a, b) => a.t - b.t).map((h) => h.id) };
+    });
+};
+/**
+ * Meshes a valid model's shell areas, spreads its line restraints and marks the members that
+ * generated nodes split (see `meshShells`).
+ * Without shells and line restraints the model is returned as it is.
+ */
+const expandModel = (model) => {
+    const shells = model.shells ?? [];
+    if (shells.length === 0) {
+        return {
+            model: { ...model, shells: [], restraints: spreadLineRestraints(model, model.nodes) },
+            mesh: undefined,
+            warnings: [],
+        };
+    }
+    const meshed = (0, mesh_shells_1.meshShells)(model.nodes, shells);
+    const nodes = [...model.nodes, ...meshed.nodes];
+    const meshNodes = new Set(meshed.elements.flatMap((e) => e.nodeIds));
+    return {
+        model: {
+            ...model,
+            nodes,
+            members: splitMembers(model, meshed.nodes),
+            shells: meshed.elements,
+            restraints: spreadLineRestraints(model, nodes),
+        },
+        mesh: {
+            nodes: meshed.nodes,
+            shells: meshed.elements.map(({ id, areaId, nodeIds }) => ({ id, areaId, nodeIds })),
+        },
+        warnings: disconnectedNodes(model, meshNodes),
+    };
+};
+exports.expandModel = expandModel;
+//# sourceMappingURL=expand-model.js.map

package/dist/fea/frame/mesh/expand-model.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"expand-model.js","sourceRoot":"","sources":["../../../../src/fea/frame/mesh/expand-model.ts"],"names":[],"mappings":";;;AAAA,uEAA6E;AAG7E,6DAAoE;AAGpE,+CAA8E;AAE9E,MAAM,QAAQ,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AAoB/D;;;GAGG;AACH,MAAM,oBAAoB,GAAG,CAAC,KAAiB,EAAE,KAA0B,EAAe,EAAE;IAC1F,MAAM,KAAK,GAAG,KAAK,CAAC,cAAc,IAAI,EAAE,CAAC;IACzC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC,UAAU,CAAC;IAChD,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,MAAM,MAAM,GAAG,IAAI,GAAG,CAAqB,KAAK,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC9F,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;QAC3D,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;YACzB,MAAM,EAAE,GACN,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;gBAC9B,IAAA,4BAAc,EAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,KAAK,IAAI,CAAC;YAChE,IAAI,CAAC,EAAE;gBAAE,SAAS;YAClB,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC;YAC7D,KAAK,MAAM,GAAG,IAAI,QAAQ;gBAAE,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,IAAI;oBAAE,SAAS,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;YAC1E,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,SAAS,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IACD,OAAO,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC;AAC9B,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,iBAAiB,GAAG,CAAC,KAAiB,EAAE,SAAuB,EAAY,EAAE;IACjF,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,EAAE,CAAC;IAClC,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,WAAW,EAAE,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IACxF,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,IAAA,sBAAS,EAAC,MAAM,CAAC,CAAC;QAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACzE,MAAM,SAAS,GAAG,wCAAuB,GAAG,IAAI,CAAC;QACjD,KAAK,MAAM,EAAE,IAAI,WAAW,EAAE,CAAC;YAC7B,IAAI,SAAS,CAAC,GAAG,CAAC,EAAE,CAAC;gBAAE,SAAS;YAChC,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC;YAC5B,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;YACtF,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;YACjC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACrD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;gBAAE,SAAS;YAC/E,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACrD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACrD,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE;gBAChD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC3C,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,SAAS,GAAG,IAAI,CAAC;YACxE,CAAC,CAAC,CAAC;YACH,IAAI,MAAM,EAAE,CAAC;gBACX,QAAQ,CAAC,IAAI,CAAC;oBACZ,SAAS,EAAE,oBAAoB;oBAC/B,OAAO,EAAE,SAAS,MAAM,CAAC,EAAE,CAAC,6BAA6B,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,yGAAyG;iBACnL,CAAC,CAAC;YACL,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,YAAY,GAAG,CAAC,KAAiB,EAAE,SAAsB,EAAiB,EAAE;IAChF,IAAI,SAAS,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC,OAAO,CAAC;IACjD,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,OAAO,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QAClC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,WAAW,CAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAE,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,KAAK,GAAG,wCAAuB,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1F,MAAM,OAAO,GAAiC,EAAE,CAAC;QACjD,KAAK,MAAM,IAAI,IAAI,SAAS,EAAE,CAAC;YAC7B,MAAM,CAAC,GAAS,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;YACzC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;gBAAE,SAAS;YACvE,MAAM,CAAC,GAAG,IAAA,4BAAc,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAClC,IAAI,CAAC,KAAK,IAAI;gBAAE,OAAO,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;QACnD,CAAC;QACD,OAAO,OAAO,CAAC,MAAM,KAAK,CAAC;YACzB,CAAC,CAAC,MAAM;YACR,CAAC,CAAC,EAAE,GAAG,MAAM,EAAE,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;IACjF,CAAC,CAAC,CAAC;AACL,CAAC,CAAC;AAEF;;;;GAIG;AACI,MAAM,WAAW,GAAG,CAAC,KAAiB,EAAiB,EAAE;IAC9D,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,EAAE,CAAC;IAClC,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACxB,OAAO;YACL,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,MAAM,EAAE,EAAE,EAAE,UAAU,EAAE,oBAAoB,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,CAAC,EAAE;YACrF,IAAI,EAAE,SAAS;YACf,QAAQ,EAAE,EAAE;SACb,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,IAAA,wBAAU,EAAC,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAChD,MAAM,SAAS,GAAG,IAAI,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;IACrE,OAAO;QACL,KAAK,EAAE;YACL,GAAG,KAAK;YACR,KAAK;YACL,OAAO,EAAE,YAAY,CAAC,KAAK,EAAE,MAAM,CAAC,KAAK,CAAC;YAC1C,MAAM,EAAE,MAAM,CAAC,QAAQ;YACvB,UAAU,EAAE,oBAAoB,CAAC,KAAK,EAAE,KAAK,CAAC;SAC/C;QACD,IAAI,EAAE;YACJ,KAAK,EAAE,MAAM,CAAC,KAAK;YACnB,MAAM,EAAE,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,CAAC;SACpF;QACD,QAAQ,EAAE,iBAAiB,CAAC,KAAK,EAAE,SAAS,CAAC;KAC9C,CAAC;AACJ,CAAC,CAAC;AA1BW,QAAA,WAAW,eA0BtB"}

package/dist/fea/frame/mesh/mesh-shells.d.ts

@@ -0,0 +1,38 @@
+import { type Vec3 } from "../elements/shell/shell-axes";
+import type { FrameId, FrameNode, ShellArea } from "../types/model.type";
+/** A shell element of the analysed model: a shell area, or one element of its mesh. */
+export interface ShellElement extends ShellArea {
+    /** The shell area it belongs to (its own id when the area is one element). */
+    areaId: FrameId;
+}
+/** The result of meshing a model's shell areas. */
+export interface ShellMeshResult {
+    /** Generated nodes, in global coordinates (mm). */
+    nodes: FrameNode[];
+    /** Every shell element, area by area in model order. */
+    elements: ShellElement[];
+    /** Areas that could not be meshed (no conforming division counts, a bad generated element, or a generated id already in use). */
+    problems: FrameId[];
+}
+/** Parameter t ∈ (0, 1) of `x` on the segment p → q when it lies strictly inside it, else null. */
+export declare const pointOnSegment: (x: Vec3, p: Vec3, q: Vec3, inclusive?: boolean) => number | null;
+/**
+ * Meshes the shell areas of a model whose shells and nodes are valid.
+ *
+ * - Each area edge is cut into segments at the model nodes lying on it (hard points), so an
+ *   edge shared whole by two areas, or partly where areas meet at a node (a T-junction), shares
+ *   its segments.
+ * - Every segment gets a division count: at least 1, at least `length / maxSize` for an area
+ *   with `mesh.maxSize`, and raised until in every area opposite edges add up to the same
+ *   count (and at least `mesh.divisions`), adding each division where the elements are longest.
+ * - Segment nodes are evenly spaced; interior nodes come from transfinite interpolation of the
+ *   four edges, so a planar area stays planar.
+ *
+ * An area that ends with 1 × 1 divisions is its own element, with its own id. Generated ids:
+ * `"<a>~<b>:<k>"` for the k-th node from model node a to b on a segment, `"<areaId>:<i>,<j>"`
+ * for interior nodes and for elements.
+ *
+ * @param nodes - Model nodes
+ * @param shells - Shell areas, each with 4 valid corner nodes
+ */
+export declare const meshShells: (nodes: FrameNode[], shells: ShellArea[]) => ShellMeshResult;

package/dist/fea/frame/mesh/mesh-shells.js

@@ -0,0 +1,244 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.meshShells = exports.pointOnSegment = void 0;
+const frame_constant_1 = require("../../../core/constants/frame.constant");
+const shell_axes_1 = require("../elements/shell/shell-axes");
+const sub = (p, q) => [p[0] - q[0], p[1] - q[1], p[2] - q[2]];
+const dot = (p, q) => p[0] * q[0] + p[1] * q[1] + p[2] * q[2];
+const lerp = (p, q, t) => [
+    p[0] + (q[0] - p[0]) * t,
+    p[1] + (q[1] - p[1]) * t,
+    p[2] + (q[2] - p[2]) * t,
+];
+/** A string that tells `1` and `"1"` apart, for ordering and map keys. */
+const idKey = (id) => `${typeof id}:${String(id)}`;
+/** Parameter t ∈ (0, 1) of `x` on the segment p → q when it lies strictly inside it, else null. */
+const pointOnSegment = (x, p, q, inclusive = false) => {
+    const d = sub(q, p);
+    const L2 = dot(d, d);
+    if (!(L2 > 0))
+        return null;
+    const L = Math.sqrt(L2);
+    const tolerance = frame_constant_1.AREA_GEOMETRY_TOLERANCE * L;
+    const t = dot(sub(x, p), d) / L2;
+    const margin = inclusive ? -tolerance / L : tolerance / L;
+    if (t < margin || t > 1 - margin)
+        return null;
+    const off = sub(x, lerp(p, q, t));
+    return Math.hypot(off[0], off[1], off[2]) <= tolerance ? t : null;
+};
+exports.pointOnSegment = pointOnSegment;
+/**
+ * Meshes the shell areas of a model whose shells and nodes are valid.
+ *
+ * - Each area edge is cut into segments at the model nodes lying on it (hard points), so an
+ *   edge shared whole by two areas, or partly where areas meet at a node (a T-junction), shares
+ *   its segments.
+ * - Every segment gets a division count: at least 1, at least `length / maxSize` for an area
+ *   with `mesh.maxSize`, and raised until in every area opposite edges add up to the same
+ *   count (and at least `mesh.divisions`), adding each division where the elements are longest.
+ * - Segment nodes are evenly spaced; interior nodes come from transfinite interpolation of the
+ *   four edges, so a planar area stays planar.
+ *
+ * An area that ends with 1 × 1 divisions is its own element, with its own id. Generated ids:
+ * `"<a>~<b>:<k>"` for the k-th node from model node a to b on a segment, `"<areaId>:<i>,<j>"`
+ * for interior nodes and for elements.
+ *
+ * @param nodes - Model nodes
+ * @param shells - Shell areas, each with 4 valid corner nodes
+ */
+const meshShells = (nodes, shells) => {
+    const position = new Map(nodes.map((n) => [n.id, [n.x, n.y, n.z]]));
+    const segments = new Map();
+    const segmentOf = (p, q) => {
+        const forward = idKey(p) < idKey(q);
+        const [a, b] = forward ? [p, q] : [q, p];
+        const key = `${idKey(a)}|${idKey(b)}`;
+        let segment = segments.get(key);
+        if (!segment) {
+            const d = sub(position.get(b), position.get(a));
+            segment = { a, b, length: Math.hypot(d[0], d[1], d[2]), count: 1 };
+            segments.set(key, segment);
+        }
+        return { segment, forward };
+    };
+    // ---- Edges cut at the model nodes lying on them -------------------------------------
+    const plans = shells.map((shell) => {
+        const corners = shell.nodeIds;
+        const edges = [0, 1, 2, 3].map((i) => {
+            const [p, q] = [corners[i], corners[(i + 1) % 4]];
+            const [P, Q] = [position.get(p), position.get(q)];
+            const inside = [];
+            for (const node of nodes) {
+                if (corners.includes(node.id))
+                    continue;
+                const t = (0, exports.pointOnSegment)([node.x, node.y, node.z], P, Q);
+                if (t !== null)
+                    inside.push({ id: node.id, t });
+            }
+            const chain = [p, ...inside.sort((x, y) => x.t - y.t).map((h) => h.id), q];
+            return chain.slice(1).map((id, k) => segmentOf(chain[k], id));
+        });
+        const d = shell.mesh?.divisions;
+        return { shell, edges, minimum: d ? [d[0], d[1]] : [1, 1] };
+    });
+    // ---- Division counts ---------------------------------------------------------------
+    for (const { shell, edges } of plans) {
+        const maxSize = shell.mesh?.maxSize;
+        if (maxSize === undefined)
+            continue;
+        for (const { segment } of edges.flat()) {
+            segment.count = Math.max(segment.count, Math.ceil(segment.length / maxSize - 1e-9));
+        }
+    }
+    const total = (edge) => edge.reduce((sum, { segment }) => sum + segment.count, 0);
+    const growTo = (edge, target) => {
+        while (total(edge) < target) {
+            // one more division where the elements are longest
+            const longest = edge.reduce((best, s) => s.segment.length / s.segment.count > best.segment.length / best.segment.count ? s : best);
+            longest.segment.count++;
+        }
+    };
+    const balanced = () => plans.every(({ edges, minimum }) => total(edges[0]) === total(edges[2]) &&
+        total(edges[1]) === total(edges[3]) &&
+        total(edges[0]) >= minimum[0] &&
+        total(edges[1]) >= minimum[1]);
+    const maxPasses = 10 + 4 * plans.length;
+    for (let pass = 0; pass < maxPasses && !balanced(); pass++) {
+        for (const { edges, minimum } of plans) {
+            for (const [p, q, least] of [
+                [0, 2, minimum[0]],
+                [1, 3, minimum[1]],
+            ]) {
+                const target = Math.max(total(edges[p]), total(edges[q]), least);
+                growTo(edges[p], target);
+                growTo(edges[q], target);
+            }
+        }
+    }
+    if (!balanced()) {
+        return {
+            nodes: [],
+            elements: [],
+            problems: plans
+                .filter(({ edges }) => total(edges[0]) !== total(edges[2]) || total(edges[1]) !== total(edges[3]))
+                .map((p) => p.shell.id),
+        };
+    }
+    // ---- Nodes along segments ---------------------------------------------------------
+    const generated = [];
+    const generate = (id, at) => {
+        generated.push({ id, x: at[0], y: at[1], z: at[2] });
+        position.set(id, at);
+    };
+    for (const segment of segments.values()) {
+        const [A, B] = [position.get(segment.a), position.get(segment.b)];
+        segment.chain = [segment.a];
+        for (let k = 1; k < segment.count; k++) {
+            const id = `${String(segment.a)}~${String(segment.b)}:${k}`;
+            generate(id, lerp(A, B, k / segment.count));
+            segment.chain.push(id);
+        }
+        segment.chain.push(segment.b);
+    }
+    // ---- Areas ------------------------------------------------------------------------
+    const elements = [];
+    const problems = [];
+    for (const { shell, edges } of plans) {
+        // node ids along each edge, corner i → corner i + 1
+        const along = edges.map((edge) => edge.flatMap(({ segment, forward }, k) => {
+            const chain = forward ? segment.chain : [...segment.chain].reverse();
+            return k === 0 ? chain : chain.slice(1);
+        }));
+        const n1 = along[0].length - 1;
+        const n2 = along[1].length - 1;
+        const { mesh: _mesh, ...properties } = shell;
+        if (n1 === 1 && n2 === 1) {
+            elements.push({ ...properties, areaId: shell.id });
+            continue;
+        }
+        const bottom = along[0];
+        const right = along[1];
+        const top = [...along[2]].reverse();
+        const left = [...along[3]].reverse();
+        const fractions = (ids) => {
+            const out = [0];
+            for (let k = 1; k < ids.length; k++) {
+                const d = sub(position.get(ids[k]), position.get(ids[k - 1]));
+                out.push(out[k - 1] + Math.hypot(d[0], d[1], d[2]));
+            }
+            return out.map((v) => v / out[out.length - 1]);
+        };
+        const [uB, uT, vL, vR] = [bottom, top, left, right].map(fractions);
+        const C = shell.nodeIds.map((id) => position.get(id));
+        const grid = [];
+        for (let i = 0; i <= n1; i++) {
+            grid.push([]);
+            for (let j = 0; j <= n2; j++) {
+                if (j === 0)
+                    grid[i].push(bottom[i]);
+                else if (j === n2)
+                    grid[i].push(top[i]);
+                else if (i === 0)
+                    grid[i].push(left[j]);
+                else if (i === n1)
+                    grid[i].push(right[j]);
+                else {
+                    // transfinite (Coons) interpolation of the four edges
+                    const v0 = (vL[j] + vR[j]) / 2;
+                    const u0 = (uB[i] + uT[i]) / 2;
+                    const u = (1 - v0) * uB[i] + v0 * uT[i];
+                    const v = (1 - u0) * vL[j] + u0 * vR[j];
+                    const [B, T, L, R] = [bottom[i], top[i], left[j], right[j]].map((id) => position.get(id));
+                    const at = [0, 1, 2].map((k) => (1 - v) * B[k] +
+                        v * T[k] +
+                        (1 - u) * L[k] +
+                        u * R[k] -
+                        ((1 - u) * (1 - v) * C[0][k] +
+                            u * (1 - v) * C[1][k] +
+                            u * v * C[2][k] +
+                            (1 - u) * v * C[3][k]));
+                    const id = `${String(shell.id)}:${i},${j}`;
+                    generate(id, at);
+                    grid[i].push(id);
+                }
+            }
+        }
+        let ok = true;
+        for (let j = 0; j < n2; j++) {
+            for (let i = 0; i < n1; i++) {
+                const nodeIds = [grid[i][j], grid[i + 1][j], grid[i + 1][j + 1], grid[i][j + 1]];
+                if (!(0, shell_axes_1.isValidShellQuad)(nodeIds.map((id) => position.get(id))))
+                    ok = false;
+                elements.push({
+                    ...properties,
+                    id: `${String(shell.id)}:${i},${j}`,
+                    nodeIds,
+                    areaId: shell.id,
+                });
+            }
+        }
+        if (!ok)
+            problems.push(shell.id);
+    }
+    // ---- Generated ids must be new ----------------------------------------------------
+    const taken = new Set(nodes.map((n) => idKey(n.id)));
+    for (const node of generated) {
+        if (taken.has(idKey(node.id))) {
+            const owner = elements.find((e) => e.nodeIds.includes(node.id));
+            if (owner && !problems.includes(owner.areaId))
+                problems.push(owner.areaId);
+        }
+        taken.add(idKey(node.id));
+    }
+    const shellIds = new Set();
+    for (const element of elements) {
+        const key = idKey(element.id);
+        if (shellIds.has(key) && !problems.includes(element.areaId))
+            problems.push(element.areaId);
+        shellIds.add(key);
+    }
+    return { nodes: generated, elements, problems };
+};
+exports.meshShells = meshShells;
+//# sourceMappingURL=mesh-shells.js.map

package/dist/fea/frame/mesh/mesh-shells.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mesh-shells.js","sourceRoot":"","sources":["../../../../src/fea/frame/mesh/mesh-shells.ts"],"names":[],"mappings":";;;AAAA,uEAA6E;AAE7E,6DAA2E;AAmC3E,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAChF,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1E,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,CAAS,EAAQ,EAAE,CAAC;IAClD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACxB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACxB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;CACzB,CAAC;AACF,0EAA0E;AAC1E,MAAM,KAAK,GAAG,CAAC,EAAW,EAAE,EAAE,CAAC,GAAG,OAAO,EAAE,IAAI,MAAM,CAAC,EAAE,CAAC,EAAE,CAAC;AAE5D,mGAAmG;AAC5F,MAAM,cAAc,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,CAAO,EAAE,SAAS,GAAG,KAAK,EAAiB,EAAE;IAC5F,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrB,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAAE,OAAO,IAAI,CAAC;IAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACxB,MAAM,SAAS,GAAG,wCAAuB,GAAG,CAAC,CAAC;IAC9C,MAAM,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IACjC,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,CAAC,CAAC;IAC1D,IAAI,CAAC,GAAG,MAAM,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM;QAAE,OAAO,IAAI,CAAC;IAC9C,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClC,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;AACpE,CAAC,CAAC;AAXW,QAAA,cAAc,kBAWzB;AAEF;;;;;;;;;;;;;;;;;;GAkBG;AACI,MAAM,UAAU,GAAG,CAAC,KAAkB,EAAE,MAAmB,EAAmB,EAAE;IACrF,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAgB,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAmB,CAAC;IAC5C,MAAM,SAAS,GAAG,CAAC,CAAU,EAAE,CAAU,EAA0C,EAAE;QACnF,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACtC,IAAI,OAAO,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChC,IAAI,CAAC,OAAO,EAAE,CAAC;YACb,MAAM,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;YAClD,OAAO,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;YACnE,QAAQ,CAAC,GAAG,CAAC,GAAG,EAAE,OAAO,CAAC,CAAC;QAC7B,CAAC;QACD,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC;IAC9B,CAAC,CAAC;IAEF,wFAAwF;IACxF,MAAM,KAAK,GAAkB,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE;QAChD,MAAM,OAAO,GAAG,KAAK,CAAC,OAAO,CAAC;QAC9B,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;YACnC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAClD,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;YACpD,MAAM,MAAM,GAAiC,EAAE,CAAC;YAChD,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;gBACzB,IAAI,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;oBAAE,SAAS;gBACxC,MAAM,CAAC,GAAG,IAAA,sBAAc,EAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBACzD,IAAI,CAAC,KAAK,IAAI;oBAAE,MAAM,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;YAClD,CAAC;YACD,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC3E,OAAO,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAChE,CAAC,CAAC,CAAC;QACH,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,SAAS,CAAC;QAChC,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;IAC9D,CAAC,CAAC,CAAC;IAEH,uFAAuF;IACvF,KAAK,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,IAAI,KAAK,EAAE,CAAC;QACrC,MAAM,OAAO,GAAG,KAAK,CAAC,IAAI,EAAE,OAAO,CAAC;QACpC,IAAI,OAAO,KAAK,SAAS;YAAE,SAAS;QACpC,KAAK,MAAM,EAAE,OAAO,EAAE,IAAI,KAAK,CAAC,IAAI,EAAE,EAAE,CAAC;YACvC,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,EAAE,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,GAAG,OAAO,GAAG,IAAI,CAAC,CAAC,CAAC;QACtF,CAAC;IACH,CAAC;IACD,MAAM,KAAK,GAAG,CAAC,IAAkC,EAAE,EAAE,CACnD,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAC5D,MAAM,MAAM,GAAG,CAAC,IAAkC,EAAE,MAAc,EAAE,EAAE;QACpE,OAAO,KAAK,CAAC,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC;YAC5B,mDAAmD;YACnD,MAAM,OAAO,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CACtC,CAAC,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CACzF,CAAC;YACF,OAAO,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC;QAC1B,CAAC;IACH,CAAC,CAAC;IACF,MAAM,QAAQ,GAAG,GAAG,EAAE,CACpB,KAAK,CAAC,KAAK,CACT,CAAC,EAAE,KAAK,EAAE,OAAO,EAAE,EAAE,EAAE,CACrB,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACnC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACnC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC;QAC7B,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,CAChC,CAAC;IACJ,MAAM,SAAS,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IACxC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,SAAS,IAAI,CAAC,QAAQ,EAAE,EAAE,IAAI,EAAE,EAAE,CAAC;QAC3D,KAAK,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,KAAK,EAAE,CAAC;YACvC,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,IAAI;gBAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;gBAClB,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;aACnB,EAAE,CAAC;gBACF,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;gBACjE,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gBACzB,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;YAC3B,CAAC;QACH,CAAC;IACH,CAAC;IACD,IAAI,CAAC,QAAQ,EAAE,EAAE,CAAC;QAChB,OAAO;YACL,KAAK,EAAE,EAAE;YACT,QAAQ,EAAE,EAAE;YACZ,QAAQ,EAAE,KAAK;iBACZ,MAAM,CACL,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAC1F;iBACA,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;SAC1B,CAAC;IACJ,CAAC;IAED,sFAAsF;IACtF,MAAM,SAAS,GAAgB,EAAE,CAAC;IAClC,MAAM,QAAQ,GAAG,CAAC,EAAW,EAAE,EAAQ,EAAE,EAAE;QACzC,SAAS,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QACrD,QAAQ,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACvB,CAAC,CAAC;IACF,KAAK,MAAM,OAAO,IAAI,QAAQ,CAAC,MAAM,EAAE,EAAE,CAAC;QACxC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAE,CAAC,CAAC;QACpE,OAAO,CAAC,KAAK,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;QAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,MAAM,EAAE,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;YAC5D,QAAQ,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;YAC5C,OAAO,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACzB,CAAC;QACD,OAAO,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;IAChC,CAAC;IAED,sFAAsF;IACtF,MAAM,QAAQ,GAAmB,EAAE,CAAC;IACpC,MAAM,QAAQ,GAAc,EAAE,CAAC;IAC/B,KAAK,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,IAAI,KAAK,EAAE,CAAC;QACrC,oDAAoD;QACpD,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAC/B,IAAI,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,CAAC,EAAE,EAAE;YACvC,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,KAAM,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAM,CAAC,CAAC,OAAO,EAAE,CAAC;YACvE,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC1C,CAAC,CAAC,CACH,CAAC;QACF,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/B,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,UAAU,EAAE,GAAG,KAAK,CAAC;QAC7C,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;YACzB,QAAQ,CAAC,IAAI,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,EAAE,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC;YACnD,SAAS;QACX,CAAC;QAED,MAAM,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC;QACpC,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,CAAC,GAAc,EAAE,EAAE;YACnC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACpC,MAAM,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;gBAChE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtD,CAAC;YACD,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;QACjD,CAAC,CAAC;QACF,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QACnE,MAAM,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;QACvD,MAAM,IAAI,GAAgB,EAAE,CAAC;QAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,IAAI,CAAC,KAAK,CAAC;oBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;qBAChC,IAAI,CAAC,KAAK,EAAE;oBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;qBACnC,IAAI,CAAC,KAAK,CAAC;oBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;qBACnC,IAAI,CAAC,KAAK,EAAE;oBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;qBACrC,CAAC;oBACJ,sDAAsD;oBACtD,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBAC/B,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBAC/B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBACxC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBACxC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAC7D,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAC1B,CAAC;oBACF,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CACtB,CAAC,CAAC,EAAE,EAAE,CACJ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBACd,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBACR,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBACd,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBACR,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;4BAC1B,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;4BACrB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;4BACf,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CACnB,CAAC;oBACV,MAAM,EAAE,GAAG,GAAG,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;oBAC3C,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;oBACjB,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;gBACnB,CAAC;YACH,CAAC;QACH,CAAC;QACD,IAAI,EAAE,GAAG,IAAI,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC5B,MAAM,OAAO,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBACjF,IAAI,CAAC,IAAA,6BAAgB,EAAC,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;oBAAE,EAAE,GAAG,KAAK,CAAC;gBAC1E,QAAQ,CAAC,IAAI,CAAC;oBACZ,GAAG,UAAU;oBACb,EAAE,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE;oBACnC,OAAO;oBACP,MAAM,EAAE,KAAK,CAAC,EAAE;iBACjB,CAAC,CAAC;YACL,CAAC;QACH,CAAC;QACD,IAAI,CAAC,EAAE;YAAE,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;IACnC,CAAC;IAED,sFAAsF;IACtF,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACrD,KAAK,MAAM,IAAI,IAAI,SAAS,EAAE,CAAC;QAC7B,IAAI,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;YAC9B,MAAM,KAAK,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;YAChE,IAAI,KAAK,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,KAAK,CAAC,MAAM,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAC7E,CAAC;QACD,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5B,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAU,CAAC;IACnC,KAAK,MAAM,OAAO,IAAI,QAAQ,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,KAAK,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC;QAC9B,IAAI,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,OAAO,CAAC,MAAM,CAAC;YAAE,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC3F,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,CAAC;IAED,OAAO,EAAE,KAAK,EAAE,SAAS,EAAE,QAAQ,EAAE,QAAQ,EAAE,CAAC;AAClD,CAAC,CAAC;AA5MW,QAAA,UAAU,cA4MrB"}