@theerapat-s28/ts-structural-eng-tools 2.1.0 → 2.2.0

package/CHANGELOG.md

@@ -0,0 +1,145 @@
+# Changelog
+
+All notable changes to `@theerapat-s28/ts-structural-eng-tools` are recorded here.
+
+The format follows [Keep a Changelog](https://keepachangelog.com/en/1.1.0/), and the package uses [semantic versioning](https://semver.org/): a **major** release can break existing code, a **minor** release adds features, a **patch** release fixes bugs.
+
+## [Unreleased]
+
+## [2.2.0] — 2026-10-24
+
+### Added
+
+- **Partial-fixity end releases**: any member end DOF in `releases.start` / `releases.end` (`u1`, `u2`, `u3`, `r1`, `r2`, `r3`) can be a spring stiffness in kN/mm or kN·m/rad instead of `true`, modelling semi-rigid connections. The spring acts in series with the member end; `0` is the same as `true` (fully released). Springs never make a member a mechanism (307).
+- **Release deformations**: every load case and combination result has `releaseDeformations`, the node-minus-member displacement across each released or spring end DOF of the members that have any — the rotation across a hinge, or the rotation of a semi-rigid connection.
+- **Automatic joint restraint**: joint rotations that no member, support or spring resists (for example every joint of a pin-jointed truss) are held automatically instead of throwing 304, with no effect on the results. They are listed in `calculationDetails.autoRestrainedDofs`. A load with a moment about such a rotation still throws 304, now naming the load case.
+- Types `DofRelease`, `ReleaseDeformation`, `MemberReleaseDeformations` and `AutoRestrainedDof`, exported from `/fea` and the package root.
+
+### Changed
+
+- `MemberEndRelease` fields are typed `DofRelease` (`boolean | number`) instead of `boolean`. Release values that are not `true`, `false` or a finite number ≥ 0 are now rejected with 314 (`details.field` such as `"releases.end.r3"`); before, any truthy value counted as a release.
+
+## [2.1.1] — 2026-10-23
+
+### Fixed
+
+- The package root now re-exports `distributeAreaLoads`, the `TrapezoidalMemberLoad`, `TorqueMemberLoad`, `NullArea`, `AreaLoad`, `AreaDistribution` and `AreaLoadDirection` types, and the `AREA_GEOMETRY_TOLERANCE` constant. Since 2.1.0 they were reachable only through the `/fea` and `/core` subpaths.
+- The quick-start example on the docs landing page used `fc` instead of `fc_`.
+
+### Documentation
+
+- README: added a quick-start example, the `core` and `utils` exports, error codes, release steps and the license.
+
+## [2.1.0] — 2026-10-23
+
+### Added
+
+- **Trapezoidal member load** (`type: "trapezoidal"`): a load varying linearly from `w1` at `start` to `w2` at `end`. Member diagrams stay exact for it.
+- **Torque member load** (`type: "torque"`): a uniform torque per unit length about the member's local 1 axis.
+- **Null areas** (`model.areas`) and **area loads** (`model.areaLoads`, kN/m²): 3- or 4-node areas that add no stiffness and pass their loads to the edge members, distributed `"twoWay"`, `"oneWay"` or `"cantilever"`.
+- `distributeAreaLoads(model)` returns the member loads generated from the areas, for checking or previewing (from `/fea`).
+- `FrameAnalysisError` codes 315 (area geometry), 316 (an area edge without members) and 317 (supports that do not fit the distribution).
+
+### Changed
+
+- Package description and README now present the library as a headless analysis engine for front-end applications.
+
+## [2.0.1] — 2026-10-11
+
+Maintenance release with no changes to the library code.
+
+## [2.0.0] — 2026-10-11
+
+### Breaking changes
+
+- `FrameModel` now requires `loadCases`. Loads move to the top-level lists `nodalLoads` and `memberLoads`, and every load names its `loadCase`.
+- `analyzeFrame` returns `{ cases, combinations, envelope, … }` keyed by id, instead of a single result.
+
+### Added
+
+- **3D frame analysis** (`/fea` subpath) by the direct stiffness method: uniform and point member loads in local or global directions, end releases, fixed and spring supports, ETABS-style local axes.
+- **Load cases and combinations**: a self-weight multiplier per load case, combinations by superposition (including nested combinations), and an `envelope` reporting the governing combination.
+- Member diagrams of forces (`N`, `V2`, `V3`, `T`, `M2`, `M3`) and deflections (`u1`–`u3` local, `ux`–`uz` global) along every member.
+- `createFrameAnalysis` to compute member diagrams on demand, and `options.diagrams: "extremes"` for smaller results.
+- `validateFrameModel` reports every problem in a model, with `details` locating it. `FrameAnalysisError` codes 301–314.
+- `memberDemandEnvelope` turns analysis results into `MuPositive`, `MuNegative`, `Vu` and `Nu` for the RC capacity functions.
+- `rectBeamToSectionProperties` (in `/rc`) converts an RC section into frame section properties, with the ACI 318-19 6.6.3.1.1 inertia modifiers and a unit weight for self-weight.
+- `beam3dLocalStiffness`, `beam3dLocalStiffnessRightHanded` and `memberAxes`.
+- Matrix and skyline solver helpers in `/utils`. Large models solve much faster (an 847-node model in about 1.5 s instead of 16.6 s).
+- ES-module build in `dist/esm/` alongside CommonJS, with `"sideEffects": false` so bundlers can drop unused code.
+
+## [1.5.0] — 2026-08-17
+
+### Added
+
+- **Side plate shear strengthening**, with two models because ACI 318-19 has no provisions for steel side plates:
+  - `sidePlateShearCapacityByWebYielding` (plate as a supplementary web, AISC 360 G2.1)
+  - `sidePlateShearCapacityByTensionTie` (plate as a tie across the crack, ACI 440.2R-17 §11.4 with steel in place of FRP)
+  - `compareSidePlateShearCapacity` runs both and reports the lower as governing.
+- `hasSidePlate` type guard.
+- `RCDesignError` codes 205 (inconsistent side plate geometry) and 206 (no shear headroom left for strengthening).
+
+## [1.4.0] — 2026-08-14
+
+### Changed
+
+- Doubly reinforced sections are solved by strain compatibility. Compression steel below the neutral axis no longer throws: it is treated as tension steel and reported with a warning citing ACI 318-19 22.2.1.2.
+- Steel plate jacketing measures every depth from the extreme compression fibre. A top plate now turns a singly reinforced section into a doubly reinforced one.
+- `calculationDetails` reports `c`, `a`, `beta1`, `d`, `As` and `ro` for every section. `ro_` uses the same `b·d` basis as `ro`.
+- `rectBeamBarLayout` takes `topBars` as optional.
+
+### Deprecated
+
+- `RCDesignError` code 101 (`COMP_STEEL_DISTANCE_OVER_AC`) is no longer thrown. It will be removed in the next major release.
+
+## [1.3.0] — 2026-08-14
+
+### Added
+
+- API reference site, published on GitHub Pages: <https://theerapat-s28.github.io/ts-structural-eng/>.
+- `homepage`, `repository` and `bugs` links in the package metadata.
+
+## [1.2.0] — 2026-08-13
+
+### Added
+
+- **Plate-to-concrete interface bolts** for steel plate jacketing:
+  - `plateInterfaceShearFlow` (shear flow `q = V·Q/I`, cracked or uncracked section)
+  - `boltShearCapacity` (ACI 318-19 17.7.1.2b)
+  - `plateInterfaceBoltRequirement` (required bolt pitch and count).
+- Top and bottom jacketing plates are each optional, with the `hasTopPlate` and `hasBottomPlate` type guards.
+- `RCDesignError` codes 201–204 for strengthening.
+
+### Fixed
+
+- The RC shear functions are now exported from the package root, not only from `/rc`.
+
+## [1.1.0] — 2026-07-19
+
+### Added
+
+- RC shear design per ACI 318-19: `concreteShearCapacity`, `stirrupShearCapacity` and `checkStirrupRequirement`, with `RCDesignError` code 105.
+- `rectBeamBarLayout`: arranges bar groups into layers that meet the ACI 318-19 spacing rules and returns `d` and `d_`.
+
+## [1.0.1] — 2026-07-19
+
+### Fixed
+
+- npm packaging: scoped package name `@theerapat-s28/ts-structural-eng-tools`, a build output that matches the `exports` map, and path aliases rewritten in the published JavaScript.
+
+## 1.0.0
+
+First version, under the unscoped name `ts-structural-eng-tools`: `rectBeamMomentCapacity` (singly and doubly reinforced rectangular beams per ACI 318-19), steel plate jacketing, and concrete property helpers.
+
+[Unreleased]: https://github.com/theerapat-s28/ts-structural-eng/compare/v2.2.0...HEAD
+[2.2.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v2.1.1...v2.2.0
+[2.1.1]: https://github.com/theerapat-s28/ts-structural-eng/compare/v2.1.0...v2.1.1
+[2.1.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v2.0.1...v2.1.0
+[2.0.1]: https://github.com/theerapat-s28/ts-structural-eng/compare/v2.0.0...v2.0.1
+[2.0.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v1.5.0...v2.0.0
+[1.5.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v1.4.0...v1.5.0
+[1.4.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v1.3.0...v1.4.0
+[1.3.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v1.2.0...v1.3.0
+[1.2.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v1.1.0...v1.2.0
+[1.1.0]: https://github.com/theerapat-s28/ts-structural-eng/compare/v1.0.1...v1.1.0
+[1.0.1]: https://github.com/theerapat-s28/ts-structural-eng/releases/tag/v1.0.1

package/README.md

@@ -59,6 +59,58 @@ Or install this repository's dependencies to work on the library itself:
 pnpm install
 ```
 
+## Quick start
+
+Analyse a simply supported RC beam under self-weight and live load, envelope the design demands, and check them against the section's flexural capacity:
+
+```ts
+import { analyzeFrame, memberDemandEnvelope } from "@theerapat-s28/ts-structural-eng-tools/fea";
+import type { FrameModel } from "@theerapat-s28/ts-structural-eng-tools/fea";
+import {
+  rectBeamMomentCapacity,
+  rectBeamToSectionProperties,
+} from "@theerapat-s28/ts-structural-eng-tools/rc";
+
+// 300×600 beam, f'c 28 MPa, fy 420 MPa, As 1800 mm²
+const beam = { Es: 200_000, fc_: 28, fy: 420, As: 1800, b: 300, h: 600, d: 540 };
+const { sectionProperties } = rectBeamToSectionProperties(beam, { member: "beam" });
+
+const model: FrameModel = {
+  nodes: [
+    { id: "A", x: 0, y: 0, z: 0 },
+    { id: "B", x: 6000, y: 0, z: 0 }, // mm
+  ],
+  members: [{ id: "AB", startNodeId: "A", endNodeId: "B", section: sectionProperties }],
+  restraints: [
+    { nodeId: "A", ux: true, uy: true, uz: true, rx: true },
+    { nodeId: "B", uy: true, uz: true },
+  ],
+  loadCases: [
+    { id: "D", type: "dead", selfWeight: 1 },
+    { id: "L", type: "live" },
+  ],
+  memberLoads: [{ loadCase: "L", type: "udl", memberId: "AB", direction: "Z", w: -15 }], // kN/m
+  combinations: [
+    {
+      id: "1.2D+1.6L",
+      terms: [
+        { case: "D", factor: 1.2 },
+        { case: "L", factor: 1.6 },
+      ],
+    },
+  ],
+};
+
+const result = analyzeFrame(model);
+const demand = memberDemandEnvelope(result, "AB");
+// demand.MuPositive → 131.328 (kN·m), demand.Vu → 87.552 (kN), governed by "1.2D+1.6L"
+
+const { phiMn, warnings } = rectBeamMomentCapacity(beam);
+// phiMn → 331.39 (kN·m), warnings → []
+```
+
+Each module is also its own subpath export (`/fea`, `/rc`, `/strengthening`, `/core`, `/utils`), so a front end can import only what it uses. The package root re-exports the most-used functions and types from every module.
+
 ## Scripts
 
 | Command                | Description                                     |
@@ -80,6 +132,7 @@ src/
 ├── rc/                   # Reinforced concrete design
 │   ├── index.ts          # Module exports
 │   ├── general.ts        # Concrete properties (beta1, Ec, unit conversion)
+│   ├── rc-bar-layout.ts  # Bar groups → layers satisfying ACI spacing, giving d and d'
 │   ├── rc-beam-design.ts # Beam moment capacity (singly/doubly reinforced) and shear capacity (Vc, Vs, stirrup requirement)
 │   └── rc-section-properties.ts # Section properties (E, G, A, I22, I33, J) for frame analysis
 ├── strengthening/        # Strengthening methods
@@ -92,17 +145,17 @@ src/
 │   └── frame/            # 3D frame analysis (direct stiffness method)
 │       ├── index.ts
 │       ├── area/             # Null areas: geometry, edge-member chains, area-load distribution
-│       ├── elements/         # Reference beam stiffness matrix, fixed-end forces, end-release condensation
+│       ├── elements/         # Reference beam stiffness matrix, fixed-end forces, end releases and partial-fixity condensation
 │       ├── transformation/   # ETABS-style local axes and the 12×12 transformation
-│       ├── solver/           # System build + skyline factorisation, per-case solve, combinations, analyzeFrame
+│       ├── solver/           # System build (with automatic joint restraint) + skyline factorisation, per-case solve, combinations, analyzeFrame
 │       ├── validation/       # validateFrameModel: every input problem, with its location
 │       ├── post/             # Result presentation, member force diagrams, envelopes, design demands
 │       └── types/            # Model (incl. null areas, load cases, combinations) and result types
 ├── core/                 # Shared infrastructure
 │   ├── index.ts
 │   ├── constants/        # RC design constants (ACI 318) and frame-analysis tolerances
-│   ├── errors/           # Custom error classes
-│   └── types/            # Shared type definitions
+│   ├── errors/           # RCDesignError and FrameAnalysisError with their coded entries
+│   └── types/            # Shared types: RC sections, plate jacketing, section properties, warnings, units
 └── utils/                # Helper utilities
     ├── index.ts
     ├── math.ts           # Quadratic solver, rounding
@@ -122,6 +175,9 @@ html-docs/                # HTML API reference, deployed to GitHub Pages
 ├── rc/index.html
 ├── strengthening/index.html
 └── fea/index.html
+ROADMAP.md                # FEA roadmap and implementation plan (not published)
+.agents/skills/update-version/ # AI agent skill: bump the version and roll the changelog
+.claude/skills/npm-publish/ # AI agent skill: release checklist for publishing to npm
 ```
 
 ## TypeScript path aliases
@@ -174,32 +230,55 @@ ACI 440.2R-17 11.4.2 (spacing, which defers to the ACI 318 limits) and 11.4.3 (`
 
 Linear-elastic 3D frame analysis by the direct stiffness method, imported from `@theerapat-s28/ts-structural-eng-tools/fea`. Global axes X, Y, Z (Z up) and ETABS-style local axes 1, 2, 3; inputs in mm, MPa, kN and kN·m, results in kN, kN·m, mm and rad.
 
-- **`analyzeFrame(model, options?)`** — Analyses every load case in `model.loadCases` (dead, superimposed dead, live, wind…, each with an optional `selfWeight` multiplier) and forms every combination in `model.combinations` by superposition. The loads are flat top-level lists, `model.nodalLoads` and `model.memberLoads` (uniform, trapezoidal and point loads in local or global directions, and torques about the member axis), each load tagged with its `loadCase`, so large models stay readable. Area loads (`model.areaLoads`, kN/m²) on null areas (`model.areas`) are first distributed to the members along each area's edges — see `distributeAreaLoads`. A combination is a list of factored terms `{ case, factor }` or `{ combination, factor }`; factors are entered by the user, and nested combinations are expanded. Supports end releases and fixed or spring supports. The stiffness matrix is factorised once for all cases. Returns `cases` and `combinations` keyed by id — each with displacements, reactions, member end forces, member diagrams along each member (forces `N`, `V2`, `V3`, `T`, `M2`, `M3` and deflections `u1`–`u3` local / `ux`–`uz` global, ready for plotting) and an equilibrium check — and an `envelope` of max/min over the combinations with the governing combination. Combination diagrams are rebuilt exactly, not summed from case extremes. The stiffness is stored in skyline (banded) form with nodes in reverse Cuthill–McKee order, so a building of ~850 nodes runs in about 1.5 s. Model and result are plain JSON, so it runs in a Web Worker. Invalid models, load cases or combinations and mechanisms throw `FrameAnalysisError` (codes 301–317) with `details` locating the problem.
+- **`analyzeFrame(model, options?)`** — Analyses every load case in `model.loadCases` (dead, superimposed dead, live, wind…, each with an optional `selfWeight` multiplier) and forms every combination in `model.combinations` by superposition. The loads are flat top-level lists, `model.nodalLoads` and `model.memberLoads` (uniform, trapezoidal and point loads in local or global directions, and torques about the member axis), each load tagged with its `loadCase`, so large models stay readable. Area loads (`model.areaLoads`, kN/m²) on null areas (`model.areas`) are first distributed to the members along each area's edges — see `distributeAreaLoads`. A combination is a list of factored terms `{ case, factor }` or `{ combination, factor }`; factors are entered by the user, and nested combinations are expanded. Supports end releases — full (`true`) or partial fixity (a spring stiffness in series with the member end) on any of the six end DOFs, with the deformation across each release reported in `releaseDeformations` — and fixed or spring supports. Joint rotations that no member resists, as at pin-jointed truss joints, are held automatically and listed in `calculationDetails.autoRestrainedDofs`. The stiffness matrix is factorised once for all cases. Returns `cases` and `combinations` keyed by id — each with displacements, reactions, member end forces, member diagrams along each member (forces `N`, `V2`, `V3`, `T`, `M2`, `M3` and deflections `u1`–`u3` local / `ux`–`uz` global, ready for plotting) and an equilibrium check — and an `envelope` of max/min over the combinations with the governing combination. Combination diagrams are rebuilt exactly, not summed from case extremes. The stiffness is stored in skyline (banded) form with nodes in reverse Cuthill–McKee order, so a building of ~850 nodes runs in about 1.5 s. Model and result are plain JSON, so it runs in a Web Worker. Invalid models, load cases or combinations and mechanisms throw `FrameAnalysisError` (codes 301–317) with `details` locating the problem.
 - **`createFrameAnalysis(model, options?)`** — Same analysis, keeping the solution: returns `{ result, memberDiagram(resultId, memberId) }` so full member diagrams can be computed on demand. Combine with `options.diagrams: "extremes"`, which keeps diagram extremes (and the envelope) but drops the station arrays, for results about 4× smaller on large models.
-- **`validateFrameModel(model)`** — Checks a model without analysing it and returns every problem, each with `details` (`nodeId`, `memberId`, `areaId`, `edge`, `loadCaseId`, `combinationId`, `field`, `dof`) for highlighting in a UI. `analyzeFrame` runs it first and throws the first problem as a `FrameAnalysisError` whose `issues` lists all of them; non-finite numbers, negative springs and unknown option values are rejected (314).
+- **`validateFrameModel(model)`** — Checks a model without analysing it and returns every problem, each with `details` (`nodeId`, `memberId`, `areaId`, `edge`, `loadCaseId`, `combinationId`, `field`, `dof`) for highlighting in a UI. `analyzeFrame` runs it first and throws the first problem as a `FrameAnalysisError` whose `issues` lists all of them; non-finite numbers, negative springs, invalid end-release values and unknown option values are rejected (314).
 - **`distributeAreaLoads(model)`** — Returns `{ memberLoads, warnings }`: the member loads that the model's null areas send to their edge members, exactly what `analyzeFrame` adds, for checking or previewing. A null area (`{ id, nodeIds, distribution, supportEdges? }`) has 3 or 4 nodes forming a planar convex polygon and adds no stiffness. `"twoWay"` sends each point to its nearest edge (the corner-bisector split, any convex shape); `"oneWay"` spans strips between the listed `supportEdges`, parallel to the unsupported edges, half to each end; `"cantilever"` (rectangles) puts `q·c` and a torque `q·c²/2` on its one support edge. Area loads (`{ loadCase, areaId, direction: "3" | "X" | "Y" | "Z", q }`) are per unit true area, `"3"` being the area's normal. Edges split by intermediate nodes are cut at them. Area problems throw 315 (geometry), 316 (an edge without members) or 317 (supports that do not fit the distribution).
 - **`memberDemandEnvelope(result, memberId, options?)`** — Reduces a member's diagrams to `MuPositive`, `MuNegative`, `Vu` (optionally at distance `d` from the supports) and `Nu`, enveloped over all combinations (or the `combinations` / `cases` you select), with the `governing` combination of each, for the RC capacity functions.
 - **`rectBeamToSectionProperties(section, options?)`** (from `rc`) — Converts a rectangular RC section into `SectionProperties` (`E`, `G`, `A`, `I22`, `I33`, `J`, `unitWeight`) with the ACI 318-19 6.6.3.1.1 inertia modifiers (0.35 beams, 0.70 columns) and a unit weight for self-weight (default 24 kN/m³, common practice rather than an ACI value).
 - **`beam3dLocalStiffness(E, G, A, I22, I33, J, L)`** — The 12×12 local stiffness matrix exactly as in the reference "Complete Beam Stiffness Matrix". Its bending rotations R2 and R3 follow the opposite sign to the right-hand rule; **`beam3dLocalStiffnessRightHanded`** returns the right-handed form that `analyzeFrame` uses.
 - **`memberAxes(start, end, angle?)`** — The ETABS-style local axes of a member as a direction-cosine matrix.
 
-### General Utilities
+### Concrete properties (from `rc`)
 
 - **`concreteBeta(fc_)`** — ACI 318 β₁ factor
 - **`concreteElasticModulus(fc_)`** — Ec = 4700√f'c (MPa)
 - **`psiToMpa(psi)`** — Unit conversion
 
+### Errors and warnings (from `core`)
+
+A design that is valid but breaks a code provision returns normally with an entry in `warnings`, citing the clause (e.g. `"ACI318-19, 9.6.1.2"`). A calculation whose assumptions break down throws instead, never returning a silently wrong result:
+
+- **`RCDesignError`** — `code` and `message` from `Errors`: 1XX for RC beam design (e.g. 102, tensile steel not yielding; 105, `Vs` above the 22.5.1.2 limit), 2XX for strengthening and plate connections (201–206).
+- **`FrameAnalysisError`** — `code` from `FrameErrors` (301–317), `details` locating the problem, and `issues` listing every problem `validateFrameModel` found.
+
+`core` also exports the shared input types, the type guards `isSinglyReinforced`, `hasTopPlate`, `hasBottomPlate` and `hasSidePlate`, and the engineering constants (φ factors, ultimate concrete strain, inertia modifiers, frame-analysis tolerances).
+
+### Utilities (from `utils`)
+
+Low-level helpers the engine is built on, exported for reuse: `solveQuadratic`, `roundToDecimalPlaces`, `mergeWarnings`, dense matrix helpers (`zeros`, `identity`, `transpose`, `matMul`, `matVec`), Cholesky and general linear solvers (`choleskyFactor`, `choleskySolveFactored`, `choleskySolve`, `solveLinearSystem`), and the skyline solver (`createSkyline`, `skylineAdd`, `skylineGet`, `skylineSize`, `skylineFactor`, `skylineSolve`, `reverseCuthillMcKee`).
+
 ## Build
 
 ```bash
 pnpm run build
-node dist/index.js
 ```
 
-The build emits a CommonJS build in `dist/` (for `require`, with the `.d.ts` types) and an ES-module build in `dist/esm/` (for `import` and bundlers such as Vite); `package.json` `exports` routes each subpath to the right one, and `"sideEffects": false` lets bundlers drop unused code.
+The build emits a CommonJS build in `dist/` (for `require`, with the `.d.ts` types) and an ES-module build in `dist/esm/` (for `import` and bundlers such as Vite); `package.json` `exports` routes each subpath to the right one, and `"sideEffects": false` lets bundlers drop unused code. `pnpm run dev` runs [`examples/basic-usage.ts`](examples/basic-usage.ts) against the sources without building.
 
 ## Tests
 
 ```bash
-pnpm run test
+pnpm run test                                   # all tests
+pnpm vitest run tests/rc/rc-beam-design.test.ts # one file
 ```
+
+Tests live in `tests/`, mirroring `src/`, and check results against hand-calculated or textbook values.
+
+## Releasing
+
+Releases follow semantic versioning and are published to npm with tags `vX.Y.Z`. What changed in each version is in [CHANGELOG.md](CHANGELOG.md). Preparing a version (changelog, version number, checks) is described in [`.agents/skills/update-version/SKILL.md`](.agents/skills/update-version/SKILL.md), and publishing (release commit and tag, dry run, publish, push tags) in [`.claude/skills/npm-publish/SKILL.md`](.claude/skills/npm-publish/SKILL.md).
+
+## License
+
+[ISC](LICENSE)

package/dist/esm/fea/frame/elements/end-releases.js

@@ -0,0 +1,20 @@
+/** Release keys in local DOF order at one member end. */
+export const RELEASE_KEYS = ["u1", "u2", "u3", "r1", "r2", "r3"];
+/**
+ * Spring stiffness unit conversion per DOF of one node or member end:
+ * kN/mm → N/mm for translations, kN·m/rad → N·mm/rad for rotations.
+ */
+export const SPRING_UNIT_FACTOR = [1e3, 1e3, 1e3, 1e6, 1e6, 1e6];
+/**
+ * The released local end DOFs of a member, in local DOF order. `true` and `0` are full releases,
+ * a positive number is a partial-fixity spring; `false` and omitted DOFs stay rigidly connected.
+ * Expects values that passed `validateFrameModel`.
+ */
+export const memberEndReleases = (releases) => [releases?.start, releases?.end].flatMap((release, end) => RELEASE_KEYS.flatMap((key, k) => {
+    const value = release?.[key];
+    if (value === undefined || value === false)
+        return [];
+    const spring = value === true ? 0 : value * SPRING_UNIT_FACTOR[k];
+    return [{ index: 6 * end + k, spring }];
+}));
+//# sourceMappingURL=end-releases.js.map

package/dist/esm/fea/frame/elements/end-releases.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"end-releases.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/end-releases.ts"],"names":[],"mappings":"AAEA,yDAAyD;AACzD,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAU,CAAC;AAE1E;;;GAGG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAU,CAAC;AAQ1E;;;;GAIG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,QAA4B,EAAmB,EAAE,CAChF,CAAC,QAAQ,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAW,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,GAAG,EAAE,EAAE,CACnE,YAAY,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,EAAmB,EAAE;IAC/C,MAAM,KAAK,GAAG,OAAO,EAAE,CAAC,GAAG,CAAC,CAAC;IAC7B,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,KAAK;QAAE,OAAO,EAAE,CAAC;IACtD,MAAM,MAAM,GAAG,KAAK,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IAClE,OAAO,CAAC,EAAE,KAAK,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;AAC1C,CAAC,CAAC,CACH,CAAC"}

package/dist/esm/fea/frame/elements/static-condensation.js

@@ -1,55 +1,106 @@
 import { FrameAnalysisError, FrameErrors } from "../../../core/errors/frame-analysis.error.js";
 import { solveLinearSystem } from "../../../utils/matrix.js";
 /**
- * Statically condenses released local DOFs out of a member, so those DOFs carry no force.
+ * Condenses released local DOFs out of a member. Each released DOF r splits into a member-side
+ * value u_m and a node-side value u_n, tied by a spring s (s = 0 for a full release):
  *
- * With released set r and retained set c:
- * `k* = k_cc − k_cr·k_rr⁻¹·k_rc`, `f* = f_c − k_cr·k_rr⁻¹·f_r`.
- * The stiffness is condensed once; the returned function condenses the fixed-end forces of
- * each load case.
+ * ```
+ *        u_c     u_m      u_n
+ *     [ k_cc    k_cr       0 ]      A = k_rr + S,  S = diag(s)
+ *     [ k_rc    A         −S ]
+ *     [  0     −S          S ]
+ * ```
  *
- * @param k - 12×12 local stiffness (right-hand)
- * @param released - Local DOF indices (0–11) to release
- * @throws {FrameAnalysisError} 307 if the released set makes the member a mechanism
+ * Eliminating u_m leaves a 12×12 matrix in the node DOFs:
+ * `k*_cc = k_cc − k_cr·A⁻¹·k_rc`, `k*_cn = k_cr·A⁻¹·S`, `k*_nn = S·A⁻¹·k_rr`
+ * (the stable form of `S − S·A⁻¹·S`), and `f*_c = f_c − k_cr·A⁻¹·f_r`, `f*_n = S·A⁻¹·f_r`.
+ * The member side follows from `u_m = A⁻¹·(S·u_n − k_rc·u_c − f_r)`.
+ * With every s = 0 this is the classical condensation with zero rows and columns at r.
+ *
+ * @param k - 12×12 local stiffness (right-hand), N and mm
+ * @param releases - Released local DOFs (0–11) with their springs, N/mm or N·mm/rad
+ * @throws {FrameAnalysisError} 307 if the full releases make the member a mechanism
  */
-export const condenseMember = (k, released) => {
-    if (released.length === 0) {
-        return { k, condenseFixedEndForces: (fixed) => [...fixed] };
+export const condenseMember = (k, releases) => {
+    if (releases.length === 0) {
+        return {
+            k,
+            condenseFixedEndForces: (fixed) => [...fixed],
+            memberEndDisplacements: (nodeDisplacements) => [...nodeDisplacements],
+        };
     }
+    const released = releases.map((r) => r.index);
+    const s = releases.map((r) => r.spring);
     const retained = Array.from({ length: 12 }, (_, i) => i).filter((i) => !released.includes(i));
+    const nr = released.length;
+    const nc = retained.length;
     const sub = (rows, cols) => rows.map((r) => cols.map((c) => k[r][c]));
     const krr = sub(released, released);
     const krc = sub(released, retained);
-    const kcr = sub(retained, released);
-    const kcc = sub(retained, retained);
-    // [k_rr⁻¹·k_rc | k_rr⁻¹]
-    const rhs = released.map((_, i) => [...krc[i], ...released.map((__, j) => (i === j ? 1 : 0))]);
-    const solved = solveLinearSystem(krr, rhs);
-    if (!solved)
+    // A = k_rr + S, scaled to D·A·D with unit diagonal (D = diag(1/√A_ii)) so the singularity test
+    // does not depend on mixing N/mm and N·mm/rad; then A⁻¹·B = D·(D·A·D)⁻¹·(D·B).
+    const d = krr.map((row, i) => 1 / Math.sqrt(row[i] + s[i]));
+    const scaledA = krr.map((row, i) => row.map((v, j) => d[i] * (i === j ? v + s[i] : v) * d[j]));
+    // [A⁻¹·k_rc | A⁻¹·k_rr | A⁻¹]
+    const rhs = released.map((_, i) => [...krc[i], ...krr[i], ...released.map((__, j) => (i === j ? 1 : 0))].map((v) => d[i] * v));
+    const scaledSolution = solveLinearSystem(scaledA, rhs);
+    if (!scaledSolution)
         throw new FrameAnalysisError(FrameErrors.INVALID_RELEASE);
-    const nc = retained.length;
+    const solved = scaledSolution.map((row, i) => row.map((v) => d[i] * v));
+    const AinvKrc = solved.map((row) => row.slice(0, nc));
+    const AinvKrr = solved.map((row) => row.slice(nc, nc + nr));
+    const Ainv = solved.map((row) => row.slice(nc + nr));
     const kStar = k.map((row) => row.map(() => 0));
-    retained.forEach((r, ri) => {
+    retained.forEach((r) => {
         retained.forEach((c, ci) => {
-            let value = kcc[ri][ci];
-            released.forEach((_, j) => {
-                value -= kcr[ri][j] * solved[j][ci];
-            });
+            let value = k[r][c];
+            for (let j = 0; j < nr; j++)
+                value -= k[r][released[j]] * AinvKrc[j][ci];
             kStar[r][c] = value;
         });
+        released.forEach((c, j) => {
+            // k*_cn = k_cr·A⁻¹·S, mirrored into k*_nc
+            let value = 0;
+            for (let m = 0; m < nr; m++)
+                value += k[r][released[m]] * Ainv[m][j];
+            kStar[r][c] = value * s[j];
+            kStar[c][r] = value * s[j];
+        });
     });
-    // k_cr·k_rr⁻¹, applied to f_r for every load case
-    const kcrKrrInv = retained.map((_, ri) => released.map((__, m) => released.reduce((s, ___, j) => s + kcr[ri][j] * solved[j][nc + m], 0)));
+    released.forEach((r, i) => {
+        released.forEach((c, j) => {
+            // k*_nn = S·A⁻¹·k_rr, symmetrised against round-off
+            kStar[r][c] = (s[i] * AinvKrr[i][j] + s[j] * AinvKrr[j][i]) / 2;
+        });
+    });
+    // k_cr·A⁻¹, applied to f_r for every load case
+    const kcrAinv = retained.map((r) => released.map((_, m) => released.reduce((sum, rj, j) => sum + k[r][rj] * Ainv[j][m], 0)));
     return {
         k: kStar,
         condenseFixedEndForces: (fixed) => {
             const fStar = new Array(12).fill(0);
             const fr = released.map((r) => fixed[r]);
             retained.forEach((r, ri) => {
-                fStar[r] = fixed[r] - kcrKrrInv[ri].reduce((s, v, m) => s + v * fr[m], 0);
+                fStar[r] = fixed[r] - kcrAinv[ri].reduce((sum, v, m) => sum + v * fr[m], 0);
+            });
+            released.forEach((r, i) => {
+                fStar[r] = s[i] * Ainv[i].reduce((sum, v, m) => sum + v * fr[m], 0);
             });
             return fStar;
         },
+        memberEndDisplacements: (nodeDisplacements, fixed) => {
+            const u = [...nodeDisplacements];
+            released.forEach((r, i) => {
+                let value = 0;
+                for (let m = 0; m < nr; m++) {
+                    value += Ainv[i][m] * (s[m] * nodeDisplacements[released[m]] - fixed[released[m]]);
+                }
+                for (let ci = 0; ci < nc; ci++)
+                    value -= AinvKrc[i][ci] * nodeDisplacements[retained[ci]];
+                u[r] = value;
+            });
+            return u;
+        },
     };
 };
 //# sourceMappingURL=static-condensation.js.map

package/dist/esm/fea/frame/elements/static-condensation.js.map

@@ -1 +1 @@
-{"version":3,"file":"static-condensation.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/static-condensation.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AACxF,OAAO,EAAE,iBAAiB,EAAe,MAAM,mBAAmB,CAAC;AAQnE;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,QAAkB,EAAmB,EAAE;IAC/E,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC1B,OAAO,EAAE,CAAC,EAAE,sBAAsB,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAC9D,CAAC;IACD,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,IAAc,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1F,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,yBAAyB;IACzB,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/F,MAAM,MAAM,GAAG,iBAAiB,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC3C,IAAI,CAAC,MAAM;QAAE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,eAAe,CAAC,CAAC;IACvE,MAAM,EAAE,GAAG,QAAQ,CAAC,MAAM,CAAC;IAE3B,MAAM,KAAK,GAAW,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;QACzB,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;YACzB,IAAI,KAAK,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACxB,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;gBACxB,KAAK,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACtC,CAAC,CAAC,CAAC;YACH,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtB,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,kDAAkD;IAClD,MAAM,SAAS,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CACvC,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAC/F,CAAC;IAEF,OAAO;QACL,CAAC,EAAE,KAAK;QACR,sBAAsB,EAAE,CAAC,KAAK,EAAE,EAAE;YAChC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC5C,MAAM,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACzC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;gBACzB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5E,CAAC,CAAC,CAAC;YACH,OAAO,KAAK,CAAC;QACf,CAAC;KACF,CAAC;AACJ,CAAC,CAAC"}
+{"version":3,"file":"static-condensation.js","sourceRoot":"","sources":["../../../../../src/fea/frame/elements/static-condensation.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,uCAAuC,CAAC;AACxF,OAAO,EAAE,iBAAiB,EAAe,MAAM,mBAAmB,CAAC;AAmBnE;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,QAAyB,EAAmB,EAAE;IACtF,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC1B,OAAO;YACL,CAAC;YACD,sBAAsB,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;YAC7C,sBAAsB,EAAE,CAAC,iBAAiB,EAAE,EAAE,CAAC,CAAC,GAAG,iBAAiB,CAAC;SACtE,CAAC;IACJ,CAAC;IACD,MAAM,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAC9C,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACxC,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,MAAM,EAAE,GAAG,QAAQ,CAAC,MAAM,CAAC;IAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,MAAM,CAAC;IAC3B,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,IAAc,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1F,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACpC,+FAA+F;IAC/F,+EAA+E;IAC/E,MAAM,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,MAAM,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/F,8BAA8B;IAC9B,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAChC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAC3F,CAAC;IACF,MAAM,cAAc,GAAG,iBAAiB,CAAC,OAAO,EAAE,GAAG,CAAC,CAAC;IACvD,IAAI,CAAC,cAAc;QAAE,MAAM,IAAI,kBAAkB,CAAC,WAAW,CAAC,eAAe,CAAC,CAAC;IAC/E,MAAM,MAAM,GAAG,cAAc,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxE,MAAM,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACtD,MAAM,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;IAC5D,MAAM,IAAI,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;IAErD,MAAM,KAAK,GAAW,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE;QACrB,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;YACzB,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACzE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtB,CAAC,CAAC,CAAC;QACH,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACxB,0CAA0C;YAC1C,IAAI,KAAK,GAAG,CAAC,CAAC;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3B,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IACH,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACxB,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACxB,oDAAoD;YACpD,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAClE,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,+CAA+C;IAC/C,MAAM,OAAO,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CACjC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CACxF,CAAC;IAEF,OAAO;QACL,CAAC,EAAE,KAAK;QACR,sBAAsB,EAAE,CAAC,KAAK,EAAE,EAAE;YAChC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC5C,MAAM,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACzC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;gBACzB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC9E,CAAC,CAAC,CAAC;YACH,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;gBACxB,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtE,CAAC,CAAC,CAAC;YACH,OAAO,KAAK,CAAC;QACf,CAAC;QACD,sBAAsB,EAAE,CAAC,iBAAiB,EAAE,KAAK,EAAE,EAAE;YACnD,MAAM,CAAC,GAAG,CAAC,GAAG,iBAAiB,CAAC,CAAC;YACjC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;gBACxB,IAAI,KAAK,GAAG,CAAC,CAAC;gBACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,KAAK,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,iBAAiB,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrF,CAAC;gBACD,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE;oBAAE,KAAK,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,iBAAiB,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC1F,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YACf,CAAC,CAAC,CAAC;YACH,OAAO,CAAC,CAAC;QACX,CAAC;KACF,CAAC;AACJ,CAAC,CAAC"}

package/dist/esm/fea/frame/post/member-diagram.js

@@ -68,9 +68,9 @@ const internalForces = (startForce, loads, x, side) => {
  * smooth piece where N is at most quadratic and M at most cubic, and one 3-point Gauss panel per
  * interval integrates the (at most quartic) integrands exactly.
  *
- * The constants come from the member's own end displacements: end translations first, an end
- * rotation only where a translation is released. A released rotation (hinge) never enters, so
- * the curve is correct at hinged ends whose rotation differs from the node's.
+ * The constants a, b come from the member's own end translations (u at x = 0 and x = L), which
+ * differ from the node displacements across a released or spring end DOF; with them the curve
+ * is correct at hinged and semi-rigid ends whose rotation differs from the node's.
  */
 const displacementsOnGrid = (member, startForce, endDisplacements, loads, grid) => {
     const { E, A, I22, I33 } = member.section;
@@ -103,39 +103,16 @@ const displacementsOnGrid = (member, startForce, endDisplacements, loads, grid)
     }
     const [strain, k3, sk3, k2, sk2] = cumulative;
     const last = n - 1;
-    /** a, b of `u = a + b·x + P(x)` from the first independent pair of end conditions. */
-    const planeConstants = (translation, rotation, rotationSign, PL, // P(L)
-    dPL) => {
-        const { released } = member;
-        // rows [ca, cb | rhs] of ca·a + cb·b = rhs
-        const conditions = [];
-        if (!released[translation])
-            conditions.push([1, 0, endDisplacements[translation]]);
-        if (!released[6 + translation]) {
-            conditions.push([1, L, endDisplacements[6 + translation] - PL]);
-        }
-        if (!released[rotation])
-            conditions.push([0, 1, rotationSign * endDisplacements[rotation]]);
-        if (!released[6 + rotation]) {
-            conditions.push([0, 1, rotationSign * endDisplacements[6 + rotation] - dPL]);
-        }
-        for (let i = 0; i < conditions.length; i++) {
-            for (let j = i + 1; j < conditions.length; j++) {
-                const [a1, b1, r1] = conditions[i];
-                const [a2, b2, r2] = conditions[j];
-                const det = a1 * b2 - a2 * b1;
-                if (Math.abs(det) > 1e-12 * Math.max(1, L)) {
-                    return { a: (r1 * b2 - r2 * b1) / det, b: (a1 * r2 - a2 * r1) / det };
-                }
-            }
-        }
-        return { a: 0, b: 0 }; // unreachable: such releases are rejected as a mechanism (307)
+    /** a, b of `u = a + b·x + P(x)` from the end translations: u(0) = a, u(L) = a + b·L + P(L). */
+    const planeConstants = (translation, PL) => {
+        const a = endDisplacements[translation];
+        return { a, b: (endDisplacements[6 + translation] - a - PL) / L };
     };
     const P3 = (i) => grid[i] * k3[i] - sk3[i];
     const P2 = (i) => grid[i] * k2[i] - sk2[i];
-    const plane2 = planeConstants(1, 5, 1, P3(last), k3[last]); // u2, with u2' = R3
-    const plane3 = planeConstants(2, 4, -1, P2(last), k2[last]); // u3, with u3' = −R2
-    const axialStart = member.released[0] ? endDisplacements[6] - strain[last] : endDisplacements[0];
+    const plane2 = planeConstants(1, P3(last)); // u2
+    const plane3 = planeConstants(2, P2(last)); // u3
+    const axialStart = endDisplacements[0];
     const onGrid = grid.map((x, i) => ({
         u1: axialStart + strain[i],
         u2: plane2.a + plane2.b * x + P3(i),
@@ -181,9 +158,10 @@ const present = (c, value) => roundToDecimalPlaces(value * SCALE[c], c === "u1"
  * Builds the force and displacement diagrams of one member from its end forces, its end
  * displacements and its loads.
  *
- * @param member - Member geometry, section, axes and end releases
+ * @param member - Member geometry, section and axes
  * @param startForce - Force the start node exerts on the member, right-hand local, N and N·mm
- * @param endDisplacements - The 12 local right-hand displacements of the member's end nodes, mm and rad
+ * @param endDisplacements - The member's own 12 local right-hand end displacements (equal to the
+ *   node displacements except across released or spring end DOFs), mm and rad
  * @param loads - Member loads in local axes, N and mm
  * @param stationCount - Number of evenly spaced stations (at least 2); load breakpoints and
  *   zero-shear points are added, with both sides reported at a point load