第34卷第2期 Vol 34 No. 2 CHINESE QUARTERLY OF MECHANICS June 2013 基于郭仲衡先生现代张量分析 及有限变形理论知识体系的相关研究 谢锡麟 (复旦大学力学与工程科学系,上海200433) 摘要:作为我国理性力学先驱之一的郭仲衡先生在其所著《张量(理论和应用)》以及《非线性弹性理论》中记述 了现代张量分析以及有限变形理论知识体系。本文按有限维 Euclid空间上微积分以及一般赋范线性空间上微 分学认识相关知识体系的理论框架,相关思想及方法,闸述了有关思想及方法的发展及其应用。本文未涉及现 代几何学在连续介质力学中的应用。 关键词:郭仲衡;现代张量分析;有限变形理论;有限维 Euclid空间上微积分;一般赋范线性空间上微分学 中图分类号:0331文献标志码:A文章编号:0254-0053(2013)02-337-15 Some Related Studies Based on Guo zhong-heng's Knowledge Systems of Modern Tensor Analysis and Finite Deformation Theory XIE XI-Lin Department of Mechanics Engineering Science, Fudan University, Shanghai 200433, China) Abstract: It is well known that Prof. Guo Zhong -heng is one of the national pioneers of the foundations and developments of rational mechanics. His two monographs one is"Tensor( theories and applications) and the other "Nonlinear Elastic Mechanics"include systemetic knowledge of modern tensor analysis and finite deformation theory. The theoretical frameworks of related know ledge systems with some ideas and methodologies have been expatiated in the point of view of calculus in finite Euclidian spaces and differen tial calculus in general normed linear spaces, including some developments of related ideas and methodolo gies with applications. The applications of modern geometry in continuum mechanics have not been in volved in the present paper. Key words: Guo Zhong-heng: modern tensor analysis; finite deformation theory; calculus on finite eu clidian spaces; differential calculus on general normed linear spaces 作为我国理性力学先驱之一的郭仲衡先生在其所著《张量(理论和应用)》以及《非线性弹性理论》 中记述了现代张量分析以及有限变形理论的知识体系 郭仲衡所著《张量(理论和应用)》,主要内容包括:①张量的代数性质(张量定义为多重线性映照)。② 仿射量的基本性质(基于外积运算)。③张量值映照微分学(含各向同性张量值映照的表示理论等)。④微 分几何中曲线论与曲面论的基本内容(主要包括局部标架及其运动方程)。⑤现代几何学中相关思想及方 收稿日期:2012-10-11 基金项目:上海市教委2011年上海高校本科重点教学改革项日“‘现代连续介质力学理论及实践’课程体系”;上海市教委2011年重点 课程立项项目“《数学分析》一年制,面对力学等技术科学专业)”;国家自然科学基金面上项目(11172069,10872051) 作者简介:谢锡麟(1974-)男,浙江鄞县人,副教授,博士,研究方向:理学力学观点下的连续介质力学理论,力学中的数学方法并将上 述理论应用于开放流场空间动力学行为等研究。Emal: xiexin@ fudan.edu,cn
!!"" !## #$%!$&% & ' ( ) '()*+,+-./01+023456+'(/*)', 789:!"*8:# ;#$%! +,noÕp]q÷Àn9_ Z¼òmtxQrUÌÁ*L4ÉÊ A !ëW*' &'S-0¤'À"áµ #$$"!!$ 23!uv ,`&'J8~ èi¤¥¦»'¨X1\A9A=!UL?@!/ *L*-01*491+ .8&.4*"%T*%+."?%0-@+#+0%4H#%#1*N*8.">01#.%T,*."@ S*/S*1)3' !J>YG@DS>=D8U6>IVG=AIKd+=BA=>>@A=B,IA>=I>"5@KADW",VG=BVGA#$$"!!"'VA=G$ ?=+1"0/1&)DAKT>99H=8T=DVGDC@8U:F=BAK8=>8UDV>=GDA8=G9YA8=>>@K8UDV>U8E>98YS>=DK8U@GDA8=G9S>IVG=AIK:(AKDT8S8=8B@GYVK8=>AK11>=K8@!DV>8@A>KG=LGYY9AIGDA8=K$2 G=LDV>8DV>@1*8=9A=>G@+9GKDAI6>IVG=AIK2A=I9KWKD>S>DAIH=8T9>LB>8US8L>@=D>=K8@G=G9WKAKG=L UA=AD>L>U8@SGDA8=DV>8@W:1V>DV>8@>DAIG9U@GS>T8@HK8U@>9GD>LH=8T9>LB>KWKD>SKTADVK8S>AL>GKG=L S>DV8L898BA>KVGE>X>>=>\YGDAGD>LA=DV>Y8A=D8UEA>T8UIG9I+KG=LLAUU>@>=M DAG9IG9I=>@G9=8@S>L9A=>G@KYGI>K"A=I9L>E>98YS>=DK8U@>9GD>LAL>GKG=LS>DV8L898M BA>KTADVGYY9AIGDA8=K:1V>GYY9AIGDA8=K8US8L>@=B>8S>D@WA=I8=DA=IVG=AIKVGE>=8DX>>=A=M E89E>LA=DV>Y@>K>=DYGY>@: !*@5."4+&F=B#S8L>@=D>=K8@G=G9WKAK#UA=AD>L>U8@SGDA8=DV>8@W#IG9I+K#LAUU>@>=DAG9IG9I=>@G9=8@S>L9A=>G@KYGI>K uv ,`&'J8~ <»gGJ¸oÚDº:}ç!,1"T5$;(%) d:×#`ç`,1;(#) qU?R<}çsVd¬Í[,1<^þÝÀ% »gGDº:}ç!,1"T5$;"ïðÞîÎÏ&8}ç<I`j!}çàkvy¥#`˲$%9 óNç<6i`j!67üÏ$%:}ç;˲ûs'!/fÊ`}ç;˲<,1{$%;û s5qÅ#1SÅï1<6iÞî!ïðÎÏ@ALÚÏ9W0$%<<5'qS}E.W
力学季刊 第34卷 法(包括基于同态映照的推前及拉回,Le导数, Hodge星算子,内导数,外微分以及相关运算之间的关系) 对此部分内容的叙述虽然未引入微分流形的概念,但所述的相关思想及方法可以几近完全地移植于流形 上的分析,且数学分析上非常清晰。⑥张量分析在连续介质(几何形态默认为 Euclid流形)中的基本应 用,包括变形刻画,输运方程;另涉及同态扩张以及Le导数等在连续介质力学中的应用,但书著中未对这 部分内容做深入阐述 郭仲衡所著《非线性弹性理论》,主要内容包括:①有限变形理论(连续介质几何形态默认为 Euclid流 形)。理论框架上分别对初始物理构形以及当前物理构形引人曲线坐标系,理论发展上按变形梯度及其基 本性质,变形刻画,输运方程,守恒律方程等。②有限变形弹性静力学、有限变形弹性动力学若干典型事例 的半解析求解。③变分原理。值得指出,基于《张量(理论和应用)》所载张量分析的知识体系,研习《非线 性弹性理论》就显得较为自然而无数学以及力学分析上的困难 本文作者自上世纪90年代初作为大学生就被郭仲衡先生的书著所吸引并持续研习至今。本文主要 基于有限维 Euclid空间上微分学以及一般赋范线性空间上微分学认识和推广相关知识,主要包括:(1) 有限维 Euclid空间中的微分同胚。(2)张量场微分学。(3)张量值映照微分学。(4)张量场多点表示形 式。(5)非完整基理论。(6)基于外积运算的相关理论 1有限维Eucd空间中的微分同胚 1.1空间曲线坐标系(空间微分同胚) 有限维 Euclid空间中的C"微分同胚,指向量值映照 ):R"→D3X(x)∈R 此处定义域D和值域Dx:=X(D2)均为R"中的开集,X(x)实现了二者之间的双射(一一对应),且X (x)与其逆映照x(X)均为CP映照;一般可记C°微分同胚为X(x)∈CP(D4,Dx)。 g3(r) curvilinear coordiante g3(x) X()ECP(D D) g2(x) g1(3 g2( (g;(x) ocal covariant bas D()=[. 图1三维 Euclid空间中曲线坐标系示意图 Fig. 1 Sketch of curvilinear coordinates in three dimensional Euclid space 如图1所示,基于微分同胚:①物理空间Dx中点的位置刻画可等价地由参数区域D中的位置进行 刻画;一般构造微分同胚时使得D的几何形态规则。②由微分同胚X(x)∈CP(D-,Dx)的 Jacobi矩阵 可定义局部协变基{g,(x)△aX/ax2(x)}1,按对偶关系可唯一确定其对偶基{g(x)|(g2,g;)=b 1,…,m}=1;称为逆变基,协变基和逆变基统称为局部基;可基于局部基展开整体形式的张量方程以 获得可实际求解的分量方程。微分同胚的上述性质实现了一般曲线坐标系的意义及功能。郭仲衡著《张 量(理论和应用)》明确将曲线坐标系定义为微分冋胚,凸显其数学严谨且几何意义清晰
B!ÎÏ67ÊÙ˲NI"(8LB>ï4"ÞNI"üûsdS}Ï~{NI{o+aj&'qÁ"î3Y9ñs# ¬K + ûsÊ_"ç;˲ S!9$&HF IT9J924S!9$0HF bàk| T9 ";|TS& ˲# QcU L> ûsÊ_vS!9$0L> !T9"TS$(")% 8! ¿D J-69"'î37û|Í>Áqr8 !"#$! 45(*67&86-+:"9".(3+6&&+'".3*()".*7+(('"0(.)"&.39J-69"')236( Z%D"67ûsÊ_&@Ü,3{ TS q¦ !T9"TS$<;GI8XA càk@AvÍ6,(3!9$K"S/"93!9$-F 3<%"9íx}Àc` àÚíx6,(2 !9$2!(2 "(3$GF <32 3" 3<%"*"F-F 2<%"rvUÍ6"vÍ6"UÍ6Arv@A6#c67@A6^Ý[x<}çW0d äFcgKDE<sçW0%ûsÊ_<á?`jg<R QÅ#xLÀ<nkÞ\%»gGº:} ç!,1"T5$;;üÅ#xLÀàkvûsÊ_"a<ÚI'®b5nk±\% 774 & ' ( ) !!""
第2期 谢锡麟:基于郭仲衡先生现代张量分析及有限变形理论知识体系的相关研究 R(,9) X 图2当前物理构型对应之曲线坐标系显含时间情形的一则示意图 Fig. 2 A sketch of curvilinear coordinates with respect to current physical configurations including time explicitly 1.2时空曲线坐标系(时空微分同胚) 如图2所示,当类马铃薯体(当前物理构型)随时间发生变化时,我们希望能建立显含时间的曲线坐标 系,使得对应的参数区域不仅几何形态规则并且不随时间发生变化。对此,可考虑如下当前物理构型对应 之显含时间的曲线坐标系 g,t)·rsin0·cos X(x, t) X(a. t) X2|(x,t) R(0,g,t)· rsing·sing D×(a,b)3 R(0,,t)·rcos 在时空坐标系(世界坐标系)中当前物理构形对应之显含时间的曲线坐标系可严格符合微分同胚的要求。 以上述映照为例,建立时空微分同胚的充分必要性条件为:①X(x,t)∈CP(Dm×R;R3),亦即X(x,t) 关于时空变量{x,t}都具有直至p阶偏导数且均在Dm×R上连续;②X(x,t)关于空间变量x在Dw上 为单射;③X(,t)关于空间变量x的 Jacob矩阵D2X(x,t)∈R在Dm×R上都非奇异。应用上可 基于上述条件,确定一个CP映照为C微分同胚。 如果我们可对当前物理构型引入显含时间的曲线坐标系X(x,t)∈C"(v,V)(现t视为参数),使 得对应的当前参数构型几何形态规则(如为方体)且不随时间变化,则连续介质运动在物理世界坐标系(时 空坐标系)中的区域V×(t,to+T)和参数世界坐标系中区域V,X(to,to+T)之间有微分同胚x(x t)∈CP(V×(to,to+T),V×(t,to+T))(现t视为独立变量),如图3所示。我们已基于一般有限 变形理论[,平行发展了当前物理构形对应之曲线坐标系显含时间的有限变形理论。籍此理论,我们可 系统获得连续介质运动在参数世界坐标系下的分量方程,在参数世界坐标系下连续介质的参数构型演化 对应一个规则的方体,现情形下物理构型的几何非规则且随时间变化的边界被严格映射为几何规则且不 随时间变化的平面或者直线,这将十分有利于数值及理论分析 1.3基于曲面构建全空间半正交系 如图4所示,为研究边界及其邻域上的流动机制,往往可基于曲面构建全空间维数曲线坐标系,如下所示 X(x,t)=X(x,123 :/→X(xx, X n(xy,t)∈R X
8" :-xÂbÒþtû|Í>ÁTÓ3 t*×qr8 !"#$" I)5(*67&86-+:"9".(3+6&&+'".3*()Á"îY9ñs# Z#D"îacdÝ!ÐÜ,Y$ãz{ä¸ÍÅz"2e\¹¿ !TA%!EH#H!$"G S!9";$ }7z3Íç,9";-Ee > [úNIÖoTA%!EHá+a#AS!9";$}73{Íç9 oTA%!á vvN#BS!9";$}73{Íç9 ˲vL> ûsÊ_% ZÙcíÐÜ,Y9: !J ; 9"J ; $! !J ; 9E!;$";$=8$"J ; E!;$";$=8$$!<;ÿv¿Íç$"Z!D%ý67 Q¬ Í[,1(#)"ûäRÐÜ,Y[íT~Å#xLÀ</z{<¬Í[,1(N)%fb,1"c ÀAäF+ajÏ9oHIylxLÀ[<sçW0"oHIylxLÀ[+aj<HIYÂÅ íT lú}<WÝ"<;[[Ü,Y<5×ú}ãz{ÍÅ<klq®?ËNv5ú}- ãz{ÍÅ<ûïÓ«e#"üxs/7I;,1sV(N)% AID +,ûÂü½î3aÁ Z"D"v©ªklÚ|á<æ9ÇO"ÁÁc67ÅïY¹ò3{KIÅ#xLÀ"Z[D S!9";$<S!9$ "9F ";$&T9J 9$ (9F )< 9% $ + 9FB% & ' ( $ ) 9 & ' ( ) F 4S!9$ "9F ";$ L S% + SFB% S & ' ( )F !9$K$!9$ ";$=9F 0'!9$ ";$0HF !## ÐY«&67»gGJ¸<}çsV¬Í[,1^þÝÀ<S}©ª 775
力学季刊 第34卷 (o,+7) >to X(5)∈C|v, x(x,)∈C 51) V×(0,4+7) 图3时空坐标系下当前物理构形对应之曲线坐标系构成微分同胚示意 设想参数构形的几何形态规则且不随时间变化 Fig. 3 Sketch of the diffeomorphism with respect to current physical configurations in spatial-temporal ordinates in which the parametric configurations are supposed as geometrical regular and invariant with time Boundary of Solid Body x-Cine Curve (x)=2(x,…,) 图4m-1维曲面上标架以及基此构建的半正交m维(全维)曲线坐标系示意图 Fig. 4 Sketches of local frames on a m-1 dimensional surface and the full dimensional semi-orthogonal curvilinear coordinates 其空间局部协变基为 t)△ aX -(xx,c",t)=g;(xy,t)+ (x,t)=(}-x"b 1)g;(xx,t) 在曲面上此局部基正好为曲面自身所诱导的局部基9(x2,t)A少(,b)以及曲面单位法向量
8% îÍ>Áê:-xÂtÒþtû|Í>ÁÂåY9ñsqr8$ Ðu#$Ât*89tGP×vðº3m³$ !"#$% 45(*67&8*7('"88(&0&+27")0 ¨QZ+wÂü*a 0 D"½D#û|Í>Áqr8 !"#$' 45(*67()&89&6398+30()&.30Y!'"0(.)"&.39)-+836(3.'*7(8-99 '"0(.)"&.39)(0"&+*7&#&.396-+:"9".(3+6&&+'".3*() Ú3{@AvÍ6v (3!9$ "9F ";$K "S "93 $ !9$ "9F ";$<(3 $ !9$ ";$=9F "' "93 $ !9$ ";$<!32 3B9F 72 3$(2 $ !9$ ";$ (F !9$ "9F ";$K "S "9F $ !9$ "9F ";$<'!9$ "; - . / $ oÅïáb@A6vÅï¢uDpN<@A6(3 $ !9$ ";$K "$ "93 $ !9$ ";$dÅ ï v w B ç 78" & ' ( ) !!""
第2期 谢锡麟:基于郭仲衡先生现代张量分析及有限变形理论知识体系的相关研究 n(xs,t)。需指出,上述基于曲面的半正交曲线坐标系为 Euclid曲线坐标系,换言之对梯度算子可有表达 形式 V仝aX (x,t)=i ve(X, t)-a(, t)=g(, t)g(, t),a, 式中{“}m=1代表典则基,{g2(x,t)}=1代表局部逆变基。在区域上开展场论分析,是否采用曲线坐标系 或者采用何种曲线坐标系不影响张量场梯度等场论微分运算的整体形式,但采用基于曲面的半正交系研 究含有曲面边界的流固耦合问题易于揭示边界的几何特征同流动行为之间的关系,因为曲面的几何性质 自然出现于曲面局部标架的运动方程之中。如果实际计算所采用的曲线坐标系不同于基于曲面的半正交 系,则可通过如下方法进行转化,以二阶张量的梯度为例 中⑧V|x=,中:(正,t)9,⑧°⑧g?(,t) ,:a{[(0n9)a"9,+(.,m)n][(,9)a“9,+(g0,n)n]⑧[(y,g)g”9.+(g0,n)an]} 等式右方即为张量梯度相对于曲面半正交系的展开式,另分析中常引人曲面梯度算子△( t),其同全梯度算子V之间有关系,以任意的r阶张量更∈T(R)为例 =(÷+n0-)0 m=V⑧-n如=V②更-n②(n·⑧更) 吴介之等籍借R3中曲面半正交系,获得了可作有限变形运动的固体边界上连续介质变形率张量的 表达形式,推广了经典的 Caswell公式。另,我们可基于曲面半正交系及内蕴形式广义 Stokes公式 推导得对任意p阶张量S∈T(R3)成立的恒等式 n·(V×(n×S)=(nxV)·(n×S)=b(n·S)+y·S=H(n·S)+·S H△b为曲面平均曲率 考虑一般可压缩粘性流动的动量方程( Navier- Stokes方程) P=vI-V×a+pfm,:=-p+(A+2g)0,0:=V·V VX V 可得可作有限变形运动的固体边界上可压缩粘性流动的法向涡量通量表达式 mn×a-×(VⅡ)-Pn×fm+(V②a)·n Pn×a-n×(VⅡ)-Pn×fm-(n×V)·[(n×a)⑧n]+V(a·n)+H(a:n) 相对于静止固体边界情形,边界有限变形运动导致增加项ⅴ(ω·n)+H(a·n)n,ω·n (ⅴ×V)·n,可见此增加项完全由曲面自身运动决定,不直接涉及流固耦合。 2张量场(映照)微分学 2.1张量赋范线性空间 张量场映照指自变量为位置刻画,因变量为张量的映照。不失一般性,三阶张量场具有如下表达式 中(x):R"D3x卜中(x)=中.(x)9(x)⑧g(x)②9(x)∈T(R”) 此处位置刻画一般可为曲线坐标系。针对张量场映照的微分学以及积分学构成了一般曲线坐标系下张量 场分析的主要内容。为此,郭仲衡著《张量(理论和应用)》明确引入了张量范数
'!9$ ";$%ö"á?67Åï99Úx(&)%¨"c67Åï4ÑÀÞô[xGk P;,K&C Úx(N) QNFímn> [}çP08> !H!$ ¿@,D8H>KW0$ &+=!0=#+$%"%&<,0J".&<,EJ cFcu¬Í[Ï9<tÝklácX`h`æ9<Bhç~ç5x + ". "' $ <&'E+B'E!, $ G$B&'E?F =+!,M.$0' <&'E+B'E!, $ G$B&'E?F B!'E,$ $ 0(!'E.$M')=, $ !.0'$=.!.0'$' Sí7f 9 t Ý k l ; [(&)"k l ¬ Í [ Ï 9 N Y 6 h , $ !.0'$=.!.0'$'".0' < , $ ! EJ$0'"crb6hxòXÅï¢uÏ9ÿà"-efætù¨% B Àn"x®#Y9` BIA ÀnyÃ|î3 }çå˲¢Íçvwk"aÍçv}ç<˲%-P Q`"~[}çåEZ[5x 9!9$&HF ITQ+AJ9249!9$<93 02K!9$(3!9$M(2!9$M(K !9$08!!HF $ bwk QcvÅ#xLÀ%í}çå˲<ûs'ds'Y R QÅ#xLÀ[}ç åsV<ïðÞî%vb"»gGº:}ç!,1"T5$;;9:R}ç6I !## ÐY«&67»gGJ¸<}çsV¬Í[,1^þÝÀ<S}©ª 78%
力学季刊 第34卷 中|r"△√④⊙中=厘,,,Ⅴ西∈T”(R") 此定义不依赖于基的选取。需指出,此范数自然诱导的距离为d(φ,y):=|φ-y|P(",并使得张量 赋范空间具有完备性。我们将一般赋范线性空间上微分学应用于张量赋范空间,可得以下系统性结果 22张量场一阶导数(可微性)以及高阶导数 按可微性的定义,有 p(x+h)=(x)+(x)(b)+0(h)∈T”(R"),此处(x)∈L(R",TP(R") 以三阶张量为例,按 中(x+h)=中.(x+h)9(x+h)⑧g(x+h)⑧g2(x+h 考虑各组成项的可微性,可有 A(x+b)=,(x)+当(x)h+0:(h)∈R,式中7=0∈R 0.(x+h)=9(x)+9 aaf(ah+o: (hER: g'(a+ h)=g()+ (a)h+o(h)∈R 结合简单张量的基本性质 ②|②Q5r”=1lk”列151”,Ⅴ,,5∈R,bs ⑧⑧3+B·⑧⑧5,Va,B∈ ,7,7,S∈R 可推导得微分的表达式 迎(x)(h)=[Vp(x)9,⑧9⑧9②g()]·[h"n()]=:(②V)(x)·[h9n(x)] =[h"g(x)]·[V中:(x)9⑧9②y9k]=:[h2(x)]·(V②更)(x) 此处h°g。(x)为参数域中位移h"q对应至物理区域中的位移;式中 Vp:}(x):()+r-P:+r: 张量分析中称为张量场分量的协变导数。类似可引入逆变导数Vφ:}(x)=9V更:(x)。协变导数 与逆变导数之间满足指标升降。 进一步,考虑到自变量h=hi;的向量结构,可以引入张量场的偏导数:以三阶张量为例 D (a)=更 (x)△lim V:(x)9,⑧g⑧9(x)∈T(Rm) 且有:2(x)(h)=(x),在明确张量场偏导数后,我们就可定义任何形式的场论微分算子 d=9 (x)或φ-V≡中 7(x) 7() 可为张量并、点积或者叉积等。 对于高阶导数,一般有 e()(h1…,hn)=D.…D,4(x)=(h2…、1、n,:;(x)9,③g②9:() 此处(x)∈T/Rm×…xR 为p阶(多重)Rm至T(R")的线性映照。对高阶导数
2928> !HF $K槡9O993% 槡 *3> " P908> !HF $ bàk-ñò76$&28> !HF $"RQF}ç ô63{ExÊ`%ü Qô6#`3{áûs'(R)T57}çô63{"cFd[ÀA`TÙ% BIB ÀnN'$"zY#QZVN'$ 9cû` !HF $"bL9 L9 !9$0)!HF "8> !HF $$ d~[}çv"9 9!9=@$[NI" Q L> 9 L9>!9$!@%"*"@>$ 9!9$ > $(,C% Q*Q,C> 930K 02 !9$)(3M(2 M(K!9$ bL> 9 L9> !9$08> HF E*EHF >B?,[O "8!!HF ! $$v > [!y¥$HF 8!!HF $[NI 786 & ' ( ) !!""
第2期 谢锡麟:基于郭仲衡先生现代张量分析及有限变形理论知识体系的相关研究 343 自然成立 c)(h d"Φ (x)(hc,…,hp),Va∈P2(P-阶置换),亦即具有作用次序无 关性 基于高阶导数,我们有无限小增量公式:如果Φ(x)在点具有直至p阶导数,则有 (0+h)=φ(x) 1dφ (|h|-) 3张量值映照微分学 3.1张量值映照的一阶导数(可微性)以及高阶导数 张量值映照的一般形式可为 f():T(Rm)2D3φ卜f(∮)∈T(Rm) 按其可微性的定义,有 f(Φ+U)=f(①)+ (φ)(U)+o(|U|r(m)∈T(Rm) 此处(中)∈L(T(R"),T?(R")。按L(TP(Rm),T(Rm))的结构,有 彐!6∈T"(R"),满足:6()U=( 亦即张量映照的导数可由唯一确定的高阶张量(q+p阶)以及p点积表示,此处q,p分别为因变量及自 变量的张量阶数。下文仍以(④)∈TP(Rm)表示r阶导数 对于高阶导数,一般有 db()(Un…,,)=Dn…D,f() d了f (U1⑧…U,)∈T(R 实际中,可能会遇见如下形式的张量映照 T‘(Rm)3卜f(φ)。-9()∈T(R"),此处∫(Φ)∈T”(R"),g(Φ)∈T(Rm) 则有()一9((CD)-[:0(0(:10-0()+(),[出(0(.)]此式左方表示对 的全导数。对此,我们只有微分的表达式,而无法获得导数的表达式[。 另,我们也需研究多变量张量值映照,如 f(,中,日):T"(R)×T(R")xT(Rm)2{,业,已}}f(,y,日)∈T(R") 则有 r(y,(P,Q,R)=5(a,y,6)(.)P+a(0,,)(.)+50,y,)(.)R 此式左方表示对自变量全体{φ,y,e}的微分 3.2张量值映照的隐映照定理 张量值映照的隐映照定理,可以表述为
¢£ ¿&L> 9 L9>!9$!@%"*"@>$ 9 L9>!9$!@6!%$"*"@6!>$$"P60Q>!>B[&$"GEu5MS= }`% 67>[NI"=¬+6çÚx&ZÙ 9!9$o9$ ¦Ee > [NI"} 9!9$=@$ HF $ D Àn0x®Y9` DIA Àn0x®*N'$"zY#QZVN'$ }ç;˲ !HF $IT9 J924?!9$08X !HF $ 9Úcû` !HF $$08X !HF $ bL? L9!9$0)!8> !HF $"8X !HF $$%9 )!8> !HF $"8X !HF $$ !HF $"]^&D > !0$4[}ç!X=> [$d> ¦"bX"> s3vaÍç¢ Íç !HF $A [NI% í7>[NI" Q LA ? L9A!9$!4%"*"4A$ ! !0$4% $> !0$4#* > !0$4A ! 0 $!4%M*M4A$08X !HF $ gKq"c\'RrZ[[x !HF $"(!9$08X !HF $ }&?!9$QB(!9$ 9 !4$"D-24?!9">"D$08X !HF $ } ?!9">"D$ 0 !Q"^"G$"D$ A !0$Q="? ">!9">"D$ C !0$^="? "D!9">"D$ ; !0$G bx@Wí¢ÍçòÝ,9">"D-<ûs% DIB Àn0x®*{x®ox }ç;˲<;˲à,"cd?v !## ÐY«&67»gGJ¸<}çsV¬Í[,1^þÝÀ<S}©ª 787
344 力学季刊 第34卷 有映照:f(中,y):T(R")×T(Rm)2D。×Dy3{φ,y}f(,y)∈T(Rm),具有正则性f(φ, y)∈C(D。XDy,T(Rm)。如有:彐(更0,)∈D。Dy,满足 ①f(④0,y0)=0∈T(R");②5(重,y)∈L(T”(R),T(Rm))具有逆算子 则有:彐B(φ0)CD,Bn(0)CDy, 对Ⅴφ∈B,(中),彐!y。∈Bn(),满足f(φ,v。)=0 由此,可构造映照(隐映照):η(φ):B2(φ。)彐φ卜η(Φ):=ψ。,满足 ①n(Φ)∈Bn(y0);②f(重,())=0;③7重)∈C(B4(中);T”(R))。 对∫(φ,(φ))=0,按一般意义的链式求导法则,可得隐映照的导数,即有 ()=-(f(φ,()。f(,中)∈L(T(R"),T(R") 进一步,可得 d 2 (Φ,7(φ) ayr (.)]a(,) (3)(,90),[(0(,(0)3(1),3()(:) 其中,由 (重,())。(,(中)=L4,可得 ay (④,7(φ) ay/(④,(④)(a4ov(,7()+9 a (,())。(④) 上述推导中利用了一般赋范线性空间上微分学的一个结论 引理:如有∫(x):X3}f(x)∈Y可微,B(x):X3xB(x)∈L(Y,Z)可微,X,Y和Z为一般赋范 线性空间,则有:X3上Bf(x)∈Z可微,且有 B(x)°f(x) de (a)(h)of(a)+B(a).(a(h)ez 证明:由f(x)和B(x)的可微性以及B(x)∈L(Y,Z),可有 B(+b),(x+h)[B2+:2(()+0([)+:2(b+b B(x),()+[2(b-x2+B)2(x)]+a∈z 式中a(b∈L(Y,2满足n10(h)12=0∈R:0,()∈y满足,i,mkR o, (hly 考虑到 log(h)f(x)lz≤|oB(h)luy.z·lf(x)ly,可有og(h)of(x)=0(h)∈Z; /°(Y)(m)P (a)(h)≤ dB(a)(h) (x)(l) dB I hI X,L(Y 2) L(X. Y [[(=0)∈,由似分析可,即得证
˲&?!9">$&8> !HF $E8X !HF $IT9 ET> J,9">-24?!9">$08A !HF $"E}`?!9" >$0L%!T9 ET> "8A !HF $$%Z&R!9$">$$0T9 ET> "]^ @ ?!9$">$$!9$">$$0)!8X !HF $"8A !HF $$EU4% }&R00!9$$ST9"0+!>$$ST> " íP9000!9$$"R=>9 00+!>$$"]^?!9">9$9"]^ @ )!9$00+!>$$#A?!9")!9$$ $ B% !9")!9$$Q"? "9!9")!9$$0)!8> !HF $"8X !HF $$ +"cF L# ) L9#!9$> !0$4 $ B% !9")!9$$ 9 > ( !0$4)Q"? "9!9")!9$$ B "?!"> $ B% !9")!9$$Q "# ? "9#!9")!9$$= "# ? ">"9!9")!9$$QL) L9 ! !9$$> ( !0$4) Úq"X "?!"> $ B% !9")!9$$Q"? ">!9")!9$$ $ B% !9")!9$$ 9 > !0$4 $ B% !9")!9$$Q "# ? "9">!9")!9$$="# ? ">#!9")!9$$QL) "9 ! !9$$> ( !0$4)Q "?!">$ B% !9")!9$$ á?QNq/5R Qô6#`3{áûs'K0- xq,0!@$0)!Z"-$]^ 9AS@4$0S 2,0!@$2)!Z"-$ 2@2S K<,!@$%F% 788 & ' ( ) !!""
第2期 谢锡麟:基于郭仲衡先生现代张量分析及有限变形理论知识体系的相关研究 345 作为隐映照定理的一则应用,考虑如下带有约束的张量函数的极值问题: 目标函数:0(④,业)∈R,式中中∈TP(R),业∈T(Rm) 约束映照:∑:={(,y)∈T(R")×T(R")f(中,y)=0∈T(Rm)} 按隐映照定理认识约束,则局部有:f(φ,7(φ))=0。则有,在局部意义下,约束上的目标函数θ(Φ,业) 等价于:θ(φ):=θ(φ,(Φ)),故可计算获得其临界点控制方程 (④,7()+(④,(φ)·(④) Φ,()-(,(更) (φ,())°(重,7(更)) Φ,()-(,(更) (,7()) 更,(更)=0∈T(Rm) 籍此,可构造 Lagrange函数 L(,y;A):=0(,y)+A⊙(,业)=0(更,y)+A”'f;(更,y)∈R 式中A∈T(R")为 Lagrange乘子,⊙表示全点积。可计算得:L(φ,y;λ)的临界点方程相同于上述基 于隐映照定理所确定的目标函数的临界值。 3.3张量值映照的逆映照定理 张量值映照的逆映照定理(局部微分同胚存在性定理),可以表述为 对任意f(中)∈C(T(Rm)2D;T(R")),如有 3∈D满足:(。)∈L(T"(R"),T(R")可逆 则有,局部存在微分同胚∫(更)∈C(B2(更);f(B4()) 对于导数计算,基于:f°f(④)=∈T(Rm),有 dv((p))=/f (φ)∈L(T(R"),T(Rm)) 进-步有F(0)(.)v=(:)(()()()(?) V|,V∈T"(R") 此处(2)()(,)=-(:5)(0)(0(,)1(:)(0),vD∈r(R”) 4张量场多点表示形式 4.1张量场多点表示形式 郭仲衡著《非线性弹性理论》,由于对初始物理构形和当前物理构形都引人曲线坐标系(二者相互独 立),故变形梯度可具有形式F:=(,t)9,(x)②G4(5)。此处构成简单张量的向量分别隶属不同 的局部基{GA()}m=1和{g1(x)}=1(此处仅列出协变基),故将此种张量称为“二点张量”。二点张量的事 例另有转移张量l:=94g⑧G4,此处g4:=(g,GA)am。 考虑到R"中任意二个基之间的相互确定关系,“两点张量”自然可以转换为“一点张量”,由此本文称 张量的“多点表示形式”:对于任意φ∈T(R"),最多可有p点表示形式,亦即其p阶简单张量的构成来 源于p个不同的基。值得指出,“多点表示”仅是同一个张量的不同表示形式,而非不同表示对应不同的 张量。例如,转移张量实际就是度量张量或者单位仿射量:
uv;˲à,$0H"xq 908> !HF $">08X !HF $ ¶¦Ë²&$&$08> !HF $E8X !HF $2?!9">$$ {Ñ7&a%!9$&!9")!9$$QL) L9!9$ !9")!9$$Q "?!"> $ B% !9")!9$$Q"? "9!9")!9$$ !9")!9$$X !0$ "?!"> $ B% !9")!9$$ A !0$"? "9!9")!9$$ !HF $ fb"cYZ2GB@G=B>\I )!9">#0$&$=0O!9">$$=03%*3A?3%*3A !9">$0H xq008A !HF $v2GB@G=B>,4"Oò¦%cêF&)!9">#0$ !HF $IT9#8X !HF $$"Z R9$0T9"]^&L? L9!9$$0)!8> !HF $"8X !HF $$cU }"@A/oûsÊ_?!9$0L%!00!9$$#?!00!9$$$$% í7NIê"67&?B% Q?!9$ !HF $" L?B% L> !?!9$$ !HF $$% +&L# ?B% L># !?!9$$X !0$J !0$ L?!L9$ B% !9$X ( !0$J)"PJ08X !HF $ b&L?!L9$ B% !9$ 9 > !0$4 ( !0$4)Q L?!L9$ B% !9$"P408> !HF $% E ÀnVKqt EIA ÀnVKqt »gGº:×#`ç`,1;"X7í\]Ü,Y["ÐÜ,Y[9:Å#xLÀ!J«SZ ¿$"2Í[ßcE[x %& !HF $"yc> ¦[x"GÚ> [¬v}ç l-Ê<6%;F"1y¦2ÛÊ l}ç<-Ê[x"×-ÊíT-Ê< }ç%Z"K*}çgK0Ûç}çÓ«vwóNç& !## ÐY«&67»gGJ¸<}çsV¬Í[,1^þÝÀ<S}©ª 789
46 力学季刊 第34卷 I=,=9⑧g=8kGA⑧GB=GAG4②GB=g,②G4 郭仲衡所著《非线性弹性力学》出现上述关系式。综上所述,本文谨认为,称“张量的多点表示形式”较“多 点张量”可能更为适宜 4.2张量场多点表示形式下的导数计算 连续介质力学研究中,我们可能同时使用初始物理构形以及当前物理构形中的局部基,由此将自然使 用张量场的二点表示形式。 为说明一般情况,我们引入三组基:{GA(4)}A=1,{9(x)}=1,{h。(y)}m=1分别由曲线坐标X= X(),X=X(x),X=X(y)所诱导;基于微分同胚的传递性,曲线坐标{4}A=1,{x2}=1,{y°}m=1之间彼 此为微分同胚。不失一般性,我们考虑四阶张量场的三点表示形式 (,,)=:A(,,y)91(x)⑧h()⑧GA()⑧G()∈T(R”) 对此,如考虑Φ(,x,y)对的偏导数,则有 g(,x,y)g(x)⑧h(y)⑧GA(5)⑧GB( 此处:VAB ①(,,y)+P:;B一F 如考虑Φ():=中(),(),y())对的全导数则有 x)⑧h(y)⑧GA()②GB( 此处:□ (,,y)=V1重:,B(,汇,y) ()·V匝:(,,y) 8)·V 式中:V④:B (,x,y)+T 5非完整基理论 5.1单点形式对应的非完整基理论 具有微分同胚的曲线坐标系可另称为完整系,完整系的局部基称为完整基;而称不具微分同胚的基为 非完整基。一般而言,我们可构造正交的完整基而不易使其为单位正交基。然而,相对单位正交基展开的 张量场可保证其所有分量都保留其原来的量纲。场论中,我们需要获得基于完整系定义的张量场场论微 分算子在非完整的单位正交基下的表示。就此,郭仲衡著《张量(理论与应用)》叙述有非完整系理论(本文 称非完整基理论),对此我们也有思想及方法上的澄清 5.2多点形式对应的非完整基理论 类比于一般非完整系理论,可以将相关思想及方法推广至张量的多点表示形式。设对完整基 {GA()}A=1,{9;(x)}=1,{h。(y)}m1,分别定义对应的非完整基{GA()}A=1,{g;(x)}m=1 {h。(y))m=1。相应的基转换系数记为GA=:CG(B,g=:Cg,h。=:Ch(;转换系数对应有其 对偶,如G(B=:C(BGD,满足 CC=8 具体分析上,首先按完整系定义张量梯度;其次按坐标变化规则确定非完整基下的分量;再次推导得 非完整基下的协变导数计算式。我们获得如下结论
*<33 2(3M(2 <(32( 3 M(2 <3" G$" M$0 <$"0$" M$0 <( 3 "(3M$" »gGDº:×#`ç`&';<á?}Àx%ªáD?"i$býv"r1}ç<y¦[x2W1y ¦}ç2c\¤v»% EIB ÀnVKqtê*'$H +aj&'©ªq"c\ÊzQ5\]Ü,Y[dÐÜ,Y[q<@A6"Xbü¢£Q 5}çå<J¦[x% v²; Q;Î"9:~k6&,$" !,$-F "<%",(3!9$-F 3<%",@5 !:$-F 5<% s3XÅ#xL S< S!,$"S<S!9$"S<S!:$DpN#67ûsÊ_<Mw`"Å#xL,," -F "<%",93 -F 3<%",:5 -F 5<%~{j bvûsÊ_%-P Q`"XYÂ[}çå<~¦[x 9!,"9":$<930"0 0/00!,"9":$(3!9$M@/ !:$M$"!,$M$0 !,$08"!HF $ íb"ZXY 9!,"9":$í,) <úNI"} "9 ",)!,"9":$<,)930"0 0/00!,"9":$(3!9$M@/ !:$M$"!,$M$0 !,$ b&,)930"0 0/00&<"930"0 0/00 ",) !,"9":$=Q" )P%930P0 0200 B#P )0930"0 020P ZXY 9!,$&<9!,$"9!,$":!,$$í,) <òNI"} L9 L,)!,$<"9 ",)!,"9!,$":!,$$="9[ ",)!,$0"9 "9[!,"9!,$":!,$$=":% ",) !,$0"9 ":%!,"9!,$":!,$$ <&T)930" 0/00!,"9";$(3!9$M@/ !:$M$"!,$M$0 !,$ b&T)930"0 0/00!,"9":$<,)930"0 0/00!,"9":$="9[ ",) !,$0,[930"0 0/00!,"9":$=":% ",) !,$0,%930"0 0/00!," 9":$ xq&,[930"0 0/00&<"930"0 0/00 "9[ !,"9":$=#3 [C9C0"0 0/00" ,%930"0 0/00&<"930"0 0/00 ":% !,"9":$B#* %/930"0 0*00 P E}Ë+xQ PIA KtÒþ*E}Ë+xQ EûsÊ_<Å#xLÀc¨rvxÀ"xÀ<@A6rvx6#r-EûsÊ_<6v ×x6% QÂ"cYZÑ<x6-QÚvvwÑ6%£"SívwÑ6^< }çåc¤ÚDsç¤;Ú®C<ç%å1q"öðäF67xÀàk<}çåå1û s4o×x<vwÑ6[<%0b"»gGº:}ç!,1ST5$;Q?×xÀ,1!i$ r×x6,1$"íbèE.WBá<®±(!)% PIB VKtÒþ*E}Ë+xQ îé7 Q × x À , 1"c d ü S } E . W B Q G } ç < y ¦ [ x%É í x 6 ,$"!,$-F "<%",(3 !9$-F 3<%",@5 !:$-F 5<%"s 3 à k í T < × x 6 ,$" Q !,$-F "<%",(3 Q !9$-F 3<%" ,@5 Q !:$-F 5<%%ST<6K&ÀIUv $" <&L!0$ " $ Q !0$"(3<&L!2$ 3 ( Q !2$"@5 <&L!/$ 5 @ Q !/$#K&ÀIíTÚ íx"Z$ Q !0$<< !0$$T"]^ L!0$ " LT !0$<3T "% EÝsVá"IJ9xÀàk}çß#ÚM9xLÍÅú}à×x6[<sç#ÛMQNF ×x6[<vÍNIêx%äFZ[T1% 78# & ' ( ) !!""
