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

有限要素法による不完全合成桁の弾塑性解析: University of the Ryukyus Repository

N/A
N/A
Protected

Academic year: 2021

シェア "有限要素法による不完全合成桁の弾塑性解析: University of the Ryukyus Repository"

Copied!
21
0
0

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

全文

(1)

Title

有限要素法による不完全合成桁の弾塑性解析

Author(s)

有住, 康則; 筑瀬, 明弘; 浜田, 純夫; 梶田, 建夫

Citation

琉球大学工学部(20): 69-88

Issue Date

1980-09

URL

http://hdl.handle.net/20.500.12000/14690

Rights

(2)

Elastic-Plastic Analysis of Composite Beams with Incomplete

Interaction

by

the Finite Element Method

by

Yasunori ARIZUMI, Akihiro CHIKUSE, Sumio HAMADA,

and Tateo

KAJITA

Summary

Recently the load factor method has been taking the place of the allowable

stress method for composite members in the United States and European

countries. It is the most important in the load factor design method to assure

the mechanical behavior such as ultimate capacity of the structure. Ultimate

strength is based on the properties of steel, concrete and shear connectors. In

the present study some discussions on ultimate states are provided for composite

beams with incomplete interaction applying elastic-plastic analysis.

Some studies on the elastic analysis for incomplete composite beams have

been presented in the bulletin by authors, and this is developed to inelastic

analysis based on bilinear elasto-plastic properties of concrete, steel and shear

connectors by means of the finite element method. An iterative initial strain

method is employed in the present analytical procedure, where inelastic strains

and deformations of materials are evaluated to the equivalent loads. This is

also applied to analysis of cracked reinforced concrete beams.

The main results of the analysis are as follows;

(1)

There is not much difference in the mechanical behavior of composite

beams with shear connectors between spaced equally and spaced proportinally

to horizontal shear.

(2)

Composite beams with partial interaction have an advantage of stress

reduction on longitudinal reinforcing bars in negative moment regions, where it

~H

:

1980~4F.l308

L*~XO)J1'J~O)-mH: -.:>~}-cfj:B*Wlmi§~~~13lID*~fiJf~~~~':i3Ii}-C~~I..,t.:.o

*

JfitER*~I~$ij!~I~f4

* *

J\=f{~.:r./'~.::.7 1) /,,y

*

* *

JfitER*~I~£f~±*I~'=f4

****

~~~*~I~~±*I~M

(3)

有限要素法による不完全合成桁の弾塑性解析:有住・筑瀬・浜田・梶111 70

maybenearthepointofcounterflectionThistendencybecomesmoreapparent

alongwithextentoftheplasticrange

KeyWords;CompositeBeamElastic-PlasticAnalysis,LoadFactorDesign,Initial StTainMethod,PartialCompositeBeam て複合構造物の弾塑性解析を行っている。それらの二 次元および三次元的な有限要素解析は,腹板の応力状 態を忠実に表わすことができるとか,FloorSystem等 への解析に適用できるなどの利点があるが,単純に- 本の合成桁を対象として,その非線形挙動を得ようと するような場合には,その手法や計算時間等に対して 般適なモデルとは思われない。 著者らはすでに,負の曲げモーメント域におけるコ ンクリートスラブの引張りの影撚を等価な力にI撒き換 える連続合成桁の解析法を示している;`)本研究では, これと同様な有限要素モデルを用いて,コンクリート, 鋼材,およびジベルの材料非線形を考慮できる単純な 有限要素解析法を示す。また,鉄筋コンクリート(以 下RCと称す)ばりの厚さ方向のひび割れの進展を考 慮した解析法を提案し,これを不完全合成桁の解析に 適艸Iした。ここで用いた非線形解析手法は,コンクリ ートスラブのひびWilれおよび塑性進展,鋼桁の鑓.性進 展,およびジベルのカーズレ関係の非線形性の影粋を すべて初期ひずみの項として取り扱う反復初期ひずみ 法である。 この解析法を)1】いて,RC単純ばり,単純合成桁およ び2スパン連続合成桁について解析し,すでに得られ ている実験結果と比較して,この解析法が実験結果の 挙動をよく示していることを確めた。また,ジベルの 配悩法の異なる単純合成桁,2スパン連続合成桁につ いて解析を行い,ジベルの配織法の違いによる合成桁 の弾塑性挙動への影響について比較検討した。 1.まえがき 近年,アメリカ,イギリス,および西ドイツ等にお いて,設計法が許容応力度設計法から荷重係数設計法 あるいは限界状態設計法に移行されつつある。!)~4)荷重 係数設計法,限界状態設計法,あるいは塑性設計法を 設計に導入する場合は,終局状態の構造物の安全性お よび安定性等力学的挙動を明らかにする必要がある。 特に,合成桁のように,コンクリートスラブ,鋼桁, およびジベルのような異種材料からなる複合描造物に おいては,それぞれの材料の性質が繊造物の終局耐力 におよぼす影響を明らかにしなければならない。この ようなことから,この研究では,コンクリートスラブ と鋼桁の接合面にズレが生じる不完全合成桁の弾塑性 解析を行い,このような構造物の終局状態について考 察を行った。 不完全合成桁については,最初Newmark5)によって 研究されて以来,弾性挙動に関しては,Hoischenl1 IIeiling,アIHomberg,`)Plum,.)山本}・Iおよび橘11)ら がNewmarkと同じ考えで解析しており,また, MatlockJ21GustofsomJ3)MoffattJ`DAnsourianJ5)お よび著者``)'7)らが有限要素法を用いて解析している。 一方,弾塑性挙動に関して,DaiJ81Hamada〉91および 前田2oI2l1らは,接合面のズレを考虚して差分法により 解析を行っている。また,Yam221231らは基本的には Newmarkと同じ方程式を用い,予測子一修正子法に より解析を行っており,Ansourian24)らは,ジベルのカ ーズレ関係を多項式で仮定し,断面での力のつり合い を求め,これにより解析を行っている。 有限要素法を用いた合成桁の弾塑性解析については, 佐藤251およびWegmuller配り27)らの研究がある。佐藤ら は,ジベルの変形による合成桁の剛度の低下をコンク リートスラブの弾性係数を低減したものに置き換えて, 平面三角形要素を用いて解析しており,また, We鷺muIlerらは,積層要素を用い,コンクリートスラ ブのひび割れを考慮して多主桁橋の解析を行っている。 しかし,WegmuIlerの解析では接合面のズレは考慮さ れていない。一方,栖原28)らは,平面三角形要素を用い YDV V1jV1 ロロ

FiglNodalDisplacementforanlncomplete

CompositeBeamE1ement

(4)

琉球大学工学部紀要第20号,1980年 71 猟塑性,およびジベルのカーズレ関係の非線形性が考 噸できる式を誘導し用いている。この論文における非 線形解析は、塑'性およびコンクリートスラブのひびW11 れの彩瀞等を初期ひずみの項として置き換え,反復初 期ひずみ法によって行った。 (1)仮定 この解析に用いた基礎的仮定は次のとおりである。 (a)ひずみは,要素内の厚さ方向に線形に変化して いるものとし、ひび割れが生じてもなお平面保持の法 則が成立するものとする。 (b)鋼材の応力一ひずみ関係は,Fig.2(a)に示すよ 2.解析法 ここで対象としたのは,コンクリートスラブと鋼桁 の接合面にズレが生じる不完全合成桁であり,Fig.1 に示すように,コンクリートスラブと鋼桁をそれぞれ 軸力と曲げを受けるはり要素で,ジベルをコンクリー トスラブと鋼桁の接合面に働く水平せん断力のみに抵 抗するスプリング要素でモデル化している。解析には, Armen291らによって用いられたはりの塑性進展を考 噛した婆索を用いた。この他に、ここではコンクリー トスラブのひび割オLによる剛性の低下の影辮,鉄筋の 0 0 0 E △s (□)STEEL(b)CONCRETEに)SHEARCONNECTOR Fig、2AssumedStress-StrainDiaRTamsforSteelandConcrete,andLoad-SIip DiaRramforShearConnectors

ルハ`''一IBM…

となり,また,応力が弾性限度を越え,塑性状態に達 した後の全ひずみごhは,弾性成分Efjと塑性成分EVj とから成り,つぎのように示される。

どルー6$+e;’2’

はり要索は一次元の要素であるので,応力およびひ ずみは機軸方向(z方向)のみを考えるとよい。ひずみ とIは,軸方向変位(Ⅸ')および曲げによる変位(【') を用いて,つぎのように与えられる。

6J-(`ねj/ヒノ2)-(。',/ヒノ22)y(3)

うに,引張および圧縮側とも完全弾塑性体とし,ひず み硬化はないものとする。 (c)コンクリートの応力一ひずみ関係は,Fig.2(1)) に示すように,圧縮側は完全弾塑性体とする。引張側 は,引機限度までは弾性とし,それ以後の荷重の墹分 に対してはひび割れが生じ,その領域のコンクリート は引振りに抵抗できないものとする。 (。)ずオし1こめのカーズレ関係は,Fig.2(c)に示すよ うにbilinearで表示し,猟性限の力Qyおよび弾性の 111大ズレ(△曽)mmxは,実験結果との比較により数値実 験的に求めるものとする。なお,コンクリートスラブ と鋼桁の間の付着および摩擦の影響は無視する。 (c)コンクリートスラブと鋼桁の間の浮き上がりは ないものとし,コンクリートスラブと鋼桁の同一鉛DII 上のたわみおよびたわみ角は等しいものとする。 (2)はり要素の塑性による初期ひずみマトリックス 仮想仕蝋の原理によれば,応力テンソルを可j,,ひず みテンソルをEIj,変位をu,,単位面積当りの力をP,と すれば,つり合い方程式は, 一方,要素の一部が塑性に達した場合,塑性領域内の 鋤性ひずみの分布は,要素内において線形に変化して いるものと仮定する。つまり,Fig.3に示すように各 節点における上下面の塑性ひずみの値E:,,e;”ど;,,

(5)

有限要素法による不完全合成桁の弾塑性解析:有住・筑瀬・浜田・梶田

72

-.に:;

●●■■ 』u▽

|・△-------‐一N‐‐‐‐‐-‐‐Cl

ラ′毛

mLLwi

_I

Vrvi

.vI

陸封

Fig3TypicalElastic-PlasticBeamElemelut ここで,9,-$はつぎに示す形状関数である。 ど;jを用い,上面の塑性領域において,塑性ひずみEH

は,常(’。富)-{e:j(1-β)+6Aβ)

(。‘一ノーy)/ヒノ“(4-α)

下面の塑性領域において,塑性ひずみGIは,

‘;(y,:)-{e;,(1-β)+§Aβ)

(/+y-d/)/(2ノーdノ)(4-6)

とする。ここで,2tは要素の厚さであり,β=z/ノで ある。また,。u,。Iはそれぞれ上面から弾塑,性の境界面 までの距離であり,つぎのように与えられる。 。、`(z)=d利1-β)+d;β(5-a) 。!(z)=d;(1-β)+djβ(5-b) 91=2β3-3β2+192=(β3-2β2+β)ノ 93=-2β』+3β29F(β3-β2)ノ(8) 一方,はり要素の変位ベストル(u)はつぎのようにな る。 (u}=〈w,,w;,vi,W,wj,w為vj,v;>『(9) 以上の結果を用いると,はり要素の内部仮想仕事が叫

憐川伊…l…

週ル…

-(盤Ⅱ畑…[…}

-伊11》…w此

却[rT…伽…}

(10-口) となる。式(10-a)をマトリックス表示するとつぎの ようになる。

dUb=(がu}’〔Kb〕{u}-{6u}T〔K:〕に,)

(10-b) また,はり要素の塑性ひずみベクトル{ep}は,つぎの ように与えられる。

(9,)-〈・か`ル,ご;i,ご内〉ア(6,

不完全合成桁においては,ジベルによって伝達され た水平せん断力によって,コンクリートスラブと鋼桁 に作用する軸力は要素内において複雑に変化し,軸方 向変位も重要な要素となる。したがって,ここでは軸 方向変位w(z)および曲げによる変位v(z)をそれぞれ 二次多項式で仮定する。Fig.1に示すように,各節点 で軸方向変位w1と軸ひずみw1,およびたわみviとた わみ角v;を未知数とすると,変位w(z)、vに)は,つぎ のように示される。 w(z)=wig,+wig2+wjg3+wjg4 (7-a) v(z)=vjg,十V錘十vjg3十vjg4 (7-b)

(6)

琉球大学工学部紀要第20号,1980年 73 ここで.〔K・〕,〔K:〕はそれぞれはり要素の剛性マト リックスおよび初期ひずみマトリックスである。〔KID を付録に示す. 一方,外力による仮想仕覗は,y方向に分布荷IRq,、 およびA1(i=1,…、)点に集中荷重pylが作用し, z方向に分布荷UIiqz,およびβh,(i=1,…、)点に災 中荷敢Pz1が作用していると, の剛性マストリックスおよび塑性による初期ひずみマ トリックスである。また,に;)は鉄筋の塑性ひずみベ クトルであり,つぎのように示される。

(cJ)=<(eハル(cZノ),,…………

(eハル,(e力)">アⅢ)

2 0α 〃 ⑤o ご α・ r1I1lw

[,伽…

6W=

〃 〃

+ZP9jO"j+ZPgjM

8=ID=1

(6邸)γ[P〕01)

となる。ここで,〔P〕は荷DII項である。 はり要素の場合と同様に鉄筋の塑性状態も考慮する ことができる。FiI9.4に示すようにRCばりのiii心軸 から距離。,(i=1,…、)の深さに断面欄AVの鉄筋が 集中的に配置されているものとする。鉄筋が降伏に入 った後の鉄筋の蛾性ひずみの分布をはり喫紫と同様に つぎのように仮定する。

(5;)j=(q;i)j(1-β)+(gル)jβ('2)

したがって,鉄筋饗素の内部仮想仕事は,

-4c「ト

Fig4CrossSectionofReinforced ConcreteBeams (3)コンクリユトスラブのひび割れによる剛性の低 下 コンクリートは,ある引張限度以上の外力には抵抗 できず,ひび剤れが生じる。コンクリートスラブにひ びW11れが生じろと,要素内に不つり合い力が生じ,ひ び割れは急激に進展すると考えられる。有限要素法を 用いたコンクリートばりのひび割れ解析については, Scordelies3I')の研究以来多くの研究報告がある。31’-34) それらの研究の多くは二次元の平面要素を用いて解析 を行っている。著者らは不完全連続合成桁の解析にお いて,負のIMIげを受ける区間のコンクリートスラブ要 素の剛性を低下させる解析法を示したが,この研究で は,ljlさ方向のひび割れを考慮できる解析法を示す。

"し…ル…

丑L・…

ラニ〔M[…-…

ル……MII…

古〔M[(…

-M`州M1a

-{624)了〔Kr】(")-(6脚)γ[〃〕(弓)(131

となる。ここで,〔Kr〕,〔KF〕は,それぞオL鉄筋要紫

玉」よ

L’二T

'--,--」

Fig5AnalyticalModelofCrackfor ReinforcedConcreteBeams

(7)

有限要素法による不完全合成桁の弾塑性解析:有住・筑瀬・浜田・梶田

74 すなわち,引張限度以上のひずみが生じている領域は すべてひび割れが発生しているものとみなし,その領 域の剛性を低下させ,この剛性の低下を初期ひずみの 項として取り扱った。なお,ここで考慮したひび割れ とは,曲げによるひび割れであり,せん断による斜め ひび割れについては考慮していない。 Fig.5に示すように,コンクリートスラブ要素の節 点i,jに深さh,,hjの曲げひび割れが生じて,要素内 の任意点におけるひび割れの深さを要素内で線形に変 化していると仮定すると,任意点のひび割れ深さh化) はつぎのようになる。 h(z)=hA1-β)+hJβ(15) ここで,A(zLG(z).およびI(Z)はつぎのように定義され る。

胴仙

叱F['二鰐

肺niw

6M『(1-β)+’β)

6{/2-(ノー〃121)2}/2

b(/3-(/-ん(2)ハハ

1181 なお,〔KC了孤ck〕は付録に示す。 (4)ジベル要素の初期ひずみマトリックス ジベルのカーズレ曲線をFig.2(c)のように完全弾塑 性体と仮定すると,ズレが弾性限度を越えた場合,ズ レ△sは,弾性成分のズレ△:と塑性成分のズレ△:を用 いてつぎのように表わせる。 △s=△:+△:(19 ひび割れの発生している領域においてコンクリートは 引張り力に抵抗できないものとすると,コンクリート スラブの内部仮想仕事6uは,

叶川:!…

Ⅲ。“…

ここで,塑性ズレ△Rを要素内において橋軸方向に線形 に変化していると仮定すると,

Ⅲコニ腿…

△:=△Kl-β)+△鯛 剛 となる。一方,合成桁をFiglのように,コンクリー トスラブと鋼桁をそれぞれ軸力と曲げを受けるはり要 素とし,ジベルを接合面の水平せん断力のみに抵抗す るスプリング要素にモデル化すると,ズレ△閾は,コン クリートスラブと鋼桁の軸方向変位IDごとい爵,および たわみ角U`を用いてつぎのように与えられるA⑥ (16)

=6tノb-6Uゲ

と得られ・ここで、dUcはコンクリートスラブが全断面 有効とした場合の内部仮想仕事であり,万U了はひび割 れ発生による低減内部仮想仕事である。コンクリート スラブ要素全域にひび割れが生じた場合のコンクリー トスラブの内部仮想仕事はdUc=0となる。コンクリ ートスラブのひび割れによる低減内部仮想仕事5U了 はつぎのようにマトリックス表示できる。

△s=zDs-zDc+ひ'y(21)

ここで,yはコンクリートスラブと鋼桁の重心間の距 離である。したがって,ジベル要素の内部仮想仕事 dusは,ジベルに作用する力をQ,Lii位長さ当りのジベ ル剛性をq圏とすると,つぎのようになる。

扉川野…

過{》…十.

……11"…}

={6〃)r[Kcmdb〕(〃)(171

麺一正…-ル。…

ル“…,

(8)

琉球大学工学部紀要第20号,1980年 75 の塑性,鉄筋の塑性,ジベルの非線形,およびコンク リートスラブのひび割れによる剛性の低下の影響を初 期ひずみの項として取り扱い,反復法によって計算を 行う。反復初期ひずみ法におけるつり合い方程式は次 式のようになる。

〔K〕{u}!=(PC}+{P)i-I㈹

式(7),(20を式似)に代入すると,ジベル要素の内部仮想 仕事はつぎのようにマトリックス素示できる。

006={M丁[Ks〕("}=(6両)γ[禰(“)

⑪ ここで,〔Ks〕はジベル要素の剛性マトリックスであ り,〔K:〕はジベル要素の塑性による初期ひずみマトリ ックスである。〔K:〕を付録に示す。また,{u),(△:) は,不完全合成桁要素の変位ベクトルおよびジベル要 紫の塑性変位ベクトルであり,次式で与えられる。 (u)=〈w,,w1,V,,v(,w2,W2,w3,w3,V3,V3, w4,wl>『卿

(△:)=〈△K△分丁(鋤

ここで.〔P〕i-1は初期ひずみによる修正荷重項であ る。式(26)において,一度前進消去を行えば,後は各反 復段階で修正荷III項の計算を行い,後退代入を行えば よい。計算手順をFig6のフローチャートに示す。 不完全合成桁において修正荷重項は次式で示される。

{P)=[〃Tl{ご;1-1+〔〃f(eZ

i-1-8-1ト1

+〔KCγαcAJ(〃}+[Kb*〕{△,)

8-1

(6J)

U-1 (5)解析手順 以上示したように,コンクリートスラブおよび鋼桁 (2m FORWARDREOUCTIDN BACKWARDSUBSTITUTION CALCULATIONOFTOTALSTRAIlVANDPLASTICSTRAIIV CALCULATIONOFTOTALSTR CALCULATIDNOFELASTIC-P CALCULATIDNOFELASTIC-PLASTICBOUNDARYDEPTH CALCULATIONOFCRACKDEPTH CALCULATIONOF N、 i<300 1< PmNT YES ̄ TIALSTRAINHATRIX CALCULATIONDFINITIALSTRAINOOATRIX CALCULATIONDFIN 1

nT=~FJ-7

lP]=|P。]゛化・]i-1[U]Ll

[K1i-1[U]Ll

丘L-f二E

~Uv1l- NO i■i+1

シ潅

<S PRINT STOP Fig6FlowChart

(9)

有限要素法による不完全合成桁の弾塑性解析:有住・筑瀬・浜田・梶朋 76 た,収束に際しては各部材力を計算し,つり合いが満 足されていることを確認したc計算における反復回数 は,桁の種題および各荷重段階によって異なるが,RC ばりの場合およそ20回前後,合成桁の場合は100回前後 であった。また,反復回数が300回を越えても収束条件 が満足されない場合は計算を打ち切った。その場合の 合成桁のコンクリートの圧縮ひずみは0.003を越えて おり,ACIで提案されているように破壊とみなした。 (1)実験結果との比較 RC単純ばり,単純合成桁および2スパン連続合成 桁について,他で行われた実験結果と著者らの解析結 果との比較を行った。比較に用いた桁の断面および材 料の諸mlをTable-Iに示す。 LoSh3s1らによって行われたRC単純ばりの曲げ試験 3.解析結果および考察 ここでは,RC単純ばり,単純合成桁および連続合成 桁について解析を行った。解析の際の分割数は,単純 ばりについては15分割,連続ばりについては20~24分 馴とし,ひび割れおよび塑`性の進展が予測される付近 を細かく分割した。なお,桁および荷重の対称性を考 慮して桁の半分について計算を行った。一方,反復計 算における収束の判定は,各反復段階のたわみについ てノルムで誤差を計算し,この誤差が許容値より小さ いかどうかで行った。誤差cは次式で与えられる。

c=||〃'一Dj-lll/||〃jll(28)

なお,誤差cの許容値としてLOx10-`を用いた。ま TablelSummaryofSectionandMaterialProperties (a)SectionProperties (b)MaterialProperties CASEA:ReinforcedConcreteBeams CASEB:SimplySuppo「tedCompositeBeams CASEC:2-SpanContinuousCompositeBeams 結果と解析結果との比較を荷重一たわみ曲線について FIR、7に示す。図には,鉄筋断面積がO4in2(p=1.25%) と0.2in2(p=0625%)の二通り示したが,著者らのひび 割れの解析結果は比較的よく実験値と一致している。 しかしながら,鉄筋麺が少ない場合は計算結果の剛性 が低い。これは,著者らの方法では,ひび割れ発生領 域の応力をすべて解除しているため,多少桁剛性の過 小評価となるものと思われる。したがって,鉄筋量が 少ない場合は,残留応力の概念を導入して除々に応力 を逓減させることも必要であると考える。しかし,鉄 筋鹸が比較的多い場合や合成桁の一部としてのコンク リートスラブのひび割れの解析としては,著者らの解 析法は有効であると思われる。 不完全合成桁の解析では,ジベルのカーズレ曲線を bHinearで仮定し,弾性限の力Qyおよび弾性最大ズレ (△s)…は,実験結果との比較により数値実験的に

求めた。Fig.8に示すように,Chapman36)らによっ

て行われた姥in径スタッドジベルの押し抜き試験の 結果より,5つのカーズレ曲線(ND1~NbL5)を仮定 し,それぞれの仮定を用いて,単純合成桁について解 SPAN CONCRETE WIDTH THICKNESS STEEL SHEARCONNECTOR NUMBER DIMENSION CASEA 10ft. 4.0 、. 8.0in. CASEB 18ft. 48.0 、. 6.0in. 12in、x6inx441b. 100 (1/2),'×2,、 CASEC 2×11ft. 24.0 、. 2.375in. 6in、x3in.×12lb. 46 (3/8)"×2pH CONCRETE YOUNG'S MODULUS EC COMPRESSIVE STRESS □cy SPLITTING TESILESTRESS ⑪ct STEEL YOUNG,S MODULUS ES YIELD STRESS ぴy

CASEA 3.Ox103(k〔 4.0(ksi) U4(ks ) 3.Ox10‘(ksi) 60(ksi)

CASEB 1330<to、/sq 、) 7250(psi) 725(ps) 13300(to、/Sq、in) 17.2(to、/sqin)

(10)

琉球大学工学部紀要第20号,1980年 77 8 (EC])ロロ迂○] 6 (咽◎三)□迂○] 4 2 5 DEFLECTION(inCheS) Fig9LoadDeflectionRelationshipsfor

SimplySupportedCompositeBeams

9 0 0020.406081.0 折を行った。Fig9に荷重一たわみ曲線について実験 値と解析結果との比較を示す。図から明らかなように, 解析結果は実験結果と比較的よく一致している。なお, たわみについては,ジベルのカーズレ曲線の仮定の相 違による差はほとんど見られなかった。 FilglOに荷重が40tonおよび45tonの場合のスパン方 向のズレ分布をNql~NO5の仮定について示す。この 図からは,5つの仮定の内M3の仮定が実験結果の性 状をよく示していると思われろ。ズレのピークが実験 値と異なるのは,カーズレ'111線をbilinearで仮定して いるためと考えられる。また,実験結果は支点上で小 さな値になっているが,これは試験桁にジベルが配置 された張り出し部がある影辨であり:`)また,コンクリ ートスラブと鋼桁の間の浮き上がりによる影粋も多少 あると思われる。しかし,コンクリートスラブと鋼桁 の間の浮き上がりに対する実験的な検証はあまりなく, その影響を計算に含めることは難しい。Nu3の仮定の みを用いた解析結果をFig.11に示す。図から明らか なように,Nu3の仮定は,極限状態に近い荷重段階に おいてはよくズレ分布を把握するものの,荷重段階が 小さい場合は多少ジベル剛性の評価が小さいものと思 われる。しかし,本解析は不完全合成桁の弾塑性解析 をおもな目的としているので,この仮定を用いてもよ いものと考える。 DEFLECTlONATMIDSPAN(inches)

Fig.7LoadDefIectiollReIationshipsfor

ReinforcedConcTeteBeams 00

…帷'w:レ

ウ ーーご_ 80 No.3 No.4 No.5 戸、 宗 一

・ ~ ● 60 〃 〃 ム0 20 0.06 0.02o04 SLIP(Inches) 0 Fi9,8AssumedLoadSlipDiagTamsfor ShearConnector 50 40 30 20 10 0 00.5101.

|’I

…・・・、EXPERIMENT -PROPOSE、 35) ANALYSIS

’●● 。 ● 、 、 ● P夕 ノ

/く

A■=0.4 A,=02

D d D C O ロ ノ P /

jiiijiiク

〃|’

|↓

L 30PC 6000 3002 12000 NO -2二0,G口xlOO (△:)二二二 12345 叩叩ね沁印 00185 0.0259 0.0143 0.0192 0010006 =

〆二

百マー

タデン

●▼ -0■

if

|||’

-.-EXPERIMENT -PRO 3句 POSED ANALYSIS P

■」 21600

(11)

有限要素法による不完全合成桁の弾塑性解析:有住・筑瀬・浜111・梶田

78

…--o-EXPERlMENT36)

---Nol ----No、2 -No.3 --●--No.4 ----No.5 0.03

(⑪①二UE-)。-J⑪ 002

:lIiiii、

、、「

緑■←

三三三:三戒H1

OLOL 0.01

二--.~~~~~矛

● 。 ■ ■■ ● P=45ton 0 0 80 40 40 80 DlSTANCEFROMMIDSPAN(inches) FiglOSlipDistribution

~…。…ExpERIMENT36)

-----PROPOSEDANALYSlS 0,03 P=45t。、 P=4Oton P=351on ,…..-q

P7,,、

P7,,、

(功①ニリ仁一)L一]⑩ 0.02 □ p p ゜

二二三=二二jIlllP二二三三I

全ニーミ、

、。、。 .◎豆 0.01 0

#二罰怠f

。●

脈へ■こ》■■覇

●■● ■bbC 巳。

〃〆P”

OG e■● 0 8040040 80 DISTANCEFROMMIDSPANmnCheS)

FigllSlipDistribution

荷1Kの70%に対応するズレの60%を用いbilinearで 表示する。この仮定を用いると,一般に多く用いられ ている%x4inスタッドジベルについては,Mai、‐ stonc37)の実験結果より,Qツー7.5ton,(△s)max=(). O14in゜となる. FiR、12にTeraszkiewicz(Yam23Dによる)による2 一方,径の違うスタッドジベルについてもこれと相 似なカーズレ曲線が得られることが知られているので, 他の径のスタッドジベルについても%in径スタッド ジベルと|司様な仮定を用いることができると思われる。 すなわち,Fig.2(C)に示すように弾・性限の力Q暖を最大 荷重の70%に取り,弾性限の般大ズレ(△S)…を,最大

(12)

琉球大学工学部紀鍵第20号,I98U年 79 21 12 (、の二UE-)z○一」U山]L四□

(いの二。Eこぐ○一xQ-]の

0.5 10 DEFLECTION 00 00 00 0, SLlP 己【し 1 11

の○一xzlくこ』の

皿血

知血緬

STRAINATBOTTOMFLANGEOFBEAM

L()nRitudinalDeflectionS1ipandStrainatLowerF1angeProfile(P=12.2ton・

Pu=14.0t()、) Figl2 スパン連続合成桁I/)実験結果と著者らの解析結果との(2)解析例 比較を示す。図には,終hij荷敢(Pu=14.0t(〕、)に近いジペル配慨 荷砿1,=12.2t(〕l1U)場合について,たわみ,ズレ,およめ単純合成桁 ぴ鋼桁下ソランジリ)ひずみ分布が示されている。ズレを行った。単 分布(3/8x2inスタツドジベル)について兇ろと,解ジベルを桁: 折結果は実験結果とよく一致しており,他の篠のスタ(MODEしA ツドジベルについても魁、侭スタッドジベルとIiil様に少なく三角 鞍カーズレ'''1線の仮定を)Iいてよいものと巷える。夫に,集KII荷![ た,たわみおよび鋼桁rフランジのひずみ分布も,解て解析を行つ 折結果は災験結果と比絞的よく一致しており,苫者らにわたってジ の解析法は,連続合成桁の弾塑・性挙動もよく把握できと!']央支点付 るものと忠」)オしる。慨しない,い ジペル配般の典なる合成桁の弾塑性拳!Iillを調べるた め単純合成桁および2スパン連続合成桁について解析 を行った。iii純合成桁に対し,FiRJ3に示すように, ジベルを桁企及にわたって等間隔に配慨した桁 (MODELA)とジベルを支点上に多く,スパン中央 に少なく三角形分布に配置した桁(MODELB)を対象 に,集KII荷111および等分布荷重を城荷した場合につい て解析を行った。また,連続合成桁では,スパン全長 にわたってジベルを等116隔に配置した桁(MODEL-C) と111央支点付近の負の曲げモーメント域でジベルを配 慨しない,いわゆる断続合成桁(MODELD)につい

(13)

80 有限要素法による不完全合成桁の弾塑性解析:有住・筑瀬・浜田・梶IH て,桁中央に集中荷重が作用した場合について解析を 行った。 TabIe2に解析に用いた桁の断面および材料の諸凪 を示す。 (a)単純合成桁 Fig.14,15に集中荷重および等分布荷血を戯荷した 場合の荷IIi-たわみ曲線を示す。これより,ジペル配 伍の違いによる桁剛性の差はさほど見られない。 MODEL-Aについて,集中荷虚および等分布荷重を受 けた場合の塑性域の進展状況をFigl6に示す。載荷 された荷iIiの迎いによって塑性域の広がりに差が見ら れる。この単純合成桁では,荷重が降伏荷jlliを越えて 増加するにともなって,桁中央の鋼桁下フランジより 塑性域が進展し,最後にコンクリートスラブが圧壊し て桁は破壊すると考えられる。なお,MODELBの破 壊形態についてもMODEL.Aと同様な結果であった。 Fig.17に典中荷重p=60,80,100tonを受けた場合の 軸力の分布を,Fig.18に等分布荷亜q=48,60kg/cm を受けた場合の軸力の分布を示す。図より,軸力につ いてモジベル配悩の違いによる差はほとんど見られな い。Fig.19およびFig.20に集中荷重および等分布荷 重を受けた場合のズレ分布を示す。これから明らかな ように.スパン中央の塑性域の進展によってズレが大 きくなり,その影響は,スパン中央より多少離れた場 所において顕著に現われる。現道路橋示方(1)30)では,ジ ベルはせん断力に比例して配置するよう規定されてい SIMPLEBEAM UNIFORMSPACING TR1ANGULARDENSITV C UNlFORMSPACING

 ̄。Ⅱ山L-u△

L DISCONTlNUOUSSPAC1NG

Figl3ConnectoTDensityforStandard

CompositeBeams Table2Summary (a)SectionProperties ofSectionandMaterialProperties WIDTHTHICKNIzsR NUMBERDIMENSION DI-I】NIRil【 (b)MaterialProperties YOUNGDS MODULUS Ec MODELAandB:SimplySupportedCompositeBeams MODELCandD:2-SpanContinuousCompositeBeams SPAN CONCRETE WIDTH THICKNESS STEEL SHEARCONNECTOR NUMBER DIMENSION SPACING MOnFLA MODELB MODELC MODELD 32m 32m 2×32m 2×32m 250cm 20cm Flg.P1.40×2.5 Web,Pll60xO、9 F19.PI.40×2.5 (c、) 378 378 2×378 2×378 l9x100 (m、) CONTINUOUS TRIANGULAR CONTINUOUS DISCONTINUOUS CONCRETE STEEL YOUNG,S MODULUS EC COMPRESSIVE STRESS 句cy SPLITTING TESILESTRFSS ⑦Ct YOUNG,S MODULUS ES YIELD STRESS cy 2.lx105 (kg/cm2) (kg/cm2)300 (kg/cm2)30 2.1xIO6 (kg/cm2) 2400 (kg/cm2)

(14)

琉球大学工学部紀要第20号,1980年 81 1 60 (Eo-)△ロベ○J 0 4

(Euへロー)ワロペ。]

20 0 48 DEFLECTlON(c、) 0 4812 DEFLECT10N(c、) 12

Figl5Load-DeflectionRelationshipfor

MODELqAandMODELBunderUni‐ formlyDistributedLoad

FiglILoad-DeflectionRelationshipsfor

MODELAandMODELBunderConcenL ratedLoad

11

P‐`…■……'・Ni卜rm-rrr帝芹伴泙卜

q=56kgに、

「F、可-F「「「「FrrrFrn-FFn

P=90t。、 q=60kgに、 q=64kgに、 P=100ton

(q)UNDERCOhにENTRATEDLOADl

Figl6YieldProgressionforSimply

(b)UNDERUNlFORMLyDBTRIBUTEDLOA、

YieldProgressionforSimPlySupportedCompositeBeamsunderCon‐

centratedLoadandUniformlyDistributedLoad 0 0 0 8 000 642

1 ̄-,5--J

 ̄-15,--」

 ̄ ̄

′IIII q  ̄「「「rrrrr「、.

‐仁=;;口

戸--

-UNIFORM SPACING ---TRIANd」LAR DENSITY P

(15)

有限婆紫法による不完全合成桁の弾塑性解析:イj住・筑瀬・浜111.梶ul

82 fIPl ▲」0【F■■■■0■ End End 006 00▲ (Eo-)山。仁。」]く員く 2 0 0 (Eu)Q-Jm 0  ̄051015 DISTANCEFROMEND(、) Filg・l9S1ipDiagramsunderC()ncentrated Load DlSTANCEFROMEND(、)

Figl7AxialForceDiaRramsfor

MODELAandMODELBunderConcent‐ ratedLoad ‘ q End End 0020 -UNIFORM SPAC1NG -・-TRIANGULAR DENSITY

一UNl1i

-._TRIA 400

q=60kgに、、

q=48kgに、ヘノ

q=60kgに、、

q=48kgに、ヘノ 0015 〈色。》)uU区。」]く一×く 300

({■U》Q|]の 0010 200 -UNIFORM SPACING ---TRIANGULAR DENS1TY -UNIFORM SPACING ---TRIANGULAR DENS1TY -UNIFORM SPACING ---TRIANGULAR DENS1TY 100 0005 q=60kgにm q=48kgノc、 q=60kgにm q=48kgノc、

q6

10 15 051015 D1STANCEFROMEND(、) DISTANCEFROMEND(、)

Figl8AxialForceDiaEramsfor

MODELAandMODELBunderUni‐

formlyDistributedLoad

FiR、20S1ipDiaRramsunderUniformly

DistributedLoad D UNIFORM ,へ 」一.--.-つ ̄ ̄へ11

三二二二己一一

(16)

琉球大学Iユ学部紀要第2Ⅱ》),1980年 83 る。本i汁in[例において6.11&人ズレ10(に荷目するなら ば,塑性域の進展にかか』ル〕ず,雌【|】荷血に対して(よ 等分布配Iivfが,霧クMlj荷'11[については三角形配搬とい うように,ジベルはせん断ノノに比例して配脱した万が 有利であると思われる。しかし,三角形配粧については, 塑性域の進展によるズレのAij部的ll1iljllについて配虚す る必喫があると考える。 (b)迎統合成桁 連続合成桁で'111趣となるのは,’'''111支点付近で負の 【M1げを受けることにより,=】ンクリートスラプにひび 割れが'Iiじることである。二Jンクリートスラプのひび ?illれの膨秤を小きぐするため,〔iの'''1げモーメント 域でジベルを配慨せずJ|:合成とする断続合成桁といわ れるW)がある。20庭!)MODELCが桁全長にわたって ジベルを配i綴した桁であり,MODILDが断続合成桁 である。解析は,スパン'11央に災'11荷iiiが作用する場 合について行った。 Fig.21にモーメントとIlIllll支点の'''1率の関係を示 す。’又Iより明らかなように,ジベルを配IITした桁の方 が断続合成桁より剛↑liが,|、liい。FiR、22にひび割オしおよ び塑性の進展状況を7Jくす。同図より,まず鹸初に''1間 支点付近の負のIM1け゛モーメント域のコンクリートスラ ブにひびWりれが生じ,荷1腕の1W川1とともにスパンⅡ'1央 および''''81支点上の鋼桁下フランジより塑`性域が進展 する。また,本解析例では,スパン[|'央において,鋼 桁下フランジの塑性域がかなり進展すると,コンクリ ートスラブ下ihiに多少ではあるがひび割れが生じた。 なお,MODELCとMOI)ILDでは,中'111支点付近に おいて,鋼桁の塑性域の進展に迎いが見られる。 FiR23に鋼桁下フランジのひずみ分布を示す。これ より,鋼桁下フランジのひずみ分布については,ジベ ル配liw〔の迎いによる錐はほとんど兄られない。Fig.24 1こ鉄筋のひずみ分布を示す。同図より,鉄筋のひずみ は,ジベルが配憧されている区IMIについては差は見ら れないが,断続区間でその違いは大きく,中間支点上 でMODELDの鉄筋のひずみはMODEL-Cの40%減 となっている。FiR、25に軸力の分布を示す。MODEL Dは非合成区間において一定1M〔となっており,その値 は,MODEL.Cのその区間のliji大11〔より小さい。Fi屑. 26にズレ分布を示す。同図よりIリIらかなように,非合

-1000

(E0Eo-)」z四三○三

800

600

400

-200

0-10-20-30-40-50-60-70

cuRvATuRE(1.9cm)

Fig21Moment-CurvatureRelationshipsattheCenterSupportfor2・SpanContinuous

CompositeBeams

-UNIFORMSPACIト$

一口 ̄

DISCONTINUOUSSF楓CING

 ̄ P ̄

岸一弓

ノ' グ グ □■■

30

30

(17)

有限要素法による不完全合成桁の弾鋼性解析:有住・筑瀬・浜[叶梶lIl

-伝 84 ■■PLASTlCREGION P E、顕CRACKREGION 鰯リリリリZ〃IilZZP P=130ton 『‘WZH UHZWZO qIUダムMn脇 P=162ton jIUWPJ"'擁 ZUiDIiIHlIy に)MODEL-C(b)HODEL-D

Fig22CrackandYieldProgressionfor2-SpanContinuousCompositeBeams

E I

suph.,,

F1 Pl ●、。 End 、のp□TI

(中OC-x)z一くこ』の

(c0.lx》z一江区』、

HX 01020。O DISTANCEFROHEND(、) DISTANCEFROMEND(、)

Fig.23StrainDiagramsinLowerF1ange

Fig24StrainDiagramsinReinforcing Bars 3000 2000 10 -1 -2 00 0 OOC 000 -3000 0102030 、“

、6

2 0 、 4

P=1621., P広1001.,  ̄UNIFORM ゆ SPACIN -.-OISCONTINUO SPAC F、弓 il CONNECTORDENS1TY CONNECTORDENSITY■ ■

(18)

琉球大学工学部紀要第20号,1980年 85 [ P 6 P I 0 -01 (仁。})四・匹。」ヨペーズく 2 3 0 0 (Eu)q-Jm -04 -0.5 0102030 DISTANCEFROMEND(、) 0102030 DISTANCEFROMEND(、) Fig25AxialForceDiagTams Fig26SlipDiagrams 成区|H1付近において,ジペル配慨の迎いによる差が現 われている。MODELDは断続点で大きなズレ,すな わちジベルに大きな力が作H]しており,その傾向は, 荷砿の1W加により顕籍に現われる。 以上の結果より,断続合成桁は,スパン全災にわた ってジベルを配置した合成桁と比較して,桁の剛性は 多少低下するが,負のllllげモーメント域の鉄筋ひずみ の低域という点では理論上有利であると忠わオしる。し かし,断続合成桁は,合成区1111と非合成区間の境界付 近でジベルに大きなブJが作用し,その傾向は塑性域の 進展によって顕著に現われる。したがって,断続合成 桁では,継続点のジベル補強および非合成区1111災等に ついて|分考慮する必要があると思われる。 鉄筋コンクリートばりのIMIげによる厚さ方向のひび割 れの進展を考慮できる解析法を示し,それを不完全合 成桁のコンクリートスラブの解析に適用した。鉄筋コ ンクリートばり,111純合成桁および2スパン連続合成 桁について,実験結果と解析結果の比較を行い,本解 析結果が実験結果とよく一致することを示した。 本解析法は,コンクリートおよび鋼材の応力一ひず み関係,およびジベルのカーズレ関係を簡単な仮定を 11】いて表わし,それらの材料非線形の影騨,およびコ ンクリートスラブのひび割れの影騨をすべて等価な力 にiiYiざ換え反復計算を行うため、他の非線形解析手法 より短い計算時'111で結果が得られるものと思われる。 また.本解析は,極々のジベル配糊された桁について も解析が可能である。 ジベル配麗法の違いによる不完全合成桁の非線形挙 1肋を調べるため,単純および連続合成桁について解析 を行い,次のような結論を得た。 (1)ili純合成桁について,ジベルを等間隔配慨およ び三角形分布配慨した桁を対象に,築中荷重および等 分布荷〕、【を'Mt荷した場合について解析した結果,ジベ 4.あとがき この研究では,コンクリートスラブと鋼桁の接合iIii にズレが生じる不完全合成桁について,コンクリート スラブおよび鋼桁の塑性進展,およびジペルのカーズ レ関係の非線形性を考慮できる解析法を示した。また, 6 4 2 00 00 00 0 -200 -400 End

Su0Inm、rt End Suppor ■ -UNIFORM SPAClNG -.-DlSCONTINUOUS SPAClNG 「、-1 EUNNECTORDENSITY ■■

--ヱーニ△

P=162ton P=100t、 -UNIFORM SPA C -.-DISCONTINUOUS SPACING 戸.司

ll

B2

(!

CONNECTOROENSITY ■。

(19)

有限要素法による不完全合成桁の弾蝿'性解析:有住・筑瀬・浜H1・梶Ⅱ1 86 ル配倣の違いによる桁剛性の差はさほど見られなかっ た。また,塑性域の進展によって接合面のズレが大き くなり,その影騨は桁中央点より多少離れた場所にお いて顕著に現われた。現行の道路橋示方書で規定して いるように,ジベルをせん断力に比例して配慨した桁 が,そうでない桁と比較して,最大ズレ麺が小さかっ た。 (2)2スパン連続合成桁について解析した結果,断 続合成桁は,負の[lIlげモーメント域において,鋼桁下 フランジの応力を墹加させることなしに鉄筋の応力を 低減させる利点を有するが,断続点でジベルに大きな 水平せん断力が作11】し,それは塑性域の進展によって 顕群に現われることが明らかとなった。今後,断続点 のジベル補強,非合成区間長,および非合成lXI1I1の鉄 筋賦等について,|分な実験的および非線形解析を含 めたFM論的研究が必幽であると思われる。 この研究における数値計算には,名古屋大学大KII計 算機センターのFACOM230-75を11lいた。 岐後に,この研究をとりまとめるにあたり,名iI「屋 大学成岡昌夫教授からiii重な助i;を賜ったことを記 し,謝意を表します。 7)Heilig,R:TheoriedesElastischen Verbunds,DerStahIbau,VOL5.pp、104-108,1953 8)Homberg,H,:BTUckemitelastischem VerbundzwischendenStahlhaupttrtigernund derBetonfahrbahntafeLDerBauingenieur,Vol 6,pp、213-216,1952 9)Plum,,.R、andM.R,H()rll:TheAnalysis ()fContinuousCompositeBeamswithPartial lnteraction,ProclCE,Part2,N059,pp、625643, 1975 10)111本稔:不光余合成桁のlllIげ理論,」:木学会論 文!↓L,N().67,1〕p」~10,1960 11)橘善雄,足立義雄:不完全合成桁について,土 木学会論文集,N()」12,pp、11~19.1964 ]2)Taylar,T・P,andILMailook:Finite E]ementAnalysisfoTCompositeBeams,Re‐ searchReport56-lO,TheUniv、Texas,Austin, 1968 13)GustofSon,W・CandRN・WriHht:Analysis 【)[SkewedCompositeGiTderBridRes,Proc. ASCE,Vol94,No.ST4,pp、919-94L1968 14)MoffatLKP・andP.T・KLim:Finite E]ementAnalysisofCompositeBoxGirder BridgesHavinRCompleteorlncompIeteInte-raction,Proc、ICE,Part2,pp、1-22,1975 15)Ansourian,P.:AnApplicationofthe MethodofFiniteElementstotheAnalysisof C()mpositeFloorSystems,PTCC、ICE,Part2, Vo1.59,pp669-726,1975 16)浜、純夫,有住康則:不完全連続合成桁の有限 要素解析,土木学会論文報告染,No.265,pp」~9,1977 17)有住康則,浜田純夫,梶IM1L夫:不完全合成桁 の有効幅,二}z木学会論文報告集,N().273,1〕p23~33. 1978

18)Dai,K、ILandC・P、Siess:Ana]ysisStudy

ofCompositeBeamswithlnelasticShear Connection,StructuralResearchSeriesNo267, UniMofl]linois,Urbana,I1linois,1963 19)Hamada,S、andJ・Long「worth:Ultimate StrenRthofContinuousCompositeBeams,

ProcASCE,vol・'02,NC,ST7,pp」463-1478,

1976 20)梶川靖治,iiiiⅡl幸雄:断続合成桁の塑性l#Iけ`性 状および曲げ耐荷力について,第32回年次学術職域会 識滅概要集,L301,,土木学会,1977 参考文献 1)AASHTO:StandardSpecificationsfor HighwayBridges,TwelfthEditionAmerica、

AssociationofStateHighwayandTransporta‐

tionOfficialSl977

2)BSI:TheDesignandSpecificationof

SteeLConcreteandCompositeBridges,Part5,

TheDesignofCompositeBridges,British Standardslnstitution,1976 3)前田幸雄:西ドイツ鋼合成桁設計施工指針(案) について,橋梁と基礎,Vol10,IVO、7,pp,1~6,1976

4)栗田章光:外国の新しい合成桁の設計法,jssc

VoLl3,No.143,pp29~39,1977 5)NewmarkN.M、,C・PSiessand1.M.Viest;

TestandAnalysisofCompositeBeamswith

lncompleteInte「action,Proc・oftheSocietyf()r

Experimenta]StressAnaIysis,VOL9.NC」.pp

75-92.1951

6)HoishenA.:Verbundtriigermitelasti‐

scherundunterbrochenerVerdUbelunR,Der

Bauingenieur,V()1.7,pp241-244,1954

(20)

琉球大学工学部紀要第20駅1980年 87 21)1W11仲次郎,前[H幸雄,梶111端治:断続合成げた の弾塑性''11げ性状および''11げ耐荷力について,第33回 年次学術織減会繊iii概要染,1-228,土木学会,1978 22)Yam,L.P.CandJ.C・Chapman:The lnelasljcBehaviouroISimplySupported CompositeBeamsofSteelandConcrete,Proc lCE,VCM5,pp651.683,1968 23)Yam,L.P.C、andJ、CChapman:The lnelasticBehaviourofContinuousComI〕osite Beamso「SteelandConcrete、ProclCE,Part2, pp、487-501,1972 24)Ams(〕urialLP、andJ,WRoderiele:AnfLlysis ofComp〔〕siIeBeams,Proc.ASCE,VoLl()`I,No. ST10,pp・'631薑1645,1978 25)奥村敏恋,佐藤政勝:有限要素法によるスタッ ドジベルを川いた合成桁の弾塑性解析,東京大学工学 部付凪綜合試験所年維,N().32,1〕p、59~68,197:1 26)WeRmLIⅡerAW、andlI.N,Amer:Nonli-nearResp()、SGO「CompositeSteel-Concrete Brid贋es,』.()fComI〕utersandStructures・VOL7, ppl61-169,1977 27)WeRmuuer,A、W、:OverloadBell【wiol・of CompositeSteelandConcrete,ProCASCE, Vol、103,NOST9.pI).1799-1819,1977 28)楠原訓綱ル西牧輿他:鋼板とコンクリートか ら描成さオしるサンドイッチ式複合織造物の強度に|10す る研究,日本造船学会論文災、Nol42,pp、312~322. 1977 29)Armen,HJ「.,APitcoandH・SLevine:A FiniteE1ementMethodforthePlasticBending AnalvsisofStructures,NASAContractors Report,CR97375,1968 30)NRO,DandA・CSC()rdelis:FiniteElemellt Analysis()「ReinforcedConcreteBeams,J・of ACI,VoL6`1,No3,pp152.163,1967 31)Nilsol1,A.H、:NonlinearAnalysisof ReiI1foreedConcretebvtheFiniteEIement,Lof ACI,VOL65,No.9,1〕p、757-766.1968 32)狐正iiIi治,滝口克己:有限要素法による鉄筋コ ンクリート部材の二次元非線形解析,建築学会論文報 告架,N()」89,pp51~57,1971 33)磯hI11i付:有限要索法によるコンクリーMII進物 の2次元illt塑性解析,雌蕊学会論文報告災,N〔)l89pp、 43~5011971 34)六11K煕,森[ll司郎,瓶H1幸次郎:Illlけ゛ひび制れ を生じた鉄筋コンクリートはり材の応力状態と変形に 関する研究,建築学会論文綴告染,No.200pp、27~34, 1972 35)Lasi.S、,.,N・MacleodandW、BlackwelI: HigllStrengthReinfOrcementinReinforced ConcreteBeams,Part2:CrackWidtll,Deflec‐ tionFatigue,ReportNo.38,Dept・ofCivilEng., Queen,sUniv.,RingstonOntario,1965 36)Chapman,』.C、andSBalaprishman: ExperimentsonCompositeBeams1TheStruc turalEnKineer,V01.42,N0.11,pp、369.382,1964 37)Mainstone,RJ、andJ.B、Menzies:Shear ConnectorsillSteelConcreteComositeBeams fol・BridRes,Concrete,V01.1,N0.9,pp,291-302, 1967 38)日本道路協会:道路橘示方ilr・同解税,1973 付録 式(1()-b)で示されたはり要索の塑性による初期 ひずみマトリックス〔K;〕は,つぎのように与えられ る。 11111011 1『と$4163荘 KKKKKKKPn 白 ワー⑰一⑪と、色⑪』⑰』⑤缶ワ』 ■0弔凸勺切口相GB点□句⑭ロロ ノ、ロンロ、』P〃、』〆、。〃、』PfL広々、』PfL TBLTOBT01TⅡ几T0LFBLT0LT0J 333333J” 囚腿隅隅瓜雌雌〃凪 0

叫鴎陶叫降恥払恥

一 〔K;〕=Eb ((1) ここで, KII=-(l/】O)Cl-(1/`l)C2 K1z=-(3/20)C1-(】/'1)Q K1:!=(1/10)C,+(】/4)C‘KlJ=(3/20)C、+(1/4)C4 K2l=-(//120)Cl+(//24)Q K亜=-(//30)C1-(ノノ2`l)C2 K2:!=(//120)C3-(//24)CJ K2m=(//30)C]+(ノ/24)Cl I<;,,=-(]/60/)C鮖一(l/6/)Cァ

(21)

88 有限要素法による不完全合成桁の弾塑`性解析:有住・筑瀬・浜田・梶田 Cl,=(7A,+7Aj)/(5/) C1z=Aj/10 cm=(3Gj-3G,)/(5/z) Cl`=-(4G,十GJ)/(5/) Cl`=A1/10 CIB=(G1+4Gj)/(5/) C22=(3A!+Aj>//30 C23=(7G+3G』)/(10/) C24=(13G,+2Gj)/30 C26=-(A1+Aj)〃60 C28=(8G1+7G」)/30 C33=(61,+61j)//3 C3F(41,+21j)/ノz C3`=-(3Ci+7Gj)/(10/) C犯=(21,+41J)〃2 C"=(31,+I,)// C46=-(7G1+8Gj)/30 Gs=(I,+Ij)/ノ C“=(Al+3A」)〃30 C`3=-(2G1+l3Gj)/30 Caa=(I,+31」)/ノ である。式伽)のジベル要素の初期ひずみマトリックス 〔K:〕は次式で与えられる。

[…命〔:|;二二二二二二二

重::!=::::|:二二二::か

K32=-(3/20/)C5+(1/い)CG+(1/6ノ)Cァ K33=(1/60/)Cs+(l/6ノ)Clo K34=(3/20J)C8-(1/6/)C,-(1/6/)OD K⑪=(1/180)C5-(1/36)C`-(1/6)Q K`2=-(1/30)Cs+(1/36)C6 K43=-(1/180)Cs+<1/36)C,+(1/6)C1o K``=(1/30)C8-(1/36)C9 K`、=_(//120)C1-(〃24)Cz KG2=<〃20)Cl+(ノ/24)C2 K63=(//120)Ca+(ノ/24)O K``=-(〃20)C3-(//24)C4 Ksu=-(1/45)Cs+(1/36)CG K82=-(7/60)Cs+(5/36)C6+(1/6)C7 K6a=(1/45)C`-(1/36)C9 KM=(7/60)C3-(5/36)Cg-(1/6)C’0 j d 2 t l j d jlj uId dl 二眼 2 9Ⅲ tCく 3ul Zld2一

.叫飢瑚釧釧却制

。u1

,吋朗叶胖慨価戯側瀦

た一一一一一一一一一一一一一一一一一一

まGQQQqQQqq

ここで, 咄則 j Ij lyly q j|yly qq

28q』〃印,何印崎刮印

j j

,口咄印

吋》“摘鐸評峅押咄》州搾輌赫

二一一一一一一一一一一一一一一一一一一一一一一一F|黍

12123344778311 Ⅱ212Ⅱ21之00 3a塾aSSSSSSSSSs である。 式(IDのコンクリートスラブのひび割れによる初期ひ ずみマトリックス〔KC…〕は, 〔Kcrack〕= 狸塑 (し(』 QQ〔し 33麺

qqqq

一CI1C16 -C12C26 -Cl3C36 -CMC46 Cl1-C16 C66 一Cl3C18 -C23C28 -C33C38 -C34C48 Cl3-C18 -C36Cc8 C33-C38 CB8 EC (b) metric で与えられる。ここで, である。

参照

関連したドキュメント

名の下に、アプリオリとアポステリオリの対を分析性と綜合性の対に解消しようとする論理実証主義の  

振動流中および一様 流中に没水 した小口径の直立 円柱周辺の3次 元流体場 に関する数値解析 を行った.円 柱高 さの違いに よる流況および底面せん断力

本研究では,繰り返し衝撃荷重載荷時における実規模 RC

2 解析手法 2.1 解析手法の概要 本研究で用いる個別要素法は計算負担が大きく,山

鋼板中央部における貫通き裂両側の先端を CFRP 板で補修 するケースを解析対象とし,対称性を考慮して全体の 1/8 を モデル化した.解析モデルの一例を図 -1

4 Ohta, H.: Analysis of deformations of soils based on the theory of plasticity and its application to settlement of embankments, Doctor Engineering Thesis, Kyoto

そこで本解説では,X線CT画像から患者別に骨の有限 要素モデルを作成することが可能な,画像処理と力学解析 の統合ソフトウェアである

Alternating-current Magnetic Field Analysis Including Magnetic Saturation by a Harmonic Balance Finite Element Method.By.. Sotashi Pamada,Member,Junwei