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

点支持を有する偏平シェルの線形解析 : 行2辺が面内変位を拘束した固定支持, 他の平行2辺に点支持を有する偏平シェルの静的, 及び固有振動数解析

N/A
N/A
Protected

Academic year: 2021

シェア "点支持を有する偏平シェルの線形解析 : 行2辺が面内変位を拘束した固定支持, 他の平行2辺に点支持を有する偏平シェルの静的, 及び固有振動数解析"

Copied!
17
0
0

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

全文

(1)

点支持を有する偏平シェルの線形解析 : 行2辺が面

内変位を拘束した固定支持, 他の平行2辺に点支持

を有する偏平シェルの静的, 及び固有振動数解析

著者

皆川 洋一, 前畑 達実

雑誌名

鹿児島大学工学部研究報告

27

ページ

103-117

別言語のタイトル

Linear analyses of shallow translational

shells with point supports : static and

natural frequency analyses of shallow

translational shelld clamped along parallel

two edges and supported with points on the

other edges

(2)

点支持を有する偏平シェルの線形解析 : 行2辺が面

内変位を拘束した固定支持, 他の平行2辺に点支持

を有する偏平シェルの静的, 及び固有振動数解析

著者

皆川 洋一, 前畑 達実

雑誌名

鹿児島大学工学部研究報告

27

ページ

103-117

別言語のタイトル

Linear analyses of shallow translational

shells with point supports : static and

natural frequency analyses of shallow

translational shelld clamped along parallel

two edges and supported with points on the

other edges

(3)

点支持を有する偏平シェルの線形解析

一平行2辺が面内変位を拘束した固定支持,他の平行2辺に点支持を 有する偏平シェルの静的,及び固有振動数解析一

皆 川 洋 一 ・ 前 畑 達 実

(受理昭和60年5月31日) LINEARANALYSESOFSHALLOWTRANSLATIONALSHELLSWITHPOINTSUPPORTS −StaticandNaturalFrequencyAnalysesofShallowTranslationalShellsClamped alongParallelTwoEdgesandSupportedwithPointsontheotherEdges− YouichiMINAKAWAandTatumiMAEHATA

ThispaperadoptstheGalerkin'smethodtoanalyzeshallowtranslationalshellswithpoint

supports,andexaminestheeffectsofGaussiancurvaturesandboundaryconditionsondis‐

p

l

a

c

e

m

e

n

t

s

m

o

m

e

n

t

s

m

e

m

b

r

a

n

e

f

o

r

c

e

s

n

a

t

u

r

a

l

f

r

e

q

u

e

n

c

i

e

s

o

f

t

h

e

s

h

e

l

l

s

序 一般のシェルは,その形状に応じた座標系を導入し て解析しなければならないが,偏平シェルは偏平理論 を導入すると直交座標で解析することが可能となる。 この理論を導入して得られる偏平シェルを支配する Vlasov式に基づいて,偏平シェルは解析されるのが

通例である。Vlasov式に基づくことにより,4隅点

が同一平面上にあり射影面が矩形であるシェルは,辺

長比γ,Gauβ曲率比入,スパンライズ比ノu,及び

Poisson比,ノをパラメータとして,それらパラメー タの変化に伴うシェルの挙動の推移を統一的に調べる ことができる。

GeneralizedLevy解を用いて各種の境界条件の

もとで,偏平シェルの静的解析が行なわれている:''2I

その中で,周辺ローラー支持された完全HPシェルの

場合,曲率の効果がなく平板に近い挙動を示し,中央点

の変位においては,周辺ピン支持の場合と比較して約25

倍もの変位を示すことが報告されており?,副周辺で面内

変位を拘束することの重要性が指摘されている。しか しながら,現実のシェルが周辺ピンや周辺ローラーの 境界条件に近い形で利用されることは少なく,すそ梁 や小柱によって端部が支持されるものと考えられる。 一方,シェルを支配するパラメータがシェルの挙動 へ及ぼす影響を的確に把握するためには,一様静水圧 のもとでの静的解析のみでは不十分であり,これらパ ラメータがシェルの剛性行列の固有値及び固有モード へ及ぼす効果を調べる必要がある。固有振動数は重要 な固有値のひとつである。しかし,Generalized

Levy解を適用して固有振動数を解析する場合,特性

方程式の係数の中に未知の固有振動数が入ってくるた め,振動数方程式を解析するのが困難になることが報

告されている4)。

ここでは,静的解析,固有振動数解析の両方に適用 できるGalerkin法を利用して,点支持を有する偏

平シェルの挙動を調べる。研究の一環として6),平行

2辺で面内変位を拘束した固定支持とし,他の平行2 辺で任意の数の点支持を有する偏平シェルを解析する。 解析手順として,まず,平行2辺が面内変位を拘 束 し た 固 定 支 持 , 他 の 平 行 2 辺 が 自 由 辺 で あ る 偏 平 シェルの解析を行ない,次に点支持を有する偏平シェ ルの解析を行なう。点支持を有する偏平シェルは付帯 条件のある変分問題であり,自由辺を有する偏平シェ ルにおける自由辺上に点支持を付加し,その点におけ る変位拘束条件を新たに追加して解析される。

(4)

/ 104 グ 1.偏平シェルの形状 ここでは中立面の形状が次式で表わされる,射影面 が矩形で,4隅点が同一平面上にある偏平シェルを扱 う。(図1)

Z

=

;

l

+

(

'

-

6

/

6

'

2

)

,

ここに,R蕊,R3,は釘,z/方向の曲率半径,α,b は。c,Z/方向の辺長を表わす。 これにより,Gauβ曲率比入=R錘/Rgを変化させ れば偏平シェルは,E・P・シェル(入>0),円筒シ ェル(入=0),H・P,シェルいく0)と推移してい く。 / グ 2.偏平シェルの基礎式 このシェルについて面要素を作る線素の長さと,そ の線素のjC-z/平面への射影の長さが等しいと仮定 する。 シェル面が有するZノーconst曲線の接線方向変位 をu’工=const曲線の接線方向変位をひ,及び法線 方向変位を”とする。 また,中立面の範方向及び@/方向の歪をEコCO,E3,0, 面内せん断歪76cyoとして次式を採用する。 X b

y 図 l シ ェ ル の 形 状 グ

=

鋤 。 = 器 - 鈴 ’ 2 1

=

面内力Mc,W,NおZノ,モーメントMjc,M3/, M”,及び面内変位u,ひは法線方向変位”,応力 関 数 ’ を 用 い て 次 式 の よ う に 表 わ さ れ る 。 / 亭曹一一一一一一一・-一 グ ジ 鹿 児 島 大 学 工 学 部 研 究 報 告 第 2 7 号 ( 1 9 8 5 )

=

-

,

=

-

,

(

'

"

,垂=-,│等十仙器}I3l

w−Dl券十仙畿|

"=八念(券-幽窯)+淵血

"

=

(

-

,

)

+

ここに,E:ヤング係数,〃:Poisson比 D:板の曲げ剛性 偏平シェルを支配する汎関数汀の第1変分を応力 関数’,モーメントM工,M1/,MJcy,換算せん断力 Vjc,W,及びu,り,”で表わし次式を得る。

伽=〃(Dw’"一志穿一古窯

一P+,伽諜)伽加叶什

窯州幾MM’筈-W伽)d叶ノ

獅一蜘 到 N 抑一一師割

恥制

伽一師

券刊

一一一一

おお

NM

い=0) / 虎 (入>0) いく0) 図 2 シ ェ ル の 形 状 変 化

(5)

皆川・前畑:点支持を有する偏平シェルの線形解析

穿

-

6

,

-

M

+

V

)

(

4

ここに,h:シェル厚,β:質量密度 P:法線方向分布外力 (2)式のu,りを消去して得られる適合条件式を応 力関数‘を用いて表すと次式となる。

古,w+古穿十方三二=。

(5)

(4)式から得られるEuler方程式と(5)式がVlasov

型偏平シェルの基礎式である。 y

自由辺

司 由 必 図 3 自 由 辺 を 有 す る シ ェ ル 105 X

3.自由辺を有する偏平シェルの解析

図1に示すシェルにおいて,Z/=±b/2が自由辺で

あり,工=0,αにおいて面内変位を拘束した固定支

持であるシェルの解析を行う(図3)。また,ここで

は,〃=a/2,z/=Oに関して対称変形する場合を扱う。

解析手法は,法線方向変位切に関して支持辺では 幾何学的境界条件を満足し,自由辺での境界条件を導 入できるように仮定する。その仮定した変位”を用 いて(5)式の適合条件式より応力関数‘の一般解を厳

密に求め,(4)式のEuler方程式に代入し,Galerkin

法を適用する。それにより,離散化された運動方程式 が得られる。 この偏平シェルの境界条件は次式となる。 〃=0,αで〃=a”/ac=0,u=U=O (6−a,b) zノー±b/2でMzノーVz/=0,M/=NJcz/=0 (6−c,。) 法線方向変位は,(6−α)式の境界条件を満足し, Jc=a/2,z/=0に関して対称な変形を表わすように仮 定する。さらに,条件M1/=Vz/=0(z/=±b/2)が導 入できるように未知パラメータを追加し,それら条件 式をFourier級数に展開した時の定数項の処理の為 の代数関数も追加する。よって,法線方向変位めい, z/)は次式のように仮定する。

"

)

=

(

B

W

c

o

s

h

¥

+

B

'

s

i

M

K

K

×

(

c

o

s

-

,

)

+

(

"

+

"

c

o

s

(

c

o

s

"

ここに,m,72は偶数をとる。 (7)式を(5)式へ代入し,応力関数‘の特解dpは次式のように求まる。

,

=

[

[

_

B

(

,

_

(

,

A

)

s

i

M

(

)

s

i

n

c

o

s

h

,

.

(

-

)

s

i

n

;

,

S

i

m

.

(

-

)

×

c

o

s

"

+

(

B

W

-

)

c

o

s

h

g

s

i

M

,

+

K

{

(

2

-

,

l

c

o

s

+

K

[

f

+

(

3

2

c

s

(7)

(6)

106 鹿 児 島 大 学 工 学 部 研 究 報 告 第 2 7 号 ( 1 9 8 5 )

+

"

.

c

o

s

"

,

c

o

s

+

l

c

o

s

ここに,γ=α/6 ゆの斉次解は次式をとる。

,

=

[

(

C

c

o

s

h

,

+

S

i

m

,

)

c

o

s

l

D

W

c

o

s

h

¥

(

-

)

+

D

(

-

)

s

i

M

¥

(

-

)

}

c

C

s

+

'

3

(

-

)

,

1

+

(

-

)

+

(

-

)

-

3

(

-

+

+

(

)

-

'

5

(

r

-

"

+

'

5

(

-

)

(9)式中の代数関数は,条件式をFourier級数のsine展開するときGibbs現象を防ぐ為, る時に生ずる定数項処理の為に導入する。 応力関数の一般解‘=ゅp+ゅn,”を(3)式へ代入し,u(範,Z/),り(jc,Z/),NZノ(範,Z/), Mz/(Jc,z/),W(。c,z/)はそれぞれ次式のように求まる。

,

)

=

[

[

B

[

+

c

o

s

h

(

3

+

I

)

s

i

n

(

A

×

(

)

c

o

s

"

-

,

-

A

S

i

m

(

"

l

s

i

n

(

,

+

,

(

-

)

c

o

s

l

l

-

B

W

,

+

,

s

i

M

×

(

,

-

,

)

s

i

n

+

(

-

)

c

o

s

"

+

+

(

,

+

"

c

o

s

h

×

l

s

i

n

(

-

)

c

o

s

"

+

(

)

{

,

s

i

n

h

,

'

3

s

i

n

(

3

+

,

(

-

)

c

o

s

}

c

o

s

h

,

(

,

_

(

,

+

,

(

)

c

o

s

+

;

l

2

c

o

s

h

,

s

i

M

,

×

(

-

)

s

i

n

"

,

S

i

m

C

W

M

+

)

.

s

i

n

C

l

.

"

{

2

c

o

s

h

(

,

+

,

,

s

i

n

K

;

_

(

s

i

n

K

:

(

_

,

+

2

A

-

A

3

-

+

)

7

:

+

l

s

i

n

"

,

+

,

)

s

i

n

-

(

-

)

-尋"緬帝│げ『'千縛斉二噸l-llsin等…判

I

D

i

M

,

+

)

s

i

n

h

¥

(

)

.

c

o

s

.

"

"

(

,

+

×

(

-

)

c

o

s

h

¥

(

-

)

-

"

)

s

i

M

¥

(

-

)

l

c

o

s

(9) 及びcosine展開す NJcZ/(範,Z/),

(7)

γ262a 107

-

+

"

)

(

-

)

$

-

(

-

)

'

(

"

,

,

)

=

{

[

B

[

;

(

,

)

(

,

+

,

,

竿

,

,

-

(

,

-

×

s

i

n

y

l

(

-

)

s

i

n

"

-

(

)

S

i

m

,

-

(

A

竿

,

_

(

,

s

i

M

,

c

o

s

l

+

B

{

(

,

+

"

(

,

+

,

c

o

s

h

K

1

(

_

)

_

(

2

(

)

_

K

_

(

-

)

(

-

61ノZ/21

-

K

'

+

2

,

(

J

-

-

(

,

-

,

1

s

i

M

,

(

-

)

s

i

n

;

(

,

c

o

s

h

-

s

i

M

,

×

(

3

-

,

)

c

o

s

-

4

-

(

[

-

;

{

(

,

+

"

c

o

s

h

_

(

,

_

)

S

i

m

,

(

-

)

c

s

_

'

1

,

+

,

s

i

M

×

c

o

s

h

'

(

-

)

s

i

n

'

2

竿

s

i

M

'

-

3

(

×

-

s

i

M

,

)

l

c

o

s

"

l

l

-

C

,

r

(

,

+

"

)

s

i

M

×

c

o

s

-

C

l

.

(

,

+

,

c

o

s

h

(

,

-

l

s

i

n

h

,

c

s

+

{

;

(

'

'

'

7

c

o

s

-

1

+

[

'

0

-

,

-

2

'

7

;

+

:

I

2

-

2

ソ

H

A

M

l

c

o

s

I

1

-

"

.

-

"

[

(

(

'

'

-

A

}

c

o

s

1

s

i

n

[

D

γ

(

,

+

,

)

c

o

s

h

¥

(

-

)

s

i

n

+

D

{

2

c

o

s

h

¥

(

)

+

(

l

+

(

-

)

s

i

n

h

¥

(

-

)

l

s

i

n

+

K

'

(

l

+

2

1

-

÷

(

K

1

2

+

(

-

)

-

-

"

+

(

-

)

-

(

-

)

+

}

,

,

)

=

{

[

竿

s

i

M

l

2

c

o

s

(

-

×

s

-

{

2

1

(

2

+

S

i

m

-

c

o

s

h

}

皆川・前畑:点支持を有する偏平シェルの線形解析 (10 (11) 〃 称一α n s α|ワ︺ 一、lノ

範y

鞘挫α

、。n

〃︾S

鯉α耐制

、m2

・劃︵し

〃+

蝶端所与

︵uOc

−トー伽Ⅵ川川

当α釧岬|耐

、、+

Ⅲ回ノートJゾ|ゲ

α’2︲可判

紳畔勝

一vlノ凡

〃α−2+

階紳亦一

×+×

(8)

鹿 児 島 大 学 工 学 部 研 究 報 告 第 2 7 号 ( 1 9 8 5 ) 108

W

)

=

-

,

-

B

(

l

c

C

s

l

l

s

i

n

h

"

{

(

,

+

1

c

o

s

,

c

-

3

1

s

i

M

,

-

{

(

,

-

)

c

o

s

,

,

(

2

-

,

c

o

s

K

¥

{

_

1

2

_

c

s

+

K

(

)

_

(

)

c

o

s

-

"

A

c

o

s

-

"

(

c

o

s

"

W

I

D

W

c

o

s

h

¥

(

-

)

+

l

2

c

o

s

h

¥

(

+

¥

(

-

)

s

i

M

¥

(

)

]

c

s

¥

K

(

-

)

-

+

K

,

(

)

-

(

-

)

=

[

[

,

-

)

(

s

i

n

h

c

o

s

h

,

)

{

s

i

n

+

(

"

-

)

c

C

s

l

-

B

3

(

l

+

)

s

i

M

(

7

-

竿

'

c

o

s

h

+

2

(

,

_

s

i

M

,

{

s

i

n

(

_

)

c

o

s

+

(

,

_

(

s

i

M

+

,

c

o

s

h

,

H

2

(

)

c

s

(

_

)

s

i

n

}

1

C

W

S

i

m

,

+

C

l

(

S

i

m

,

+

c

o

s

h

)

s

i

n

( 12)

-

K

s

i

n

K

o

l

_

(

)

l

s

i

n

+尋"徳融け蒜鶏辿Sm等Ⅷn判

(14) (13

"

W

I

D

W

s

i

M

¥

(

-

)

+

D

:

{

s

i

n

h

¥

(

-

)

+

(

"

-

)

c

o

s

h

(

"

-

)

×

s

¥

,

+

(

-

)

(

-

)

K

,

{

(

-

(

-

)

M

'

)

=

-

,

[

[

B

'

c

s

-

1

1

c

"

{

(

,

-

,

,

)

c

o

s

"

-

,

,

S

i

m

,

+

2

(

c

o

s

"

-

,

)

c

o

s

,

応一α

婿些.

〃S

’O

陪恥

K伽

十−1叩

判別”

一一

応一α称一α

ゾーゲ叩ⅦⅡ皿

派恥砿荊

鴫判伽

十一一

(9)

皆川・前畑:点支持を有する偏平シェルの線形解析

"

{

"

:

-

1

c

o

s

"

_

"

r

1

s

i

n

(O−b)式の境界条件より(10,(11)式に範=0を代入して次式を得る。

"

(

=

[

[

-

;

(

)

c

o

s

h

(

,

)

(

,

+

,

s

i

M

-

[

;

1

1

+

A

I

c

o

s

h

'

+

'

3

(

3

-

'

s

i

M

,

竿

c

o

s

h

,

l

+

;

"

l

+

"

"

,

{

+

×

s

i

M

竿

D

(

,

+

竿

c

o

s

h

-

(

,

-

I

s

i

n

h

l

l

c

o

s

¥

+

K

l

-

;

+

(

2

+

+

K

+

K

+

I

H

,

_

2γ262

-

K

(

'

5z/2

+

=

[

{

(

,

,

-

I

l

-

A

-

4

×

s

i

M

,

[

1

,

+

A

+

4

,

1

(

c

o

s

h

_

s

i

M

,

+

[

(

2

竿

s

i

M

,

-

3

,

,

+

3

s

i

M

,

(

,

+

,

c

o

s

h

,

-

(

,

-

,

l

-

C

×

s

i

M

-

C

l

.

(

,

+

,

c

o

s

h

,

-

(

,

-

l

s

i

M

+

{

l

H

I

-

;

+

'

+

(

2

-

3

2

+

'

1

-

'

2

+

,

'

-

_

"

_

"

I

"

'

+

,

l

s

i

n

{

D

W

l

+

c

C

s

D

'

'

2

c

o

s

h

竿

(

l

+

s

i

M

l

s

i

n

K

l

2

I

l

+

2

1

+

K

K

(

:

1

2

M

-

2

一M号+M州│器-器+淵l

(0−c,。)の境界条件より,(12)∼(15)式にZ/=b/2を代入して次式を得る。

MM/21=器[尋[一端1-入祭Sim等l2cos等”

_

(

_

)

s

i

n

{

2

1

1

s

i

M

0

-

A

×

l

2

c

o

s

(

-

)

s

i

n

-

(

s

i

M

'

3

(

-

)

s

i

n

+

(

-

)

c

o

s

'

1

-

(

s

h

α

S

i

m

)

c

o

s

109 (16) ( 17)

(10)

110 鹿 児 島 大 学 工 学 部 研 究 報 告 第 2 7 号 ( 1 9 8 5 )

K

_

(

+

)

c

o

s

K

(

;

+

)

_

(

+

)

c

o

s

s

-

"

(

c

o

s

s

+

"

W

I

D

W

c

o

s

h

¥

(

-

)

+

l

2

c

o

s

h

¥

(

-

)

+

(

-

)

s

i

M

¥

(

-

}

l

c

o

s

K

,

K

(

-

)

-

}

5

1

(

)

-

(

-

)

'

N加川)=-号{尋[-B増H1(siM筈十帯cosh等)

×

s

i

n

(

"

-

)

c

C

s

"

l

-

B

[

l

3

I

l

+

I

s

i

M

(

7

-

×

c

o

s

h

2

(

,

_

S

i

m

l

l

s

i

n

(

_

)

c

o

s

+

(

(

S

i

m

c

o

s

h

)

'

2

(

_

)

c

o

s

_

(

_

)

s

i

n

-

,

,

r

l

C

l

s

i

M

C

l

(

s

i

n

h

c

o

s

h

)

l

s

i

n

K

s

i

n

K

'

_

(

;

+

)

l

s

i

n

+

(

-

)

-

(

-

)

+

-

÷

(

-

)

3

+

(

-

)

'

2

)

=

-

,

{

{

B

c

o

s

h

'

0

-

"

)

c

C

s

"

{

(

,

-

)

c

o

s

-

,

s

i

M

2

(

c

o

s

-

,

)

+

(

_

)

c

o

s

}

+

K

o

l

(

)

c

o

s

_

-

"

s

-

"

(

"

γ

2

+

c

o

s

z

-

"

'

c

o

s

'

2

1

=

-

,

{

I

-

B

I

W

s

i

M

(

1

-

'

'

c

C

s

(

(

,

+

"

)

c

o

s

-

3

}

s

i

M

-

(

,

-

)

c

o

s

,

-

K

(

2

-

c

C

s

K

¥

{

-

(

2

-

l

c

o

s

-

]

}

(16)∼(21)式をゼロにするために次式のようにFourier級数展開した式を用いる。

"

(

0

,

,

1

=

+

(

)

c

o

s

(

)

=

(

,

j

s

i

n

=

+

(

)

c

o

s

伽(M/21=-筈尋停伽)sin等”

(18) ( 19) (20) (21) (22−a) (22−b) (22−c) (22−.)

(11)

皆川・前畑:点支持を有する偏平シェルの線形解析 111

M山川)=-,号(州。+零M緬卿cCs等麺)(22-eI

WIM'21=-,号い"+靭卿cos等璽)(22-f)

(22)式が恒等的に成立するために次式を満足しなければならない。

拠 嗣 0 = 0 ( 2 3 − a ) 拠 加 " = 0 ( 7 2 = 2 , 4 , 6 … ) ( 2 3 − b )

ひ " " = 0 ( 、 = 2 , 4 , 6 … ) ( 2 3 − C ) g N , 、 o = 0 ( 2 3 − . )

gNn"=0(、=2,4,6…)(23−e)エッNm"=0(、=2,4,6…)(23−f)

シ M m o = 0 ( 2 3 − 9 ) g M 耐 " = 0 ( 〃 = 2 , 4 , 6 … ) ( 2 3 − h )

y V m o = 0 ( 2 3 − i ) 型 V 耐 " = 0 ( 7 Z = 2 , 4 , 6 … ) ( 2 3 − j )

さらに(22−6),及び(22−.)式のSine展開が閉区間で一様収束するために,端部で次式を満足させなけれ

ばならない。

"

(

/

2

)

=

[

[

-

B

(

1

-

A

c

o

s

h

-

(

1

-

s

i

M

B

,

[

(

,

4

"

)

(

c

o

s

h

s

i

M

)

+

(

,

-

[

×

(

2

s

i

M

-

3

c

o

s

h

+

3

s

i

M

)

_

'

1

,

+

"

c

o

s

h

(

,

-

)

s

i

M

l

l

l

-

C

W

(

,

+

s

i

M

(

,

+

,

×

-

I

h

I

s

i

M

+

K

§

(

+

K

(

2

-

3

-

2

+

1

1

-

2

l

-

"

.

K

+

2

1

-

K

+

K

,

(

2

+

"

1

-

−K‘号冊,器-長十赤)l=O

6

'

2

1

=

-

[

[

;

(

"

(

s

i

M

c

o

s

h

)

B

'

{

;

0

+

"

(

s

i

M

c

o

s

h

)

+

(

)

'

3

c

o

s

h

+

,

)

s

i

M

]

-

K

÷

+

K

÷

+

K

(

,

_

=

また,面内力Wが端部で大きく乱れるのを防ぐために次式を追加した。

N

b

/

=

[

I

-

B

W

+

(

1

-

s

i

M

B

(

,

_

)

-

(

)

S

i

m

-

1

,

-

}

_

(

C

F

c

o

s

h

C

,

s

i

M

)

K

_

(

:

+

+

K

(

+淵-人(た+器+器)'一零"耐祭器'Lcos等

-

"

"

W

l

D

W

c

o

s

h

竿

(

2

C

O

S

A

竿

s

i

M

)

×

c

o

s

K

K

'

3

-

+

(

,

-

)

]

=

(24) 卿 (26)

(12)

112 鹿 児 島 大 学 工 学 部 研 究 報 告 第 2 7 号 ( 1 9 8 5 )

m=2,4,6,……,2N,刀=2,4,6,……2Nとすれば,(23)式より6N+4本の条件式が得られ,さらに(24),

(25),(26)を加えて6N+7本の条件式となる。それら条件式は未定係数Bf,…BfN,B;,…B;N,Cf,…C;N,C;,…,

C;N,Df,…DfN,D;…D;N,K,,K2,K3,K4,K5,K6,K7(6N+7個)と未定係数”"n(NN+N個)の関係を与え,

C;N,Ⅸ,…,DfN,D;,…D;N,K,,…K71 次式のようにmajγな表示される。

{

S

}

=

I

F

l

l

d

l

[

S

2

l

l

d

=

[

F

J

l

d

1

ここに,|dllT=│z()20,”22,…,”2N,2Nl ld21T=lBf,…BfN,B;,…B;N,Cf,…,CfN,C;,…,C;N,Ⅸ,…,DfN,D;,・・ ld31T=│B:,…,BfN,B;,…,B;N,K5,K61 (27),剛式の'0121,|d31をldllで表わすと次式となる。

d

2

1

=

[

S

r

1

F

1

]

l

d

l

l

d

=

{

(4)式から求まるEuler方程式に法線方向変位”,及び‘を代入して次式を得る。

L

(

"

,

)

=

,

+

_

"

[

I

-

B

W

"

'

r

c

o

s

h

B

,

(

4

c

o

s

h

,

+

s

i

M

,

)

+

"

(

"

"

)

c

C

s

K

C

l

(

"

-

2

4

'

r

2

4

)

c

C

s

-

2

4

γ

1

+

"

c

o

s

"

(

'

+

"

)

c

o

s

"

W

l

c

o

s

l

+

I

H

I

(

A

竿

,

s

i

n

h

+2cosh等,│等(近-号)sin等叶;入(H等,Sim等…等エ

,

[

[

4

(

2

+

,

s

i

M

,

+

(

,

_

c

o

s

h

l

c

o

s

'

2

(

5

+

)

c

o

s

h

H

X

1

1

-

A

s

i

n

h

2

U

-

,

(

_

)

s

i

n

,

_

l

2

c

o

s

h

(

,

s

i

M

,

(

-

)

c

o

s

}

_

,

s

i

M

,

"

"

1

C

,

-

l

c

o

s

h

a

l

2

c

o

s

h

+

(

,

_

,

}

×

c

o

s

-

K

_

c

o

s

(

,

c

o

s

K

{

(

H

l

c

o

s

(

,

_

s

×

]

+

"

s

"

(

c

o

s

_

!

c

s

n

{

"

¥

(

-

)

+

D

,

(

,

-

(

-

)

s

i

n

¥

(

(27) (28) (29) (30)

(13)

皆川・前畑:点支持を有する偏平シェルの線形解析 113

-

2

s

¥

(

-

)

l

1

c

o

s

-

K

,

6

1

(

-

)

-

1

2

-

1

一M+M'''-2緋-号)鰯+淵一脈

K

,

-

A

'

(

-

)

-

(

-

)

+帯│卵cosh等州等,siM等叶K等十K‘筈

一血楓小s等顧−1)+零恥(cos等麺−1)cos判-蒜剛

エに,鰹=急

(4)式にGalerkin法を適用する。Galerkin法は境界条件が変位で与えられる場合に,この条件を満足するような

変位モードを仮定して適用するのが通例である。ここでは変位仮定において5種類の変位モードを用いているので,

Galerkin法を用いると,次式が得られる。

;

"

L

I

"

,

,

)

(

c

o

s

_

,

)

c

o

s

d

(32)

器念〃鯉L(",‘,(cos等jH)cosh等’伽㈱

器六〃'驚L(",,)(cos半工_!)半,Sim等‘畑帥

器 〃 奥 L I " 岬 筈 ( C 、 s 半 r _ ! ) 伽 開

器 〃 2 肌 ’ 1 淵 c o s 等 璽 − 1 ) 伽 ‘ ‘ ’

ここに,771=2,4,6…,2N,両=2,4,6,…,2N(32)∼(30)式を次式のようにmamJc表示する。

,¥I器淵│差│+隠監11差│=│芸’帥

ここで,(29),(30)式の関係を用いて(37)式を未定係数がldll,Idlだけの式で表わせ,次式を得る。

帯 I M l l ‘ , │ + { K } │ d 1 l = │ P } ‘ ,

]

[

,

]

[

M

l

2

1

S

r

1

F

1

]

[

S

1

2

]

[

,

+[S5lF2]T[M22ISTlF,]

K

]

[

,

,

]

[

]

[

,

]

[

]

T

[

K

2

,

]

[

S

]

T

[

K

2

2

]

[

lPl=│P,}+[S可lF2]│P21

(38)式に基づいて静的,及び固有振動数解析を行なう。ここで用いた解法は,自然境界条件をすべて満足するよう

な解法なので[K]は実対称であり,固有振動数は標準的な固有値問題の解法を利用して求められる。

4.点支持を有する偏平シェルの解析

点支持を有する偏平シェルの解析は,これまで束縛条件として境界条件だけであったのに対し,新たに束縛条件

が課せられた変分問題となる。すなわち,点支持した点の変位拘束条件を付帯条件とした変分問題である。

ここで扱う解析モデルは図4に示すように。c=0,αは従来通り面内変位を拘束した固定支持辺であり,

"=±6/2においては任意の点に固定点支持を有する偏平シェルである。支持点では変位u,ひ,”,及びZ/方向

(14)

114 鹿 児 島 大 学 工 学 部 研 究 報 告 第 2 7 号 ( 1 9 8 5 ) 支持点 X X y 1 点 支 持 3 点 支 持 5 点 支 持y y 図4点支持を有するシ土ノレ X

これらの付帯条件とLagrange乗数の積を累加することにより,

(4)式を導くのに用いた汎関数汀に,これらの付帯条件とLagral

関数万*を次式のように定義する。

*

'

r

W

(

j

r

/

2

)

:

(

J

c

i

±

/

2

)

W

'

u

(

±

/

2

±

6

'

2

ここに,ノW,ノW,入?,入::Lagrange乗数

jc‘:支持点のJr座標 (洲式の第1変分をとると次式となる。

’施薙=-〃(D両'"一志N鍾一古W-P+岬"仙

十/(N伽十N柳-M等十V伽)吻

U

N

N

,

±

b

,

±

6

/

2

l

l

-

U

N

N

b

/

+

;

±

b

/

U

V

W

(

"

,

±

b

,

6

/

2

)

{

'

l

M

6

'

,

6

'

2

)

+ヴルu(工i,±6/2)+6ノV・ひ(工‘,±6/2)+6,蝦.u)(jr‘,±6/2)

b

'

MOI式より,点支持を有する辺における応力境界条件は次式となる。

MZ/(工,±6/2)=0ただしMZ/い,±6/2)=一ノV

W(Jc,±6/2)=OただしVz/い,±6/2)=−ノリU NZ/(工,±6/2)=0ただしNZ/(Jci,±6/2)=−ノW

NJcy(gr,±6/2)=0ただしNJcZ/(Jc‘,±6/2)=一入:

また,変位境界条件として次式を得る。 uい,±b/2)=0,ひ伽,±6/2)=0

±

6

'

2

1

=

0

,

±

b

'

2

1

=

(41)式の境界条件を満足させるために,点支持したところだけでイ

Fourier級数展開すると次式を得る。

点支持したところだけで値を表現できるデルタ関数を用いる。

新たな汎 (洲 (40) (41) 伽 こ れ を

(15)

皆川・前畑:点支持を有する偏平シェルの線形解析 115

肌cos等鍾十ルーノvいcos等叶‘w)

肌cos等叶卿Vい=-鳩(零淵cos等叶叶)

M3)

(

)

c

o

s

,

N

"

=

-

;

(

c

C

s

;

(

)

s

i

n

=

-

、 α

s

i

n

(43)式を!恒等的に成立させるために次式を満足させなければならない。 9M派=一入s・脇(m=2,4,6,…) シM・=−入;・碗 9V耐=−ノW'・淵(、=2,4,6,…) 9V。=−ノW・がW ZソN、"=一入‘・脇(、=2,4,6,…) 44 9N。=−入り・6; 囚麺"N碗"=−ノW・脇伽=2,4,6,…) 伽,(28)式を用いれば,(44)式は次のようにmatγな表示される。

I S 化 │ = [ F ! l l d l + { D , l l l l M ,

[Sル│=[F211d!│+[DJ│'’‘‘’

ここに,I2lT=│ノW,入;,…,ANP,ノWj,…,聯,ノW…,服p,ノW,…,入#P│,NP:点支持数

また㈱式を77zamJc表示して次式を得る。

[ S ‘ l l d 班 │ + [ F ‘ l l d 1 } = 0 ㈹

㈱式を(47)式へ代入してMをldllのみで表すと次式となる。

│'│=_IMrⅦ!ⅡM「1F!+F‘lldl㈹

さらに,㈱,㈱式へ(48)式を代入する事により,’d2',’d3'は'。,'のみで表わされ,これを次式のように書き表わ

す。

│ d H R 1 l l d l ㈹

│ d ‘ │ = [ R ル } 剛 ’

ここに,[R,]=[S「1F,]−[SrlDlIS3・ST1D,]'[S3SrlFl+F3] [R2]=[SrlF2]−[S厨lD21S3SrlD,]![S3S「'F1+F3]

以降の計算は,自由辺を有する偏平シェルの解析において,関係式(29),(30)の換わりにMリ,(50)式を用いるだけで良

く,(3,)式以降同様に演算を行ない,静的,及び固有振動数解析を行なう。

参照

関連したドキュメント

 仮定2.癌の進行が信頼を持ってモニターできる

ベクトル計算と解析幾何 移動,移動の加法 移動と実数との乗法 ベクトル空間の概念 平面における基底と座標系

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

Series of numerical analysis to estimate structural frequency and modal damping were conducted for a two-dof model using the simulated external forces induced by impulse force and

A., Miller, J., 1981 : Dynamically consistent nonlinear dynamos driven by convection in a rotating spherical shell.. the structure of the convection and the magnetic field without

地盤の破壊の進行性を無視することによる解析結果の誤差は、すべり面の総回転角度が大きいほ

平均的な消費者像の概念について、 欧州裁判所 ( EuGH ) は、 「平均的に情報を得た、 注意力と理解力を有する平均的な消費者 ( durchschnittlich informierter,

管の穴(bore)として不可欠な部分を形成しないもの(例えば、壁の管を単に固定し又は支持す