くぴんのブログ

プロフィール

くぴん７４

フォローする

フリーページ

家電

(1)

バックナンバー

2025年11月
2025年10月
2025年09月

2025年08月
2025年07月

お気に入りブログ

5.　バリュー投資対… New! みきまるファンドさん

DOW 47716.42 +289.3… New! どらりん0206さん

楽天ラッキーくじ更…

じゃっかすさん

まさかの Maryu21さん

自分を愛することは… まりあのじいじさん

ニューストピックス

キーワードサーチ

▼キーワード検索

< 新しい記事

新着記事一覧(全5115件)

過去の記事 >

2015年07月05日

python sfepy 円孔の応力集中

テーマ：プログラミングあれこれ(231)

カテゴリ： CAE

単位系: mm,N,MPa
半径30の円孔がある200x200の平板を1/4対称としてモデル化。
これをY方向に20MPaで引っ張る。
円孔の応力集中係数は3なので、円の縁に60MPaの最大応力が発生する。

↓計算実行時のメッセージ
pyontaku14@pyontaku14-virtual-machine:~/sfepy/work_dir/plate_hole$ ../../simple.py plate_hole.py
sfepy: left over: ['verbose', '__builtins__', '__file__', '__name__', 'stiffness_from_youngpoisson', 'steel', 'data_dir', 'refinement_level', 'young', 'lame_from_youngpoisson', '__package__', 'stress_strain', 'refine_mesh', 'output_dir', 'poisson', '_filename', '__doc__']
sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (/home/pyontaku14/sfepy/work_dir/plate_hole/plate_hole.mesh)...
sfepy: ...done in 0.00 s

sfepy: Load
sfepy: Omega
sfepy: Y-sym
sfepy: X-sym
sfepy: ...done in 0.01 s
sfepy: equation "balance_of_forces":
sfepy: dw_lin_elastic_iso.2.Omega(steel.lam, steel.mu, v, u )
= dw_point_load.0.Load(Load.val, v)
sfepy: using solvers:
ts: no ts
nls: newton

sfepy: updating variables...
sfepy: ...done
sfepy: setting up dof connectivities...
sfepy: ...done in 0.00 s
sfepy: matrix shape: (644, 644)

sfepy: ...done in 0.00 s
sfepy: matrix structural nonzeros: 10640 (2.57e-02% fill)
sfepy: updating materials...
sfepy: steel
sfepy: Load
sfepy: ...done in 0.00 s
sfepy: nls: iter: 0, residual: 4.358899e+02 (rel: 1.000000e+00)
sfepy: rezidual: 0.00 [s]
sfepy: solve: 0.01 [s]
sfepy: matrix: 0.00 [s]
sfepy: nls: iter: 1, residual: 1.465871e-11 (rel: 3.362938e-14)
sfepy: equation "tmp":
sfepy: ev_cauchy_strain.2.Omega(u)
sfepy: updating materials...
sfepy: ...done in 0.00 s
sfepy: equation "tmp":
sfepy: ev_cauchy_stress.2.Omega(steel.D, u)
sfepy: updating materials...
sfepy: steel
sfepy: ...done in 0.00 s
pyontaku14@pyontaku14-virtual-machine:~/sfepy/work_dir/plate_hole$