动荷作用下沥青路面结构的变形响应分析
【技术】车轮荷载对半刚性基层沥青路面结构的力学响应分析

【技术】车轮荷载对半刚性基层沥青路面结构的力学响应分析摘要:基于我国的沥青路面设计理论及标准,选取典型半刚性基层路面结构及材料参数,采用BISAR3.0软件对不同荷载作用下的路面结构应力、应变和位移进行计算,并分析了各力学指标对道路结构的影响。
结果显示:不同荷载模式对路面结构的影响相当大,这对进一步解释路面面层的一些破坏现象提供了有益的参考。
关键词:道路工程,半刚性基层,沥青路面,车轮荷载,力学响应引言道路是车轮荷载的承受者,随着我国交通运输业的高速发展,重载超载现象普遍存在,导致路面结构的受力发生了很大变化,我国广泛采用的半刚性基层沥青路面在设计年限内出现了较为严重的早期破坏现象[1-2]。
许多学者对其破坏的机理进行了研究,并建立了路面结构的力学模型及结构响应[3-4],也分析了重载车辆对路面结构层间行为的力学研究[5-6]。
研究表明:重载交通的反复作用导致半刚性基层与沥青面层或面层与面层之间不连续,是路面结构出现破坏的主要原因[7]。
因此,本文结合我国典型的半刚性基层沥青路面结构对其应力、应变和位移进行力学计算,并分析各个力学指标沿路面结构深度方向的变化规律,进一步解释路面结构在重载作用下破坏的原因,为运营中的道路在车辆荷载控制方面提供一定的参考。
1 路面结构、计算参数及荷载的选择1.1路面结构本文选取的是我国典型的半刚性基层路面结构,沥青面层分为三层,基层分为半刚性基层和半刚性底基层,最下面的层位为土基,其具体路面结构参数如表1所示。
表1路面结构层材料参数层位材料抗压模量(MPa)(弯沉计算用)抗压模量(MPa)(应力计算用)厚度(cm)泊松比上面层细粒式沥青混凝土1400200040.25中面层中粒式沥青混凝土1200180060.25下面层粗粒式沥青混凝土1000120080.25基层水泥稳定碎石130********.3底基层水泥稳定碎石130********.3土基天然土3535-0.351.2计算参数为了了解弯拉应力、应变和弯沉在不同荷载作用下沿道路深度的影响,用bisar3.0软件进行力学计算,计算点位选取在单轮荷载中心。
沥青路面力学响应分析及其研究方法综述

沥青路面力学响应分析及其研究方法综述发布时间:2022-07-11T02:28:48.391Z 来源:《工程管理前沿》2022年5期3月作者:黄勇维[导读] 沥青路面具有比较复杂的力学特性,为了研究沥青路面的破坏机理以及路面应力应变变化规律,黄勇维重庆交通大学土木工程学院重庆 400074摘要:沥青路面具有比较复杂的力学特性,为了研究沥青路面的破坏机理以及路面应力应变变化规律,本文系统阐述了荷载、温度、路面结构类型以及层间接触状态对沥青路面结构力学的响应机理。
并且鉴于以往对路面进行力学研究不能够准确、真实、细致的反映其力学行为的问题,本文简述了对沥青路面细观力学行为的研究,使沥青路面力学的研究能够宏、细观相结合。
研究发现,细观力学分析能对内部材料变化进行量化处理,全面分析沥青路面力学响应,对改善路面性能有重要意义。
关键词:路面力学响应;荷载;温度;层间接触状态;细观力学研究0 引言我国沥青路面损坏影响因素主要有材料、荷载和温度。
因此,解决沥青路面这些问题,就要从因素出发,有必要对沥青路面力学响应因素进行分析研究。
沥青路面长期处于不同的自然环境中,并非单一不利因素影响沥青路面,在恶劣的气候条件和车辆荷载共同作用下,沥青路面材料内部逐步发生变化,路面出现宏观的损坏现象。
以往对路面进行力学研究,通常将沥青路面通过假设条件进行了不同程度的简化,与实际情况存在差别,不能够准确、真实、细致的反映其力学行为,因此,有必要对沥青路面细观力学行为进行研究,达到宏、细观相结合的目的,全面分析沥青路面力学响应,对改善路面性能有重要意义。
1沥青路面力学响应分析综述沥青路面是多层路面结构,具有比较复杂的力学特性。
国内外大量研究表明,对沥青路面力学响应有显著影响的因素主要有荷载、温度、路面结构类型和层间接触状态等。
研究不同因素影响下的路面力学响应,可以为更科学合理的路面设计方案提供必要的参考。
1.1荷载在沥青路面的力学性能分析时,通常把轮胎与路面的接触面作为路面受力分析的影响区域。
车辆动荷载作用下路面结构动响应分析最后祥解

交通荷载作用下沥青路面动应变
垂直动应变时程曲线图
在沥青面层顶部,容易有滑 移破坏问题出现,主要是因为在 顶部存在垂直动应变的拉应变与 压应变的交替,极易与基层错位。 路面结构的位移弯沉,主要是由 于存在垂直压应变的缘故,除面 层外,垂直动应变均为压应变, 且最大值出现在在基层与底基层 结合部。沥青路面结构各层在同 一层内水平应变也是交变的,既 有拉应变也有压应变,与横向应 变一样,最大拉应力均在基层与 底基层结合部产生,然横向应变 各层都是拉应变,除了沥青路面 层顶部的压应变以外。
研究现状
现有的沥青路面结构设计理论中,通常采用 静态加载模式设计路面各结构层厚度。这种方法 适用于车速较低、车载较小的情况,但随着汽车 工业的快速发展,车辆的速度越来越快,重型汽 车也越来越多,因此研究行驶的车辆对路面产生 动荷载作用下路面结构的力学响应具有重要的理 论意义和潜在的应用前景。
路面平整度
3、速度对水平拉应力的影响
荷载在沥青路面层底基层与路基结合部产生的水平及横 向拉应力最大。
不同车速时底基层最大水平拉应力
4ห้องสมุดไป่ตู้速度对剪应力的影响
不同车速时沥青路面层顶部最大剪应力
结论
分析路面结构在车辆动态荷载作用下的竖向位移、动应力、 动应变、剪应力响应,得到以下结论:
(1)在车辆动荷载作用下,随着路面深度的增加,沥青路面 的动态响应峰值都在减小。
模拟路面平整度的结论
按照频谱表示法模拟路面平整度的理论,编制生成随机 路面的matlab程序,分别对A、B和C级路面进行时域仿真, 假设行车速度分别为5m/s、10m/s、15m/s和20m/s,仿真时间 取200s,即模拟1000m、2000m、 3000m和4000m长度的路 面的平整度时域模型,得到不同速度下的路面平整度时域模 型,通过分析我们可以得出这种结论:
动、静荷载下不同沥青路面结构力学响应分析

动、静荷载下不同沥青路面结构力学响应分析作者:何基雷罗资清傅松来源:《西部交通科技》2024年第03期作者简介:何基雷(1988—),工程师,主要从事道路工程、路面养护方面的研究工作。
为探究动、静荷载下沥青路面结构的应力响应,获取不同影响因素对路面的实际作用效果,文章利用ABAQUS软件构建了沥青路面结构应力响应模型,分析荷载形式、车辆轴载、行驶速度等因素对力学响应的影响。
研究表明:路面结构的应力应变与车辆轴载存在着一定的线性关系;相较于静荷载,动荷载在相同轴载下所产生的应力应变值较低,且存在最佳行驶速度使荷载对路面产生的力学响应最小。
由此证明,在道路使用时,控制车辆的行驶速度及车辆超载可减缓路面纵向位移及路表弯沉的产生,延长道路的使用寿命。
沥青路面结构;移动荷载;力学响应;使用寿命;应力应变U416.217A1906850引言随着我国机动车保有量及道路交通量的逐年上升,道路重载及超载现象的持续增长,使得已建道路在使用过程中暴露出使用寿命不足[1-2],裂缝、坑槽、松散、剥落、车辙等病害出现频率较高的现象。
道路养护时运营成本增加,而且还影响了交通事业的发展[3]。
因此,为更好地了解路面结构在不同因素下的力学响应,需探究不同影响因素对路面的力学响应。
国内外专家学者针对沥青路面的应力响应从多方面展开了研究。
Assogba、Hu、李江等[4-6]通过建立三维有限元模型,研究了车辆速度、车辆超载对沥青路面的影响,证明较低车速会引起结构受载时间增加,扩大了载荷的冲击效应。
严战友、Ogoubi等[7-12]通过建立车辆模型和有限元道路模型,证明路面结构的动态应变应力峰值受分析点位、行車速度、沥青层厚度、车轴荷载、制动工况和道路粗糙度等因素的影响。
Liu[13]通过提出了一种将全尺度加速路面试验(accelerated pavement test,APT)、室内试验和有限元(finite element,FE)模拟相结合的方法,分析了车轮范围、温度及轴重对于沥青路面的动态响应。
移动荷载作用下Winkler地基的沥青路面动力响应

ቤተ መጻሕፍቲ ባይዱ
1 路面计算 力学模 型
为模 拟车辆荷载作 用下路面的动力响应 , 将
收稿 日期:20 —l0. 090 一6 宁 波 大学 学报 ( 工版 )网址 : t :3 b b . u n 理 ht / x . ue . p/ n d c 基 金项 目:浙 江 省 自然 科学 基金 ( 0 67); 江 省科 技计 划项 目 ( 0 8 3 00). Y173 浙 2 0C 13 第 一作 者:刘 干斌 ( 9 6 ), , 西 吉安 人 , 士眉0 17 一 男 江 博 教授 ,主要 研究 方 向 : 基路 面工 程 . — i l g 7@ 13 o 路 Emal i b 6 6 r :u cn
在研 究不平整度对路面动 力响应 影响时, 将位移
函数和 荷 载取 近 似模式 ,展 开为 正 弦级 数 . m 等 Ki
换得到数值结果以分析荷载速度 、 粘弹性地基阻尼 及双轮问距对路面竖向位移 、弯矩的影响.
人l 究了前后双轮变幅值荷载条件下, n l 粘 Wi e kr 弹性地基上无 限大板 的动力响应. 周华飞等人 采 j
( 式和( 式可得到路面动力响应的控制方程为: 1 ) 2 )
载大小为 qe o ,长 2 宽 2 ,以速度 c x方 厶, 沿 z 向运动 [ 考虑双轮荷载工 况条件下, 4 】 . 荷载 的大小
相等. 地基模型采用粘弹性 Wi l 模型, ne kr 路面的 竖向位移为 ( x, . , f z)
第2 3卷第 2 期 21 0 0年 4月
宁 波 大 学 学 报 (理 工 版 )
移动荷载下粘弹性层状沥青路面动力响应模型

A MODEL TO STUDY THE DYNAMIC RESPONSE OF VISCO-ELASTIC LAYERED SYSTEM UNDER MOVING LOAD
*
DONG Zhong-hong , LU Peng-min
(Key Laboratory for Highway Construction Technology and Equipment of the Ministry of Education, Chang’an University, Xi’an 710064, China)
q( x, y ) ,利用 Fourier 展开技术,车轮荷载可表
示为:
U nm ( z ) 、 Vnm ( z ) 、 Wnm ( z ) 仅为深度 z 的函数。
由式(6)~式(8)及几何方程和本构关系可知: ij i nV ij (9)
2 2 ij n V ij
q ( x, y ) Anm ei n x eim y
n 1 m 1
N
M
(4)
(10)
对于以速度 V 沿轴 x 移动的车轮荷载,则式(4) 转化为:
其中: ij 代表 x 、 y 、 z 、 xy 、 xz 、 yz 等 6 个应变分量。
q( x, y, t ) Anm e
———————————————
收稿日期:2010-04-12;修改日期:2011-04-29 基金项目:国家自然科学基金项目(51008030);中央高校基本科研业务费专项项目(CHD2010JC063,CHD2011ZD001) 作者简介:*董忠红(1975―),男,河南开封人,副教授,博士,从事车辆-道路系统动力学研究(E-mail: dzhong@); 吕彭民(1957―),男,陕西渭南人,教授,博士,主要从事工程力学研究(E-mail: lpmin@).
行车荷载下沥青路面响应分析

关 键 词 :公路 ; 动 荷 栽 ;动 力 有 限 元 ;阻尼 矩 阵 ; 面 弯沉 移 表
中 图 分 类 号 : 1 . 1 U4 6 2 7 文 献标 识 码 : A 文 章 编 号 : 6 1 2 6 ( 0 6 0 —0 5 —0 1 7 — 6 8 2 0 )3 0 5 3
{ ) N] ) r 一[ {
将 式 () 人 ()得 : 3代 2,
() 3
一
D]
- f D] B] l [  ̄
( 1 1)
式 中 : 系数 。 为
[ N]
() 4
作用 于单 元结 点上 的 阻尼力
维普资讯
公
56
动 的结构 中取 出一 微 小 部 分 , 它单 位 体积 上作 用 在
的 阻尼力
) 一
( )一 r
N]
() 8
[幻 { +[ ] ) [ { 一 { T) ^ ) C { + K] ) F( )
() 1
式 中 : 比例 常数 ; 为形 函数 。 a是 N 利用 式 ( ) 5 可求 得单 元 结点 上 的阻 尼力
近年来 , 随着经 济 的快 速发 展 , 界 各 国修 建 了 世 大量 的高 速 公 路 。公 路 交 通 越 来 越 表 现 为 交 通 量
大、 汽车载 重增 加 、 速加 快 等 特 点 , 有 这些 不仅 车 所
荷 载移 置 的一般 公式 为
(;肌 捌 P一 ) d z
{ F) IN] { )V— 一 d I
[ 一 IN]pN] V m] [ d I
1 2 阻 尼 矩 阵 .
如 果 假定 阻 尼 力 正 比于 质 点 的 运 动 速度 , 运 从
FWD荷载作用下沥青路面动力响应有限元分析

FWD荷载作用下沥青路面动力响应有限元分析论文
本文旨在探讨FWD(falling weight deflectometer)荷载作用下
沥青路面的动力响应情况,通过有限元分析的方式,分析其结构响应特性以及受力行为。
利用实验数据优化有限元模型,并将其应用于汽车对沥青路面进行路面质量评定。
针对FWD荷载作用下沥青路面,开展有限元分析。
根据有限
元理论,建立一个均匀的有限元模型,并运用经典的梁单元进行模拟,如Young-Von Karman模型。
同时,根据实验数据,
优化模型,使其最大程度反映真实情况。
此外,考虑地面材料的拉伸模量、剪切模量和泊松比,以及基础土的应力应变。
最后,基于不同的FWD荷载作用,计算路面的响应力,以及每
一段路面的形变,其中包括剪切变形、水平和纵向变形等。
结果表明,FWD荷载作用下沥青路面的动力响应随荷载的增
大而增大,荷载强度与响应之间呈线性关系,最终得出路面承载能力的最佳估计值。
此外,FWD荷载作用下沥青路面的形
变情况也随着荷载的增大而增大,且与不同部位的位移及形变有关。
经过有限元分析的研究,我们不仅可以更好地了解沥青路面的动力响应行为,而且还可以将最优预测值应用于汽车对路面进行质量评估中。
然而,路面在实际情况下还存在一些复杂情况,也需要进一步的研究和实验支持,更好地预测路面的响应性能。
总而言之,本文通过有限元分析的方式,研究FWD荷载作用
下沥青路面的动力响应现象,并优化有限元模型,更好地预测沥青路面的响应性能。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
A s at D nmirso s o l ee p vm ns j td dnmil d ia a zd ui bt c: a c pne a rd e et sbe e t ya c s nl e b s g r y e f y a u c o o s y a y n te e-i es nl f i -l et to i cnu ci wt N WMA K t rt n h tr dm ni aOD) t e m n m h d o jnt n h E he o ie e n e n o i R i e ai ng o shm . peet u s cm ae wt te ib terta rsl t vlae ce eT e sn r l ae pr i h aaal hoei l ut adt te h r e t r o s d h v l e c e o i h
动荷作用下沥青路面结构 的变形响应分析
侯 芸‘郭忠印’田 波“杨 , , ,
(.同济 大学 交通运输 工程学 院 上梅 1
众’
1 08 00 8)
2 0 9 :.交通部公路科学研究所 , 京 0022 北
摘 要: 应用三维有限元动力学的基本方法, 结合 N WMA K积分方法逐步求解运动方程, E R 对动 载( WD荷载) 用下 多层弹 t体 系的响应进行 了分析 并将计算结果和理论 结果进行对比, F 作 i 表明 计算模型合理 , 计算结果准确。利用上述计算模型, 对于多种路面结构研究其在不同荷载周期下的 动态反应规律。 关键词: 三维有限元;WD荷载; F 动态弯沉; E N WMA K积分 R 中图分类号 : 1- 1 U4 6 2 7 文献标识码 : A
式中:i} i 为整体坐标和 自然坐标 系统之 间的雅克 比转换 ;N」 [ 由下式表示
本次网格 自动剖分的原则为:
() 1根据道路层状结构的特点, 以每一层作为单 元划分的基本单位, 即每一层单元材料性质必须是
相同的, 一个单元将不能穿过本层的分界线 。
内 饰 卜 四
0 0 . V, 0 N
卜:
4 计算分析
‘1 路面结构的动态反应规律
通常路面结构可视为半无限地基上的弹性多层 体 系, 但在有 限元 的计算 中, 上基只能取有限尺寸。 为确定半无限体厚度取值问题 文中选用不同深度 的均值体有限元解和 L AMB问题理论计算结果进 行 比较 , 结果发现 当厚度取值大于 l m 时 , 两者的 挠度接近。同时根据文献〔 们的分析, 对于F WD的 冲击荷载, 沿深度扩展到 1 m 可满足计算要求。本 2
3 用 N WMA K 积分方法求解运 E R
动方程
可以用振型叠加法和直接积分法求解式 ( ) 1所
N (o:+(,。 , 1l。: +: = +( )。 。 音 F十() + ,,; N ‘'p+ ; 卜lo} = ‘+() 专 (1n 11 N 〔XE+ 。 卜1a} 一 :+(o 十 ( )' 1 ll N (l;+ ; 卜)op = }十(a : ( ) 1 1)
第 3期
侯
芸, : 等 动荷作用一沥青路面结构的变形响应分析 「
到弹性系统动力问题的有限元基本方程如下
[ { + ( {} K {} P t 1 [〕沙 十[ ]8 = ; ( }() M口s 一 ) 式 中:M][ ][ 分别为系统的质量矩阵、 [ ,C ,K] 阻尼
矩 阵和刚度矩阵;S .8 .创 分别为位移、 { }{ }{ 速度和 加速度向量 ;P() 为动荷载 { t全
f .1 i 1 =9
G G1 =1, 2 0 31 =1, 4
玩一+ + 1a+ 。 Ss [-I ] ,, } Sa “ l - A 2 .‘
得 到 以 S S S 表示 的 S 、 , . 的计 算 式 +
式中:. aR为按精度和稳定性要求确定的参数。
I ( k 1 0 A )
(0 推导 饥恤 并将其代 入式 ( ) 经过整理和 一 N」
0
0
: \ ,
() 2 以规整的矩形单元作为一级剖分单元, 通过 输人其最下层的结点和单元信息, 逐层推算其上部 各结点和单元的信息 () 3 对一级单 元进行 二级剖分, 获得二级 剖分 单元。它们就是进行有 限
元计算的基本单元 。
由于阻尼的影响通常比惯性和弹性的影响小, 可以利用一个近似的阻尼矩阵来描述阻尼特性的影 响。 假定阻尼力正 比于质点运动速度时, 可以得到如 下的单元阻尼矩阵
应用弹性力学三大定理和虚功原理可推得应变
转换矩阵「 ] B 和弹性矩阵L , DP具体表达形式可参考 文献1] 3 。下面介绍单元的刚度矩阵和质量矩阵, 并
对阻尼矩阵作简要说明。 单元刚度矩阵「 K月为
八T -卞a 川 。 “ 什 卞0 a 一 } } e t 汁 一厂 山 ` m= -L
CC {,:〕C一jd( CS J ’ } N ,d 8 "6 丁〔a ,J p ) I。 6一 ・T cs [ 一 d} N f
式中: 为密度托为阻尼系数。 P
22 三维有限元计算式 ・
从保 证计算精度 的角 度出发 采用 2 0结点空间 等参元( 1 , 形 函数 图 )其
为
田 t 典型 的 2 0结点等参 元
划分 。
单元质量矩阵[ ] M,为
A.一S 1J 一 Id5 l。,, [o [,d ( ‘ 丁 Nx 3・ g ) 丁 P 3 , } [ d N f
对于第 一个 时 间步, 到 均 布荷 载 尸() 受 t 的作 用, 等效结点荷载列阵为
、。一 丁N一 ( , 尸。 丁 [;、) , x ] P }A td
收稠 日期 20-80 0 1 -3 0
使计算结果可靠。 仁 述两方面的需要推动了路面动 态响应的研究。国内关 于这方面的研究多采用二维 算法, 由于实际的路面为三维条带状结构, 在向二维 简化过程中存在不少的内在假定川, 因而计算结果 有不少的局限性 。为 了更加真实地模拟道路在动态
荷载作用下的响应问题, 笔者采用 2 结点的二维块 0 状单元。 配合 N WMAR E K积分方法的使用, 对运动
第I 5卷 第 3 期 20 0 2年 7月
中 国 公 路 学 报 Ci Junl i w y Tas r h a r oHg a ad np t n o a f h n r o
V t1 \ o. ( 5 Jl 02 2 0 uy
文幸编号 10-3220 )300-5 0177(020-060
分解 。
(2 1)
次计算模型的尺寸为 1 mxl n又1 n, 2 Z l 2l 底面为固
[ 〕 。一 K 。 。
一 ・: : : ‘
如
气
() 3
一 、
走 ,
[ \丁} ']jl7 k 3’f [[1eg 〕一J ,] BJd [ n dd B .
(, , , 2 ) ej 一12 一。0
( 理)
2 三维有限元算法
2 1 网格自动划分 .
在进 行一般 的有 限元计算 中, 初始数据的输人 在技术上要求不多。而对 于单元数量较多的三维问 题, 则因其复杂和抽 象容易将数据输错或漏输 。因 此, 有必要在关键数据输入的基础 _ r进行 网格 自动
式中:o , ? ; ; ¥=E w=, }=y 色 . o 氛 n
2 ( -8 ) ‘ 二1 )
示的运动方程。其中直接积分法中的 N WM八 K E R 积分方法是应用较广的求解方法。N WMA K积 E R 分方法实质上是线性加速度法的一种推广, 采用 如
下假 设
5} = . ( 一 ) , S + [ 1 8 S十 口咨 ,山 ,. ,八 〕 () 9
( ”在形成刚度矩 阵 K、 质量矩阵 叼 和阻尼矩
阵 C的基础上, 给定 占、。占, 。子、。并选择时间步长 匆,
参数 a 使得 月 .,妻。2( 5 ), 渭, 笋。5a ,50 十召,并计算 .
积 分 常数
_ 夕
口色 妞
a, 六 1一 一 益・ 一一誊 1 3 , , 一
・ { 2 “一,加 。 子 ) “夕?卢 一 一, 警 一 , “ 一
acr y te eet oi m. s d i cnut t ea a te et o vr u cua o h p sn a rh A uy od c d vl t h ef s ai s c f r l t g t s e o u e f c f o prm t s te nmi rso s o l ee pvm ns j t t dnmilas aa ee o h d a c p ne a rd e et sb ce o a c d. r n y e f y a u e d y o K y rs tr -ies nl 3 f i- e et to ; l d dnm c fci ; e w d : e d ni a (D) t e m n m hd F o he m o i el n e WD a ; a i d l t n o y ee o
NE WMAR itgain K e rt n o
0 引
言
现行路面设计方法采用竖向静载作用下弹性多 层体系的理论模型 对于实际道路上行驶的车辆, 由 于路面不平整而产生振动, 路面受到的车辆荷载随 时间和路表特性而变化 , 是一个典型的动力荷载 为 了深人研究路面在动态荷载作用下的破损机理和改 进路面设计的力学模型, 需要关心动力荷载下路面 的响应问题。 此外, 路面结构评价中采用无破损检测 设备, F 如 WD, W 等, SS A 施加在路面结构上 的荷 载为动力荷载, 因而在进行材料特征参数反算时, 需 采用动力荷载作用下弹性多层体系响应的程序才能
方程进行求解 , 并分析 了荷载参数对路面响应的影 响。
1 运动方程
应用弹性力学的 HAMIT L ON 变分原理 , 可得