Hijink, H.

Veenstra, P.C.

English

Pure OAI Repository

 Analysis of both kinematically and statically admissiblevelocity fields in plane strain compressionCitation for published version (APA):Veenstra, P. C., & Hijink, H. (1978). Analysis of both kinematically and statically admissible velocity fields inplane strain compression. (TH Eindhoven. Afd. Werktuigbouwkunde, Laboratorium voor mechanischetechnologie en werkplaatstechniek : WT rapporten; Vol. WT0439). Technische Hogeschool Eindhoven.Document status and date:Published: 01/01/1978Document Version:Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers)Please check the document version of this publication:• A submitted manuscript is the version of the article upon submission and before peer-review. There can beimportant differences between the submitted version and the official published version of record. Peopleinterested in the research are advised to contact the author for the final version of the publication, or visit theDOI to the publisher's website.• The final author version and the galley proof are versions of the publication after peer review.• The final published version features the final layout of the paper including the volume, issue and pagenumbers.Link to publicationGeneral rightsCopyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright ownersand it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.            • Users may download and print one copy of any publication from the public portal for the purpose of private study or research.            • You may not further distribute the material or use it for any profit-making activity or commercial gain            • You may freely distribute the URL identifying the publication in the public portal.If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, pleasefollow below link for the End User Agreement:www.tue.nl/taverneTake down policyIf you believe that this document breaches copyright please contact us at:openaccess@tue.nlproviding details and we will investigate your claim.Download date: 17. Nov. 20231 WT 0439 ~ANALYSIS OF BOTH KINEMATICALLY AND STATlCALLY ADt1tSSlBLE VELOCITY FIELDS IN PLANE STRAIN COMPRESSION. P.C. Veenstra H. Hij ink October 1978 1. Introduction PT-Rapport Nr. 0439 The co-ordinate system used to describe a situation of plane strain compression (upsetting) of a specimen between parallel flat dies is def i ned In fig. 1. 1 • Fig. 1.1. The upper die moves with the velocity - t 00 with respect to the plane of symmetry z = 0, while the lower die has the velocity + tOO' In the frictionless case the velocity field is given by 0 =~O x h 0 0y = 0 (1.1) 0 = - ~ 0 z h 0 - 2 -The effect of friction, which presents itself by bulging of the specimen, can be accounted for by introducing ( 1.2) The function must be symmetrical with respect to z = 0 and assume there a maximum value. Moreover the function must be anti symmetrical with respect to x = O. The present study aims to trace those velocity fields which both stabilize or minimize the system with respect to power and satisfy equilibrium conditions in case of non strain-hardening material. 2. The calculus In the analysis reduced (dimensionless) quantities are used lit X lit Z bllt b x =11 z =1;" ; =1) lit (J •• • lit h o .. = ....!..l E: •• = €: •• 00 IJ (J I J I J O~ O· I =-I 0 0 For the sake of simplicity the asterisks will be omitted. Introduce the velocity field o = x {PO + f (x,z,P.)} x I where P. are free parameters. I The relation must satisfy the conditions 0 (x,o,P.) = max x I Q (x,z,P.) = U (x,-z,P.) x I x ' I U (x,z,P j ) = -u (-x,z,P.) x x I (2.1) (2.2) (2.3) - 3 -The continuity of the material flow requires ;.. x ::: t 0 dz o x from which follows that t Po = 1 - 2 f f (x,z,P.) dz o I It is denoted t J = f f (x,z,P.) dz o I and hence Po = 1 - 2J From eqs. 2.2 and 2.7 it is derived Because of the situation of plane strain it is found -0 ::: J E dz = (1-2J)z + f f(x,z,P.)dz + x f ~f dz z xx I oX which must satisfy the boundary conditions o ::: 0 Z for z = 0 0=- t z for z = t Next it can be calculated and in plane strain ~ = 2. (e 2 + E: 2)· t 13 . xx xz (2.4) (2.5) (2.6) (2.7) (2.8) (2.10) (2. 11) (2. 12) - 4 -from which according to Levy-von Mises is derived £ xz £ (2.13) Through this the frictional stress in the contact plane is known as E (x,LP.) xz I • (2.14) ~ (x,·k,P.) I From eqs. 2.2. and 2.7 the relative sliding velocity in the contact plane is found to be o = x {1 - 2J + f (x,t,P.)} xo I Since physically it must hold that o ~ 0 for x ~ 0 xo the kinematical constraint is f (O,!,P.) - 2J ~ - 1 I The reduced press force is defined as F~ = __ F __ b.w.a (2. 15) (2. 16) (2. 17) where w is the width of the specimen and hence the product b.w is the surface contact area. Because of equil ibrium it must hold for any value of z where (J 1 b/2 {oAVE}z = b/2 l dx On the other hand the power dissipated in the process equals (2. 18) (2. 19) (2.20) where Fd represents the contribution of deformation to the press force, whereas Ff does so with respect to friction. - 5 -Through this it can readily be shown to hold in dimensionless notat i on and 4 b/2 Fd = b f x=O 1/2 !. J£ dx dz z=O (2.21) (2.22) which through eqs. 2.12, 2.13 and 2.15 can be expressed in terms of the co-ordinates and the parameters P .• Thus it becomes possible to I investigate whether the sum Fd + Ff = F minimizes of stabilizes power as a function of the parameters. The equations for equilibrium in plane strain are dT xz ---ax (2.23) Because of symmetry it must hold that £ = 0 for z = 0 and hence T = O. xz xz From the plasticity condition it follows that 2 - -13 The boundary condition is for x = b/2 and thus b/2 {dT } . xz J ---az x z=O dx 2 = - - + 13 (2.24) and next {} 2 1 b/2 b/2 {aT xz} (} - - - + - f f -- dx dx zAVE z=O - 13 b/2 0 x az z=o (2.25) The second equation of equil ibrium renders {o } _1. = - j )'lTxz l dz + {oz}z=O Z Z-2 1) 1 ax )x (2.26) -.T') -A-s-t-e-r-j -5 k-s-a-re om i t ted. - 6 -Body equil ibrium requires that and hence 1 b/2 b/2 J o 1/2lClT t 6 d~Z~X dz dx = 0 (2.27) Now all conditions for a kinematically admissible velocity field which both minimizes power and satisfies equil ibrium are defined in the eqs. 2.21, 2.22, 2.29 and 2.31. It can be shown that all velocity fields which are developable in power series, like the parabolic field or the cosine field, the exponential field, etc., neither satisfy the condition of minimum power nor equil ibrium. The first field encountered which satisfies power conditions is the hyper elliptic field However from this class of velocity fields only the elliptic field proves to satisfy the conditions of equilibrium with respect to z = o. In order to satisfy body equilibrium the field has to be modified, as will be discussed later. 3. The ell iptic velocity field (3.1) From eq. 2.6 it follows J - t~1-(2PZ)2 dz - {lV + * arc sin p! 0.2) and thus through eqs. 2.7 and 2.1 - 7 -Ox = x [1 + B {~1-{2PZ)2 - 2J}] (3.3) £ = 1 + B {~1-{2Pz)2 - 2J} xx (3.4) -Oz == z [1 + B {t~(1-2Pz)2 - 2J}] + fr arc sin 2Pz (3.5) which satisfies the boundary conditions at z = 0 and at z == t. 30 z = 0 x (3.6) (3.7) T = -.1... Bp2 -=-_______ ::...:.xz=--________ --=o;-xz I'J [{ 1- (Zpz) Z}{1 +B (" 1- (Zpz) t - ZJ) r +{ ZPBZXZ}Z Jt 0.8) Thus the frictional stress is 1 2 x TO = - - BP -::-- -----------------:::;"";'" I'J [(1-PZ){I +f (0 -t arc sin P)}Z +{BPZX}Z]' 0.9) The relative sl iding velocity follows from eq. 3.3 ° = x [1 +! (Q - 1. arc sin P)] xo 2 P (3.10) Hence the kinematical constraint is 2 B ~ = BMAX 1. arc sin P -0 P (3.11) From eqs. 3.9 and 3.10 follows the contribution of friction to the reduced press force . fz} - 8 -The contribution of deformation is 1 b/2 -J bIZ x=O 112 [{ } 2 2]t J 1+B("t-(2Pz)2_2J 2 + {ZBP xz} 2 z=O 1-(2Pz) Numerical variational analysis of F = Ff + Fd shows that this functional minimizes with respect to the parameters Band P. A typical example is shown in fig. 3.1. Fig. 3.1 dx dz From this it is clear that the minimum is extremely flat, which implies that it is impossible to determine the value B t sufficiently op accurate. The same holds even more for P t' op However, when now considering the conditions for equilibrium, it follows from eq. 3.8 that 2 - - + B (l-2J) and hence according to eq. 2.24 (3. 14) Next it is found through eq. 2.25 {} Z {1 Bp2 'b 2} GAVE z=o = - Ij 1 +3' 1 + B (1-2J) (-2) (3.15) - 9 -The condition Fmin= {OZAVE}Z=O yields 3n13F. -1) B = min 2 P ( ~) - 3 (t 13 F min -1) (1 - 2J) (3. 15) This is visualized in fig. 3.1 where the curve according to eq. 3.15 intersects the curve corresponding to F in a well defined way. From the numerical analysis of F it appears that P ~ 1 and that variation of this parameter has only a minor influence on Fmin . Using estimated values for P as obtained from the numerical analysis and applying eq. 3.15 the table 3.1 results which is plotted in fig. 3.2. Fig. 3.2. Table 3. 1 b PEST FMIN Ff Fd B bB 1 0.98 1.332 0.084 1 .248 3.126 3.13 2 0.98 1.585 0.228 1.358 1.501 3.00· L. 0.99 2.135 0.687 1.448 0.744 2.98 6 0.99 2.702 1.206 1.496 0.496 2.98 8 0.99 3.227 1. 767 1 .510 0.373 2.98 10 0.99 3.848 2.303 1.545 0.298 2.98 12 0.99 4.423 2.766 1.657 0.249 2.99 16 0.99 5.586 3.695 1.891 0.187 2.99 It is concluded that 1. Though neither the quantity Ff nor Fd is a linear function of the compression ratio b, the minimum reduced press force F . min and hence minimum power in the system behaves virtually Ii nearl y. - 10 -This agrees with the slip line solution after Prandtl [1 J and the solution proposed by Unksow [2 J, though the latter refers to rotational symmetry. 2. As the compression ratio increases the effect of friction becomes increasingly dominant. 3. Most probably there exists a hyperbol ic relationship between the shape factor B of the velocity field and the compression ratio b B • b = 3 (3.16) Now, when introducing this relation in eq. 3.15 more precise values for the parameter P can be found as a function of b. This is not important for power, but it proves to be relevant for the stress distribution. Performing the calculations the results as 1 isted in table 3.2 are found. Table 3.2 b P FMIN B BMAX l-P b(l-P) 1 n 0 s o 1 u t i o n 2 0.9830 1.585 1.5000 1.5657 0.0170 0.0340 4 0.9907 2.135 0.7500 1.5247 0.0093 0.0372 6 0.9936 2.702 0.5000 1 .4770 0.0064 0.0384 8 0.9955 3.277 0.3750 1 .4306 0.0045 0.0360 10 0.9963 3.848 0.3000 1.4263 0.0037 0.0370 12 0.9968 4.423 0.2500 1 .4128 0.0032 0.0384 16 0.9975 5.586 0.1875 1.3487 0.0025 0.0400 0.0373 From this it is concluded: 1. The hyperbol ic relation 3.16 impl ies that the shape factor P in case of minimum power and equil ibrium is also controlled by a hyperbol ic relation 1 b (l-P) = 33 = 0.037 (3. 17) - 11 -2, In case that b < 1.8 no solution is found because the kinematical constraint eq. 3.11 is violated. The conditions for minimum power and equilibrium cannot simultaneously be satisfied. Probably a kind of mechanical instability is present. The results obtained up to now make it possible to formulate the ell iptic velocity field explicitely in terms of the compression ratio b. Consequently the stress distribution according to eqs. 2.13, 2.14, 2.24 and 2.26 resp. can be calculated for any state of compression. However, since T according to eq. 3.8 is a steady function of x, xz {cr} 1 turns out to be less in absolute value than the corresponding Z Z="2" value {cr} • For this reason {cr } ~ {cr } ,which violates z z=O zAVE z=! zAVE z=O the conditions for body equil ibrium. 4. The modified elliptic velocity field The problem that the elliptic field does not satisfy body equilibrium can be met by modifying the field in such a way that the shear stress is a non steady function of x. This is achieved if the frictional shear stress vanishes at the edge ( x = b/2 , Z = 1/2) as well as in the plane of symmetry ( x = 0 ). The modified elliptic velocity field which satisfies this requirement is defined by Ux • x [ 1 +B{(tY- ("x{ r {v'1 - (2Pz) ~ - 2J }] (4.1) where The quantity ~ is introduced in order to avoid instability of the com-putation at the very edge x = b/2 • Its physical meaning is that at this edge some frictional stress is assumed to be present, for instance of a Coulumb nature. Following the same procedure as in the previous sections it is found . -u z - 12 -2}{1 ;- 2 (2n + l)(lJX) 2' z\1-(2Pz) + 1 +1fP arc sin 2PZ - 2JZ} (4.3) au 1 X "2"-az = - 2j (b)2 2}n 2BP t '2 - (lJx) XZ { • Z ~ (21'Z) 2 + 41P arc sin 21'Z - 2Jz} From these results obtained and through eqs. 2.11,2.12 and 2.13 the quantities ~ and T can be calculated. £. xz For z = ~ the frictional .. o = - -...L BP2 J3 I r- 1 \1-pL p arc sin (4.4) (4.5) (4.6) .As a matter of fact .. = 0 for x = 0, whereas it assumes a minor value for x = ~ if lJ is suffi~iently close to one, as will be shown later in table 5.2. The relative sliding velocity in the contact plane follows from eq. 4.1. The kinematical constraint is 2 - 13 -= B max (4.8) Firstly it has been investigated whether the parameter n has a significant influence on the power dissipated in the system as well as on the optimum values of Band P. As shown in fig. 4.1. this proves to be not the case. If the parameter n is varied over a wide range as shown, minimum power varies about only 5%. However, the parameter is quite relevant for the ratio of the contribution of friction and the contribution of deformation to power. Physically this is the very reason that the parameter n takes care of body equilibrium. In order to find the n-values which satisfy the condition for body equi-IIbrium In the integral eq. 2.27 ~o z AVE 1 = -b/2 the roots for n are solved. (4.9) The procedure is fully numerical, i.e. first T is calculated from the xz relations as derived before, next the function is numerically derivated with respect to x and subsequently integrated with respect to z and x. It is assumed that the hyperbolic relations 3.16 and 3.17 may be applied, whereas ~ = 0.999 is introduced. The results are listed in table 4.1. TABLE 4.1 b n~104 3b2 ~(J *10 2 zAVE 2 12 12 -2.280 4 48 48 -0.420 6 112 108 -0.106 8 195 102 -0.012 10 300 300 +0.018 12 440 432 +0.022 14 598 588 +0.010 16 775 768 -0.014 18 990 972 -0.049 • - 14 -It appears that the calculated values as shown in the second column can be very closely approximated by the relationship -4 n = 3b 2 10 (4.10) In order to check its validity the relation is substituted into eq. 4.9 and the deviation from body equilibrium ~~AVE thus obtained is listed in the last column. It is clear that if equation 4.10 is applied the modified elliptic velocity field satisfies body equilibrium quite well. 5. Resu Its It is concluded that the modified elliptic velocity field satisfies re-quirements of minimum power, local equilibrium and body equilibrium if the relevant parameters are given by B :;: 3/b p = 1 - (5. 1) 33b n = 3b 2 • 10 -It Since T now exp 1 i cite I y can be ca I cu I ated as a funct i on of the co-ord i na tes xz for a given value of the compression ratio b, through eqs. 2.24 and 2.26 the stress distributions in the plane of symmetry as well as in the contact surface can be computed. The results are listed in table 5.1., whereas figs. 5.1,5.2 and 5.3 visualize some stress distributions. Finally in fig, 5.4 the stress at the point ix = 0, z = 1/2} is represented as a function of the compression ration. In most cases this value corresponds to the maximum tool load and for this reason ,n··is o·f technological importance. Conclusion Apart from the fact that the modified elliptic velocity field is of theoretical interest because it explains the stress peaks at the edge of the specimen as observed in experiments (3], the introduction of this particular field does not greatly affect the stress situation as derived from the elliptic velocity field, if the compression ratio ranges in 2 ~ b* ~ 10. - 15 -In th i s case and in order to calculate the maximum stress on the tool the eq. 3. 14 at the point x = 0 can be applied safely, as it overestimates s 1 i gh t 1 Y the maximum load in the contact plane. Hence BP2 2 a 2 1 (;) 2 s b s 10 Fl:l -- -zMAX V3 J3 1 + B(l - 2J) The relation is visualized by the dotted line in fig. 5.4. However, when the compression ration increases and thus according to eq. 4.10 the parameter n increases, the influence of modifying the elliptic field on the stress situation becomes more and more significant as is clear from fig. 5.4 with respect to the maximum stress on the tool surface. .. TABLE 5.2 b x 0 O. 1 O.~ 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 b/2 . {az}z=O 1. 79 1. 78 1. 76 1. 73 1.69 1.63 1.56 1.47 1. 38 1.27 1. t 5 6a 0.26 0.25 z 0.23 0.17 O. 14 O. 12 0.10 0.09 0.08 0.06 +7.23 2 {a z} z=t 1.53 1.53 1.53 1.56 1.55 1. 51 1.46 1.38 1. 30 1. 21 8.38 T 0.00 0.573 0.576 0.577 0.577 0 0.577 0.577 ' 0.577 0.577 0.577 0.004 2.61 2.60 2.56 2.48 2.38 2.25 2.08 1.89 1.67 1.42 1. 15 O. 18 0.16 0.12 0.10 0.08 0.07 0.06 0.05 0.04 0.03 +3.97 4 10 2.43 2.44 2.44 2.38 2.36 2.18 2.62 1.84 1.63 1.39 5.12 0.00 0.524 0.562 0.571 0.574 0.575 0.576 0.576 0.576 0.577 0.003 3.49 3.47 3.39 3.28 3. 11 2.90 2.64 2.33 1.98 1. 58 1. 15 0.13 O. 11 0.08 0.07 0.05 0.04 0.04 0.03 0.03 0.02 +2.75 6 1D 3.34 3.36 3.31 3.21 3.06 2.86 2.60 2.30 1.95 1.56 3.90 0.00 0.521 0.562 0.570 0.573 0.575 0.575 0.576 0.576 0.576 0.002 4.43 4.39 4.29 4. 13 3.89 3.59 3.23 2.80 2.30 1. 75 1. 15 0.12 0.09 0.06 0.05 0.04 0.03 . 0.03 0.02 0.02 0.02 +2.11 8 10 4.31 4.30 4.23 4.08 3.85 3.56 3.20 2.78 2.28 1.73 3.26 0.00 0.527 0.564 0.571 0.574 0.575 0.576 0.576 0.576 0.576 0.002 5.46 5.41 5.28 5.06 4.75 4.36 3.87 3.31 2.66 1.93 1. 15 0.10 0.07 0.05 0.04 0.03 0.03 0.02 0.02 0.02 0.01 +1.72 10 1D 5.36 5.34 5.23 5.02 4.72 4.33 3.85 3.29 2.64 1.92 2.87 0.00 0.536 0.566 0.572 0.574 0.575 0.576 0.576 0.577 0.577 0.002 6.63 6.57 6.40 6.12 5.72 5.22 4.60 3.88 3.05 2. 13 1. 15 0.09 0.06 0.05 0.04 0.03 0.02 0.02 0.02 0.01 0.01 +1.45 12 10 6.54 6.51 6.45 6.08 5.69 5.20 4.58 3.86 3.04 2.12 2.60 0.00 0.544 0.568 0.573 0.575 0.576 0.576 0.577 0.577 0.577 0.001 ... TABLE 5.2 CNTD v ~ "" b x' , 0 O. 1 0.2 0.3 0.4 0.5 0'.6, tw 0.8 0.9 1.0 -' , 0.7 b/2 .:,',',.'.' ' 8.00 7.93 7.71 7.36 6.86 6.22 5.45 4.54 3.50 2.35 1.15 0.08 0.06 0.04 0.03 0.02 0.02 0.02 0.01 0.01 0.006 +1.26 14 1D 7.92 7.87 7.67 7.33 6.84 6.20 5.43 4.53 3.49 2.34 2.41 0.00 0.551 0.570 0.574 0.576 0.576 0.576 0.577 0.577 0.577 . 0.00 1 9.64 9.55 9.29 8.84 8.22 7.42 6.45 5.32 4.04 2.62 1. 15 0.08 0.05 0.04 0.03 0.02 0.02 0.01 0.01 O. 01 0.004 +1.12 16 10 9.56 9.50 9.25 8.81 8.20 7.1,.0 6.44 5.31 4.03 2.62 2.27 0.00 0.557 0.572 0.575 0.576 0.576 0.577 0.577 0.577 0.577 0.001 11.67 11.56 11.22 10.66 9.88 8.88 7.68 6.27 4.69 2.94 1.15 0.08 0.05 0.03 0.02 0i-'02 0.01 0.01 0.01 0.006 0.003 +1.00 18 1D 11 .59 11. 51 11. 19 10.64 9.86 8.87 7.67 6.26 4.63 2.94 2.15 0.00 0.563 0.574 0.576 0.576 0.577 0.577 0.577 0.577 O.S77 0.001 14.21 14.07 13.64 12.94 11.96 10.71 9.20 7.46 S.49 3.34 1. 15 20 1D 0.08 o.OS 0.03 0.02 0.02 0.01 0.01 0.007 0.005 0.001 ' +0.90 14. 13 14.02 13.61 12.92 11.94 10.70 9.19 7.39 5.48 3.34 2.05 0.00 0.567 0.575 0.576 0.577 0.577 0.577 0.577 0.577 0.577 0.001 17.46 17.29 16.75 15.86 14.62 13.04 11. 15 8.96 6.50 3.84 1. 15 0.08 0.05 0.03 0.02 0.02 0.01 0.008 0.006 0.004 0.000 +0.81 22 10 17.38 17.24 16.72 15.84 14.60 13.03 11. 14 8.95 6.50 3.84 1.96 0.00 0.571 0.576 0.577 0.577 0.577 0.577 O.S77 0.577 0.577 0.001 21. 71 21.48 20.80 19.66 18.08 16.08 13.67 10.90 7.81 4.49 1. 15 24 1D 0.09 0.05 0.03 0.02 0.01 0.01 0.007 0.005 0.003 0.000 +0.72 21.62 21.43 20.77 19.64 18.07 16.07 13.66 10.89 7.81 4.49 1.87 0.00 0.573 0.576 0.577 0.577 0.577 0.577 0.577 0.577 0.577 0.001 

Analysis of both kinematically and statically admissible velocity fields in plane strain compression

Abstract

Similar works

Full text

Available Versions

Pure OAI Repository

Repository TU/e

NARCIS