1. No category

Muroran-IT Academic Resources Archive
Title
Author(s)
Citation
Issue Date
URL
部分分布荷重を受ける有限円柱の3次元応力解析について
松岡, 健一; 能町, 純雄
室蘭工業大学研究報告.理工編 Vol.9 No.1, pp.169-179, 1976
1976-12-18
http://hdl.handle.net/10258/3642
Rights
Type
Journal Article
See also Muroran-IT Academic Resources Archive Copyright Policy
Muroran Institute of Technology
部分分布荷重を受ける有限円柱の 3次元応力解析について
松岡健一・能町純雄*
O日
aThree DimensionalStressAnalysisofa Finite Cylinder
SubjectedtoPartiallyDistributedLoads
KenichiMatsuokaandSumioNomachi*
A
b
s
t
r
a
c
t
巴msc
o
n
c
e
r
n
i
n
gaf
i
n
i
t
ec
y
l
i
n
d
e
r旦r
es
o
l
v
巴dt
r
e
a
t
i
n
g
Thet
h
r
e
ed
i
m
e
n
s
i
o
n
a
ln
o
n
a
x
i
s
y
m
m
e
t
r
i
c
a
ls
t
r
e
s
sp
r
o
b
l
ff
i
n
i
t
eF
o
u
r
i
e
r
H
a
n
k
e
lt
r
a
n
s
f
o
r
m
s
t
h
ee
q
u
a
t
i
o
n
so
fe
q
u
i
l
i
b
r
i
u
mo
ff
o
r
c
e
si
nc
y
l
i
n
d
r
i
c
a
lc
o
o
r
d
i
n
a
t
e
s，bymeanso
巴c
y
l
i
n
d
e
ri
ssubmit
t
e
dt
oap
a
i
ro
ft
r
a
c
t
i
o
n
s
Then
u
m
e
r
i
c
a
lc
a
l
c
u
l
a
t
i
o
nwasp
e
r
f
o
r
m
e
df
o
rt
h
ec
a
s
ewhent
h
w
i
t
hr
e
c
t
a
n
g
u
l
a
rd
i
s
t
r
i
b
u
t
i
o
nonb
o
t
he
n
d
so
fad
i
a
m
e
t
er
.Ther
e
s
u
l
t
sa
r
ecomparedw
i
t
ht
h
o
s
巴 i
nc
a
s
eo
f
i
n
f
i
n
i
t
l
l
yl
o
n
gc
y
l
i
d
e
rc
a
r
r
i
n
gt
h
esamel
o
a
d
sandt
h
er
e
s
u
l
t
sc
o
r
r
e
s
p
o
n
d
i
n
gt
oas
i
m
i
l
a
rl
o
a
d
i
n
gc
o
n
d
i
t
i
o
no
ftwo
ea
l
s
oshowi
nt
h
ef
i
g
u
r
e
sf
o
rr
e
f
e
r
e
n
c
es
a
k
巴.
d
i
m
e
n
s
i
o
n
a
ls
t
r
e
s
s丘r
1.まえカずき
円柱座標で表わされる次元応力問題の一種として円柱の問題がある。円柱の軸対称問題はI.
N. Sneddonが軸対称側圧が作用する場合の無限円柱の解析を Fourier変換を用いて行なって
いる 1)。また有限円柱の軸対称問題としてはL.N.G.Filonの研究がありへその後今井，斉
藤が行なっているが，特に斉藤は側圧および端面に部分的圧力が作用する場合の解析を行なっ
ている 3)。非軸対称問題は現在でもあまリ多く研究されてはいないが，無限円柱に対しては G
Mathewが側面に任意の外力が作用する場合の理論を発表しているが数値計算は行なっていな
，。この数値計算は吉野・宇津木によって部分的な側圧をうける問題として行なわている
.2)
。
4
)
さらに非軸対称変形をする有限円柱の解は秦によって導かれ 5)，これを発展させた奥村の研究
があリ数値計算も行なわれている 6)。
一方能町は，有限 Fourier-Hankel交換による円柱座標の 3次芯応力問題の解 j
去を示し 7) 著者
等はこれを有限円筒の非軸対称問題にまで拡張し，若干の解析例も示した 8)。ここではこれを有
限円柱に適用し，部分分布荷重を受ける場合の解析を行なった。数値計算は荷重の位置を 2種
変え，無限円柱とした場合の結果および非軸対称 2次元応力問題としてえられる結果と比較検
討した。
*北海道大学工学部教授
工博
(
1
6
9
)
松
1
7
0
健
岡
田I 和
l雄
一・能
2，有隈円柱の F
o
u
r
i
e
r
-Hankel変換による解
円柱座標で表わされる
3次元応力問題の有限 Fo
w
i
e
r_
H
a
n
k
e
l変換による解法にすでに発表し
ている 8)のでここではその結果のみを示す。
r
(u)
1のように円柱の半径を a
， 長さ C， 半径方向， 円周方
図
(v)
e
向，軸方向の座標を r
，，
zとし，各方向の変位をそれぞ、れ U，
V，
W
とすれば，
u 二 五L
_(Aorz+Borz)十 12(Amrz+Bm)cosM(1)
ムπ
π
u=12(Amγz-B
m
r
z
)s
i
nm
e
冗
w
こ
図 -1
m
(2)
m
τLWoTz+12W
L π π
m
明間
c
o
sm
e
(3 )
ここで
A mrz
十
二
1_ ， ，
(1
1
~， m +1
3
/
1
ー
ト A A
~fmp(r)l古川十 2(2μ 十入)向。 +-J-Z石 Amo
m-1_
1
l
:
_
士
ム
口mo一 :
1 一一
一
、
l
一
一
一
一
つ
1
2μ+λ 【
2(
2
μ +A)
¥W
mcwmo
'j
I
1J
r
1
1
3 m+1 A
十 汗 ∞sNZl_
NXmp(N川三五 αmn十27β mnサ - 7 A
つ
Tr
/
>T
¥
(
， /1+A n
I
1 m-1 ~
2
μ 十人
)
μ+A 1ω ρ (
I H _ ¥f 1 n
，mn1
1
十 ← 一一.J.. Bmn- ~ ~ I~ H Cm
Nr)
i
βmn-，
C
~1
nr
2A vmnj
←
2A NWmp¥nf'
l2 μ j
J
{
五
乍 DmhhALiA)(3γ mik+omik)}
十三平 {Jれ己決;会Jφ(k)(~iZ )
出い(k)(~iZ) - ljr(k)(~iZ)}{ ルー去( 3
γm
i
k十
九
品
)
}
]
(4)
α(r¥m-1I
1
¥
，
/1+
A 1
， 2J
n
μ 十μ十人 1
(~J
I C)/，)~ 1
1¥(α
m
o
，
:土
一一
A
m
o+
μ
μ
B
m
r
z一←
C ¥α)
L
2(3
μ 十人)¥
um
o一向。)+
jJm
o， 3
μ
十
A a"
<
m
u
'ー 3
μ +A
n
A
I
μ十
よ
￨
+
十一一一一一
。l
-E;"ω)一 μ(μ3
品
。
1- O
，
"
'
0
2
)
一 一 (Wm
c-W
)
I
o
3
μ 十 A (E品
μ 十入)(
¥
V
m
O
l
v
m
o
z
，I 一2一
(
3
μ一十人)¥
m
c
u
.
rm
m
o
!
J
(1
!
m
s
(r
)
1
αm
l
2
μ o
1
2
(
2
μ+A)
n
，
β問。+
m+
1 /1↓~
， m-1 3
υ↓ 1
一ー:: ，
/
1
.
， A問。十一一一一色一'
+
B
α 2
μ +A"<mu' a 2μ十人
A
1
1
J
十 人 /
、
2
(
2
μ 十人 ¥WmC-Wmoノ
f
u/
(
17
0
)
1
7
1
部分分布荷重を受ける有限円柱の 3次元応力解析について
I1
r
υ
ー
f
-A
1 m+1
+~
~Txms(Nr)~
一一一
α mn2
十丘一一九
n~~~A
c~cosNzl N
A
.m
s
¥
l
2μ
μA n
2 a
I >T
~\....vV .L't' Nl
.J.'f'
n
¥
!
A
"~mn
fJTTt
3 m-1B
2
a L
J
m
n
ー ←
一
当1Cmn} 日
主 川){Lsmn一Cmn}]
A
2~
ωm
r
]m
l
(
t
i
r) ~ ~(剖
r
A μ 十λ )
)
}
2
g
zE
I
[
φ(
k
)(
t
;i
Z)
l2
J
1
'
平
工 Emik μ(2μ 十入 f(γmik+3Omik)f
+ ;2手{jm+l(印
品士(
φ
((tiZ)}{Emik-L
(
γ
m
川い)}]
吊
+
m
(
f
Y
1
(ん(Z)(wmo-2ゐC2)九
(c-z)(wm
1
)}
c-2ら 0
(5)
Wmrz= hl(Z)
g
m
(r
)Wmc+h(c-z)Wmo
I1 ~ '
， ，
(1 ~
m+1
m-1~
2J1ー
ト .
¥~
1
+-2mlb17VG
山)~
n 一α
一 Amn十 一
Bmnー 」
ト CmnJ
n
'--，~ L
N -"
" ~ .， '
l入 βm
)
+
ヰ
生 F山
けは
τ
A
y
βm
n一 Cmn}]
記
r
{
Q
(
k
)
(μ)+
点土よ ρ
(
k
)
(
t
i
Z
)~(Dmik -Em
i
k)
{
]
"
山(
t
i
a
)}
2t
i!
:
l
l
l
"
" U，
-'N/' 2 μ + λ )
+~~~m( か)
~
( μ
t
Q
(
k
)
(
t
i
Z
+ λ L
，，¥
I
(6)
ここで， μ，λは L
a
m
e
'の弾性定数であり ，N = nπ
/c
(n=
m 次の第 l種 B
e
s
s
e
l関数である。 t
i(
i二
土
1，2，………)， ]m(;
t
i
r
) l
1，2，・・・…・・・)は ]m(;
t
i
a
)= 0の根を小さいものか
)
頃に並べたものである。また式中の関数は，
ら1
1
m
l(Nr)
;
;
:
N
G
)
'
X(Nr) τ
(N
a)
_ Im(Nr)
I >.T..¥ _
1 ( N r)
Gm(Nr) 五 両 ，
xmp(Nr)
二
I>T
¥
二
問
Fm(Nr)= Nrxms(Nr)-Naxms(Na)Gm(Nr)
ωmp(Nr)= NrGm(Nr)-Naxms(Na)xmp(Nr)
ωms(Nr)= NrGm(Nr)-Naxmp(Na)xms(Nr)
か (
r
)= (: )m，
fmp(r)二 京 北 百
r
玄品川:
(7)
1
(
1
7
1)
イ(:
+
r
1
削
(
r
)
松岡健一・能町純雄
1
7
2
Q(1)
(
f
i
Z
)
lc
oshfi(C-Z)土 coshfiZ
Q
(
2
)
(
f必)J
c
o
s
hf
ic:
t1
φ(1)(EiZ)l_ sinhfi(C-Z)平 s
i
n
hfiZ
一 一 一 一 一 一 一
f
i
Z
)J
c
o
s
hf
ic:
t1
φ(2)
ψ
(
2
)
(fiZ)1 fizcoshf
i
(C-Z)
干f
i
(
C-Z)
co
s
hfiZ
cοshfiC:
t1
ψ
(川 fiZ)J
ρ(1)(fiZ)1 主z
s
i
n
江主包二三注 fi(C-z)si凶 fiZ
c
o
s
hf
iC士 1
ρ(2)(fiZ)J
hl(Z)=
三 h2(Z)=三2c¥ 立6
これらの関数の主なものの導関数は
l
t
Gm(Nr)±Nxms(NY)一 町 m
(N r)
= NXm
p
(Nr)十竺 Gm(Ny)，
r
l l uN
..¥
a
ar Xmρ(Nr)=NGm(Nr)一 一m一+
r)，
l
f
-Ixms(Nr)=NGm(Nr)+
一一~ X
ms(Nr)，
rー か ゆjJ(
¥...."， a
r
l
I
l
. fIL
jrFm(Nr)=Nhs(Nr)一 咋 川 ) 二 Nω mp(Nr)十 町 川r
，
)
f
m十 1
jLω mp(Nr)= NFm(Nr) 一 一 ωm
p
(Nr)+2NG，
川 r
，
)
r
m+1
j L山 ( ル )= NFm(Nr)十 一 一 山 (Nr)十 問 問 (Nr)，
r
j
E
φ
(臼 ) ニ ー ω ( μ )， θ，_θz 似
Q(μω)二一イg
州
ω
必
φ
糾州(
v_
a
θ
否Zψ (fiZ)= f
i
{Q(fiZ)-P(fiZ)
，
} 万三~ P(fiZ)= f
i
{
φ (fiZ)一 ψ(fiZ)，
}
式中の係数 α間 n
，βm
n，Amn …等は各境界の変位および応力で与えられる積分定数であり，
ooo
これらは境界条件から決定される。
従って各応力成分は， σr
，σe
，σzをそれぞ、れ r，e
，Z 方向の直応力， τT(}，T
f
)Z ，T
z
rを せ
ん断応力とすれば， H
ookeの法則を用いて
r
σ二
去[
2
μ
(
j
h
o+aar Borz)
十 A.
1o
]
ロ
十士事 cosmθ
[
2
μ
(ご
す Amrz+ å~
1，.
"，.
μ
2
1ー
(Ao +Bo
0
1I
~
σ。
=τ一
~7T
I
r
ロ
，~
Bmrz)十A.
1m
]
1，
l
μ
1 ~ __
_ n I~ ( m+]
c
o
smBI2 ヱ午ナ土 Amrz
問)十人 .1 1 十土~
I
(7)
7
T m
(
17
2
)
L
¥
r
1
7
3
部分分布荷重を受ける有限円柱の 3次元応力解析について
ペ
ユ
l
B
m
r
z
)+ぬ
(8)
μ
一
.
_
1I
~
a
，. 1 ， 1'<0
n I
~
a
1
白ご一一
12μ←~Worz
，L:1
l+ --=
'
:
_~ c
osm
B
W
十人
L
lm
ol
r
z
2πL
-f-< az .V'~ 十，，
，
v
π
....，......~ .
.
.
.
.
'
"I
[2
'
J
;
v
.
- a
z，
. m
ff
tJ，
i
G
1¥..... lI
l
m
TTT
士号 T日inmθ
ai
寸
(Amrz-Bmrz)-~Wmrz]=
114
)
+
(
1
m
μ山
士
号 μsinme[ å~
「三五
(9)
唱E 4
z
1
)ハリ
日
" .
J
(
柑
=
士
号
r
TTT
Ia
a
1 1
rコ
r
z+Bれ ) 十 万 W
μ
l
万 (Ao
o
r
z1+才 μ
∞smBI å~ (Amrz十Bmrz)
，~
1.
"
，
TTT
~
n
d
LWmrzl=l_~ Tmrzc。
…
。
r
I
1[
m+1 ¥.
/θ
(
12
)
m
m-1¥ 月
，( a
(
13
)
I~一十一一一一 lAm 間十 I~一一←一一一 lBm 間十て子-W m ロ
¥ 0γ γ /
¥ ο γ γ
/
OZ
3，境界条件
各境界に任意の外力が作用する場合の一般的な境界条件は
i) z=0または
z=cにおいて
2'
" (:m\7~'~\121
]m
(C
ir
) r
1 T rt r~. 1 ，'" T rt r 1
n1
z q
l
(
r，
e
)= て二2~
~ JmCo[ql]+~JmCm[ql]cosmBI
{
]
m
+
l
(
f
i
a
)
}
2l2
m
"
'
O
L':ll
J ~tI m
v
m
l
l
f
l
JLJ(14)
r
二
t
t
I
…
I
ω dB
JmC [ql]=fL231(?， 川 刈 け c
ただし
加
VCI fftU
τZ=q3(r
，
θ
)
τr
z= q
2
(r，
e)， θ
(
1
5
)
式 (
1
5
)は
Ju
r
Ahb
r
+
、
m
IJ
lry
Au
r
QU
m
n
qd
σ
Aγ勺
q
n
u
+
Q
J
V
m
FL
O
QU
M
γ円
ptt'fo
q
，
。
πil
ptt'FO
a
EEEJ
q
一
一
、
q
u
“
ヮ
G
A
+
m
+
J
Jm山
とすれば
y
1'
" ]
m
+
l
(
t
i
r
)
T
m
r
z十 T
;
"
r
z=
-;
;ヲ万
戸
高
司
Jm+
l
[
Q
Z
+
q
3
]
T
m
r
z
-Tんz
(
1
6
)
手 川 1J11q2Cosh3sinmovd
肋
(
17
3
)
1
7
4
松岡健一・能町純雄
1 ~ ]m-l(tir)
2lJIq-q3l
:
; Um l(tiaW
(
1
7
)
十
また式の誘導過程からの条件として
(
1
8
)
Dmik+Emik二 γm
i
k十 O
m
i
k
を満足しなければならない。
i
i
)
r
二
αにおいて
1
ニ
一
一
σγ=のθ
(，
z
)
~ ~
r
"
1~
~(
o[
CoCO
[
q
4
]
+_
_
J
，
_
_~ c
osmθ{
で CmC
q
4
]
π
m
1C
ム1[C;
十与~
cosNzCmCn[q4H
L m
(
1
9
)
J
]
=
f
f
V
4
(
θ，z)cosmθ
cosNzd
紘
CmC
[
の
n
，
1..，-，. ~(1~ ~r 1
2~
>r~ _ r ，
i
τγθ ニ q5(B ， z) = ー~ ~sinmB~
~ SmC
[
q
5
]+C
~ n
~COSNZSmCn[q5] (
O
m
lC
..~-
SmCn[
Q
5
]
二
'.'-~J
(
2
0
)
f
γQ5(Bz
，)
s
i
n
m
B
c
o
sNzd
紘
，
(l~~rl ，
τ「仇 (θ ， z) 二汗 sinNzl2~ C
O
S
n
[
Q
6
]十
l ..，-，
1
∞smB~_~rlì
CmSn[Q6]r
(
2
1
)
またここでも式の誘導過程から
A
m
r
z
)
日 二
1Amo+22Amncoslvz
し
L
(
2
2
)
n
Bmrz) γ二G ご 1Bmo 十三~
BmncosNz
L n
(
2
3
)
し
を満足しなければならない。
以上の条件より Ql~Q6 が与えられれば積分定数を
全て決定することができる。
4，数値解析
有限円柱の数値解析例として図
2に示すような部
分分布荷重を側面に受ける場合を考える。このときは
1
4
)， (
1
6
)， (
17
)， (19)~(21) において
式(
q
l= Qz二 Q
3二 Q
5二 Q
6二 O
図
(
1
7
4
)
2
1
7
5
部分分布荷重を受ける有限円柱の 3次元応力解析について
であるから，式(16
)， (
17
)から
γmik- Omik =0
(
2
4
)
(25)
. Dmik士 一 Emik
また式 (
2
0
)から
αmn- 0
(
2
6
)
式から(19)から
βmn=ffq4(r，
θ
)∞sm
(
Jc
osNzm
倣
キスι ( 2 7 )
= 4{1一(一 1)
n
}
{1-(-1)吋」れi
n
m
B
o
c
o
sNClsin
Hμv
残った積分定数 ，Amn，Bmn， Cmn，Dmik等は式 (
1
4
)，(
21
)(
2
3
)から決定する。
0.25，c
a=O.2
数値計算で、得られた結果を凶 -3-11に示す。計算は全てポアソン比 ν=
o/
の場合で，級数項は幻二 6
4， 1ニ 32， m ニ
4
0としたものである。
3，4は軸対称荷重で cl/a=0.9すなわち円柱の中央部に荷重が作用する場合の σ円 の
図
2はめの分布であるが図中に同
の z方向の分布を rの 各 断 面 ご と に 示 し た も の で あ る 。 図
仁条件で無限円柱とした場合の値を r=0.8aと r= 0のものについて破線で示した。この計
算は既に斉藤によって計算されていてはほ同様の結果がえられているが，著者等の値は多少大
きくなっている。これは級数項のとり方による差と思われる
無限円柱との比較では r/aの
r/a=1.0
1
.0
r/a=1.0
0.8
0.8
h
〆 u
0
L
ぺ
QJ
甘
lLnopL
001
8
戸
nuqLnu
，/一守
zにd=
a
・
一
一リν
ヴ¥ h b
O / 一司/
0.6
•
.•
0.4
o.
2
nu
。
ー0
.2
1
.0
0.8
0.6
0.4
0.2
。
。
0.0
1
.0
0.8
0.4
0.6
一一一一一- z/a
z
/a
図 -3
図
(
17
5
)
4
o.
2
0.0
1
7
6
松岡健一・能町純雄
大きなところではほとんど差はなしまた荷形に近いところでは r/aが小さく伝ってもほぼ
等しい分布をしている。
1
.0
0.8
0.4
LIJ
0.2
¥)=0.25
co/a=0.2
日
。 =π/2
ミ0.6 ~W …
。
』
cl/a=O.l
σo.つ
、
。
い
ハハ
0.4 ト
1
' -0.2
o.
2 .
，
O
.4 ~
¥)=0.25
c，
/a=0.2
由
。 =π/2
cl/a=O.l
-0.2
-0.8
0.0
0.2
0.4
0.6
z/a
0.8
1
.0
0.0
0.2
0.4
0.6
z/a
0.8
1
.0
0.2
0.0
図 -6
図 -5
r/a=1.0
1
.0
0.8
¥)=0.25
co/a=0.2
。
0=0.2
『
、
0.6
。
0.6
cl/a=0.9
、¥
日=
0
b
亡
ア
N
与
司
I 0.4
0.4
‘
。2
0.2
0.2
0.0
0.0
-0.2
-0.2
1
.0
0.8
0.6
0.4
0.2
0.0
1
.0
0.8
0.6
0.4
z/a
z/a
図
7
図 -8
(
17
6
)
1
7
7
部分分布荷重を受ける有限円柱の 3次元応力解析について
さくなるに i
也、急激に減少し ， rニ 0
.
8aでは r=αの分布と比較し，ほぽ全体的に符号が逆転
している。 r=Oでは荷重とは逆符号の割合大きな応力を生じている。無限円柱の結果はあほ
ど ~i文せずー当然のことながら z/a 二 O ではかな 1) 差がある。
i
勾
5，6は同様に軸対称荷重が作用した場合であるが
部分荷重が作用するときの結果ーである。めは
ず，特に
fニ
0
.
8aでは逆に荷重強度よ
Zニ
Cl/α
=0.1，すなわち円柱の両端部に
Oでは rが 小 さ く な っ て も そ れ ほ ど 減 少 せ
1
)大きくなっている。 Zが大きくなると図 -3の相当位
置の値と近くなっている。
同
4の分布とは全体としてかなけ異なっている。 z/α=0
.
2
6の の の 分 布 は や は り 図
で応力は不連続となるが，この場合は z/a>0.2において荷重とは逆符号に大きな応力を生じ
ており ， z/
α<0.2では境界の条件を満足させるように応力が減少している。ここでは無限円
村との比較は行なっていない。
図-7，
8は非軸対称荷重 (80=0.2) が C，
/α=0.9すなわち円柱中央部に作用した場合の
のの分布を示した。図一 7はめの分布であるが合力の大きさが小きいため特に r方的jへの
h
減少が大きく中心部では括以下になっている。図
7には同じ条件の無限円柱の備も示してい
、
るがほとんど一致しているため図にはその差が表われていな L。
の は r=αではかなり大きい値であるが
f方向にも
z方向にも急激に減少しているがその減
少の程度は軸対称の場合よりかなり著しい。この場合も無限円柱の値を破線で示したが，ほと
r/a=1
.0
1
.0
0.4
0.8
v=0.2S
0.2
co/a=0.2
σ0.6
、
。
日
。
相
=0.2
cl/a=O.l
0.4
0'
、
N
b
0.0
日=
0
0.4
-0.2
0.2
co/a=O.2
0.4
。
。
白
。
=0.2
cl/a=O.l
日=
0
-0.6
-0.2
0.0
o.
2
o.
6
0.4
0.8
1
.0
一一一一一一~ z
/
a
図
。
。
0.2
0.4
0.6
一一一一~ z
/a
図 1
0
9
(
17
7
)
0.8
1
.0
1
7
8
松岡健一・能町純雄
んど一致しており rニ G でもそれ程大きな差で
1.0
はない。
9，1
0は同じ非軸対称荷重で c，
/α=0
図
0.8
lすなわち円柱端部に作用する場合である。め
の 分 布 は z=oで は そ れ 程 急 激 に は 減 少 し な
℃
ず
i
o口 6
0
いc 従って図 -7とはかなり異なり大きな値であ
るが ， zの大きなところでは図
7と一致して
0.4
くる。
0.2
のはやはリかなり特徴的な分布をしている o
yニ G では
z=oからはほ、直線的に増加し z/
。
。
α=0.2で、のキ 0.42q と な リ z/α>0
.
2で は
z/a 0
.
2の の キ
二
~一一一一一 r/ a
0
.
5
8qか ら か な り 急 に 減
より小きな rで は の は 急 に 減 少
少する。 y=a
し，減少の程度は軸対称の場合よリ大きい
1.0
図-11には，以上で計算したいろいろの場合
からめと σ
。について θニ Oにおける荷重直下の
f
直，c
，
/aニ0
.
1に対しては
Zニ
co/a==O.2
え0.8
0， c，
/aニ0
。
.
(
>
9に対しては
Z二
G の値の 7
方向の変化を示し
た。悶には無限円柱の値も示しているがめ，ぬ
0.4
一十
，(
j
。
ニ π
/2， z/α=1
.0の線
ともの /α=0.9
に添っている破線と c
，/α=0.9，80ニ0.2，z/
α=1.0 の 線 に 添 っ て い る 破 線 が そ れ ぞ れ 軸
対称および非軸対称荷重に対する値である。ま
た ， そ れ ぞ れ の 図 に 2次 元 非 軸 対 称 問 題 の
1.0
0.8
4
(
]
.6
日 4
0.2
0日
一一一一_r/a
図 -11
値 9)を合せて一点鎖線で示した。
めの分布ではすでに示したように，軸対称，非軸対称とも c
，
/aニ0
.
9の場合には無限円 f
iの
債とほとんど一致しており
2次元問題としてゐ =0.2で計算した値も 3次元問題では軸刈科、
とした場合とほぼ等しくなっている。
c
，/a =0.1 の場合には軸対称，非軸対称とも 0.8<
r/α<0.9で最大値をとり荷重強度より大きくなる。特に軸対称では r/a与 0.85で γ
σ キ1.1
1
q位とかな 1
)大きく r方向にゆるやかに減少する。
ぬの分布は図から明らかに θ
。=0.2の場合は r
が小さくなるに従い急激に減少している。ま
たc
，
/aニ0
.
1の方が c
，/α=0.9の場合より小さくなっているのはあの分布とは逆である。
無 限 円 柱 と し て の 値 は や は り 破 線 で 示 し た が 。'
0=
=π/2で 多 少 差 が あ る が 8
0.2ではほ
0 =
(
1
78
)
部分分布荷重を受ける有限円柱の 3次元応力解析について
1
7
9
とんど一致している。また 2次 元 問 題 と し て の 解 は ぬ の 場 合 は 80 =0.2の場合に近いが全体
としてはあまり一致しない。
数値計算は先に記したように z
=32，刀 工 6
4，m=40として級数和を求めた結果であるが精
度を検討するため種々の項数の和をとり精度を検討した。結果は記さないが，上の項数で卜分
であることは確認している。
5
.む す び
Fourier-Hauhel変換による有限円柱の解を用いて若干の数値計算を行ない，無限円柱として
の解， 2次元非軸対称問題における解と比較も行なった。計算例は多くないが， c/σ=2.0の
場合でも長きが載荷巾の 5倍もあれが無限円柱としてもそれ程の差を生じない(中央載荷の場
合) 2次元問題としての解は単純には有限円柱には適用出来ないが， σ
。についてはある程度適
0
用出来ることがわかった。
以上の計算は北海道大学大型電子計算機センターの FACOM-230/75で行なったものであ
(昭和5
1年 5月2
2日受理)
る
。
参考文献
1)Sn
巴d
don1
.N
.:F
o
u
r
i
e
rT
r
a
n
s
f
o
r
m
s，1
s
t，E
d
.p
p
.504-510，1
9
5
1
， 1
9
6
7，
2) 宮 本 博 3次元弾性論，第 4章
3)斉藤秀雄・機械学会論文集，第 1
8巻 6
8号 pp.21-28
，1
9
5
2，
4)吉野利男，宇津木諭:第 2
3回応用力学連合講演会講演論文抄録集. pp.303-304， 1
9
7
3，
5)秦 謹一:北海道大学工学部研究報告，第 1
3号
， p
p
.13-44， 1
9
5
5，
6)奥村 勇:第 2
4回応用力学連合講演会講演論文抄録集， p
p
.183-184， 1
9
5
3，
7) NomachiS
.G
.:t
h
eMem.MuroranI
n
s
.T
e
c
h
.，Vo
.
l3
，
N
o
.
3，pp.477-50
l
， 1
9
6
0
.G
.andS
.G
.Nomachi:THEORETICALANDDPPL
lEDMECHANICS，Vo
.
1
2
2，s
p
.199-209，
8)
M
a
t
s
u
o
k
aK
1
9
7
4
(
1
7
9
)