=== 2026-10-20 21:54:41  verify_solver.py

  solve: dof=2994 it=1 rel=9.37e-11 1.4s
cantilever h=5.0 tet=433 tip=0.19477 (an 0.19664, -0.95%) sxx=29.805 (an 30.000, -0.65%)

  solve: dof=11211 it=1 rel=2.35e-10 0.8s
cantilever h=3.3 tet=1924 tip=0.19504 (an 0.19664, -0.82%) sxx=29.967 (an 30.000, -0.11%)

  solve: dof=19812 it=1 rel=2.60e-10 1.0s
cantilever h=2.5 tet=3574 tip=0.19508 (an 0.19664, -0.79%) sxx=29.999 (an 30.000, -0.00%)

  solve: dof=9219 it=1 rel=9.35e-13 0.8s
torsion h=2.5 tet=1660 theta=0.007749 (an 0.007751, -0.02%) vm=88.49 (an 88.21, 0.31%)

  solve: dof=29724 it=1 rel=1.93e-12 1.6s
torsion h=1.6 tet=5971 theta=0.007751 (an 0.007751, -0.00%) vm=88.32 (an 88.21, 0.12%)

  solve: dof=103653 it=1 rel=4.04e-12 9.5s
torsion h=1.0 tet=22438 theta=0.007751 (an 0.007751, -0.00%) vm=88.25 (an 88.21, 0.05%)
--- 46 s  exit 0
=== 2026-10-20 21:55:26  hand_calc.py
Ra 7.653 Rb 12.524 (report 7.653 12.524)  M_gear 137.76 N·mm
y_gear bend -3.687e-05 total 4.969e-04 / end 9.048e-04 ; twist end 1.219e-04 rad u 6.097e-04
歯車のキー溝(4×4・Ø12)              M=68.9 T=570.3 sb=0.406 tau=1.681 Kt=2.55 Kts=2.91 vm_nom=2.939 vm_peak=8.536 s1=5.928
軸端のキー溝(3×3・Ø10)              M=0.0 T=570.3 sb=0.000 tau=2.904 Kt=2.38 Kts=2.77 vm_nom=5.030 vm_peak=13.920 s1=8.037
段 z=4(Ø14.5/Ø12・R0.3)        M=19.1 T=0.0 sb=0.113 tau=0.000 Kt=2.43 Kts=1.00 vm_nom=0.113 vm_peak=0.274 s1=0.274
段 z=15.5(Ø14.5/Ø12・R0.3)     M=107.1 T=0.0 sb=0.632 tau=0.000 Kt=2.43 Kts=1.00 vm_nom=0.632 vm_peak=1.532 s1=1.532
段 z=33(Ø12/Ø10・R0.5)         M=0.0 T=570.3 sb=0.000 tau=2.904 Kt=1.00 Kts=1.56 vm_nom=5.030 vm_peak=7.848 s1=4.531
平らな所 z=35(Ø10)               M=0.0 T=570.3 sb=0.000 tau=2.904 Kt=1.00 Kts=1.00 vm_nom=5.030 vm_peak=5.030 s1=2.904
--- 1 s  exit 0
=== 2026-10-20 21:55:27  mesh_shaft.py fillet coarse
shaft_fillet_coarse {'n_tet': 28966, 'n_node': 48787, 'min_scaled_jacobian': 0.20967327774439512, 'variant': 'fillet', 'level': 'coarse', 'volume': 5733.308313694591, 'volume_step': 5729.829555492926, 'sizes': {'hmax': 2.0, 'h_step': 0.3, 'h_key': 0.16}, 'path': 'mesh/shaft_fillet_coarse.msh'}
--- 37 s  exit 0
=== 2026-10-20 21:56:04  mesh_shaft.py fillet medium
shaft_fillet_medium {'n_tet': 62765, 'n_node': 103685, 'min_scaled_jacobian': 0.5487344082084985, 'variant': 'fillet', 'level': 'medium', 'volume': 5733.308313694591, 'volume_step': 5729.829555492926, 'sizes': {'hmax': 1.5, 'h_step': 0.17, 'h_key': 0.09}, 'path': 'mesh/shaft_fillet_medium.msh'}
--- 186 s  exit 0
=== 2026-10-20 21:59:10  mesh_shaft.py fillet fine
shaft_fillet_fine {'n_tet': 160198, 'n_node': 259423, 'min_scaled_jacobian': 0.10343613702821082, 'variant': 'fillet', 'level': 'fine', 'volume': 5733.308313694591, 'volume_step': 5729.829555492926, 'sizes': {'hmax': 1.2, 'h_step': 0.1, 'h_key': 0.05}, 'path': 'mesh/shaft_fillet_fine.msh'}
--- 853 s  exit 0
=== 2026-10-20 22:13:23  mesh_shaft.py r0 coarse
shaft_r0_coarse {'n_tet': 17168, 'n_node': 28826, 'min_scaled_jacobian': 0.21463802952795827, 'variant': 'r0', 'level': 'coarse', 'volume': 5729.829555492926, 'volume_step': 5729.829555492926, 'sizes': {'hmax': 2.0, 'h_step': 0.3, 'h_key': 0.16}, 'path': 'mesh/shaft_r0_coarse.msh'}
--- 11 s  exit 0
=== 2026-10-20 22:13:34  mesh_shaft.py r0 medium
shaft_r0_medium {'n_tet': 32736, 'n_node': 53673, 'min_scaled_jacobian': 0.530225649440445, 'variant': 'r0', 'level': 'medium', 'volume': 5729.829555492926, 'volume_step': 5729.829555492926, 'sizes': {'hmax': 1.5, 'h_step': 0.17, 'h_key': 0.09}, 'path': 'mesh/shaft_r0_medium.msh'}
--- 23 s  exit 0
=== 2026-10-20 22:13:57  mesh_shaft.py r0 fine
shaft_r0_fine {'n_tet': 55459, 'n_node': 89561, 'min_scaled_jacobian': 0.6158476852323155, 'variant': 'r0', 'level': 'fine', 'volume': 5729.829555492926, 'volume_step': 5729.829555492926, 'sizes': {'hmax': 1.2, 'h_step': 0.1, 'h_key': 0.05}, 'path': 'mesh/shaft_r0_fine.msh'}
--- 87 s  exit 0
=== 2026-10-20 22:15:24  mesh_housing.py coarse
housing_coarse {'n_tet': 18761, 'n_node': 36304, 'min_scaled_jacobian': 0.13708230427384327, 'level': 'coarse', 'volume': 124676.11588011465, 'sizes': {'hmax': 4.0, 'h_local': 1.5}, 'path': 'mesh/housing_coarse.msh'}
--- 7 s  exit 0
=== 2026-10-20 22:15:31  mesh_housing.py medium
housing_medium {'n_tet': 42256, 'n_node': 77372, 'min_scaled_jacobian': 0.0963233722699222, 'level': 'medium', 'volume': 124676.11588011465, 'sizes': {'hmax': 2.8, 'h_local': 1.0}, 'path': 'mesh/housing_medium.msh'}
--- 14 s  exit 0
=== 2026-10-20 22:15:46  mesh_housing.py fine
housing_fine {'n_tet': 101647, 'n_node': 173927, 'min_scaled_jacobian': 0.1087317840108388, 'level': 'fine', 'volume': 124676.11588011465, 'sizes': {'hmax': 2.0, 'h_local': 0.7}, 'path': 'mesh/housing_fine.msh'}
--- 95 s  exit 0
=== 2026-10-20 22:17:20  run_shaft.py fillet coarse

shaft_fillet_coarse nel 28966 dof 146361
  total load {'Fx': 6.900977, 'Fy': 18.960278, 'Fz': 0.0, 'Mz': 0.0, 'Mx': -369.725417, 'My': 134.569047}
  solve: dof=146361 it=1 rel=9.09e-13 8.7s
  vm_max 52.515 MPa at [77.99  3.6  21.78] / u_max 1.439e-03 mm at [80.62 -1.91 51.7 ] / 28s
   歯車のキー溝(4×4)              52.515 MPa at [77.99  3.6  21.78]
   軸端のキー溝(3×3)              52.194 MPa at [74.51  3.3  43.96]
   段 z=4(Ø12/Ø14.5)         0.416 MPa at [73.87 -5.63  3.82]
   段 z=15.5(Ø14.5/Ø12)      4.543 MPa at [80.26  4.25 15.69]
   段 z=33(Ø12/Ø10)          8.036 MPa at [75.48  4.99 33.37]
  reactions {'front': {'Fx': 2.7111928986506877, 'Fy': 7.422632134106432, 'Fz': -9.391705654544866e-13, 'Mz': -3.1596947280831955e-12, 'Fr': 7.90228035012447}, 'rear': {'Fx': 4.18978384588581, 'Fy': 11.53764564369136, 'Fz': 0.0, 'Mz': -2.7755575615628914e-15, 'Fr': 12.274834242246854}}
--- 30 s  exit 0
=== 2026-10-20 22:17:50  run_shaft.py fillet medium

shaft_fillet_medium nel 62765 dof 311055
  total load {'Fx': 6.900977, 'Fy': 18.960278, 'Fz': 0.0, 'Mz': 0.0, 'Mx': -369.725417, 'My': 134.569047}
  solve: dof=311055 it=1 rel=1.10e-12 22.1s
  vm_max 59.815 MPa at [77.95  3.55 21.59] / u_max 1.405e-03 mm at [80.78 -1.47 51.7 ] / 73s
   歯車のキー溝(4×4)              59.815 MPa at [77.95  3.55 21.59]
   軸端のキー溝(3×3)              57.983 MPa at [74.54  3.25 43.01]
   段 z=4(Ø12/Ø14.5)         0.424 MPa at [74.46 -5.81  3.77]
   段 z=15.5(Ø14.5/Ø12)      4.664 MPa at [80.05  4.44 15.73]
   段 z=33(Ø12/Ø10)          8.044 MPa at [76.06  5.02 33.35]
  reactions {'front': {'Fx': 2.701292884604858, 'Fy': 7.410012225292168, 'Fz': -3.1073786715263415e-12, 'Mz': -9.986303450837397e-11, 'Fr': 7.8870314077855825}, 'rear': {'Fx': 4.199683859925251, 'Fy': 11.55026555247899, 'Fz': 0.0, 'Mz': 8.43769498715119e-15, 'Fr': 12.290076438171548}}
--- 76 s  exit 0
=== 2026-10-20 22:19:06  run_shaft.py fillet fine

shaft_fillet_fine nel 160198 dof 778269
  total load {'Fx': 6.900977, 'Fy': 18.960278, 'Fz': 0.0, 'Mz': 0.0, 'Mx': -369.725417, 'My': 134.569047}
  solve: dof=778269 it=1 rel=1.37e-12 73.5s
  vm_max 62.648 MPa at [74.56  3.24 44.77] / u_max 1.429e-03 mm at [80.7  -1.71 51.7 ] / 234s
   歯車のキー溝(4×4)              62.095 MPa at [77.96  3.55 21.8 ]
   軸端のキー溝(3×3)              62.648 MPa at [74.56  3.24 44.77]
   段 z=4(Ø12/Ø14.5)         0.464 MPa at [73.3  -5.42  3.87]
   段 z=15.5(Ø14.5/Ø12)      5.011 MPa at [80.12  4.36 15.76]
   段 z=33(Ø12/Ø10)          8.143 MPa at [75.73  5.01 33.36]
  reactions {'front': {'Fx': 2.6855932520631645, 'Fy': 7.385906202842538, 'Fz': -4.146488606301697e-12, 'Mz': 2.6662005936373134e-11, 'Fr': 7.8590089421450005}, 'rear': {'Fx': 4.2153834924760405, 'Fy': 11.57437157493119, 'Fz': 0.0, 'Mz': -1.2212453270876722e-15, 'Fr': 12.318097878455692}}
--- 239 s  exit 0
=== 2026-10-20 22:23:05  run_shaft.py r0 coarse

shaft_r0_coarse nel 17168 dof 86478
  total load {'Fx': 6.900977, 'Fy': 18.960278, 'Fz': 0.0, 'Mz': 0.0, 'Mx': -369.725417, 'My': 134.569047}
  solve: dof=86478 it=1 rel=7.96e-13 5.5s
  vm_max 43.866 MPa at [74.5   3.26 45.31] / u_max 1.439e-03 mm at [80.62 -1.91 51.7 ] / 17s
   歯車のキー溝(4×4)              40.449 MPa at [78.    3.57 20.92]
   軸端のキー溝(3×3)              43.866 MPa at [74.5   3.26 45.31]
   段 z=4(Ø12/Ø14.5)         0.336 MPa at [73.09 -5.25  3.84]
   段 z=15.5(Ø14.5/Ø12)      4.128 MPa at [79.62  4.79 16.47]
   段 z=33(Ø12/Ø10)          8.262 MPa at [73.25 -4.17 33.15]
  reactions {'front': {'Fx': 2.6922110280116396, 'Fy': 7.450430315298384, 'Fz': -3.054588103724304e-12, 'Mz': 2.1674217975942156e-11, 'Fr': 7.9219260349024125}, 'rear': {'Fx': 4.208765716518542, 'Fy': 11.509847462435284, 'Fz': 0.0, 'Mz': 8.43769498715119e-15, 'Fr': 12.255215112965981}}
--- 19 s  exit 0
=== 2026-10-20 22:23:24  run_shaft.py r0 medium

shaft_r0_medium nel 32736 dof 161019
  total load {'Fx': 6.900977, 'Fy': 18.960278, 'Fz': 0.0, 'Mz': 0.0, 'Mx': -369.725417, 'My': 134.569047}
  solve: dof=161019 it=1 rel=9.82e-13 10.9s
  vm_max 60.976 MPa at [74.5   3.23 43.64] / u_max 1.420e-03 mm at [80.78 -1.47 51.7 ] / 33s
   歯車のキー溝(4×4)              54.629 MPa at [78.    3.54 21.47]
   軸端のキー溝(3×3)              60.976 MPa at [74.5   3.23 43.64]
   段 z=4(Ø12/Ø14.5)         0.383 MPa at [73.12 -5.27  3.91]
   段 z=15.5(Ø14.5/Ø12)      4.247 MPa at [80.04  4.43 16.3 ]
   段 z=33(Ø12/Ø10)          9.630 MPa at [76.95  4.91 33.09]
  reactions {'front': {'Fx': 2.663960810712275, 'Fy': 7.354559698011016, 'Fz': 1.9686622564605866e-12, 'Mz': 1.0826206597869259e-10, 'Fr': 7.822163099336314}, 'rear': {'Fx': 4.23701593382091, 'Fy': 11.605718079742328, 'Fz': 0.0, 'Mz': -1.7763568394002505e-15, 'Fr': 12.354958363746526}}
--- 36 s  exit 0
=== 2026-10-20 22:24:00  run_shaft.py r0 fine

shaft_r0_fine nel 55459 dof 268683
  total load {'Fx': 6.900977, 'Fy': 18.960278, 'Fz': 0.0, 'Mz': 0.0, 'Mx': -369.725417, 'My': 134.569047}
  solve: dof=268683 it=1 rel=1.22e-12 22.3s
  vm_max 69.557 MPa at [78.    3.52 20.31] / u_max 1.436e-03 mm at [80.7  -1.71 51.7 ] / 65s
   歯車のキー溝(4×4)              69.557 MPa at [78.    3.52 20.31]
   軸端のキー溝(3×3)              68.295 MPa at [74.5   3.22 42.27]
   段 z=4(Ø12/Ø14.5)         0.441 MPa at [73.06 -5.23  3.95]
   段 z=15.5(Ø14.5/Ø12)      4.072 MPa at [80.23  4.25 15.55]
   段 z=33(Ø12/Ø10)          11.419 MPa at [75.99  5.   33.06]
  reactions {'front': {'Fx': 2.679159736516605, 'Fy': 7.371047418954528, 'Fz': 9.231829446135643e-12, 'Mz': 3.0309088572266774e-11, 'Fr': 7.842846227374851}, 'rear': {'Fx': 4.221817008006291, 'Fy': 11.589230358834426, 'Fz': 0.0, 'Mz': 3.3306690738754696e-16, 'Fr': 12.334261192273363}}
--- 68 s  exit 0
=== 2026-10-20 22:25:08  run_housing.py coarse

housing_coarse nel 18761 dof 108912
  solve: dof=108912 it=1 rel=9.12e-13 5.6s
  load [ 4.28334744 11.76845462  0.        ] reac {'Fx': -4.283347435842391, 'Fy': -11.768454623529507, 'Fz': -8.611548973913585e-13}
  vm_max 0.2388 at [117.86  44.7   35.  ] (hole dist 0.00) / away 0.2248 at [81.68  8.83 28.  ] / seat 0.2248 / u_max 6.760e-05 at [77.49 10.39 28.  ]  16s
--- 19 s  exit 0
=== 2026-10-20 22:25:27  run_housing.py medium

housing_medium nel 42256 dof 232116
  solve: dof=232116 it=1 rel=1.55e-12 14.9s
  load [ 4.28351079 11.76839517  0.        ] reac {'Fx': -4.283510787633374, 'Fy': -11.768395168341408, 'Fz': 1.7242370725645273e-13}
  vm_max 0.2556 at [117.8   44.47  35.  ] (hole dist -0.00) / away 0.2247 at [82.88  7.94 28.  ] / seat 0.2247 / u_max 6.777e-05 at [77.49 10.39 28.  ]  37s
--- 40 s  exit 0
=== 2026-10-20 22:26:07  run_housing.py fine

housing_fine nel 101647 dof 521781
  solve: dof=521781 it=1 rel=3.06e-12 52.2s
  load [ 4.28340078 11.76843521  0.        ] reac {'Fx': -4.283400780760787, 'Fy': -11.768435207566291, 'Fz': -7.249652267393714e-13}
  vm_max 0.2717 at [117.85  44.66  35.  ] (hole dist 0.00) / away 0.2256 at [82.31  8.39 28.  ] / seat 0.2256 / u_max 6.786e-05 at [77.56 10.38 28.  ]  123s
--- 127 s  exit 0
=== 2026-10-20 22:45:42  export_glb.py
/files/lp-engineering/fem/web/shaft_vonmises.glb 0.003560318687615759 93.97258525398894 91402 182800
/files/lp-engineering/fem/web/housing_vonmises.glb 3.896286189473091e-05 0.4076132626578346 78012 156080
--- 2 s  exit 0
=== 2026-10-20 22:45:44  figures.py
figures done
--- 458 s  exit 0
=== 2026-10-20 22:53:22  report_pdf.py
/files/lp-engineering/fem/report/LP-GB-01_FEM_report.pdf
--- 7 s  exit 0
