formulab 0.1.0 → 0.5.0

package/CHANGELOG.md

@@ -5,7 +5,288 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
-## [0.1.0] - 2025-02-05
+## [0.5.0] - 2026-02-07
+
+### Added
+
+- **Machining** (new domain, 12 functions):
+  - `truePosition()` — GD&T True Position with MMC bonus tolerance (ASME Y14.5)
+  - `boltCircle()` — Bolt hole pattern coordinate calculation
+  - `sineBarHeight()` — Sine bar gauge block height with rounding error analysis
+  - `radialChipThinning()` — Radial chip thinning compensation for light cuts
+  - `toolDeflection()` — End mill cantilever deflection (carbide/HSS)
+  - `cuspHeight()` — Ball end mill scallop height and surface roughness
+  - `effectiveDiameter()` — Ball end mill effective cutting diameter at depth
+  - `boringBarDeflection()` — Boring bar deflection with L/D ratio guidance
+  - `threadOverWires()` — 3-wire thread measurement (60°/55°/29°/30° angles)
+  - `gaugeBlockStack()` — Gauge block combination using successive subtraction (47/88/81-pc sets)
+  - `triangleSolver()` — Triangle solver (SSS/SAS/ASA/AAS/SSA with law of cosines/sines)
+  - `cycleTimeEstimator()` — CNC cycle time estimation with operation breakdown
+
+## [0.4.0] - 2026-02-07
+
+### Added
+
+- **Battery** (new domain, 10 functions):
+  - `energyDensity()` — Wh/kg and Wh/L energy density calculation
+  - `cRate()` — C-rate ↔ current/time bidirectional conversion
+  - `stateOfHealth()` — SOH % with degradation status thresholds (IEEE 1188)
+  - `batteryPackConfig()` — Series/parallel cell configuration calculator
+  - `cycleLife()` — Cycle life estimation with chemistry (7 types), DOD, and temperature factors
+  - `internalResistance()` — DCIR calculation from OCV and load voltage (IEC 61960)
+  - `selfDischarge()` — Self-discharge rate (daily and monthly %)
+  - `thermalRunaway()` — Thermal safety margin analysis (SAE J2464)
+  - `bmsBalancing()` — BMS passive cell balancing time estimation
+  - `chargingProfile()` — CC-CV charging profile timing with phase breakdown
+
+- **Environmental** (new domain, 10 functions):
+  - `scope1Emissions()` — Direct combustion emissions with 6 fuel types (EPA/IPCC 2006 factors)
+  - `scope2Emissions()` — Purchased electricity emissions for 12 grid regions (IEA 2023 factors)
+  - `scope3Emissions()` — Supply chain spend-based emissions for 8 Scope 3 categories (EPA EEIO)
+  - `vocEmissions()` — VOC emissions with capture/destruction efficiency
+  - `productCarbonFootprint()` — Product lifecycle carbon footprint with stage breakdown (ISO 14067)
+  - `gwpCalculator()` — GWP conversion for 8 GHG gases × 3 time horizons (IPCC AR6)
+  - `energyIntensity()` — Energy intensity per unit/revenue (ISO 50001)
+  - `waterFootprint()` — Water footprint blue/green/grey breakdown (ISO 14046)
+  - `emissionsIntensity()` — Emissions intensity per unit/revenue/employee
+  - `esgSummary()` — ESG reduction tracking with projections (TCFD/CDP/SBTi)
+
+- **Metal** (3 new functions, 22 → 25):
+  - `materialGradeConverter()` — ASTM/EN/JIS/GB/KS grade cross-reference (20 grades × 5 standards)
+  - `pipeSpec()` — ANSI/ASME B36.10 pipe dimensions lookup (15 sizes × 6 schedules)
+  - `flangeSpec()` — ASME B16.5 WN flange dimensions lookup (8 sizes × 3 pressure classes)
+
+### Changed
+
+- Total domains: 11 → 13
+- Total function count: 111 → 134 (+23)
+- Total test count: 1705 → 1887 (+182)
+
+## [0.3.0] - 2026-02-06
+
+### Added
+
+- **Quality**: `controlChart()` — SPC X-bar/R and X-bar/S control charts with AIAG/ASTM E2587 constants (n=2..25)
+- **Chemical**: `pipeFlow()` — Darcy-Weisbach pipe flow calculator with Swamee-Jain friction factor (laminar/transitional/turbulent)
+- **Chemical**: `heatTransfer()` — Three modes: conduction (Fourier), convection (Newton), radiation (Stefan-Boltzmann)
+- **Construction**: `momentOfInertia()` — Section properties (A, Ix, Iy, Sx, Sy, rx, ry) for 7 cross-section shapes (rectangle, circle, hollow rectangle, hollow circle, I-beam, T-section, C-channel)
+- **Electronics**: `ohmsLaw()` — V/I/R/P calculator with discriminated union input
+- **Energy**: `solarOutput()` — Solar panel output estimation based on PVWatts methodology
+- **Safety**: `ventilationRate()` — Required ventilation rate (ACH/CFM/L/s) per ASHRAE 62.1 / OSHA guidelines
+- **Automotive**: `brakingDistance()` — Stopping distance calculator per AASHTO method with grade adjustment
+
+### Changed
+
+- Total function count: 103 → 111
+- Total test count: 1630 → 1705
+
+## [0.2.0] - 2026-02-06
+
+### Breaking Changes
+- **`MetalWeightInput`**: Now a discriminated union by `shape`. Each shape variant requires its specific fields (e.g., `shape: 'plate'` requires `width` and `thickness`). The `materialName` field is now a `MaterialName` literal union type instead of `string`.
+- **`DilutionInput`**: Now a discriminated union by `solveFor`. Each variant provides exactly the three known values, eliminating optional fields and `!` assertions.
+- **`ReactorInput`**: Now a discriminated union by `shape`. `shape: 'cylindrical'` requires `height`; `shape: 'spherical'` does not accept `height`.
+- **`BoltInput`**: Now a discriminated union by `mode`. `mode: 'torqueToPreload'` requires `torque`; `mode: 'preloadToTorque'` requires `preload`.
+- **QR code types removed**: `QrcodeInput`, `QrcodeResult`, `ErrorCorrectionLevel` types and `qrcode.ts` module removed from `formulab/utility` (was empty shell with no implementation).
+
+### Enhanced
+- **`roundTo()`**: Fixed negative number rounding (`-2.555` now correctly rounds to `-2.56` instead of `-2.55`). Uses sign-aware epsilon correction.
+- **`oee()`**: Results now consistently rounded via `roundTo()` (factors to 4 decimals, percentages to 1 decimal).
+
+### Fixed
+- **Duplicate code removal**: Consolidated `normalCDF()`, `normalInvCDF()`, and `clamp()` from 6 files into shared `math.ts` module.
+- **`Math.round` inconsistency**: Replaced manual `Math.round(x * N) / N` patterns with `roundTo()` in 8 files (`metalWeight`, `fallClearance`, `welding`, `carbonFootprint`, `powerCost`, `bendAllowance`, `weldHeat`, `roughness`).
+- **Duplicate test file**: Removed `src/metal/pressTonnage.test.ts` (14 duplicate tests); canonical version at `src/metal/__tests__/pressTonnage.test.ts` (19 tests) retained.
+
+### Changed
+- **Coverage thresholds**: Raised from baseline (5% lines) to realistic levels (90% lines, 95% functions, 85% branches, 90% statements).
+- Total test count: 1630 tests across 102 test files.
+
+## [0.1.11] - 2026-02-06
+
+### Enhanced
+- **`fallClearance()`**: Complete rewrite with ANSI Z359.1/Z359.4 compliance
+  - Added `rescueClearance` parameter (default 0.9m per ANSI Z359.4)
+  - Added `obstacleHeight` parameter for elevated obstacle detection
+  - Added `freeSpaceRequired` and `clearanceAboveObstacle` outputs
+  - Added `warnings` array with validation against OSHA/ANSI limits
+  - Separated physical fall distance from safety factor (correct per standards)
+
+- **`weldHeat()`**: Major enhancement with industry standards
+  - Added AWS D1.1:2020 Table 5.8 preheat requirements with CE×thickness interaction
+  - Added Rosenthal cooling time equations (t8/5) for thick and thin plate heat flow
+  - Added Yurioka HAZ hardness formula with cooling rate factor
+  - Added Carbon Equivalent Pcm formula for low-alloy steels (C < 0.18%)
+  - Added hydrogen control level determination
+  - Added source attribution for preheat recommendations
+
+- **`pressTonnage()`**: Deep drawing enhancements with DIN 8584 basis
+  - Added Siebel formula for drawing force with friction and bending effects
+  - Added blank holder force calculation: F_bh = π/4 × (D² - d²) × p_bh
+  - Added multi-draw detection and estimation for deep draws (d/D < 0.55)
+  - Added `frictionCoefficient`, `blankHolderPressure`, `dieRadius` parameters
+  - Added `warnings` array for process risk assessment
+
+### Added
+- **Test suite**: `pressTonnage` (19 test cases covering blanking, bending, enhanced drawing)
+- Total test count: 1648 tests across 103 test files
+
+### Changed
+- Formula review completed for high-risk domains (safety, metal)
+- Three functions upgraded from 70-85% to 95%+ theoretical completeness
+
+## [0.1.10] - 2026-02-06
+
+### Added
+- **Test suites (Cycle 14)**: 28개 함수 완전 테스트 추가 (437 test cases)
+  - metal: `bendAllowance`, `bolt`, `cuttingStock`, `flatPattern`, `gear`, `kFactorReverse`, `material`, `pressFit`, `pressTonnage`, `roughness`, `screw`, `spring`, `tap`, `thread`, `tolerance`, `vibration`, `weldHeat`, `welding`
+  - logistics: `containerFit`, `fillRate`, `freightClass`, `pallet3d`, `palletStack`, `pickTime`, `shipping`, `tsp`
+  - construction: `roof`, `stair`
+- Total test count: 1622 tests across 102 test files
+
+### Changed
+- Coverage dramatically improved: 62.09% → 98.11% lines (+36.02%)
+- Branch coverage: 88.72% → 90.55%
+- Function coverage: 89.63% → 100%
+- All domains now have 95%+ line coverage
+- **Milestone**: 100% function coverage achieved, 98%+ line coverage, library fully tested
+
+## [0.1.9] - 2026-02-05
+
+### Added
+- **Test suites (Cycle 13)**: 25개 함수 대규모 테스트 추가 (488 test cases)
+  - quality: `aql`, `cycle`, `downtime`, `lineBalancing`, `ppk`
+  - chemical: `concentration`, `ph`, `reactor`, `shelfLife`, `injectionCycle`
+  - electronics: `reflow`, `resistor`, `smt-takt`, `solder`, `trace`, `capacitor`, `stencil`, `via`
+  - construction: `beamLoad`, `concreteMix`, `formwork`, `rebar`, `aggregate`, `brick`, `pert`
+- Total test count: 1185 tests across 74 test files
+
+### Changed
+- Coverage improved: 36.46% → 62.09% lines (+25.63%)
+- quality 도메인: 98.68% lines
+- chemical 도메인: 99.40% lines
+- electronics 도메인: 98.37% lines
+- construction 도메인: 83.30% lines
+- **Milestone**: Line coverage 60% 돌파, 6개 도메인 90%+ 커버리지 달성
+
+## [0.1.8] - 2026-02-05
+
+### Added
+- **Test suites (Cycle 11)**: `yieldCalc`, `rpn`, `ppm` (quality), `safetyStock`, `kanban` (logistics)
+- **Test suites (Cycle 12)**: 17개 함수 대규모 테스트 추가
+  - safety: `fallClearance`, `noiseExposure`, `wbgtCalculate`, `havsCalculate`, `respiratorCalculate`
+  - food: `haccp`, `expiry`
+  - automotive: `batteryRuntime`, `evCharging`, `fuelEconomy`, `tireCompare`, `torque`, `power`
+  - energy: `compressedAirCost`, `motorEfficiency`, `pfCorrection`, `vfdSavings`
+- Total test count: 697 tests across 49 test files
+
+### Changed
+- Coverage improved: 20.74% → 36.46% lines (+15.72%)
+- safety 도메인: 39.63% → 100% lines
+- food 도메인: 45.45% → 100% lines
+- automotive 도메인: 11.89% → 68.96% lines
+- energy 도메인: 26.45% → 100% lines
+- **Milestone**: Line coverage 35% 돌파, 4개 도메인 100% 커버리지 달성
+
+## [0.1.7] - 2026-02-05
+
+### Added
+- **Test suites**: `cutting`, `bearing`, `hardness` (metal domain, 45 test cases)
+- Total test count: 371 tests across 27 test files
+
+### Changed
+- Coverage improved: 18.71% → 20.74% lines (+2.03%)
+- metal 도메인: 4.78% → 11.54% lines
+- **Milestone**: Line coverage 20% 돌파
+
+## [0.1.6] - 2026-02-05
+
+### Added
+- **Test suites**: `calculateUnit`, `getUnitCategories`, `getUnitsForCategory` (unit conversion, 29 test cases)
+- **Test suites**: `solveAssignment` (Hungarian algorithm, 12 test cases)
+- Total test count: 326 tests across 24 test files
+
+### Changed
+- Coverage improved: 15.21% → 18.71% lines (+3.5%)
+- utility 도메인: 0% → 78.83% lines
+- Branches coverage: 76.89% → 80.96%
+
+## [0.1.5] - 2026-02-05
+
+### Added
+- **Test suites**: `ledResistor`, `awgProperties` (electronics domain, 33 test cases)
+- **Test suites**: `calorie`, `nutrition` (food domain, 23 test cases)
+- Total test count: 285 tests across 22 test files
+
+### Changed
+- Coverage improved: 12.08% → 15.21% lines
+- electronics 도메인: 0% → 18.77% lines
+- food 도메인: 0% → 45.45% lines
+
+## [0.1.4] - 2026-02-05
+
+### Added
+- **Test suites**: `slope`, `earthwork` (construction domain, 27 test cases)
+- **Test suites**: `powerCost`, `carbonFootprint` (energy domain, 24 test cases)
+- Total test count: 229 tests across 18 test files
+
+### Changed
+- Coverage improved: 10.38% → 12.08% lines
+- construction 도메인: 0% → 7.73% lines
+- energy 도메인: 0% → 26.45% lines
+
+## [0.1.3] - 2026-02-05
+
+### Added
+- **Test suites**: `dilution`, `batch` (chemical domain, 25 test cases)
+- **Input validation tests**: `eoq` validation edge cases (5 test cases)
+- Total test count: 178 tests across 14 test files
+
+### Fixed
+- **`eoq()` input validation**: Returns zeros for invalid inputs (negative/zero demand, cost)
+- **`gearRatio()` NaN bug**: `mechanicalAdvantage` now returns 0 when `inputTorque` is 0
+
+### Changed
+- Coverage improved: 8.8% → 10.38% lines
+
+## [0.1.2] - 2026-02-05
+
+### Added
+- **Test suites**: `eoq`, `dimWeight`, `gearRatio` (40 additional test cases)
+- Total test count: 148 tests across 12 test files
+- Coverage measurement infrastructure fully operational
+
+### Fixed
+- **Coverage tool version mismatch**: Downgraded `@vitest/coverage-v8` from 4.x to 3.x to match vitest 3.x
+
+### Changed
+- **Coverage thresholds**: Adjusted to realistic baseline (lines: 5%, functions: 50%, branches: 50%)
+- Current coverage metrics: 8.8% lines, 62.38% functions, 72.25% branches
+
+## [0.1.1] - 2026-02-05
+
+### Fixed
+- **Main entry point**: Added missing exports for 5 domains (metal, chemical, electronics, construction, energy)
+- **`roundTo()` precision bug**: Fixed IEEE 754 floating-point rounding issues (e.g., `roundTo(0.615, 2)` now correctly returns `0.62`)
+- **Massive duplicate code removal**: Consolidated `roundTo()` from 36 files to single source in `utils.ts`
+  - Removed from: quality (8 files), metal (15 files), electronics (8 files), chemical (3 files)
+- **README documentation**: Updated all API examples to match actual implementation
+
+### Changed
+- **BREAKING**: `GearInput` → `GearRatioInput` (automotive domain) to avoid collision with metal domain
+- **BREAKING**: `GearResult` → `GearRatioResult` (automotive domain)
+- **BREAKING**: `getCategories()` → `getUnitCategories()` (utility domain) to avoid collision with metal domain
+- Removed `--passWithNoTests` flag from test script (tests are now mandatory)
+
+### Added
+- Input validation for `oee()`: throws error when `goodCount > totalCount`
+- Input validation for `oee()`: throws error when `goodCount < 0`
+- Test suites: `oee`, `cpk`, `taktTime`, `dpmo`, `mtbf`, `metalWeight`, `cbm`, `nioshLifting`, `roundTo` (108 test cases)
+- NaN/Infinity handling in `roundTo()`
+- Vitest coverage configuration with thresholds
+- `test:coverage` npm script
+
+## [0.1.0] - 2026-02-05
 
 ### Added
 
@@ -134,7 +415,7 @@ Initial release with 103 calculation functions across 11 domains.
 #### Utility (3 functions)
 - `solveAssignment()` - Hungarian algorithm optimization
 - `calculateUnit()` - Unit conversion
-- `getCategories()` - Get unit categories
+- `getUnitCategories()` - Get unit categories
 
 ### Technical Features
 - Zero dependencies

package/README.md

@@ -10,14 +10,36 @@ A comprehensive collection of engineering formulas and calculations for manufact
 
 ## Features
 
-- **100+ industrial calculations** — OEE, Cpk, metal weight, CBM, NIOSH lifting, and more
-- **11 specialized domains** — Quality, Metal, Chemical, Electronics, Construction, Automotive, Logistics, Energy, Safety, Food, Utility
+- **146 industrial calculations** — OEE, Cpk, SPC control charts, metal weight, CNC machining, GD&T, pipe flow, CBM, NIOSH lifting, battery SOH, GHG emissions, and more
+- **14 specialized domains** — Quality, Metal, Chemical, Electronics, Construction, Automotive, Logistics, Energy, Safety, Food, Utility, Battery, Environmental, Machining
 - **Zero dependencies** — Lightweight and fast
 - **TypeScript first** — Full type definitions included
 - **Tree-shakeable** — Import only what you need
 - **Well tested** — Comprehensive test coverage
 - **Research-based** — Built on industry standards (ISO, OSHA, NIOSH, IPC)
 
+## Verification Status
+
+| Domain | Functions | Golden Tests | Key References |
+|--------|-----------|-------------|----------------|
+| Quality | 14 | oee, cpk, controlChart | ISO 22400-2, AIAG/ASTM E2587, JIPM |
+| Metal | 25 | metalWeight | Machinery's Handbook, ASME B36.10/B16.5 |
+| Logistics | 14 | cbm | Physical formula |
+| Safety | 7 | nioshLifting | NIOSH 94-110 |
+| Chemical | 9 | — | Darcy-Weisbach, Fourier |
+| Electronics | 11 | — | IPC-2221 |
+| Construction | 12 | — | AISC, Timoshenko |
+| Automotive | 8 | — | AASHTO |
+| Energy | 7 | — | NREL PVWatts |
+| Food | 4 | — | HACCP |
+| Utility | 3 | — | — |
+| Battery | 10 | — | IEEE 1188, IEC 62620, Battery University |
+| Environmental | 10 | — | GHG Protocol, IPCC AR6, IEA 2023 |
+| Machining | 12 | — | Machinery's Handbook, ASME Y14.5, Sandvik Coromant |
+
+> Functions with golden reference tests have been verified against authoritative engineering sources.
+> See each function's JSDoc for specific references.
+
 ## Installation
 
 ```bash
@@ -35,42 +57,50 @@ import { oee, metalWeight, cbm } from 'formulab';
 
 // Calculate OEE (Overall Equipment Effectiveness)
 const result = oee({
-  availability: 0.90,
-  performance: 0.95,
-  quality: 0.99
+  rawData: {
+    plannedTime: 480,      // minutes (8 hours)
+    runTime: 432,          // minutes (90% availability)
+    totalCount: 1000,
+    goodCount: 990,        // 99% quality
+    idealCycleTime: 0.456, // minutes per piece (95% performance)
+  },
 });
-console.log(result.oee); // 0.846 (84.6%)
+console.log(result.percentages.oee); // 84.6%
 
 // Calculate metal weight
 const weight = metalWeight({
-  material: 'steel',
   shape: 'plate',
-  dimensions: { length: 1000, width: 500, thickness: 10 },
+  materialName: 'steel',
+  length: 1000,    // mm
+  width: 500,      // mm
+  thickness: 10,   // mm
 });
 console.log(weight.weight); // 39.25 kg
 
-// Calculate CBM
+// Calculate CBM (Cubic Meter)
 const volume = cbm({
   length: 120,
   width: 80,
   height: 100,
-  unit: 'cm'
+  quantity: 1,
+  unit: 'cm',
 });
-console.log(volume.cbm); // 0.96 m³
+console.log(volume.totalCbm); // 0.96 m³
 ```
 
 ## Domains
 
-### Quality & Production (13 functions)
+### Quality & Production (14 functions)
 
 ```typescript
-import { oee, cpk, taktTime, dpmo } from 'formulab/quality';
+import { oee, cpk, taktTime, dpmo, controlChart } from 'formulab/quality';
 ```
 
 | Function | Description |
 |----------|-------------|
 | `oee()` | Overall Equipment Effectiveness |
 | `cpk()` | Process Capability Index |
+| `controlChart()` | SPC X-bar/R and X-bar/S charts |
 | `cycleTime()` | Cycle Time analysis |
 | `taktTime()` | Takt Time calculation |
 | `aql()` | AQL sampling inspection |
@@ -83,7 +113,7 @@ import { oee, cpk, taktTime, dpmo } from 'formulab/quality';
 | `rpn()` | Risk Priority Number (FMEA) |
 | `yieldCalc()` | First Pass Yield / RTY |
 
-### Metal & Machining (22 functions)
+### Metal & Machining (25 functions)
 
 ```typescript
 import { metalWeight, bendAllowance, cutting, bearing } from 'formulab/metal';
@@ -113,11 +143,14 @@ import { metalWeight, bendAllowance, cutting, bearing } from 'formulab/metal';
 | `vibration()` | Natural frequency analysis |
 | `weldHeat()` | Weld heat input calculation |
 | `welding()` | Welding parameters |
+| `materialGradeConverter()` | ASTM/EN/JIS/GB/KS grade cross-reference |
+| `pipeSpec()` | ANSI/ASME pipe dimensions lookup |
+| `flangeSpec()` | ASME B16.5 flange dimensions lookup |
 
-### Chemical & Process (7 functions)
+### Chemical & Process (9 functions)
 
 ```typescript
-import { dilution, concentration, ph, reactor } from 'formulab/chemical';
+import { dilution, concentration, ph, reactor, pipeFlow, heatTransfer } from 'formulab/chemical';
 ```
 
 | Function | Description |
@@ -125,19 +158,22 @@ import { dilution, concentration, ph, reactor } from 'formulab/chemical';
 | `batch()` | Batch scaling calculation |
 | `concentration()` | Concentration conversion |
 | `dilution()` | Dilution (C1V1 = C2V2) |
+| `heatTransfer()` | Conduction/convection/radiation heat transfer |
 | `ph()` | pH and buffer calculations |
+| `pipeFlow()` | Darcy-Weisbach pipe flow pressure drop |
 | `reactor()` | Reactor sizing |
 | `shelfLife()` | Shelf life prediction (Arrhenius) |
 | `injectionCycle()` | Injection molding cycle time |
 
-### Electronics & SMT (10 functions)
+### Electronics & SMT (11 functions)
 
 ```typescript
-import { traceWidth, solderPaste, resistorDecode } from 'formulab/electronics';
+import { traceWidth, solderPaste, resistorDecode, ohmsLaw } from 'formulab/electronics';
 ```
 
 | Function | Description |
 |----------|-------------|
+| `ohmsLaw()` | Ohm's Law V/I/R/P calculator |
 | `reflowProfile()` | Reflow temperature profile |
 | `resistorDecode()` | Resistor color code decoder |
 | `smtTakt()` | SMT line takt time |
@@ -149,14 +185,15 @@ import { traceWidth, solderPaste, resistorDecode } from 'formulab/electronics';
 | `stencilAperture()` | Stencil aperture design |
 | `viaCurrent()` | Via current capacity |
 
-### Construction (11 functions)
+### Construction (12 functions)
 
 ```typescript
-import { concreteMix, rebarWeight, slope, stair } from 'formulab/construction';
+import { concreteMix, rebarWeight, slope, stair, momentOfInertia } from 'formulab/construction';
 ```
 
 | Function | Description |
 |----------|-------------|
+| `momentOfInertia()` | Section properties (Ix, Sx, rx) for 7 shapes |
 | `beamLoad()` | Beam load calculation |
 | `concreteMix()` | Concrete mix ratio |
 | `earthwork()` | Earthwork volume |
@@ -169,14 +206,15 @@ import { concreteMix, rebarWeight, slope, stair } from 'formulab/construction';
 | `roof()` | Roof calculation |
 | `stair()` | Stair dimension calculation |
 
-### Automotive (7 functions)
+### Automotive (8 functions)
 
 ```typescript
-import { batteryRuntime, evCharging, torque } from 'formulab/automotive';
+import { batteryRuntime, evCharging, torque, brakingDistance } from 'formulab/automotive';
 ```
 
 | Function | Description |
 |----------|-------------|
+| `brakingDistance()` | Stopping distance (AASHTO method) |
 | `batteryRuntime()` | Battery capacity/runtime |
 | `evCharging()` | EV charging time estimation |
 | `fuelEconomy()` | Fuel economy conversion |
@@ -208,14 +246,15 @@ import { cbm, eoq, safetyStock, kanban } from 'formulab/logistics';
 | `shipping()` | Shipping cost estimation |
 | `tsp()` | Traveling salesman problem |
 
-### Energy & Utilities (6 functions)
+### Energy & Utilities (7 functions)
 
 ```typescript
-import { powerCost, motorEfficiency, carbonFootprint } from 'formulab/energy';
+import { powerCost, motorEfficiency, carbonFootprint, solarOutput } from 'formulab/energy';
 ```
 
 | Function | Description |
 |----------|-------------|
+| `solarOutput()` | Solar panel output estimation (PVWatts-based) |
 | `carbonFootprint()` | Scope 2 emissions |
 | `compressedAirCost()` | Compressed air cost |
 | `motorEfficiency()` | Motor upgrade ROI |
@@ -223,14 +262,15 @@ import { powerCost, motorEfficiency, carbonFootprint } from 'formulab/energy';
 | `powerCost()` | Electricity cost with demand |
 | `vfdSavings()` | VFD energy savings |
 
-### Safety & Ergonomics (6 functions)
+### Safety & Ergonomics (7 functions)
 
 ```typescript
-import { nioshLifting, noiseExposure, wbgtCalculate } from 'formulab/safety';
+import { nioshLifting, noiseExposure, wbgtCalculate, ventilationRate } from 'formulab/safety';
 ```
 
 | Function | Description |
 |----------|-------------|
+| `ventilationRate()` | Required ventilation ACH/CFM (ASHRAE/OSHA) |
 | `fallClearance()` | Fall protection clearance |
 | `nioshLifting()` | NIOSH lifting equation |
 | `noiseExposure()` | TWA/Dose calculation |
@@ -261,7 +301,66 @@ import { solveAssignment, calculateUnit } from 'formulab/utility';
 |----------|-------------|
 | `solveAssignment()` | Hungarian algorithm optimization |
 | `calculateUnit()` | Unit conversion |
-| `getCategories()` | Get unit categories |
+| `getUnitCategories()` | Get unit categories |
+
+### Battery (10 functions)
+
+```typescript
+import { energyDensity, cRate, stateOfHealth, cycleLife } from 'formulab/battery';
+```
+
+| Function | Description |
+|----------|-------------|
+| `energyDensity()` | Wh/kg and Wh/L energy density |
+| `cRate()` | C-rate ↔ current/time conversion |
+| `stateOfHealth()` | SOH % with degradation status |
+| `batteryPackConfig()` | Series/parallel cell configuration |
+| `cycleLife()` | Cycle life estimation (chemistry/DOD/temp) |
+| `internalResistance()` | DCIR from OCV and load voltage |
+| `selfDischarge()` | Self-discharge rate calculation |
+| `thermalRunaway()` | Thermal safety margin analysis |
+| `bmsBalancing()` | BMS cell balancing time estimation |
+| `chargingProfile()` | CC-CV charging profile timing |
+
+### Environmental (10 functions)
+
+```typescript
+import { scope1Emissions, scope2Emissions, gwpCalculator, esgSummary } from 'formulab/environmental';
+```
+
+| Function | Description |
+|----------|-------------|
+| `scope1Emissions()` | Fuel combustion direct emissions (6 fuels) |
+| `scope2Emissions()` | Purchased electricity emissions (12 regions) |
+| `scope3Emissions()` | Supply chain spend-based emissions (8 categories) |
+| `vocEmissions()` | VOC emissions with capture/destruction |
+| `productCarbonFootprint()` | Product lifecycle carbon footprint |
+| `gwpCalculator()` | GWP conversion (8 gases × 3 time horizons) |
+| `energyIntensity()` | Energy intensity (MJ/unit, kWh/unit) |
+| `waterFootprint()` | Water footprint (blue/green/grey) |
+| `emissionsIntensity()` | Emissions intensity per unit/revenue/employee |
+| `esgSummary()` | ESG reduction tracking and projections |
+
+### Machining & CNC (12 functions)
+
+```typescript
+import { truePosition, boltCircle, toolDeflection, threadOverWires } from 'formulab/machining';
+```
+
+| Function | Description |
+|----------|-------------|
+| `truePosition()` | GD&T True Position (diametral, MMC bonus) |
+| `boltCircle()` | Bolt hole pattern coordinates |
+| `sineBarHeight()` | Sine bar gauge block height |
+| `radialChipThinning()` | Chip load compensation for light radial cuts |
+| `toolDeflection()` | End mill cantilever deflection |
+| `cuspHeight()` | Ball mill scallop height |
+| `effectiveDiameter()` | Ball mill effective cutting diameter |
+| `boringBarDeflection()` | Boring bar deflection with L/D guidance |
+| `threadOverWires()` | 3-wire thread measurement |
+| `gaugeBlockStack()` | Gauge block combination (47/88/81-pc sets) |
+| `triangleSolver()` | Triangle solver (SSS/SAS/ASA/SSA) |
+| `cycleTimeEstimator()` | CNC cycle time estimation |
 
 ## API Examples
 
@@ -271,19 +370,19 @@ import { solveAssignment, calculateUnit } from 'formulab/utility';
 import { oee } from 'formulab/quality';
 
 const result = oee({
-  plannedTime: 480,      // minutes
-  operatingTime: 420,
-  totalPieces: 1000,
-  goodPieces: 990,
-  idealCycleTime: 0.5,
+  rawData: {
+    plannedTime: 480,      // minutes
+    runTime: 420,          // actual running time
+    totalCount: 1000,      // total pieces produced
+    goodCount: 990,        // good pieces
+    idealCycleTime: 0.4,   // minutes per piece
+  },
 });
 
 console.log(result);
 // {
-//   oee: 0.846,
-//   availability: 0.875,
-//   performance: 0.952,
-//   quality: 0.99,
+//   factors: { availability: 0.875, performance: 0.952, quality: 0.99, oee: 0.825 },
+//   percentages: { availability: 87.5, performance: 95.2, quality: 99, oee: 82.5 }
 // }
 ```
 
@@ -293,22 +392,22 @@ console.log(result);
 import { nioshLifting } from 'formulab/safety';
 
 const result = nioshLifting({
-  loadWeight: 23,
-  horizontalDistance: 25,
-  verticalLocation: 75,
-  verticalTravel: 25,
-  asymmetryAngle: 0,
-  couplingQuality: 'good',
-  frequency: 1,
-  duration: 'moderate',
+  loadWeight: 23,           // kg
+  horizontalDistance: 25,   // cm
+  verticalDistance: 75,     // cm (height at lift origin)
+  verticalTravel: 25,       // cm (vertical travel distance)
+  asymmetryAngle: 0,        // degrees
+  coupling: 'good',         // 'good' | 'fair' | 'poor'
+  frequency: 1,             // lifts per minute
+  duration: 'short',        // 'short' | 'medium' | 'long'
 });
 
 console.log(result);
 // {
-//   rwl: 23.0,              // Recommended Weight Limit (kg)
-//   liftingIndex: 1.0,      // LI = Load / RWL
-//   riskLevel: 'low',
-//   multipliers: { hm: 1.0, vm: 1.0, dm: 0.93, am: 1.0, fm: 0.94, cm: 1.0 }
+//   rwl: 21.62,             // Recommended Weight Limit (kg)
+//   liftingIndex: 1.06,     // LI = Load / RWL
+//   riskLevel: 'moderate',
+//   hm: 1.0, vm: 1.0, dm: 1.0, am: 1.0, fm: 0.94, cm: 1.0
 // }
 ```
 

package/dist/automotive/brakingDistance.d.ts

@@ -0,0 +1,16 @@
+import type { BrakingDistanceInput, BrakingDistanceResult } from './types.js';
+/**
+ * Braking Distance Calculator (AASHTO method)
+ *
+ * Calculates stopping distance based on:
+ * - Reaction distance: d_r = v × t_reaction
+ * - Braking distance: d_b = v² / (2 × g × (f ± G))
+ *   where G = grade/100 (positive uphill, negative downhill)
+ *
+ * Reference: AASHTO "A Policy on Geometric Design of Highways and Streets"
+ *
+ * @param input - speed, friction, reaction time, grade
+ * @returns stopping distance components
+ */
+export declare function brakingDistance(input: BrakingDistanceInput): BrakingDistanceResult;
+//# sourceMappingURL=brakingDistance.d.ts.map

package/dist/automotive/brakingDistance.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"brakingDistance.d.ts","sourceRoot":"","sources":["../../src/automotive/brakingDistance.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,oBAAoB,EAAE,qBAAqB,EAAE,MAAM,YAAY,CAAC;AAE9E;;;;;;;;;;;;GAYG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,oBAAoB,GAAG,qBAAqB,CAqClF"}