56 lines
3.3 KiB
Python
56 lines
3.3 KiB
Python
#!/usr/bin/env python3
|
|
"""Render the confirmed station table (지반고 = A열, 계획고 = B열, 20 m station
|
|
numbering per user decision 2026-10-21) into DOC/03_설계 and check the design
|
|
grade against the drawn label S=2.4913% (정상 구배 검산). No value is adjusted."""
|
|
import csv
|
|
import sys
|
|
|
|
SRC, OUT = sys.argv[1], sys.argv[2]
|
|
rows = []
|
|
with open(SRC, encoding="utf-8") as fh:
|
|
rd = list(csv.DictReader(fh))
|
|
gl, dg = {}, {}
|
|
for r in rd:
|
|
sta = int(r["station_anchor"])
|
|
if float(r["across_m"]) < -25:
|
|
gl[sta] = float(r["value"])
|
|
else:
|
|
dg[sta] = float(r["value"])
|
|
stas = sorted(set(gl) & set(dg))
|
|
n = len(stas)
|
|
xs = [s for s in stas]
|
|
mx = sum(xs) / n
|
|
my = sum(dg[s] for s in stas) / n
|
|
slope = sum((s - mx) * (dg[s] - my) for s in stas) / sum((s - mx) ** 2 for s in stas)
|
|
grade_pct = slope * 100
|
|
resid = [dg[s] - (my + slope * (s - mx)) for s in stas]
|
|
max_res = max(abs(r) for r in resid)
|
|
lines = ["# 본선 종단 측점값 — 지반고·계획고 (초안, 확정 해석 반영)\n",
|
|
"| 버전 | 일자 | 변경 내용 |", "|---|---|---|",
|
|
"| v0.1.00 | 2026-10-21 | 최초 작성. 02-002 종단 눈금 4+740~5+480, 지반고(A열)·계획고(B열) 확정 해석 반영 |\n",
|
|
"## 1. 확정 사항 (사용자 답변 2026-10-21)",
|
|
"- 소수 A열(중심선 측 약 -33m) = **지반고**, B열(약 -16m) = **계획고** [출처: C0060302-002, CR-TABL-TEX2]",
|
|
"- 측점은 **측점 번호 20m 기준**. 도면상 눈금 거리 21.27m는 사용하지 않음 (축척 1.0635배 불일치 기록)",
|
|
"- 공백 구간(4+700~4+740, 5+480~5+500)은 **[미확정]**, 확보 구간만 모델링\n",
|
|
"## 2. 측점별 지반고·계획고",
|
|
"| 측점(m) | 표기 | 지반고 GL(m) | 계획고(m) |", "|---|---|---|---|"]
|
|
for s in stas:
|
|
lines.append(f"| {s} | {s // 1000}+{s % 1000:03d} | {gl[s]:.3f} | {dg[s]:.3f} |")
|
|
lines += ["",
|
|
"## 3. 검산",
|
|
f"- 전구간 직선 회귀 구배 {grade_pct:.4f}%는 **단일 구배가 아님**을 뜻하므로 채택하지 않음 (잔차 최대 {max_res:.3f} m)",
|
|
"- 구간별 구배 (계획고 2점 기울기, 도면 표기 S=2.4913%와 비교):"]
|
|
for a, b in [(4740, 4800), (4800, 4900), (4900, 5000), (5000, 5100), (5100, 5200), (5200, 5300), (5300, 5400), (5400, 5480)]:
|
|
seg = [k for k in stas if a <= k <= b]
|
|
if len(seg) >= 2:
|
|
g = (dg[seg[-1]] - dg[seg[0]]) / (seg[-1] - seg[0]) * 100
|
|
lines.append(f" - {a // 1000}+{a % 1000:03d} ~ {b // 1000}+{b % 1000:03d}: {g:.4f}%")
|
|
lines += ["- 해석: 4+740~5+300은 2.491%로 S=2.4913%와 일치. 5+300 이후 구배 감소는 종곡선 구간으로 보이나 [추론], 종곡선 제원은 도면에서 확인 못함 [확인 못함]",
|
|
"- 지반고는 구배 검산 대상이 아님 (지형 값)\n",
|
|
"## 4. 미확정·후속 (OPEN 등재 대상)",
|
|
"- 4+700~4+740, 5+480~5+500 종단 값 없음 → [미확정]",
|
|
"- 허용오차 ±50mm 판정: 본 초안은 도면 표기 값 그대로이며, 3D 대조 판정은 Task 7 이후 수행",
|
|
"- 라이닝 두께 수치: 04-004~009 콘크리트라이닝 구조도 확인 예정 (사용자 지시 2026-10-21)\n"]
|
|
open(OUT, "w", encoding="utf-8").write("\n".join(lines) + "\n")
|
|
print(f"stations={n} grade={grade_pct:.4f}% max_resid={max_res:.3f}")
|