• 検索結果がありません。

機械工学教室和田知之   〃    盛  中  清  和

N/A
N/A
Protected

Academic year: 2021

シェア "機械工学教室和田知之   〃    盛  中  清  和"

Copied!
5
0
0

読み込み中.... (全文を見る)

全文

(1)

切欠きをもつ異方性鋼板の降伏域の

         進展について

(第2報:異方性値を上降伏応力値と下降伏応力値        ともに考慮した場合)

(昭和48年10月20日 原稿受理)

トヨタ自工KK本田正則

大学院機械工学科阿部清輝

機械工学教室和田知之

   〃    盛  中  清  和

The Progress of the Yielding Zone of the Anisotropic        Steel Plate that has Notches

 (2nd report:Considering both upper and lower yield       strength as the anisotropic values)

by Masanori HONDA

  Kiyoteru ABE

  Tomoyuki WADA

  Kiyokazu MORINAKA

  By the Finite Element Method, we analyzed the elastic−plastic stress−strain of rolled steel plates under the plane stress axial tension. The steel plates were anisotropy and had notches. And we considered the progress of the yielding zone.

The result under considering both upper and lower yield strength agreed with

phenomena bettr than that under collsidering only upper yield strength.

1.まえがき       2.解 析

 すでに第1報1)において述べたように,有限要   表1に記号,図1に計算に用いた応力・ひずみ

素法を用いて,切欠きをもつ圧延鋼板を単軸引張  関係を図示する。図1のE は,後に述べる修整弾

り下で,異方塑性体として解析を行なった。    性係数である。

 前報では,単一の降伏応力のみを有する材料   計算は直接反復法2)によった。

(アルミニウムなど)についての解析を行なった        表1記   号

が,一般に用いられている低・中炭素鋼は,上降

伏応力,下降伏応力とも…る聴上・下降ミEN㌔欝 晋AX警⊇警

伏応力を考慮に入れて解析を行なった。      D  節点変位   YU 上降伏応力

鯉性の応力・ひずみ関係を縦等方鰍び1:Q繋鵬応力㌃憂麟;嵯縦

異方性についての〔D〕マトリクスは前報1)を見   SMAX最大相当応力     弾性係数・)

られたし。       SA  要素のひずみ  〆  塑性域での修整ボ

      アソン比1)

(2)

相 当

応 力

A=YU/SMAX 1〜NEM

SE=SE*A rEQ=SEQ*

1〜JMAX

RニR*A t=U*A

v=W*A

       降伏してない

      降伏?

SE=SE*A     前ステップで

SEQ=SEQ*   1降伏している このステップ        1       で降伏した要素        l   B_YL/YU        [       一        ・      SE=SE*B       

      ・。U・A l S㌫§1Q*B        U=U*B        [

       L _       相当ひずみ      ε

       図3 フローチャート  図4 フローチャート

図1 計算に用いた応カーひずみ関係図      (1)の部分       (H)の部分  降伏状件はミゼスの降伏状件を用いて,要素の  力ε Cに乗じ,最大の相当応力を有する要素を上 相当応力がデータで与えられた上降伏応力に達す  降伏応力状態にする。

ると,その要素を降伏させた。      図2の(n)は,図4に示すように,その荷重

 図2に計算のフローチャートを示す。図2の  で降伏させた要素のみにβ一γL/γσを乗じて,

(1)は,図3に示すように,降伏していない要  下降伏状態にする計算である。

素の中で最大の相当応力を有する要素の応力を,   上述の計算を行なって,次に示す修整縦弾性係 上降伏応力γσに達するまで,荷重を増加させる  数E ,修整ボアソン比〆を計算する。詳細は前 計算である。       報1)を見られたし。

 すなわち,メーγ0/SMλXを全ての要素の応        E       E

       E・ =1+E〃1E・ =1+Eゐ2

START

」1礎データの q込と印刷 ステップ数だ ッ繰返す 剛性マトリッ

Nスのr1卜算

戊本有重を加える R,D,SE,SEQの計算

「久伏していない要 fで最大の相当応 ヘの要素を探す

(1)

最大の相当応力を

Lする要素が嗣s

キるように荷重を

¥げる

START      ※     ※※

(II);㌶勇鵠鮭

 下降伏点まで下げ

r イ犬しプいる要素 の叩性域でのEIV をr,1嘉する

}げる      1

卒も

※     ※※

    図2 フローチャート       図5 解析対象材の寸法

    (1)(皿)は図3,図4参照       (応力集中係数3)α=3.5)

(3)

direction

P

R

由rec亘◇澄

Young,s modulus  (kg/mm2)

2、1×104

PoisO】α,sratio

◎.30

Strainhardenlng

 Coe伍¢ie砿

0.18

Stτe登gth C◎e麦董《cie茸t

  (kg/mm2)

Upρe∫Yieldρ◇《培   (kg/mm2)

10wer yieldρ◎i滋   (kg/mm2)

25.7  ]

Q6.1

24.6 Q5.0

P 圧延方向と平行に引張る場合  R 圧延方向と直角に引張る場合

y

剛体

荷 重       機械的性質を表2に示す。計算は図6のように・

      対称性から全体の四分の一の部分を解析した。す

      なわち,y軸上では, X方向への変位は零y方

      向へのみ変位し,X軸上では,γ方向への変位は       零,X方向へのみ変位するとした。

S 結果および考察

      図7,図8は荷重の増加に対して降伏域が進展

       していく様を試験片の四分の一について表わした        ものである。両図において,P方向, R方向は,

       表2に示すように,P方向は圧延方向と平行に引

       張る場合であり,R方向は,圧延方向と直角に引        張る場合である。図中の各曲線に付された引張り        荷重下での降伏域の広がりをあらわしたものであ        る◇

      図7は,前報1)のV切欠付試験片の降伏域進展       x    図である。これは,上降伏応力のみを考慮した場        合である。図より横軸方向への進展の程度が少な        いことがわかる。

      反力      図8は,上降伏応力・下降伏応力ともに考慮し      図6 計算対象の模式図 ・      た場合である。この図より,x軸方向へ向う降伏        域進展が,よくあらわれているのがわかる。

  ,_ρ+励、/2  , ッ+万島/2

 ガ1− 1十疏正 物= 1十Eカ2       異方性について・図8より次のことがわかる。

       X軸方向への降伏域進展はP方向に比べて,R ただし

       方向は同一荷重に対して,21%(2612kg),28%

 〃、=壁  〃、_肇       (2746kg)程度速く降伏して行くことがわかる。

   σ1       σ2

       また荷重の増加にともなってれR両方向の降伏

   6♪:相当塑性ひずみ      域進展の差が増加していくことがわかる。

   ∂ :相当応力       次に,切欠底から左斜上方への進展(軸に対し

解析の対象を図5に示す。計算に用いた材料の  て45°をなす方向一勇断応力が最大の面に相当

(4)

y    ・」  法      y    寸  }去

0  10   20   30

3・]

寸20   −一一一 法1。!

30・

寸20

1・!

        一  >x       O

   ぷ 寸《NΦ再oo◎うトト       ぷ 寸F司Φ〔寸oo うトト     トト◎り.う◎oΦoうo       卜卜 り りΦoo㊤oうo

    NNNN再〔〔r       《NNNN−H〔H←4

試験片に加わる引張り荷重      試験片に加わる引張り荷重

    P方向       R方向

  図7  上降伏応力のみを考慮した場合の降伏域進展図 0  10  20  30一

}一一一一一一      ∠    ._.一_.一

リ)o■うΦΦ【力Φoめ◎o

30

y    ・j  法      y     ・1  法

0 10 20 30    /0 10 20 30

  −一一」  一一.一一一一一    一一.」一      _一___.___L____一__⊥__

寸 法20、

10

30三

・1@i

法  i

  エ       し     

0}},、÷。.,。,N。一一・・   oL一皿西N。 一一・

  メ課品訳閑ε      ぷさ闘88閤ミ

     へふロ       ヘバヘト ぺ

試験片に加わる引張り荷市      試験片に加わる引張り荷重

    P方向      R方向

 図8  上,下両降伏応力を考慮した場合の降伏域進展図

電 kg 3000

2000

1000

       している)は,異方性の影響が顕著にみられない        ,、餐莞1{      ことがわかる。わずかにR方向の方がX軸へ傾

        !

     /      いていることがわかる。このことは前報1)の結果

     ノ

    ,・/       においても一致している。

    ノノ

   ./       図9は,切欠底要素の荷重とひずみの関係(注)を

 /      示したものである。これによると,荷重の増加と

,,

@       とも,P方向, R方向との差が増えていき,同一

・      荷重に対してP方向に引張った場合の方が伸び        やすいことがわかる。このことは,上に述べた荷

       重の増加にともなう,PとR両方向の降伏域進        展の差の増加と対応している。

・ α2 α4 α6 ・・8%  4結 論

      ひ ず み

図9 切欠底要素の荷重一ひずみ線図       本解析において,次のことがわかった。

(5)

  方を考慮すると,図8のように,それらを有  三角形要素1個で代表させているため・切欠底の荷重一   する材料(鉄鋼)の降伏域の進展を,上降伏  ひずみの関係をみるには・近似が荒くなっている。しか 応力のみを端する場合よりも実際の蜘こ ケ罐麗欝㌶漂豊㌶蕊1;

  近く解析できた。       かわらず,よい近似が得られていると思われる。

 2) 異方性は,x軸方向への降伏域進展に対し

  て顕著にみられ,また進展が進むとともに・      文    献

  異方性の影響が大きくなっていくことがわか   1)阿部,稲田他 九州工業大学研究報告(工学)第   った。       26号,P・38, S48・3

計算は九州大学大型計算機セ〃一で行なつ 2)法ρ富、鵠露 cz他マトリックス有限螺

た。      3)西田正孝 応力集中,p.ユ92,森北出版

参照

関連したドキュメント

Let F be a simple smooth closed curve and denote its exterior by Aco.. From here our plan is to approximate the solution of the problem P using the finite element method. The

Recently, Velin [44, 45], employing the fibering method, proved the existence of multiple positive solutions for a class of (p, q)-gradient elliptic systems including systems

Analysis and numerical results are presented for three model inverse problems: (i) recovery of the nonlinear parameter in the stress-strain relation for a homogeneous elastic rod,

Keywords: compressible Navier-Stokes equations, nonlinear convection-diffusion equa- tion, finite volume schemes, finite element method, numerical integration, apriori esti-

The finite element method is used to simulate the variation of cavity pressure, cavity volume, mass flow rate, and the actuator velocity.. The finite element analysis is extended

Based on these results, we first prove superconvergence at the collocation points for an in- tegral equation based on a single layer formulation that solves the exterior Neumann

In Section 2, we discuss existence and uniqueness of a solution to problem (1.1). Section 3 deals with its discretization by the standard finite element method where, under a

To derive a weak formulation of (1.1)–(1.8), we first assume that the functions v, p, θ and c are a classical solution of our problem. 33]) and substitute the Neumann boundary