995
591
INCLUSIONS IN ELASTO-PLASTIC SOLIDS UNDER THE INFLUENCE OF INTERNAL AND EXTERNAL PRESSURES
TAPAN KUMAR DAS
IndianInstituteofMechanicsof Continua 201 ManicktolaMainRoad, SuiteNo.42
Calcutta- 700054,India P.R. SENGUPTA
Department
ofMathematics University of Kalyani Kalyani,West Bengal,
Indiaand LOKENATHDEBNATH
Department
of Mathematics Universityof CentralFlorida Orlando,FL
32816,U.
S. A.(Received Febrt,arv 3, 1992 :rd it revised form April 10, 1993)
ABSTRACT.
An
attemptismadetostudytheproblemsofsphericaland circular inclusions in elasto- plasticsolidsunderthe action of internalincreasingpressure
and externalconstantpressure,
takinginto considerationofthework-hardening effect. Particular attention isgiventothe linear workhardening effectonbothproblems.
Itisshownthatresultsofthisanalysisare ingood agreementwiththose of idealplasticsolids.KEY WORDS AND PHRASE.
Sphericaland circular inclusions,elasto-plasticsolids.1991
AMS
subject classification codes: 73E.1.
IN’IODUCTION.
Severalauthorsincluding
Mott
and Nabarro ],Eshelby[2
5],Jaswon
andBhargava[6]
Willis [78]
have made systematic investigations oftheproblems of sphericalandcylindric’,dinclusions in an elastic medium. BothBhargava[9]
andSengupta [10]
haveusedtheprincipleof minimumenergy
tostudy problems
of symmetric anisotropicinclusions in anelastic medium.Tresca
11 has shownthat a metal yields plastically whenthe maximum shearstress attains acritical value.Bhargava [12]
has investigatedtheproblemsof inclusions inelasto-plasticsolidsunder theassumptions ofthe infinitesimal theory ofstrainandperfectplasticity conditionssatisfiedbyanelastic-plasticsolid. Fromaphysical point of view, thework-hardeningeffectisimportantin aplasticmaterialeventhoughtheconsequences
of this effect arequite complex. Hopkins 13]has included thework-hardeningeffectof the material in592 ’1. K. I)AS, P. R. SI..NGIlP’I,\ ,\NI 1.. I)ICBNA[’lt
his
problem
of dynamical expansion of sphericalcavities in metals.Sengupta
and his associates[14 15]
have studiedtheproblems
ofinclusions inelastic-plasticsolidsof work-hardeningmaterialoffinite and infiniteextent.On
the otherhand, Tokuoka[16]
has investigated theplastic
deformations and instability ofspherical shells under internalpressure. In
spiteof thisprogress,furtherinvestigationof the inclusionproblemsin anelastic-plastic mediumisneeded.Thispaperisconcernedwiththeproblemsofsphericaland circular inclusionsinelasto-plastic solids under the action of internal increasing
pressure
and external constantpressure,
taking into considerationof thework-hardening effect. Special attention is givento the linearwork-hardeningeffect
onbothproblems.
2.
SPHERICAL INCLUSION IN ELASTO-PLASTIC SOLIDS.
We
consider a spherical region ofradius a in a finiteelasto-plastic material which tends to undergoadimensionalchangetoasphereof radius a(1 +i)
where lies inthe region of plasticstrain.We
designate the spherical region asinclusion,and theoutermaterialasmatrix,where its external boundary isasphericalsurface of radius b. Thesphereisunderconstantexternalpressure,
because of theconstraintsofmatrix, stressesappear
both in the inclusion and the matrix. Hill 17]considered the problemofsphericalshellunderuniformpressure
on itscavity surfaceundertheassumptionsof finite plasticstrains.We
assume that the internalboundary ofthe matrix(a
sphere of radiusa)
isunderagradual
increasingpressure p. At
firstatmoderatepressure p,
thematrixbehaves likeanelastic material.Introducing spherial polarco-ordinates(r, 0,
)
and thecorresponding displacement components (u,v, w), for radial symmetry of deformation wesuppose
u=u(r), v=0 and w=0
(2.1)
In
absence ofbodyforce, theonly non-vanishing equation of equilibriumsatisfiedbythedisplacement component uis+-- =o. (2.
The solutionofthe above differentialequationis u
Ar + B
r2, (2.3)
where
A
andB
areconstants. The stress components are,therefore,o"
(3;1. + 2)
4p-y, B (2.4ab)
fro=o, =(3) + 2/)+ 2/ 7. B
AccordingtoHenckyandyonMises,the
yielding
commenceswhenthemaximumofIo’ o’,1
reachesacritical value
v,
wherev
isa material constant.Now o’ rr, 6t B
ismaximum when r a. Therefore,theyielding commencesatr a and thecorresponding pressure
P0isdeterminedby
the followingconditionsI ]
[ar],=a
=--P0’[r]__ =-p,, 6P7 7 (2.5abc)
where P0
>>
P,.Thus,
yielding
begins whenIN;IL:51t)NS IN t’.I.,\ST()-II,\SI IC b)ll!)S 593
Po 7’
+
p(2.6)
andthecorrespondingdisplacementatthe innerboundaryisgivenby
2( a_ /
aU (a)
- ?’a
3,+ 2,u
-7b+ 3& +
2m(2.7)
Withincreasing
pressure
aplastic region spreads intotheshell. The plasticboundary
beaspherical
surface, its radiusatany
moment is denotedby
c.Thestressesanddisplacementin the elasticregion
(b 3-1 / (a3
r-/
o’,
=-A +B --1,
c
< ?’ <
bare stillofthe fomacr
o=O’=A
+1-B
+1 (2.8abc)b
-
u=At 3r+2/ r+4// -B r+3z+2// 4//
intheregion c<r
<
b.We
mustconsidertheplasticsolid in theregion a<r<c0o’,
2+ 7(err tro)
0subjecttotheyieldscondition
tro-
o’,?’.
Solvingthe differentialequationandusingthe condition ofcontinuityof thenormalstress atthe elasto- plasticinterface r c, we have
or, _2?’logC r- 2( c3)
1--p,
o=,=?’_2rlogC_
r--?’
2(c33)
1--p, (2.10ab)
If
o-
-p, at r a, thenp
2?’logc+
1-+ p[
a
(2.11)
From
the condition intheplastic regione + 2eoo -E2 v (O’ + 20"o)
where
v
andE
arethe Poisson’s ratio andYoung’s
modulusrespectively,and the conditionof continuity of displacementatthe interface, the radialdisplacementuintheplastic region a<
r<
cis_rr 27"c .41.tc3_+3_Kb )
c27’(1 2v) [
log-c-lc3 1]
--3-+
ru=3K
b, 12,u -3"+ -3"+
r E r3 -P" (2.12)
Thedisplacementatinnerboundaryis
I( ) ] [c-1c31]
a27’c3 41.tc3+3Kb
lc27’(1-2V)
log--ff-+ a-p,
u(a)=- Lt i -fi ’-’’+
aE
a7
We
know thattheradialdisplacement ofthe inclusionisa(8- e),
anda and therefore pl
3K’(- e) u’(a) p
3K-
(2.13)
Now by usingtheconditionofcontinuityof normalstressatthe inclusion andmatrix,we get
3K’(8 e) 2?’logC +
1-+ (2.14)
Thedisplacementof the innerboundaryof themau’xisawhichby equation (2.13)is
ae c-
3K
]2kt b75-3 -5
+ a E log +33
-5- (2.15)Therefore
= (1-
YC3 v)
(2.16)The relation
tween
eand isgivenby-+--- P (2.17)
a 3 3
3K"
The material in the region a g r gc is under elasto-plastic strainand satisifes the work-hardening conditions. Ifwedefine
e;, e;
as the radial and tangentialelastic sains,e, e
asthe radial andtangential plastic strainsand
e,eo
asthe coespondingquantitiesof total sain and if u the radial displacement, then we have thefollowingrelationse e + e,
ee; + e (2.18)
Ou
ue , e =-r (2.19)
E e
2Vo,Ee (1 V) V (2.20)
It
iscustoma
toassume thattheplasticsain satisfies theincompressibilityconditione + 2e
0.(2.21)
Then the following compressibility can derivedbyusingthe equations(2.18) (2.21)
+ 2 0r + (2.22)
where
K
isthebulk modulus.The elasto-plasticmaterial safifiesthe following
equilibrium
equationO
2a, =7 (a a’) (2.23)
The
sess-sn
ce forawork-hdening materialinuni-axialcompressionisof the fo(sHill [J])r + H(e) (2.24)
where
,
earecompressivesess
andstrain(bothten
aspositive),T
the initialyieldsss, H
ise
degree of hardening expressed as a function of total strain. Evidently, in radially symmeic deformations,
y
elementof material issubject
to a uni-axi, radialcompressive sss
state- together
withahy&ostatictensile stressa.
Th lattersess, by
itself,produces
a posiveisotopic,elastic sain ofamount
{(1- 2v)a}
/E.
Rememberingthesignconventionfor8
de,e
appropfiageneralyieldcriterionis
o-, y+H{-e,+ (1- 2v)}E (2.25)
If
ere
isnoBauschinger effect, thenH(e)=-H(-e),i.e.,H
(e) isan odd functionofstrn.
Thuse gener
yieldcriterion fora work-hdeningmatefi isfro ff + H -e + E
1(2.26)
In
caseoflinework-hardening,H
is a linear functionoftotalsain andan analytical discussionis possible.In
such acase therateof work-hdeningis constantand wesupposethe yield cri:eNen as+ffo
(2.27)where
H’(e)= E,
gradient of thestress-strain curveinthe plasticrange, supposedconstant. Solving (2.22)and(2.27),wefindcr
o’
(1- E,
3(2.28ab)
Eliminatingthestresses
o’
and0" between theequations(2.23)
and(2.28ab)we obtain thefollowing ordinarydifferentialequationfor u.d
(d--r + 27u ) 6y(1--/1 (2.29)
dr
(3K + E,)
rThe solutionof the differentialequation(2.29)is
U
A,r ++r ?’
3K+ E, r(31ogr- 1) (2.30)
where
A
andB
arearbitraryconstants. Therefore, with thehelpofboundaryconditions(2.5abc)
and the continuity condition, the stresscomponents and displacementcan be obtained from(2.30)
and(2.28ab)
inthe plastic region a_<r<
cand theyare2y(1---)log(_)3 4?’ E,(1-v)+4’ E,(1-v)
O’= 3(1+3__)
3E(I+-)
3E(I+E__)()
-Pl,1-v
](C3r3 c3a c3 Pl
+)
log2, +J 5 )+
2-_v
3KE(b-c)-E’(Ec+ 3Kb) +(1-@)log
r+E(I+)
"=-( E(ll + )
3KEb::[3bc(4-3K)(1-5){4E, r+9Ka(1-5)} ],
and
(2.31)
where590 I. Ix. I)AS, P. R. SENG[;I)[.\ AND l,. I)EBN,\’It]
q
=9Ka1--- 12Kb3(E,-31a)+121uc’(E, +3K)]-4E,[3Kb(_K+411)]. (2.32)
Therefore, thedisplacementattheinternalboundary ofthematrixis
L 3 +
1+logTa
)’(1E(1 -v)(1- + _) E--)c3
aP" q [ab3c3(4 3K)(1- ){4E,
a+ 9Ka(I (2.33)
It
isnotedthatthisdisplacementinduces the initial elasticpart of thedisplacements.We
considertheinclusion which isunder a normalpressure p
onitsexternal
boundary andis therefore< -p,, a; -p,, r;o o.
Itisnotverydifficulttomark with thehelp oftheyieldcriteria of
Tresca
11 that the materialof the inclusion neveryieldsandit isalwaysinastateof elastic deformation. Thedisplacementatthesurface of theinclusion isgivenbyu’(a)= ap( (2.34)
3K’
By
usingthe conditionof continuity ofnormalstress atthe inclusionand matrixboundary, we getE
(_) 4?’ E,(I-v)c
3K’(d;-e)= 21’ 1----
l+31og+
E,
aE,
3E(1
+_3__)
3 3K
-2’cA+-
b(2.35)
Thedisplacementof theinnerboundaryofthe matrix will be ae,whichby equation
(2.33)
is=-(--3 E
1--2V[3KE(b3-c3)-E,(Ec3+3Kb3) (+ 1--/logC3]a
t 1+- )
P[g663c3(4-3K)(1-){4E, a+9Ka(1-)}]
E 1+
(2.36)
Solving ford; with the aid ofequations
(2.35)
and(2.36),
weobtain-v) c3[1+E,(
4(2.37)
Ifthebulk modulus of thematrixand the inclusion are same, thent y(1- v)c’ p, bc3(4u 3K’)(E,- 9K’)
E
a3K’ q (2.38)
where
INCI.USIt)NS IN t.I.ASIt)-I’I.\SI 1(; SOI,IDS _597
q
12a3(9K ’- E,)[K’b’(E,- 3)+ lIc’(E, + 3K’)]- 12E, K’b(3. + 2//)
This result is essentially the same ifthe work-hardening effect of the plastic material is altogether absent.
Now for theelasto-plasticsolidsof work-hardeningmaterial,the relation between e and
t
is gven by27
1-E,
() --(l+31og.) 41,’ E,(1-v) (.):+2Y(c_)
e=- I+--3KE’ 9KE(I+)
9K,b)q
Forapicularly plasticmaterial,
e
corresponding resultisgivenbye log-
+ (2.40)
3Kk
a 3 3Comping theabove results itfollowsthat the relation betweene and for thecase of sameinclusion andthe
maNx
materialdepends(i)onPoisson’s ratio
v, Young’s
modulusE,
bulkmodulusK’
of the inclusionandontheyieldsss
in case of perfectly plastic solid, and (ii) ontherate ofconstant work-hdeningfactorE,
sidesqutitiesalreadymentionedin(i).
Theequilibriumpressureof the inclusioninthe present case isgivenby
(2.31),
whichmaycomp
touation
(2.11)for theperfectly plasticcase.Another impoant pointinthepresentcase isthejumpinthe
hoop sess
one surface
inclusiongivenby31-@ -+ Yc + 3p, (2.41)
+
aThisresult is independent of ande butdependsontheconstant teofwork-hdeningfactor
In e
case cfperfectly
plasticsolids thejumpin thehoopsess
is%theyieldstss,whichcan ducedom
theaboveresult by puttingE,
0.All
oer
impoant results involved in(2.31), (2.33)
and(2.39)
are found toagree
withcoesponding
resultsforpefftly
plastic solids given byBhgava 12],
if weput b, p,
0 dE,
0,inthe aboveuations.
3.
CIRCULAR INCLUSION UNDER CONDIONS OF PLA SAIN
Bhargava 12]discussedtheproblemofcircular inclusion inanelasto-plasticmaterial of infinite extent.
We
consdierhere circular cylindrical inclusion in a finiteelasto-plasticmaterialwitheffects of work-hardening. The cylinderisunderconstantexternalpressure p,. As
theproblem
ofcircular cylindrical inclusioneven under theconditionsoftheplane
strain is muchmoredifficultthan that of sphericalinclusion, wesuppose v for compressiblematerial. Thisassumption greatly simplifies the2 solutionoftheproblem.
Introducingthecylindricalcoordinates
(r,O,z)
and the corresponding displacements (u,v, w), weassume for the presentproblemu=u(r),v w
O"
Proceeding in asimilarmanneras the case of sphericalinclusion, wesupposethattheinternalcircular boundaryof the matrixis underapressurep within the elasticlimit. Ourdiscussionisconfinedtothe scopeof the infinitetheoryof strain.
The
only
equationofequilibriumunder nobodyforce is+
=0(3.1
dr k,dr and hence the solutionis
u=Cr+-- D
(3.2)t"
whereCandDarctwoarbitraryconstants.
Theyieldingcommences when the maximum of
Io’o o’rl
attains acritical value ?, where],isa materialconstant.Now o’o o, 4// D
ismaximumwhen ," a. Thereforethe yielding beginsatr a and the correspondingpressurePo
is determinedby thefollowingconditions[(Tr]r=a
=--Po,[(Tr]r=o=--Pe
and
7 4//-y D
(3.4)where P0
>> Pe
Therefore yieldings commenceswhen
p0
- + p (3.5)
andatthisstage thedisplacement of theinternalboundaryis
Uo(a)=[b (A. +//)
a+ -] 2( a2 +//)
pe(3.6)
For
anincreasingpressure beyond 1--- + p,
theplasticzone isdeveloped
inthe matrix,andif c be theradiusof elasto-plastic interface, then proceeding as inthe case of spherical inclusion under externalpressure,
theelastic stressesand displacementinthematrixbeyond
c aregivenby
(r,,
c - p,,
a-c 7’ +
p, ov(o- +
o-o)
v2p, (3.7abe)
and
"=
4
b( +//)
r+ 7 2(. +//) p (3.8)
The material in the region a_<r_< c iselasto-plastic and satisfies the work-hardeningcondition.
Presentingananalysissimilartothatofasphericalinclusion, theequilibrium equationis
--- 0o’r +
r(o’ ’)
0(3.9)
and the yield criterion of linear work-hardeningmaterial is
Theplasticmaterial satisfiesthe compressibility equation
Ou
uOr
r(3.10)
(3.11)
Sincethe displacementcomponent u is continuousonthe elasto-plasticinterface r=c,thenit isgivenby
4
b"(& +p)c+
r2(; + p)
(3.12)Eliminatinguando" fromequations (3.9),
(3.10)
and(3.12),weobtainrr -r-
1-+ E, -c b2()+p)
c+
c jrE, 2(A. +p),
"3(3.13)
Solvingthe above differentialequationfor the normal stressandusingthecontinuitycondition onthe elasto-plasticinterface r=c we have
o’, 7’log-
+
r
b2(2 +
2
=
(3.14ab)
2
7’
1- 1-2log+cZ E,
-C b2(+,f/)c+v + V- -P" 4(;/.t)+1
given by
-Pl o -Pl z
-2VPN1 (3.17ab)
u’ (1 + v’)(1- 2v’)
-p (3.18)
r
E’
where
E"
andv"
aretheYoung’smodulus and Poisson’sratiorespectively.According to Bhargava [9] the inclusion never yieldsfor infinitesimal strains. The circular inclusion of radius of radius a spontaneously undergoes dimensional change to a circle of radius
a(1 + )
inthe absenceof matrix and attains the radiusa(1 + e),
when it is inequilibriuminthepresence
of thematrix. Presentingananalysissimilartothatof asphericalinclusion, we obtain]1 Pe
:=4 kb2(;I,+//)c+-
a2(+//)
(1 + V’)(1-
2v’) E"
(3.19) (3.20)
where Plisgivenby (3.16).Theseresults
(3.15)
(3.20)for elasto-plastic solidswithwork-hardeningeffects are ingoodagreement withthe correspondingresults of aperfectly plastic solid if we put 1), 0 andE,
0 in the aboveresults.
The
pressure
Plattheinternalboundary
ofthe matrix is" @_) 7’ a2[ff_ { l}(1 1)(1 cl__/]
p,:il- 7’log- b2(X+/---c’+-- + F 5---- E,c + -
As
regdsthe inclusion, if itisunder unifopressurep,
thepfincipNstressanddisplacementfield isThe jumpinthehoopstressaswecross thecommonboundaryof the inclusion and the matrix is
60
o0
1-+
9( + flc + p, (3.21)
If
E, O,
thisresult reducestothatofperfectly
elasticsolid.It
isimportanttopointoutthatin case ofperfectlyplasticsolidthehoop suessisindependent of the external boundary of thematrix.Finally, all the results obtained in the above analysis depend on the constant rate of work- hardening factor
E,
whichplaysanimportant roleinplasticity.REFERENCES
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
Mott,
N.
F.,and Nabarro,F.
R. N.,Proc. Roy. Soc. A
52(1940)
p. 25.Eshelby,
J.
D.,Phil,Trans. Roy
Soc,A. 244 (1951)
p. 87.Eshelby, J. D.,
Prof. Roy. Soc.
(London).A.
241,376(1957) p. 376.EshelbyJ. D.,
Proc. Roy. Soc.
(London) A.252
(1959b) p.561.Eshelby, J. D.,
Progress
insolidMech:ni;Vol.II (Ed"Sneddon and Hill)(1961).
Jawson,M. A.,andBhargava, R. D., PrQc.
Carnie.
Phil..Soc.
57 (1961) p. 669.Willis,J. R., Jo0r,Mech.Ph.y,Solids
V91. 1
(1965) p. 377.Willis,
J.
R.,Ouart. Jour.
Mech. Appld. Math..Vol.17(1964) p.
156.Bhargava, R.D., Appl.
Sci.Res A,
Vol.II, (1962) p.
80.Sengupta, P. R.,
Ind.J.
Theo. Phys. Vol.XI, No. (1963) pp.
27-37.Tresca,H., Compl,Rend.59
(1864)
p. 754.Bhargava, R. D.,Ind.
Jour.
M0.1;h, Vol. 5,No.2(1963)Hopkins,H. G.,
Progress
insolidMechanics,Vol. (Ed:Sneddon and Hill)(1961) p. 83.Sengupta,
P. R.,Ind.J.
Mech. Math.Part II,
Specialissue,(1969)
pp. 80 89.Sengupta,
P.R.,
and Kundu,H. N., InO, J,
The0,Phys.Vol. 18,No.
4(1970) pp.
109 121.Toku ka,Tatsuo,Indiana Univ.,l,20 (1970/71) pp. 193 206.
Hill,