轉(zhuǎn)向節(jié)文獻(xiàn)翻譯-一種無(wú)摩擦接觸問(wèn)題的有限元方法

1、一種無(wú)摩擦接觸問(wèn)題的有限元方法 摘要文章提出了一種新的解決包括整體可能經(jīng)歷一定的運(yùn)動(dòng)和變形的接觸性問(wèn)題的有限元方法。這個(gè)方法是基于兩結(jié)構(gòu)接觸性問(wèn)題變成兩個(gè)同時(shí)發(fā)生的問(wèn)題的分解，最終結(jié)果是在幾何學(xué)上而不是在離散的接觸表面上。一個(gè)表面接觸單元是特定設(shè)計(jì)在無(wú)條件下，允許在兩接觸表面間進(jìn)行一直的牽引動(dòng)力的傳輸。關(guān)鍵詞：無(wú)摩擦性接觸，大的變形，有限元1.引言有限元方法廣泛應(yīng)用在解決接觸性問(wèn)題上。從完全的計(jì)算角度來(lái)看，接觸性的檢測(cè)和隨后的強(qiáng)制性約束的實(shí)現(xiàn)是有待解決一般算法結(jié)構(gòu)的發(fā)展的兩個(gè)重要問(wèn)題。從康里和瑟瑞吉，禪和圖巴的早期作品以后，很多方法論都在接觸機(jī)械學(xué)的文獻(xiàn)中被提出了。一個(gè)相當(dāng)廣泛的在題目中的調(diào)查

2、在參考書(shū)目中發(fā)現(xiàn)。目前的工作是參與適合的解決在大的運(yùn)動(dòng)和變形的兩機(jī)構(gòu)的接觸性問(wèn)題有限元分析方法的發(fā)展。這一類(lèi)問(wèn)題在許多實(shí)際應(yīng)用中有特別重要的意義，如金屬成形過(guò)程和車(chē)輛事故分析。各種商業(yè)和科研計(jì)算機(jī)編碼采用的算法，以解決這種問(wèn)題。拉格朗日算子方法和它們的規(guī)格化（罰數(shù)和增廣的拉格朗日方法）是典型地被用在執(zhí)行不可測(cè)知中。 工作性條件的集成方法的選擇在線性動(dòng)力平衡方面的微弱形式下與接觸牽引動(dòng)力有聯(lián)系并在接觸原理建設(shè)中起著關(guān)鍵的作用。節(jié)點(diǎn)正交的使用包括一個(gè)（二個(gè)）點(diǎn)或一個(gè)（兩個(gè)）基座或點(diǎn)和面的接觸。假設(shè)存在一些連續(xù)的接觸表面，另外的整合法則也都是合適的。根據(jù)兩機(jī)構(gòu)問(wèn)題接近于連續(xù)的兩同時(shí)發(fā)生問(wèn)題，目前的文

3、章的重要貢獻(xiàn)是鑒定一般的程序。像傳統(tǒng)的雙行程運(yùn)算法則，兩相互作用機(jī)構(gòu)的表面被用來(lái)分析而不需要采用中間的接觸表面(任意的被選擇)。提出的方法的重要優(yōu)勢(shì)與雙程節(jié)點(diǎn)上表面的運(yùn)算法則相比，如果允許一個(gè)一體化規(guī)則的簡(jiǎn)單解釋用于接觸面，并為可容許領(lǐng)域的適當(dāng)選擇，允許一個(gè)物體到另一個(gè)物體的連續(xù)壓力的精確傳播。根據(jù)源于鐐銬的工作中的斑貼試驗(yàn)，及其以后的概括，恒壓（在大小和方向）的精確代表性的能力被看做是一個(gè)健全的必要條件和總體接觸算法的收斂。接觸力學(xué)的一個(gè)簡(jiǎn)潔的闡述在第2節(jié)中出現(xiàn)，特別強(qiáng)調(diào)在接連發(fā)生的算法的發(fā)展中常常用公式表示。在第3節(jié)中提出并分析了一個(gè)二維接觸單元，而在第4節(jié)中，提出和討論了正在使用的這個(gè)單

4、元的選擇數(shù)值模擬的結(jié)果。在第5節(jié)中，給出了結(jié)束的評(píng)論。2. 兩物體接觸問(wèn)題考慮物體認(rèn)為等同于線性空間的開(kāi)連通集，配有標(biāo)準(zhǔn)基(,)和歐幾里德準(zhǔn)則。至少有一個(gè)物體是假定可變的。在參考位形中的一個(gè)典型的重要點(diǎn)是在線性構(gòu)造中重要點(diǎn)用向量數(shù)學(xué)地說(shuō)明，對(duì)于每一個(gè)t給出，移植向量是通過(guò)來(lái)規(guī)定的。假定至少在內(nèi)，映射在定義域內(nèi)是連續(xù)和可轉(zhuǎn)置的。向量，在線性結(jié)構(gòu)（以t變化）中它的范圍是用和來(lái)確定的，因此，則有和 。并且的外在的單位標(biāo)準(zhǔn)記作。多個(gè)物體（包括一個(gè)物體）的任何系統(tǒng)的提議都受制于物質(zhì)的不可測(cè)知性，在引文11（224頁(yè)）被Truesdell和Toypin規(guī)定。這意味著這兩機(jī)構(gòu)問(wèn)題一直有 （1）在任何給定的時(shí)

5、間上，所說(shuō)的這兩個(gè)機(jī)構(gòu)都與它們的邊界子集C有聯(lián)系，當(dāng)且僅當(dāng) (2)根據(jù)上面的定義，每個(gè)物體的定義可以被唯一的分為三個(gè)相互排斥的區(qū)域，根據(jù)下面的式子當(dāng)?shù)依锟巳R和諾艾曼邊界條件分別地被強(qiáng)加在 上 。雖然沒(méi)有明白的之處，但是一般依照時(shí)間來(lái)說(shuō)，它應(yīng)該不包含, 和 C。不連續(xù)函數(shù)，可能是多值的，每個(gè)物體的邊界按照如下的定義：對(duì)每一個(gè)，是給定的為 （2a）當(dāng)是在這樣的式子中看做指數(shù)是1。的凸面性使的值是惟一的，盡管有如此的一個(gè)幾何條件的限制也不能強(qiáng)加在開(kāi)端。的一個(gè)完整的相似的定義生成下面的結(jié)果： (2b)又有對(duì)每一,滿(mǎn)足。定義方程式(2a) 和(2b)即指不連續(xù)函數(shù) 和 在C上是相等的，都等于0.也就是

6、(3)因此不可測(cè)知條件（1）根據(jù)上面的不連續(xù)函數(shù)被重寫(xiě)作。在不存在慣性效應(yīng)時(shí)，方程式的局部形式控制著如下給出的每個(gè)物體的運(yùn)動(dòng)，div in (4a) on （4b） on (4c) on (4d)當(dāng)是克西的應(yīng)力張量，是質(zhì)量密度，是每單位質(zhì)量的質(zhì)量力，是法定的界壁位移，是在上的法定牽引向量。標(biāo)準(zhǔn)加權(quán)殘值法的應(yīng)用，與拉格朗日乘數(shù)的引入p 0對(duì)費(fèi)解的限制相結(jié)合，結(jié)果在運(yùn)動(dòng)方程的若形式中指出，位移解方程（ 4 ）和拉格朗日乘數(shù)外地p滿(mǎn)足 (5a) (5b)和q是對(duì)所有的允許函數(shù)。在大部分的缺失不被用在方程(5a)中時(shí)，則不連續(xù)函數(shù)定義在C上。位移區(qū)域?qū)儆诳臻g即質(zhì)量函數(shù)也屬于空間，定義為允許函數(shù)p 0是

7、分段連續(xù)的。為了證明方程式(5a)的正確性，把這項(xiàng)工作看做是允許函數(shù)和上沿C的觸點(diǎn)壓力，即為 (6)在C上，在沒(méi)有摩擦和回顧時(shí)，柯西引理應(yīng)力矢量意味著 (7)借助于方程(7)，方程(6)可寫(xiě)為這表明拉格朗日乘數(shù)區(qū)域是自然等同于正常的在接觸范圍內(nèi)的牽動(dòng)引力（壓力）。壓力場(chǎng)p一般假定僅僅是分段光滑的，因此認(rèn)為每個(gè)物體的特征材料的接觸面在C的附近。而且由于 on (8)不等式(5b)是從(4d)的形式中得到的，而且假定q是非負(fù)的。為了進(jìn)一步證明拉格朗日乘數(shù)在兩物體接觸問(wèn)題中的作用，為接著發(fā)生的近似數(shù)值提供了一些動(dòng)機(jī)，用(2b) 和(8)改寫(xiě)(5b)中在接觸面C上構(gòu)成整體所必需的，即 (9)符號(hào)僅僅被

8、用來(lái)強(qiáng)調(diào)接觸面C不普遍的看做的一部分，就像(3)中所表明的。整式(9)表明，在每個(gè)物體的范圍內(nèi)，適當(dāng)?shù)慕o定一個(gè)拉格朗日區(qū)域，(5b)可以在每個(gè)接觸面上分別作用。然而，很顯然區(qū)域應(yīng)該滿(mǎn)足在（一般地）接觸面上的線動(dòng)量的平衡，即為以上資料數(shù)據(jù)將被應(yīng)用在(5a) 和 (5b)的近似解決方案中。.(9)式的推導(dǎo)用下面的程序，等式(5a) 和(5b)被記作 (10a) (10b)當(dāng)，對(duì)所用可容許的和都適用。(5a) 和 (5b)的典型的罰數(shù)調(diào)整是在的條件下獲得。 當(dāng) 是一個(gè)Macauley支架和時(shí)，上面的（5a）式的積分關(guān)系是一個(gè)連續(xù)的損失參數(shù)。約束條件(4d)是無(wú)約束的（因此允許受約束的滲透來(lái)替代）替代

9、區(qū)域是確定的，則有 (11)在約束條件下，原始的無(wú)約束問(wèn)題等同與一個(gè)凸面的最小化問(wèn)題，序列的解決方案當(dāng)時(shí)表明收斂是約束問(wèn)題的解決方案(5).實(shí)踐中，罰數(shù)方法成功的用在平衡狀態(tài)下與最小總勢(shì)能不相關(guān)的情況下。由于罰數(shù)參數(shù)的有限值的普及率，接觸面C的一個(gè)獨(dú)特的定義和缺陷函數(shù)并不是容易地應(yīng)用在(11)中的。在提出的明確表達(dá)中關(guān)鍵的一步是(10a) 和 (10b),易于接受的規(guī)則化的罰數(shù)的引入，例如由(12)，等式(10a)可以寫(xiě)為 (13)上面的(13)中的任意兩個(gè)積分在形式上是相同的，即積分邊界來(lái)源于一個(gè)新問(wèn)題的規(guī)則化罰數(shù)。盡管在目前的背景下這兩種積分在缺陷函數(shù)的兩個(gè)定義(2a) 和(2b)的基礎(chǔ)上

10、是有聯(lián)系的，但是這種分解要設(shè)計(jì)要涉及到重要的計(jì)算，這在下一節(jié)中再討論。這里提出的事態(tài)發(fā)展可以很容易地靠歸納兩體問(wèn)題擴(kuò)展到多體接觸問(wèn)題。3. 三維接觸的兩個(gè)應(yīng)用這篇文章的其余部分專(zhuān)門(mén)討論在明確表達(dá)罰數(shù)的基礎(chǔ)上的離散二維接觸問(wèn)題和兩表面上的缺失函數(shù)的鑒定，像在第二節(jié)中所述。這里待確定的兩個(gè)關(guān)鍵問(wèn)題是接觸單元的幾何作圖和有限元近似值的可允許區(qū)域的選擇。3.1 接觸面的離散化如果有一個(gè)接觸面的離散化是連續(xù)的，那么接觸面 和 是分別在 和的條件下是唯一確定的。撇開(kāi)推行細(xì)節(jié)，接觸面是被一系列正常的從 到的投影（不定靠近點(diǎn)）所離散，如圖表2中所示。的離散化采用相同的步驟。因此每個(gè)接觸單元在接觸面上的立體空間

11、元素與其對(duì)立面相聯(lián)系。顯然地，非光滑邊界的離散會(huì)導(dǎo)致單位法線和缺失函數(shù)的不連續(xù)。這里目的不是設(shè)法回避這個(gè)問(wèn)題，盡管它的結(jié)果不可忽略，特別是在滾動(dòng)問(wèn)題中。對(duì)于一個(gè)特殊的離散來(lái)說(shuō)會(huì)導(dǎo)出代表兩個(gè)層面的光滑邊界，見(jiàn)參考資料13 。3.2 有限元區(qū)域在這節(jié)中假設(shè)一個(gè)具體的接觸有限元是建立在方程式(10a) 的(12)的基礎(chǔ)上。選擇有限維領(lǐng)域的指導(dǎo)方針是以前的作品特別是在青年受理領(lǐng)域位移和壓力導(dǎo)致單方面接觸單元，能夠穩(wěn)定地復(fù)制費(fèi)解的條件（即滿(mǎn)足的LBB條件）是準(zhǔn)確的。這種收斂性分析，目前盡管并不適用于這里的兩體接觸問(wèn)題的動(dòng)非線性這方面，但是它為可允許區(qū)域的選擇提供了一個(gè)準(zhǔn)則。數(shù)值積分典型地把所有邊界術(shù)語(yǔ)用

12、在接觸面上。不考慮離散化使用空間的元素（例如三角形的三個(gè)定點(diǎn)和四邊形的四個(gè)點(diǎn)）的線性規(guī)則。在單個(gè)接觸單元中的缺失函數(shù)是關(guān)于軌跡邊界坐標(biāo)系的非標(biāo)準(zhǔn)函數(shù)。因此，一般情況下在上的數(shù)值積分是不準(zhǔn)確的，所有的引入錯(cuò)誤的積分規(guī)則直接影響到用公式的確切表示。所提出的接觸單元是建立在從發(fā)源到標(biāo)準(zhǔn)Q9元素（等參的九個(gè)節(jié)點(diǎn)）的移置區(qū)域和引用的辛普森標(biāo)準(zhǔn)上的。A NOVEL FINITE ELEMENT FORMULATION FOR FRICTIONLESS CONTACT PROBLEMSSUMMARYThis article advocates a new methodology for the finite

13、 element solution of contact problems involving bodies that may undergo finite motions and deformations. The analysis is based on a decomposition of the two-body contact problem into two simultaneous sub-problems, and results naturally in geometrically unbiased discretization of the contacting surfa

14、ces. A proposed two-dimensional contact element is specifically designed to unconditionally allow for exact transmission of constant normal traction through interacting surfaces.KEY WORDS frictionless contact; large deformations; finite elements1. INTRODUCTIONFinite element methods are used extensiv

15、ely in the solution of contact problems. From a purely computational standpoint, detection of contact and subsequent satisfaction of the impenetrability constraint are the two key issues to be addressed in the development of a general algorithmic framework. Numerous methodologies have been proposed

16、in the literature of compulational contact mechanics since the early works of Conry and Seircg,and Chan and Tuba.A fairly comprehensive survey on the topic is found in Reference.The present work is concerned with the development of finite element methods suitable for the solution of two-body contact

17、 problems in the presence of large motions and deformations. This class of problems is of particular significance in numerous practical applications, such as metal forming processes and vehicular crash analyses. Various commercial and research computer codes employ algorithms for the solution of suc

18、h problems. Lagrange multiplier method,and their regularizations (penalty and augmented Lagrangian methods) are typically used in enforcing impenetrability.The choice of integration method for the work-like term associated with the contact tractions in the weak form of linear momentum-balance plays

19、a pivotal role in the construction of contact elements. Use of nodal quadrature involving the contacting nodes of one of (resp. both) surfaces yields the standard one-pass (resp. two-pass) node-on-surface algorithms.Other integration rules are also applicable, provided there exists a continuous disc

20、retization of the contact interface.The main contribution of the present paper is in the identification of a genneral procedure according to which the two-body problem is approached as a sequence of two simultaneous sub-problems. As in the traditional two-pass algorithms, the surfaces of both intera

21、cting bodies are used in the analysis without need for introduction of an (often arbitrarily chosen) intermediate contact surface. The main advantage of the proposed approach over the two-pass node-on-surface algorithms if that it allows for a straightforward interpretation of the integration rules

22、used on the contacting surfaces and, for appropriate choices of admissible fields, permits the exact transmission of constant pressure from one body to another. In the spirit of the patch test originating in the work of Irons,and its subsequent generalizations, capability for exact representation of

23、 constant pressure (in both magnitude and direction) is viewed as a necessary condition for robustness and convergence of the overall contact algorithm.A brief exposition to contact mechanics is presented in Section 2, with particular emphasis on formulations to be used in the ensuing algorithmic de

24、velopments. A two-dimensional contact element is proposed and analysed in Section 3,while the results of selected numerical simulations using this element are presented and discussed in Section 4.Concluding remarks are given in Section 5.2. THE TWO-BODY CONTACT PROBLEMConsider bodiesidentified with

25、open and connected sets in linear space , equipped with canonical basis (,) and the usual Euclidean norm. At least one of the bodies is assumed to be deformable. A typical material point of in the reference configuration is algebraically specified by vector of material point in the current configura

26、tion is given byat each time t, and the displacement vector is defined according toThe mapping is assumed smooth throughout its domain and invertible at least on .Body and its boundary in the current configuration (at time t) are identified with and ,respectively, henceand .Also, the outer unit norm

27、al to is denoted by .The motion of any system of bodies (including a single body) is subject to the principle of impenetrability of matter, as stated by Truesdell and Toypin in Reference11 (p.244).This implies for the two-body problem that at all times （1）At any given time, the two bodies are said t

28、o be in contact along a subset C of their boundaries if, and only if, (2)It follows from the above definetion that the boundary of each body can be uniquely decomposed into three mutually exclusive regions according toWhere Dirichlet and Neumann boundary conditions are enforced on and , respectively

29、. Although not explicitly noted, it should be clear from the above that , and C generally depend on time.Gap functions , possibly multi-valued, can be defined on the boundary of each body as follows: for each , is given by （2a）Where is such that see Fifure 1.Convexity of renders single-valued, altho

30、ugh such a restrictive geometric condition will not be imposed at the outset . A completely analogous definition for yields (2b)Where,again, for each , satisfies .Defining equations (2a) and (2b) imply that gap functions and are identically equal to zero on C, namely that (3)Consequently, impenetrab

31、ility condition (1) can be rewritten in terms of the above gap functions asAt the absence of inertial effects, the local form of the equations governing the motion of each body are given bydiv in (4a) on (4b) on (4c) on (4d)Where is the Cauchy stress tensor, the mass density, the body force per unit

32、 mass,the prescribed boundary displacement, and the prescribed traction vector on .Application of the standard weighted-residual method,in conjunction with the introduction of Lagrange multiplier p0 for the impenetrability constraint,results in the weak form of the equations of motion, which states

33、that the displacement solution of equations (4) and the Lagrange multiplier field p satisfy (5a) (5b)For all admissible function and q. Without loss of generality is used in equation (5a) for the definition of gap function on C. Displacement fields belong to spaces withand weight functions belong to

34、 spaces defined asThe admissible functions q0 are piecewise continuous.To prove that equation (5a) holds, note that the work done by contact forces along C on admissible functions and is given by (6)At the absence of friction and recalling that on C, Cauchys lemma on the stress vector implies that (

35、7)With the aid of equation (7), equation (6) is written as,Which shows that the Lagrange multiplier field is naturally identified with normal traction (pressure) along the contact region. The pressure field p is generally assumed to be only piecewise smooth, thus allowing each body to feature materi

36、al interfaces in the neighbourhood of C. Moreover, since on (8)inequalities (5b) follow from (4d) and the assumed non-negativeness of q.In order to further clarify the role of the Lagrange multipliers in the two-body contact problem and provide some motivation for the ensuing numerical approximation

37、s, use (2b) and (8) to rewrite the integral in (5b) on surface C as (9)The notation is employed to merely emphasize that the contact surface C can be selectively viewed as part of ,as indicated in (3).The integral expression (9) suggests that, given a properly defined Lagrange multiplier field on th

38、e boundary of each body, integration of (5b) can be performed separately on each of the contacting surfaces. However, it is clear that fields should satisfy balance of linear momentum on the (common) contact surface, namely that The above observation will be exploited in the approximate solution of

39、(5a) and (5b). Folowing the procedure used in the derivation of (9), equations (5a) and (5b) are rewritten as (10a) (10b)for all admissible and ,where .The classical penalty regularization of (5a) and (5b) is obtained by settingin the last integral term of (5a),where is the Macauley bracket and is a

40、 constant penalty parameter. Constraint conditions (4d) are relaxed (thus permitting controlled penetration to take place) and displacement fields are determined so that (11)for all .Under restrictions such as formal equivalence of the original unconstrained problem to a convex minimzation problem，t

41、he sequence of solutions as can be shown to converge to the solution of the constrained problem (5).In practice, penalty methods are used successfully even when the equilibrium state is not associated with minimization of the total potential energy.Due to the occurrence of penetration for finite val

42、ues of the penalty parameter ,a unique definition of the contact surface C and gap function g is not readily available for use in (11).A crucial step in the proposed formulation is the introduction of a well-defined penalty regularization for (10a) and (10b),such that (12)With the aid of (12),equati

43、on (10a) becomes (13)Each of the last two integrals of (13) identical in form to the boundary integral emanating from penalty regularization of a Signorini problem. Although in the present context these two integrals are clearly coupled by the definitions (2a) and (2b) of gap functions , this decomp

44、osition has important computational implications, as will be discussed in the next section.The developments presented here can be easily extended to encompass the muli-body contact problem by reducing it to a series of coupled two-body problems.3. APPLICATION TO TWO-DIMENSIONAL CONTACTThe remainder

45、of this article is devoted to the discretization of two-dimensional contact problems based upon the penalty formulation and the identification of gap functions on both surfaces, as suggested in Section 2.The two key issues to be addressed here are the geometric construction of contact elements and t

46、he choice of admissible fields for the finite element approximation.3.1. Discretization of the contact surfacesA continuous discretization of the contact surfaces is advocated. Surfaces and are uniquely determined as those on which and, respectively. Setting aside implementational details, contact s

47、urface is discretized by a series of normal projections (not necessarily closest-point) from to ,as shown in Figure 2. An analogous procedure is followed for the discretization of .Consequently, each contact element relates a single spatial element edge on surface to the opposite surface.The apparent non-smoothness of the discrete boundaries results in