连续介质力学2-2
2-2 转动定律
刚体转动定律
❖ 扳手拧螺丝 ❖ 四两拨千斤
力矩
❖ 力的作用效果不仅取决于力的大小,而且与力的 方向和作用点的位置有关.只有与转轴既不相交、 也不平行的作用力才能使物体转动
且一在刚转力体动矩绕平面O内z ,轴r旋转为由, 点力OF到作力用的在作刚用体点上P点的P径,
矢.
M
F
对转轴Z M
上式两边同时乘以 ri
riFi m iri2
上式两边求和得
riFi ( miri2) J
O ri
r Fi
mi
上式左边等于刚体受到的合力矩 riFi
所以有
J
刚体对定轴的转动惯量乘以角加速度等于刚体所受的外力矩之 和。这就是刚体定轴转动的转动定律。
转动定律:刚体定轴转动的角加速度与它所受的合外力
解 在两个力的作用线中间任意 选取一个转轴O
例2.4
解 在两个力的作用线外部任意 选取一个转轴P 两个力的力矩转向相反
例2.5 总重量为mg 总长度为L 的均匀细棒,绕端点O 转 动.当棒与水平方向的角度为θ 时,重力的力矩多大?
解 质量微 元的重力
质量微元 的力矩
总力矩
例题2.6 轻绳跨过定滑轮,其两端分别悬挂着质量为 m1和m2 的物体,且m1 < m2 .滑轮的转动惯量为J , 半径为r,绳与滑轮之间无相对滑动.试求物体的加速 度a和绳的张力
矩成正比,与刚体的转动惯量成反比 M J J d
dt
➢ 转动惯量物理意义:转动惯性的量度.
➢ 质量连续分布刚体的转动惯量
J miri2 r2dm 质量元: dm
注意
i
转动惯量的大小取决于刚体的密度、几何
形状及转轴的位置.
连续介质力学
矢量与张量
连续介质力学基础
克罗内克符号(Kronecker delta)
ij
1 0
i i
j j
1mam 1a 1 11a 2 21a 3 3 a 1
a a 2mam 2a 12a 2 22a 3 a 2 im m
i
3mam3a 13a 2 23a 3 3 a 3
如果e 1 ,e 2,e 3 是相互正交的单位矢量,则有 ei e j ij
变到 x1, x2 , x3 时,它们又是如何变换的。
如果变量系在变量xi中只有一个分量Φ,在变量 x 中i 只有一
个分量 ,并且在对应点,Φ 和 相等,则称为数量场。
(x1,x2,x3) (x1,x2,x3)
矢量与张量
连续介质力学基础
数量、向量和张量的解析定义
如果变量系在变量xi中有三个分量 i ,在变量x i 中有三个分
/
ei )
xi/ ijxj
矢量与张量
连续介质力学基础
一般坐标变换
一组独立的变量x1, x2, x3 可以一点在某一参考标架中的坐标。
通过方程 xi fi(x1,x2,x3) 把变量x1, x2, x3 变成一组新的变
量 x1, x2 , x3 这就规定了一个坐标变换。
逆变换 xi gi(x1,x2,x3)
绪论
连续介质力学基础
连续介质力学中的“基元”
时空系:
时间和空间是运动物体的客观存在形式。空间 表示物体的形状、大小和相互位置关系;时间 表示物体运动过程的顺序。
为描述物体的运动,需要在时间和空间中选取一特 定的标架,作为描述物体运动的的基准,这种标架 称为时空系。
绪论
连续介质力学基础
连续介质力学中的“基元”
连续介质力学-2d
度
23
量
H
=
1 2
ln(I
+
2E)
=
ln(I
+
E工 )
=
1 2
(2E
−
1 2
(2E)2
+
1 (2E)3 3
− ...)
=
E工
−
1 2
E工2
+
1 3
E工3
−
...
Green应变,工程应变和对 数应变张量之间的关系
18
第 二
同理, 对参考构形中随体坐标系{ xi, t0 }变形状态的描
章 述还可定义工程应变张量 e工和对数应变张量 h :
(C AB
− GAB )
(2.69)
应
称为Green应变张量
变
度 量
Green应变张量 E 是对当前构形中随体坐标系
{ XA, t }的变形状态描述,称为空间应变张量。
—
E的主方向:Lα ,α = 1,2,3
E的主值:
1 2
( λ α2
− 1), α
= 1,2,3
14
第 二 章
任取变形前在参考构形R中的单位矢量 L ,记
L ⋅ E ⋅ L 为Green应变张量在 L方向的法分量:
变 形 和 运
∴L ⋅ E ⋅ L = 1 (L ⋅C ⋅ L − L ⋅ I ⋅ L) 2
=
1 2
(λ2L
−1)
=
ds2 − dS 2dS 2
2
(2.70)
动
2.考虑{ xi, t0 }
{xi} 线元长度的改变(E描述):
—
应 ds2 − dS 2 = dx ⋅ i ⋅ dx − dx ⋅ c ⋅ dx = dx ⋅ (i − c) ⋅ dx
连续介质力学
连续介质力学的应用领域包括:工 程力学、流体力学、固体力学、生 物力学等。
连续性假设:假设介质是连续的没 有空隙或裂缝
各向同性假设:假设介质在各个方 向上都是相同的
添加标题
添加标题
添加标题
添加标题
均匀性假设:假设介质在各个方向 上都是均匀的
小变形假设:假设介质的变形很小 不会影响其物理性质
流体:不可压缩、连续、无固定形状的 物质如空气、水等
多尺度连续介质力学:研究不同尺度下的连续介质力学问题如分子动力学、介观力学等
跨学科连续介质力学:与其他学科交叉如生物力学、环境力学等
计算连续介质力学:发展高效的计算方法和软件解决复杂问题如流体动力学、固体力学 等
PRT SIX
连续介质力学是研究流体和固体力学 的重要学科
连续介质力学的特点包括:连续性、 守恒性、对称性等
研究方法:数学模型、数值 模拟、实验验证等
研究对象:连续介质如液体、 气体、固体等
基本概念:应力、应变、位 移、速度、加速度等
应用领域:工程力学、流体 力学、固体力学等
PRT THREE
弹性力学的定义:研究弹性体在外力作用下的变形和应力分布的学科 弹性力学的基本假设:连续性假设、小变形假设、均匀性假设、各向同性假设 弹性力学的基本方程:平衡方程、几何方程、物理方程 弹性力学的应用:工程结构设计、地震工程、材料科学等
,
汇报人:
CONTENTS
PRT ONE
PRT TWO
连续介质力学是研究连续介质(如 液体、气体、固体等)在力作用下 的变形、流动和应力分布的学科。
连续介质力学的研究内容包括:应 力、应变、变形、流动、热传导等。
添加标题
添加标题
添加标题
连续介质力学第二章.
即得( i ),将( i )作相应的指标替换, 展开化简,将得其余三式。
二维置换符号 e (, 1, 2)
从三维退化得到
e ei j3 e 3
其中
e11 e22 0, e12 e21 1
有下列恒等式
e e
又如,方程
12
2 2
32
111
2 22
333
用指标法表示,可写成
i i i ii i ii i ii
i 不参与求和,只在数值上等于 i
1.2 Kronecker 符号
在卡氏直角坐标系下,Kronecker 符号定义为:
ij
表示
e1 A11e1 A12e2 A13e3 e2 A21e1 A22e2 A23e3 e3 A31e1 A32e2 A33e3
ei Aije j i 为自由指标,j 为哑标
表示
e1 A11e1 A12e2 A13e3 e2 A21e1 A22e2 A23e3 e3 A31e1 A32e2 A33e3
新旧基矢量夹角的方向余弦:
ei e j | ei || e j | cos(ei , e j ) cos(ei , e j ) ij
1.5.1 坐标系的变换关系
ij cos(ei , e j ) ei e j
旧 新
e1
e 2 e 3
e1
11 21 31
ai xi a1x1 a2 x2 a3x3 bjj b11 b22 b33
cmem c1e1 c2e2 c3e3
双重求和
33
S
力学讲义第六章连续介质力学
第六章 连续介质力学连续介质模型:物质(气,液,固)连续地分布在它们所占有的区域内连续介质质元: 宏观小, 微观大物质讨论宏观力: 包括外力以及外力作用下形变or 运动引起内部的弹性恢复力 讨论内力的一般方法:假想将其切开,切下部分的作用由内力代表;由平衡条件求力.例: (不计重力)连续介质是比质点、刚体更普遍的经典力学模型,应用也最普遍。
物理状态量在连续介质模型下成为点函数. 不计微观内力 §6.1 应力和应变6.1.1 应力固体为例截面π , 方位 n ; P 处邻域 ∆S 上 张力∆TP 处应力σ = lim ∆∆ TS = d T /dS =σ(P, n ) =σt +σn正应力(法向应力, 张力) σn 单位:P a (压强)(>0为拉应力 ; <0为压应力) 剪应力 (or 切应力) σt应力状态:对同一点P 处,方位不同的截面上应力σ不同。
函数关系σ=σP ( n)叫P 处的应力状态. 由平衡方程可以证明,互相垂直的三个截面上的6个应力(正,切应力)就可以完全决定一点处的应力状态 (由此6个应力可以计算出该处任意方位截面上的应力)应力主面: 该面上只有正应力, 称为主应力. 一点处必有三个互相垂直的应力主面6.1.2 应变固体有两种基本的应变形式:线(拉,压)应变 ;剪应变1. 线应变 ε均匀形变 : 长度l , 总形变∆l (截面法向x ) 则 εx = ∆l / l形变不均匀:一点处位移uAB 段形变=∆u x =u x (x+∆x) -u x (x)=∂∂u xx∆x A 处x 方向线应变εx = lim (∆u x /∆x) = ∂u x / ∂x类似: y 方向线应变 εz =∂u y / ∂y z 方向线应变 εz =∂u z / ∂z 一般情况下应变也是点函数, 不均匀形变时各处应变也不相同.应变是位移的空间变化率(位移的偏导数)2. 剪应变以xy 平面为例, 矩形 → 菱形定义:A 点剪应变(xy 平面上,小变形)为 εt = lim (δ1+δ2)= ∂u x /∂x + ∂u y /∂y δ1 ≈tan δ1=B’B’’/A’B’’=[u y (x+∆x) -u y (x)]/∆x → ∂u y /∂x 类似, 当 ∆x →0 , ∆y →0时 , δ2 → ∂u x /∂y3. 体应变均匀形变时, 体应变 εV = 体积增量/体积 =∆V / V不均匀形变时, 讨论一点处体应变一点附近小长方体(∆x,∆y,∆z) 小形变后为[(1+εx )∆x ,(1+εy )∆y, (1+εz )∆z] V=∆x ∆y ∆z ∆V ≈(εx +εy +εz )∆x ∆y ∆z 小变形 εV =εx +εy +εz 剪应变引起的体应变为高阶小量.自然状态无内力内力与外力平衡F F 内∆S →0 ∆x →0∆x →0∆y →0 y+∆侧平面)∆ll x∆x)6.1.3 胡克定律——应力和应变的关系 1678年胡克提出单向拉伸时 ε ∝ σ , 后来推广到三维 (实验定律) 1. 单一正应力引起的线应变 σx 引起 纵向线应变 εx = σx /Y 横向线应变εy =εz = -μεx = -μσx /Y Y —杨氏模量(压强量纲)μ ——泊松比(无量纲) 0≤ μ ≤ 0.5 σy , σz 的贡献类似 2. 总线应变与正应力的关系——广义胡克定律(在一定的形变范围内—比例极限) εx =1Y [σx -μ(σy +σz )] εy =1Y [σy -μ(σx +σz )] εz =1Y [σz -μ(σx +σy )] 3. 体应变与正应力εV =εx +εy +εz =(1-2μ)(εx +εy +εz )/Y ≡ σ0/K σ0≡(σx +σy +σz )/3 K=Y/[3(1-2μ)] K —体弹性模量 由4. 剪应变与剪应力εt =σt /G G —剪切弹性模量5. 各向同性固体只有两个独立的弹性模量, Y 、G 、K 、μ中只有两个独立K= Y / [3(1-2μ)] G=Y /2(1+μ) < Y一般 μ ≈ 0.35 G 、K 、Y 的量级为1010 —1011 P a , 差别不太大部分材料的弹性模量材料 铝 铜 金 电解铁 铅 铂 银 熔融石英 聚苯乙烯 K 7.8 16.1 16.9 16.7 3.6 14.2 10.4 3.7 0.41 G 2.5 4.6 2.85 8.2 0.54 6.4 2.7 3.12 0.133 Y 6.8 12.6 8.1 21 1.51 16.8 7.5 7.3 0.36 μ 0.355 0.37 0.42 0.29 0.43 0.30 0.38 0.17 0.353 说明: K 、G 、Y 的单位 为1010P a补充题4. 矩形截面杆在轴向拉应力σz =2.0⨯105 P a作用下变形,已知Y=19.6⨯1010 P a , μ=0.3 .求:εV 补充题5. 矩形悬臂梁的一端有作用力P.已知l =2 m, h=20cm,梁宽b=5 cm ,P=1000kg 力, 求梁内最大正应力§6.2 固体拉伸.弯曲.扭转讨论三种情况下的应力状态,计算应力与应变 6.2.1等截面直杆的拉压 圆形截面直杆;两端均匀压强p (拉>0;压<0)横截面 σz =p σt =0 应力状态: 与z 轴互垂两面上 σR =σφ=0 ——单向应力状态 ∴ σz =p= Y εz = Y ∆l / l 均匀形变 弹性形变势能: E P = ⎰ F 外du = ⎰0∆lSY u ldu=YS ∆l 2 / 2l u 为z 方向位移, S 为横截面积(近似不变) 弹性形变势能密度 e P =E P /V=12Y εz 2 =12σz εz (也适于不均匀形变) 说明:其他均匀截面直杆σR ≈0 σφ≈0 可以近似按圆杆处理6.2.2 矩形梁纯弯曲矩形梁(高h,宽b) 力偶矩M纵向画线弯曲:上短—压; 中不变—中性面; 下长—拉横截面上 σx , σt =0应力状态: σy =σz =0——单向应力状态M ⇒ 应力σx , 形变θ0P 处:εx= lim (PP’-oo’)/oo’= lim[(ρ+y)∆θ-ρ ∆θ]/ρ ∆θ=y/ρ σx =Y εx =Yy / ρ ∝ y 下面求ρ 横截面上:∑F =0 (∴中性面正在中点)∆θ→0 ∆θ→0 p z φM 内= ⎰y σx dS = Y ⎰ y 2 dS /ρ ≡YρI z =(应该)= M ——柏努力. 欧勒定律∴ Y/ρ = M/I z σx =M I z y σx max =M I z 2h ρ=YI z /M θ0 = l /ρ(θ0 为转角,代表形变;l 为中性面的长度) 定义对z 轴惯性矩 I z ≡ ⎰y 2 dS 对矩形截面 I z =2b ⎰02h /y 2dy =112bh 3 为节约材料:h ↑ , b ↓ ; 减少中性层还有鸟骨、麦杆…说明:(1)其他形状截面的梁在力偶矩作用下弯曲时,σy ≠ 0 σz ≠0, 非单向应力状态,但σy ≈0 σz ≈0 ,与单向应力状态偏差不大,可以近似按单向应力状态计算(2)非力偶矩作用时,一般可以忽略剪应力,近似按纯弯曲处理:(不计重力) 悬臂梁M 内=M(x)=P(l -x)简支梁 x ∈(0,l /2) M 内=M(x)= P x/2仍有: σx (x)=M(x) y/I z ρ(x) =YI z / M(x) 注意:σx (x),ρ(x),M(x)不再是常数 (3)仍有:e P =12Y εz 2 =12σz εz6.2.3 圆柱扭转表面画上圆周和母线圆周线不变, 横截面保持平面——横截面上 σtR =0应力状态: 横截面上 σt =σt φ σz =0 (只有M) σR =σφ=0 横截面上形变:圆周处εt (R)=R φ /h r 处εt (r)=r φ /h ∴ σt (r)=Gr φ /h ∝ r下面求φ M 内= ⎰ σt r dS = ⎰0R σt r 2πrdr=12h πGR 4φ ≡D φ =(应该)=M ∴G φ/h=2M/(πR 4) σt (r)= G φr/h M=D φ ∴ σt (r)=24M R πr σt max (r)=2M /πR 3 φ=M/D 扭转弹性系数 D=πGR 4/2h (悬丝扭矩 M=D φ D ∝ R 4/h ) 扭转弹性势能E P = ⎰0φM d φ=D φ2 /2 可证e P =12G εt 2 =12σt εt6.2.4 允许应力.强度计算1. 只有正应力or 剪应力材料极限应力(正or 剪)σj , 许可应力[σ]=σj /K 安全系数=1.4—3.0 — 14材料 屈服极限σs 强度极限σb 许可应力 [σ] (kg/cm 2)A 3 2200—2400 3800—4700 1700 16Mn 2900—3500 4800—5200 2300 300#水泥 拉21,压210 拉6,压105 红松(顺纹) 拉981,压328 拉65, 压100 注:A 3—普通低碳钢 16 Mn —低合金钢 常温、静态、一般工作条件材料中最大应力(正or 剪) 应满足 σmax ≤ [σ] 2. 复杂应力情况——按相应的强度理论计算§6.3 流体静力学——流体力平衡下内应力的分布 流体:液,气; 具流动性; 主要讨论液体; 设: 连续、均匀6.3.1 静止流体内应力δσt1. 一点处应力状态σt≡0 只有正应力σ , 且正应力大小与截面无关σ( n)≡σ证: 因为可流动流体静摩擦力=0 ∴σt≡0如图四面体受力平衡设S面上正应力为σ ,x向Sσ⋅x -σx S x=0σ=σ n S=S n S x=S ⋅ x∴σx S x=Sσ⋅x =σS⋅x= σS xσx=σ类似σy=σ=σzx,y,z任选, ∴任意截面上的正应力的大小皆为σ由四面体受力平衡, 从三个坐标平面的应力⇒任意截面S上的应力. 注意:忽略了体积力2. 流体内压强定义:流体内压强为P= -σ(流体中一般没有拉应力,∴σ<0 P>0)说明:(1)压强为标量,严格定义P= -σ0 = (σx+σy+σz) /3(2) 由一点处应力状态, σ与方位无关∴P与方位无关(3) 从证明知,关键σt=0 . 所以对理想流体(无内摩擦)在流动(包括加速流动)时结论也对(4)对粘滞性流体流动时有剪应力,各截面σ不相同.但若σt较小可以忽略,各截面正应力近似相等为σ , P ≈-σ(5) 流体中负压强(拉应力).特定条件(稳定,缓慢过程)下,流体中可出现负压. 水的负压可以达到300atm6.3.2 静止流体平衡方程——临近点处压强关系取小段柱状流体f—单位质量..上的体积外力x向: [P(x) - P(x+∆x)] ∆S + ρ∆S ∆x f x =0∴∂P /∂x = ρf x类似: ∂P /∂y = ρf y ∂P /∂z = ρf z合起来:∇P = (∂P/∂x) x +(∂P/∂y) y +(∂P/∂z) z = ρf 6.3.3 重力场中静流体1. 流体中压强随高度分布小范围g为常矢量f = (∆m g) /∆m =g = g y ∂P/∂x =∂P/∂z = 0 ⇒P与x,z无关, 在同一高度上P相等∂P/∂y = ρg若ρ为常数(液体or高度差不大的气体)积分得:P(y)=P0+ρgy P0=P(0)不同密度液体(鸡尾酒)的稳定分界面为水平面2. 帕斯卡定律定律:加在密闭液体中的压强等值地传到液体中各处以及壁上.解释: 设压强加在o处,使P0等值地改变,但ρgy 保持不变,所以P(y)随P0同样增加.3. 阿基米德定律定律:浸在流体中物体所受浮力等于物体排开的流体的重量证明:设物体外表面为S .流体对物体作用通过压强体现.∴浮力=⎰-Pd S保持S不变,则浮力不变. 将物体换成流体,该流体应处于平衡,即外界对S的压力之和等于流体重量:⎰-Pd S +m g =0∴浮力= -m g 浮力作用点即该流体重心(一般情况下不是物体的重心)附: 等温理想气体压强随高度的分布已知其密度ρ=cP (c为常数)解: dP/dy = -ρg = -cgP ⎰PPdPP= ⎰y-cg dy 得:P(y)=P0e-cgy又例: 以ω匀速转动的水平试管,内部充满流体. 以试管为参考系, 则惯性离心力为体积力,产生径向压强差.§6.4 流体的定常流动6.4.1 描述流体运动的两种方法1. 两种方法拉格郎日法: 认准各个质元,分别描述其运动状态(r i,v i,a i)及其变化规律r i,v i,a i只是t的函数, v=d r/dt , a=d v/dt ; 应用牛顿定律必须用拉格郎日法. 困难:如何认准?如何跟踪?描述不便欧拉法: 讨论流体场(流体性质场)的场分布∆x)主要是流速场v=v(r,t) . 还有a=a(r,t)P=P(r,t) 压强场……2. 欧拉法中质元的加速度质元加速度a = d v/dt (速度全导数or实质导数)是对一个确定质元速度v(即拉格郎日法中的速度v)的导数.流速场v(r,t)在地点不变下对t的偏导数∂v/∂t ≠a (流速场中同一地点不同时刻的v是不同质点的速度)认准m i :a=d v(x,y,z,t)/dt=∂v/∂t+[∂∂vxdx +∂∂vydy+∂∂vzdz]/dt=∂∂vt+v x∂∂vx+v y∂∂vy+v z∂∂vz=∂∂vt+ v ⋅∇v3. 流体流动的图象表示拉格郎日法: 流体质元的实际运动轨迹——迹线流管——流线围成的细管;流束——流管中流体6.4.2定常流动: v与t无关,v=v(r) ;不定常流动: v与t有关定常流动特点:∂v/∂t =0 a = v⋅∇v≠ 0流线不变,与迹线重和∴迹线也不变P,ρ与t无关是否为定常流与参考系有关设迹线如图. V1,2,3为t1,2,3时刻同一质点的速度.若v与t无关,则v也是速度场中1,2,3点的速度,迹线也是流线. 迹线不变则场中质元数不变,∴ρ不变圆柱在理想流体在匀速直线运动. 在静系中流体为非定常流动,在圆柱参考系中为定常流动§6.6 粘滞流体的流体长时间、长距离、相对速度很大时,粘滞性不可忽略主要讨论层流. 层流:流体分层流动,彼此不混淆流体粘滞性的体现:固、液相对运动时出现摩擦力;液体内部流速不同,各层之间出现摩擦力6.6.1流体的粘滞性板A匀速直线运动引起层流,各层之间粘滞力fz层假想剖面∆S, 两侧粘滞力∆f牛顿摩擦定律:(实验定律) ∆f ∝ (dv/dz) ∆S 即∆f = ηdvdz∆Sdv/dz : z方向速度(空间)变化率(速度梯度)η: 粘滞系数(黏度)温度T↑⇒η↓ (液体) η↑(气体)(f本质: 液体主要来自层之间分子力;气体是通过该层交换宏观定向动量)[η]=ML-1T -1SI(MKS)制为Pa ⋅s CGS制为“泊”1泊=0.1 Pa⋅s η/ρ——运动黏度(比黏度)满足牛顿摩擦定律的流体——牛顿流体(否则叫非牛顿流体—少数如血液)6.6.2 粘滞流体的运动规律1. 动力学方程(介绍) 纳维—斯托克斯方程(Nevier,M. , Stokes,G.G.)-∇P+ρf+η∇2 v = ρ (d v/dt)2. 修改后的伯努力方程定常流动,不可压缩,沿流管(有粘滞性) 由功能原理dW粘1→2 +(P1-P2)dV = dE= (dm v22/2+dm gz2)-(dm v12/2+dm gz1)dm=ρdV∴ P1+ρv12/2+ρgz1=P2+ρv22/2+ρgz2 +w12——修正后的伯努力方程∆t)∆t)m i运动轨迹m质点t2t时刻:3流线w 12 = -w 粘1→2 = dW 粘1→2 /dV >0 为单位体积..流体克服..粘滞阻力做的功水平均匀细管中: v,z 相同, P 1 -P 2=w 12=P 2 -P 3=…=P 0’-P 1=ρg(H 1-H 2)=…=ρg ∆H=ρg(H 0’-H 1) ∴P 0’-P B =P 0’-P 0=ρgH 0’=w 细管 将液面A 与出口B 联系:P 0+ρgH 0+0=P 0+0+ρv 2/2+w 细管+w 粗管∴ρv 2/2=ρg(H 0-H 0’) -w 粗管=ρgh 0-w 粗管≈ρgh 0 v ≈(2gh 0)1/2w 细管, w 粗管分别是单位体积流体在细管和粗管中流动克服阻力做的功∴粘滞流体水平均匀流动必有压强差——流水水面不水平 , 熔岩流动高度差很大3. 哈根—泊肃叶(Hagen,G. , Poiseuille, J.L.M.)方程——水平圆管层流哈—泊定律由哈根1839年实验证实, 后为泊肃叶1842年独立发现水平圆管, 定常流动柱坐标(r,φ,z)v z 与r,φ无关v =v z (r)z d v /dt=0忽略体积力f =0 , 流线平行直线, ∴同一横截面上P 相同对小圆柱, 1、2两横截面上对应处速度相同 ∴合外力为零 即 (P 1-P 2)πr 2 + ηdv drz⋅2πr l =0 (f 粘为-z 方向, dv z /dr<0 ∴取 “+”)⎰0v r z ()dv z = ⎰R r -12ηl(P 1-P 2)r drv z (r)= (P 1-P 2)(R 2 -r 2) / (4ηl ) Q V = ⎰ v ⋅ d S = ⎰0Rv z 2πr dr = π(P 1 -P 2)R 4 / (8ηl ) ——哈—泊公式由此可以讨论石油、天然气、水输送问题(管径、压差与流量);隧道、河流的流量…平均流速 v =Q V /S= (P 1 -P 2)R 2 / (8ηl ) P 1 -P 2=8ηv l R -2 ∝ l R -2,l光滑金属管光滑同心环缝滑阀口Re C2000—2300 1100 260例. 日常生活. 水管d=0.025m Re C =2000 1atm 20︒C时η=1.0⨯10 -3Pa⋅ s 则临界水流速v C = ηRe C /ρd = 0.079 m/s∴一般管流为湍流。
3.连续介质力学
加速度
vX, t 2 uX, t aX, t v t t 2
2 变形和运动
运动描述:
独立变量是空间坐标 x 和时间 t,称为空间或Eulerian描述
v v( x, t ) vΦX, t , t
通过链规则得到材料时间导数 (全导数)
空间时间导数
2 变形和运动
运动条件
连续可微,一对一(F可逆),J > 0 第一个条件,变形梯度通常在材料的界面上是非连续的。在 某些现象中,例如扩展裂纹,运动本身也是非连续的。要求在运 动及其导数中非连续的数量是有限的。实际发现,有些非线性解 答可能拥有无限数量的非连续。然而,这些解答非常罕见,不能 被有限元有效地处理,所以不关注这些解答。 第二个条件,即运动为一对一的,要求在参考构形上的每一 点,在当前构形上有唯一的点与之对应,反之亦然。这是 F 规则 的必要充分条件,即F是可逆的。当变形梯度F是正常的, 则 J 0 ,因为当且仅当 J 0 时F的逆才存在。因此,第二个和 第三个条件是有联系的,后者更强。 第三个条件,更强的条件, J 必须为正而不仅非零,在第 3.5.4节看到这遵循了质量守恒。这个条件在零尺度集合上也可能 违背,例如在一个裂纹的表面上,每一个点都成为了两个点。
如何建立x与X的关系? 需要借助单元坐标和初始构形
2 变形和运动
解:三角形3节点线性位移单元的构形(见附录3)
1 2 3 1
xξ , t xI t I x1 t 1 x2 t 2 x3 t 3
运动条件
除了在有限数量的零度量集合上,假设描述运动和物体变形的映射
ΦX, t
满足以下连续性条件: 连续可微,一对一(F可逆),J > 0
《连续介质力学》课件
动量矩守恒定律
描述物质系统动量矩变化规律的定律。
动量矩守恒定律也是连续介质力学中的基本定律之一。它指出在一个没有外力矩作用的封闭系统中,系统的总动量矩保持不 变。动量矩是系统动量和位置矢量的乘积,因此这个定律说明系统的旋转运动状态只与系统的初始状态有关,而与时间无关 。
能量守恒定律
描述物质系统能量变化规律的定律。
金属材料的疲劳和断裂 研究
01
02
03
复合材料的细观结构和 力学行为分析
04
无损检测和结构健康监 测技术
环境科学
01
土壤和岩石的力学性质研究
02
地质工程和地震工程中的稳定性分析
03
生态系统和自然资源的可持续性发展研究
04
环境流体力学的模拟和分析
06
连续介质力学的未来发展
新材料与新结构的挑战
新材料特性
能量守恒定律是物理学中的基本定律之一,它在连续介质力学中也有重要应用。这个定律指出在一个 封闭系统中,系统的总能量保持不变。能量的形式可以包括动能、势能、内能等,但不论能量的形式 如何转化,总量始终保持不变。
熵增原理
描述系统无序程度变化规律的定律。
熵增原理是热力学中的基本定律之一,它指出在一个 封闭的热力学系统中,系统的熵(表示系统无序程度 的物理量)总是趋向于增加。也就是说,系统总是倾 向于向更加混乱和无序的状态发展,而不是向更加有 序和有组织的状态发展。这个原理在连续介质力学中 也有重要的应用,例如在研究流体和热传导等问题时 需要考虑熵增原理的影响。
THANKS
感谢观看
《连续介质力学》ppt课 件
• 连续介质力学概述 • 连续介质力学的基本概念 • 连续介质力学的物理定律 • 连续介质力学的数学模型 • 连续介质力学的应用领域 • 连续介质力学的未来发展
连续介质力学
x
c
xdm dm
xRd
Rd
2 R cos d
α
O
dθ
dl
Rsin
d R
α
θ
x
§1 刚体的动量和质心运动定理
[例题3] 求质量均匀,半径为R的半球的质心位置. [解]如图 ,质心一定在x轴上。 设半球的密度为,将半球分割成许多厚为dx的圆片
的连线总是平行于它们的初始位置间的
连线.
刚体平动
质点运动
§1 刚体的动量和质心运动定理
(2)刚体的转动(对点,对轴) 刚体中所有的点都绕同一直线做圆周运动,这条直线叫作转轴。 固定转轴:转轴不随时间变化—— 刚体定轴转动 瞬时转轴:转轴随时间变化 —— 一般转动
定轴转动的特点: •各质点都作圆周运动; •各质点圆周运动的平面垂直于轴线,圆心在轴线上; •各质点的矢径在相同的时间内转过的角度相同。
§1 刚体的动量和质心运动定理
1、刚体:在任何情况下,其形状和大小都保持不变的物体。(任意两质
点间距离保持不变的特殊质点组) 对“刚体”概念的说明 •“刚体”和质点一样是一个理想化的力学模型; •任何两点之间的距离在运动过程中保持不变; •刚体可以看成是无数质点组成的质点系 2、刚体的运动类型 (1)刚体的平动: 若刚体中所有点的运动轨迹都保 持完全相同,或者说刚体内任意两点间
mi xi
n
xC 6.8 1012 m
12 rC 6.8 10 mi
y
d
O
H C
52.3
o
o
d
H
x
52.3
o
§1 刚体的动量和质心运动定理
连续介质力学
一、引论连续介质力学研究物体的宏观力学微观粒子性质.在宏观现象中,物体变化的最小特征尺度远大于原子的尺度,虽然物理上物体是物质点的集合,质量连续性假设对物休的宏观力学过程的研究却是合理的,在连续介质力学中可以对物体进行无限的分割,也就是说,可以用场的观点来描述物体的内部变化和作用过程.质量连续性假设要求物体连续地充满它所占据的空间,即可以用三维欧氏空间的一个开集表示物体的客观存在、指示其位置.开集中的一点表征占据该位置点的一个微小介质团,这样的介质团我们称之为物体单元,开集中所有点表征的物体单元组成了物体.若要用严格的数学理性推演连续介质力学,必须知道物体单元在数学上的确切涵义,即要回答: 表征物体单元的点是开集还是闭集?若是闭集,则物体单元表现为数学上离散的点,物体是连续点的集合,可以用构形(物体在空间所占的区域)表示;若是开集,则物体单元表现为数学上点的无穷小邻域,物体是作为拓扑基的所有点邻域的并集,可以用微分流形(容许拓扑结构改变的物体表示空间)表示.从逻辑上看,目前的连续介质力学是从经典质点力学类推得出的,它一方面把物体看作连续的质点系,物体单元具有离散特征,一方面又以场的观点看待物体的内部变化和受力,物体单元变化特征要求是连续的.在质量连续性假设下,物体单元虽然宏观意义上可以看作无穷小但总还是有尺度内涵的,即具有连续性适用的典型尺度,而经典力学中的质点却没有尺度内涵德冈辰雄指出“, 连续介质无论怎样分割也不会成为质点,质点无论怎样连续也不是连续介质”我们知道,经典力学中的质点在数学上表现为三维欧氏家间中的一点(闭集),把表征物体单元的数学上的点看作闭集,无异于沿用质点力学的观点,抹杀连续介质与质点系的区别,这样导出的连续介质力学(简称为质点观点的连续介质力学)是质点观点和场观点的大杂烩,这样的一种结合虽然使连续介质力学在其发展过程中可以同时借鉴经典力学和场论的一些成果,却妨碍了连续介质力学的现代发展,比如运用场论的现代发展—规范理论于连续介质时就显得不伦不类.实际上,质点观点在赋予物体变化连续性的同! 讨,对物体的表示空间强加了过分的约束.限制了场的观点的发挥,使连续介质力学在描述物体复杂宏观力学过程时困难重重.为了使连续介质力学摆脱质点观点的限制,.采用与现代场论一致的基本观点,物体单元用数学上的开集表示是必须的,这时连续性可以用邻域而不是距离定义从而与拓扑学的概念一致,称之为拓扑观点.我们知道,拓扑学是现代微分几何的概念基础,现代微分几何是规范场论的数学基础,因此,拓扑观点的连续介质力学是连续介质的纯粹的场理论,它可以容许物体空间拓扑结构的改变,能够刻划物休的复杂变化过程.可见,物体单元的开集表示与场的现代观点是同气共枝的,由此导出的理论保证了数学概念上的连贯、逻辑上的统一,并且能接纳耗散结构作为物体复杂变化的物理基础.二、流动与变形物体的流动由物沐单元的运动组合而成,物体的变形由物件单元的变形组合而成.物体单元不同于质点: 物体单元的开集表达隐含着单元具有尺度内涵,作为开集的点不仅有平移特征还有方向特征和尺度特征,从而可以独立地体现介质的变形和转动.物体单元的这些特征预示着单元的变形和单元的运动是两个不同的变化过程,物体单元的变形表现为点(及其邻域)的特征的改变,包括尺度的改变和方向的改变,物体单元的运动则表现为点(及其邻域)的平移(空问位置的改变)和转动(方向的改变),可见,单元的变形与其空间位置无关,单元的运动与其尺度特征无关.与此不同,作为闭集的点不具备尺度特征和方向特征,不能独立地体现介质的变形和转动,介质的变形是通过介质点之间距离及相对方位的改变体现的,介质的转动也是通过不同介质点之间的方位关系体现的,这就客观上对物体表示空间提出了要求,难以刻划复杂的变形过程,而单元的运动由于缺乏方向性,对物休单元具有曲线运动的流运过程就无法准确把握.三、局部与整体物体的局部变化是指组成物休的各个单元的变化,物体的整体变化是指物体整体特征或性质的变化.物体单元的变化除了运动和变形外,还有该单元的相邻其它单元的物质交换,这种交换可能是微观的(分子级的),也可能是细观的(源于结构的变化并具有耗散结构尺度的),一般物体单元的转动不均匀性会严重影响这种交换过程;物体的整体变化不仅包括组成物体的各单元的变化,还包括物体表示空间的拓扑结构的变化,后者可以用单元问的变化联络关系表达.一般来说,物体的整体变化不能用其局部变化的直和表示.质收观点的连续介质力学限制了物体空间性质的改变,各个变化阶段的物体的表示空问要求是拓扑等价的,物体单元变化的直和等价于物体的整体变化,因此客观上要求:l)单元间的物质交换与方一向无关;2)单元的尺度变化与方向无关,也就是说,物体单元的变化是各向同性的,这相当于平直层流和均匀变形或者转动影响可忽略的微小变形的情况.在大多数宏观现象中,物体实际变化状态不满足上述要求,质点观点的连续介质力学不再适用,必须用拓扑观点考察物体单元间的变化联络关系的影响,全面研究物体的整体变化过程.四、内应力物体的变形使物体的各部分之间存在相互作用,物体这种反抗变形的内部作用称为内应力,包括应力和应力偶.具体而言,在各物件单元的表面作用有应力和应力偶,这种作用不仅与该单元的纯变形有关,还与该单元的相对转动(净转动)有关,这样,质点观点的连续介质力学中的应力原理必须修正,而非极性物体内应力偶的存在成为可能的了.拓扑观点的连续介质力学给出的非均匀有限变形理论更合理和先进,可统一壳体等转动(方向性)占优的变形理论,并且在这一新观点下,加深了对物体塑性的理解。
[精品]连续介质力学习题二.doc
连续介质力学习题二二.变形与运动2-1如果物体在运动过程中保持任意两点间的距离不变,则称这样的运动为刚体运动,试证:物体的运动若为刚体运动,则参考构形中的物质点反变换到当前构形中的空间位置玉时,必满足:je = <2(0-(X-A) + a(Z),其中如)为正常正交仿射量。
2-2现取物质坐标系{X%和空间坐标系{铲}为同一个直角坐标系,其单位基向量为01,药,码),有一物体的变形为:x1 =X1+k0X2,x2 =X2,x3 =X3,试写出以下各量:1)变形梯度张量户和变形梯度张量之逆F-1;2)右,左Cauchy-Green张量C,B;并计算。
和万的三个主不变量;3)写出。
和万的特征方程,并求出三个特征值化和相应的特征方向兀和八0 =1,2,3) o4)试给出极分解F=V R中的左伸长张量V和正交张量在的矩阵表示。
2-3现取物质坐标系{X”为直角坐标系{X, Y, Z},空间坐标系{?}为圆柱坐标系{r,0,z},令z轴与Z轴重合,0 = 0与X轴重合,图示长方体发生纯弯曲, 变形满足r = r(X),。
= 0(丫),z = z(Z),且存在逆关系:X = X(r), , Y = Y(3), Z = Z(z),试写出以下各量:1)变形梯度张量户和变形梯度张量之逆F-1;2)右,左Cauchy-Green张量C,B;并计算。
和万的三个主不变量;3)写出。
和万的特征方程,并求出三个特征值化和相应的特征方向兀和n(a =1,2,3)。
2-4现取空间坐标系任}为直角坐标系,其单位基向量为(外弓,弓),有一物体的小变形位移场为u = (%1 - x3)(x' - x3)^ + (%2 + x3)(x2 + x3)e2 - ,试求:(1)P (0, 2, -1)点的小应变张量小转动张量。
及其反偶矢量污;(2)求P点在7 = (8召一弓+4弓)/9方向上的线应变;(3)求P点在7 = (8甬-弓+4弓)/9和万=(4弓+4弓-7弓)/9二方向上的直角的变化量。
连续介质力学
目录1简介2基本假设3研究对象4古典连续介质力学5近代连续介质力学6主要分支学科简介研究连续介质宏观力学性状的分支学科。
宏观力学性状是指在三维欧氏空间和均匀流逝时间下受牛顿力学支配的物质性状。
连续介质力学对物质的结构不作任何假设。
它与物质结构理论并不矛盾,而是相辅相成的。
物质结构理论研究特殊结构的物质性状,而连续介质力学则研究具有不同结构的许多物质的共同性状。
连续介质力学的主要目的在于建立各种物质的力学模型和把各种物质的本构关系用数学形式确定下来,并在给定的初始条件和边界条件下求出问题的解答。
它通常包括下述基本内容:①变形几何学,研究连续介质变形的几何性质,确定变形所引起物体各部分空间位置和方向的变化以及各邻近点相互距离的变化,这里包括诸如运动,构形、变形梯度、应变张量、变形的基本定理、极分解定理等重要概念。
②运动学,主要研究连续介质力学中各种量的时间率,这里包括诸如速度梯度,变形速率和旋转速率,里夫林-埃里克森张量等重要概念。
③基本方程,根据适用于所有物质的守恒定律建立的方程,例如,热力连续介质力学中包括连续性方程、运动方程、能量方程、熵不等式等。
④本构关系。
⑤特殊理论,例如弹性理论、粘性流体理论、塑性理论、粘弹性理论、热弹性固体理论、热粘性流体理论等。
⑥问题的求解。
根据发展过程和研究内容,客观上连续介质力学已分为古典连续介质力学和近代连续介质力学。
基本假设连续介质力学的最基本假设是“连续介质假设”:即认为真实的流体和固体可以近似看作连续的,充满全空间的介质组成,物质的宏观性质依然受牛顿力学的支配。
这一假设忽略物质的具体微观结构(对固体和液体微观结构研究属于凝聚态物理学的范畴),而用一组偏微分方程来表达宏观物理量(如质量,数度,压力等)。
这些方程包括描述介质性质的方程(constitutive equations)和基本的物理定律,如质量守恒定律,动量守恒定律等。
研究对象固体:固体不受外力时,具有确定的形状。
连续介质力学引论
流体运动学
断裂力学
流体动力学
用哲学上的一个比喻: 连续介质力学是“共性”,它的研究具有 一般性;弹塑性力学等是“个性”。
连续介质的运动学
1.物质坐标
P(t) P(t0)
X
x
0
那么有如下形式:
x x( X, t )
(3.1) 就描述了在t=t0是位于X的每个质点的轨迹 (对于不同的质点,X不同)。
二连续介质力学关于上述连续体的力学即为连续介质力学1连续介质理论物质是由许多微小的粒子组成的所以物质并不是连续的但是在很多情况下为了描述宏观现象间的关系而不考虑微观尺度上物质结构的理论称为连续介质理论
连续介质力学引论
绪论
一、连续介质
1、质量分布密度
设有一定质量的物质充满一定的空间,并 设P为空间内的一点。取一系列子空间,使 其收敛于P。以 S n 表示第n子空间,其体积 为 Vn ,其中的物体质量为M n,则如果
(1)当某个采用物质描述时, θ= θ(X1,X2,X3,t) 那么,
D Dt t
Xi — —固定的
(3.3)
(2)当这个量采用空间描述时, θ= θ(x1,x2,x3,t) 又,xi=xi(X1,X2,X3,t) 于是:
D x1 x2 x3 Dt t X i 固定 x1 t x2 t x3 t t
Mn lim n V n V 0
n
存在,则将此定义为在P点处物体的质量分 布密度。
如果在所设空间内各点处都有这样的密度, 则质量被认为是连续分布的。同样,可以 定义动量密度、能量密度等等。所谓连续 介质,即指这样一类物体,它的质量密度、 能量密度、动量密度等,从上述意义上说, 都是存在的。简单地说,即认为真实的流 体和固体是由连续的,充满全部空间的介 质组成 。
演示文稿连续介质力学第二讲
所以:
J J X A vi J vi Jdiv v xi X A xi
div v 0
2. 动量方程 (Balance of linear momentum )
2.1 以前的推导
在即时构形中,任意取一个域V ,体积元记为dV
对此域运用动量定理:
σ nda fdV aˆdV
d dt
vdV
f
dV
ቤተ መጻሕፍቲ ባይዱ
σ
nda
dv vdivv σ f
dt
div v 0
3. 角动量方程 (Balance of angular momentum ) 所以:
4. 守恒率的一般形式 如果采用欧拉描述,上述三个守恒率可表达为:
固体力学常采用拉格朗日描述:
其中: 拉格朗日描述中,体元体积不变:
可以推广于多个二阶张量点积的情况,例如 tr(a b c d)
w Jσ : D τ : D 的其它表达形式
由于: τ P FT F T FT
有: w P FT : L tr P FT T L tr F PT L
引理1:设a与b为二阶张量, 则:
a : b tr(a bT ) tr(aT b) aT : bT
引理2:
即: aijbij aijbTji aTjibij aTjibTji
tr(a b c) tr(b c a) tr(c a b)
即: aijbjk cki bjk ckiaij ckiaijbjk
对物质坐标求散度
5. 能量平衡律 在即时构型中任意v域内的总能量P由动能K与内能E组成,即
PKE
E edV
V
根据热力学第一定律,总能量P的物质导数,即对时间的 变化率等于作用于v域的外力功率与每单位时间从v域外部 所加的热:
连续介质力学
所以: J 0 其中: J dV dV 0
对方程两边求物质导数:
J
0
可证明: J div v
J
J
所以: div v 0 率形式的质量守恒律
证明: J div v J
引理: 设矩阵a的行列式为: a , 元素 amn 的代数余子式记作 Amn
将行列式 a 看作它的9个元素的函数,则有:
a
V
式中h表示热流矢量(或称热通量),即每单位时间每单位 面积的热流,k表示每单位质量接受外部的热(称为热源)
而 P K E
其中K为动能.
动能 K 1v2dv
v2
其中 v2 v v
由质量守恒知: dv的物质导数为零
所以: K 1 d v2 dV V 2 dt
又 1 d v2 v dv v a
三、应力理论
1. Cauchy应力
ijeie j
定义在即时构形中的应力张量 又称真应力.
变形后斜截面上的应力矢量:
作用于 da上的力:
σn σn
pnda σ nda σ da
Cauchy应力是以即时构形 中的面积为基准来度量的。
由微六面体的力矩平衡,可知经典连续介质学理论中 σ
为对称张量,即:
d dt
vdV
f
dV
σ
nda
dv vdivv σ f
dt
div v 0
3. 角动量方程 (Balance of angular momentum ) 所以:
4. 守恒率的一般形式 如果采用欧拉描述,上述三个守恒率可表达为:
固体力学常采用拉格朗日描述:
其中: 拉格朗日描述中,体元体积不变:
对物质坐标求散度
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
[
λ2
ˆ λ3 • P
]
1 ˆ = P • diag λ1
T
1
λ2
1 ˆ •P λ3
ˆ ˆ ˆ R = F • U −1
例1. 图示二维线元变形 v v v v ˆ ˆ • dX + F • dX ˆ dx = F • dX = F I II dx1 = F1I dX I + F1II dX II = x1,I dX I + x1,II dX II x2 v dx 2 e 2 dx 2 = F2 I dX I + F2 II dX II = x 2,I dX I + x 2,II dX II
L形式 形式: 形式
B XⅢ
v ∂r v v dr = dX K R ∂X K O v ∂ ( xi ei ) XⅠ dX K = ∂X K v v v ∂ xi ei E K = • d X M EM ∂X K
v dR P
XⅡ
v dr
p
v r
b(t)
(
) (
)
∂x i FiK = ∂X K
v v v v dr = ( x∇ X ) • dR
同理: 同理
2 E IJ = U I , J + U J , I + U M , I U M , J 2e ij = ui , j + u j ,i − um ,i um , j
v v 是同一个矢量。 注:U与u是同一个矢量。
讨论: 讨论:
v v (dr )2 = dr • dr = δ ij dxi dx j
1ˆ ˆ ˆ e= I − F −1 2
( )
T
ˆ • F −1
1 ˆ ˆ = I −c 2
[
]
ˆ 应变张量, ˆ 变形张量。 称e为Euler应变张量, c为Cauchy变形张量。
应变张 量与位 移关系: 移关系
v v ˆ ˆ ˆ v = x∇ v − I = F − I U∇ X X v ˆ = U∇ v + I ˆ F X
∂x i ∂x i v v dr( 1) • dr( 2 ) = dX 1 dX 2 = C12 dX 1 dX 2 ∂X 1 ∂X 2
∴ dr( 1) dr( 2 ) cos θ = C12 dX 1 dX 2
而dr( 1) = C11 dX 1 , dr( 2 ) = C 22 dX 2 ∴ cos θ = C12 C11 C 22 = 2 E12 2 E11 + 1 2 E 22 + 1
v dX I I E I I
XⅡ v dX XⅠ
v dx
v ˆ ⋅ dX F II
v dX I E I
v ˆ ⋅ dX F I
x1
v dx 1 e 1
分析剪切变形(两坐标架重叠 两坐标架重叠) 例2. 分析剪切变形 两坐标架重叠 x1 = X I + bt , x 2 = X II , x 3 = X III
v 使得dX沿三个方向分别伸长 µ i 倍。
ˆ 是正交张量, 旋转。 而R是正交张量,只是使之 旋转。
分解方法: 分解方法: ˆ ˆ ˆ ˆ ˆ ˆ F = R • U , F T = U • RT
ˆ ˆ ˆ ˆ C = FT • F = U2
ˆ 是对称正定的, C是对称正定的,可正交 变换为对角张量
一 般 情 况 下
∴ Λ ⋅ (Λ + 2) = 2 N I N J E IJ
v ˆ 、dR,即可计算 Λ。 给定E
2.主应变和主方向 主应变和主方向 下面寻找Λ的极值 下面寻找 的极值
Λ ⋅ (Λ + 2) = 2 N I N J E IJ
注:Λ(Λ+2)=(Λ+1)2-1,其极值即 极值 ,其极值即Λ极值
v v dX K Q NK = 为dR与E K 的方向余弦 dR
∴ N K N Lδ KL = 1
取泛函 φ ( N K ) = N K N L E KL + E (1 − δ KL N K N L )
其中E为Lagrange乘子 其中 为 乘子 极值条件为
dφ =0 dN K
即( E KL − Eδ KL )N L = 0
EKL为二阶实对称张量,故有三个主值与主方向 为二阶实对称张量,
3.变形梯度的极分解 变形梯度的极分解
ˆ ˆ ˆ ˆ ˆ F = R •U = V • R v v v v ˆ ˆ • U • dX ˆ ˆ • U • dX ˆ dx = F • dX = R =R
(
)
(
)
v U是正定张量,如果取坐 标轴为其主轴方向,则 是正定张量, 标轴为其主轴方向, v v v v v v v U = µ 1 E1 E 1 + µ 2 E 2 E 2 + µ 3 E 3 E 3 v v v v v U • dX = µ 1 dX 1 E1 + µ 2 dX 2 E 2 + µ 3 dX 3 E 3
其中b = b ( X II )
db 另记 2 tan β = ⋅t dX II
2 tan β 1 + 4 tan 2 β 0 0 0 1
解:
1 2 tan β {FiK } = 0 1 0 0
0 0 1
1 {C IJ } = 2 tan β 0
(dR ) 2 (dr )
2
v v = dR • dR = δ IJ dX I dX J v v = dr • dr = δ ij dx i dx j
(dr )
2
2
∂x i ∂x j = δ ij dX I dX J ∂X I ∂X J
2
(dr )
− (dR )
∂x i ∂x j = δ ij − δ IJ dX I dX J ∂X I ∂X J
v 2 ∀dX ≠ 0, (dr ) > 0
ˆ ∴ C对称正定
如果仅仅刚体旋转,不变形, 如果仅仅刚体旋转,不变形,则dR=dr,对应的 , C必为单位张量。故C只与变形有关,与刚体转动无 必为单位张量。 只与变形有关, 而变形梯度与变形和刚体转动都有关。 关。而变形梯度与变形和刚体转动都有关。
§2-3 有限变形分析 1.线元长度的变化及微线元间夹角的变化 线元长度的变化及微线元间夹角的变化 v v v 设原E1方向的微线段为 dR = dX 1 E1
ˆ = 1 FT • F − I ˆ ˆ ˆ E 2 v ˆ T = ∇ vU + I ˆ F X
[
]
ˆ 1 ˆ ˆ ˆ E = FT • F − I 2
[
]
v v v v ˆ = 1 [∇U + U∇ + ∇U • U∇ ] E 2 ˆ e= 1 v v v v [∇u + u∇ − ∇u • u∇] 2
2 2
1 + sin β 1 − sin β λ1 = , λ2 = , λ3 = 1 cos β cos β
主方向
{P }
T IJ
cos β 2(1 + sin β ) − cos β = 2(1 − sin β ) 0
E形式 形式: 形式
(dR )
2
= δ IJ
∂X I ∂X J dx i dx j ∂x i ∂x j
(dr )
2
− (dR )
2
∂X I ∂X J = δ ij − δ IJ dx i dx j ∂x i ∂x j
1 ∂X I ∂X J e ij = δ ij − δ IJ 2 ∂x i ∂x j
∂x i ∂x i dX I dX J = ∂X I ∂X J
= C IJ dX I dX J
C12 C 22 C 32
v v ˆ • dX = dX • C
= (dX 1
dX 2
C11 dX 3 ) C 21 C 31
C13 dX 1 C 23 dX 2 C 33 dX 3
diag[λ1
ˆ ˆ ˆ λ2 λ3 ] = P • C • P T
ˆ λ2 λ3 ] • P
ˆ 的主值, ˆ λi为C的主值, P的第i行是对应的(单位)主 方向。 行是对应的(单位) 方向。 ˆ ˆ C = P T • diag [λ1
ˆ ˆ ∴ U = P T • diag λ1 ˆ ∴U
−1
dr − dR 记Λ = dR
v v dX K NK = 为dR与E K 的方向余弦 dR
dr − dR dr − dR + 2dR = ⋅ dR dR
= Λ ⋅ (Λ + 2)
(dr )2 − (dR )2 由 (dR )2
(dr )2 − (dR )2 且 (dR )2
2 E IJ dX I dX J = dR ⋅ dR
记E IJ
1 ∂x i ∂x j = δ ij − δ IJ 2 ∂X I ∂X J
E IJ
1 ∂x i ∂x j = δ ij − δ IJ 2 ∂X I ∂X J
v v ˆ = 1 δ ∂x i ∂x j − δ E E E IJ I J ij 2 ∂X I ∂X J
而E IJ
1 = [C IJ − δ IJ ] 2
小变形时即E 小变形时即 11
v v dR( 1) = dX 1 E1
v v dR( 2 ) = dX 2 E 2
∂x i v v v dr( 1) = dx i e i = dX 1 e i ∂X 1
∂x j v v v dr( 2 ) = dx j e j = dX 2 e j ∂X 2 微 线 元 间 夹 角