基于统一硬化模型的结构性松砂应力_应变关系.docx

基于统一硬化模型的结构性松砂应力 应 变 关 系 *姚 仰 平 余 亚 妮( 北京航空航天大学交通科学与工程学院 ，北 京100191 )摘 要 : 在临界状态理论的基础上 ，通 过 引 入 修 正 屈服函数的结构 性 参 数 ，可更好地描述松砂的塑性变 形 特 性 。 通过引入状态参数来调整硬化参数和剪胀方程 ，从而更好地反映松砂的剪缩特性和对密度 、有 效 主 应力的双重依赖性 。 同 时 ，模 型 采 用 砂 土等压固结线和临界状态线在有效应力幂函数平面内呈线性 关 系 的 特 性 ，更准确地反映了砂土的应力 应 变 特 性 。 通过数值模拟结果与试验结果对比 ，验证了模型的有效性 。关 键 词 : 松 砂 ; 临 界 状 态 ; 结 构 性 参 数 ; 硬 化 参 数 ; 本 构 模 型stress-strain relationship for structural loosesand based on unified hardening modelyao yangping yu yani( school of transportation science and engineering，beihang university，beijing 100191 ，china)abstract : based on the critical state theory，a constitutive model was proposed for structural loose sand thestructural parameter was used to revise the yield function，thereby predicting the plastic deformation of loose sand more appropriately another parameter was used to adjust the hardening parameter and dilatancy equation， thus giving a better description of contraction behavior for loose sand the introduction of the parameters emphasized the dependency on both density and effective principal stress of sand linear isotropic consolidation and critical state lines were applied in void ratio versus the power of effective mean stress space，and more accurately reflected the stress- strain behavior of sand a comparison between model simulations and experimental results showed the validity of the proposed modelkeywords : loose sand; critical state; structural parameter; hardening parameter; constitutive model实 际 工 程 中 ，受 自 然 环 境 的 影 响 ，砂 土 由 于 受 扰动 、冻 结 、原 填 土 中 的 冰 块 或 杂 质 等 的 溶 解 、冲 刷 的 影 响 以 及 地 震 液 化 后 砂 土 颗 粒 的 重 排 列 等 ，孔 隙 比 相 应 变 大 ，形 成 松 砂 。 根 据 been 和 jefferies1 的 定 义 ，状 态 参 量 为 当 前 有 效 平 均 正 应 力 下 的 孔 隙 比 e 与 临 界 孔 隙 比 ec 之 差 ，当 大 于 0 时 ，砂 土 的 初 始 状 态 点 都 位 于 临 界 状 态 线 ( csl ) 以 上 ，砂 土 剪 切 时就 会 发 生 剪 缩 ，处于 此 状态的砂土即为结构性松砂 。从 工 程 问 题 考 虑 ，进 行 数值分析需要分析松砂时 ，研 究 结 构 性 松 砂 的 应 力 应 变关系具有重要的指导意 义 。大 量 的 试 验 研 究 表 明 ，砂土的变形行为不仅依 赖于其相对密 度 也 依 赖于所处的固结压 力1 2 。本文针对结 构 性 松 砂，基 于 临 界 状 态 概 念，通 过引入结构性参数 来修正屈 服 函 数，以 调 整 松 砂 的 屈 服曲面形状，又通过 引 入 状 态 参 量 ，使 得 硬 化 参 量 h 9 11 和剪胀方程都 依 赖 于 状 态 参 量 而 变 化。同 时 模型采 用 正 常 固 结 线 ( ncl) 和 临 界 状 态 线 ( csl) 在e ( p / pa ) 平面内呈线性关系的假设，强化了模型中砂土对密度和平均有效应力的双 重 依 赖 性 。最 后，对 模型进行预测并与试验结果进行对比 。1结 构 性 松砂及其等压固结特性 1. 1临 界 状 态 线临 界 状 态 线 ( csl) 是 指 剪 切 中 砂 土 的 体 应 变 增been 和 jefferies1提 出 状 态 参 量 之 后 ，* 国家自然科学基 金 项 目 ( 50879001 ，10872016 ，11072016 ) ; 国 家 重点基础研究发展规划 ( 2007 cb714103 ) 资 助 项 目 。第 一 作 者 : 姚 仰 平 ，男 ，1960 年 出 生 ，博 士 ，教 授 ，博 士 生 导 师 。e mail: ypyao buaa edu cn收 稿 日 期 : 2010 06 06345 6，lingjefferies、manzari和dafalias、li等等7 以 及 constantine 8 都 对 砂 土 的 变 形 性 状 开 展了 卓 有 成 效 的 研 究 ，提出了很多适 用于密砂和中密 砂 的 本 构 模 型 。30industrial construction vol. 41 ，no. 9 ，2011工 业 建 筑 2011 年 第 41 卷 第 9 期量 为 0 ，偏 应 变 增 量 为 无 穷 大 的 一 个 状 态 。 li( p / pa )和在 e 平 面 内 弹 性 卸 载 线 的 斜 率 假 设wang12 对 早 期 文 献 研 究 后 提 出 ，砂 土 的为 :csl 在e ( p / pa )平 面 内 近 似 为 直 线 ，表 达 式 为 : ( 6 ) = ed dp c ( p)( 1 )ec =ec01 + e其 中 =a2( 2. 97 e)1 + ed式 中 : ec0 为 p = 0 时 对 应 于 csl 上 的 孔 隙 比 ; c 是csl 的 斜 率 ; pa 是 标 准 大 气 压 ，通 常 取 100 kpa， 是材 料 常 数 。1. 2砂 土的等压固结 特 性砂 土 的 固 结 特 性 与 黏 土 截 然 不 同 ，在 =d2( 2. 97 ed )13式 中 : 和 d 参 照 richart 等 提 出 的 经 验 公 式 ，是关 于孔隙比的修正系数 。e ( p / pa )平 面 内 不 存 在 唯 一 的 等 压 固 结 线 ，而 是2松 砂 本 构 模 型 描 述有 无 数 条 等 压 固 结 线 ( icl ) ，彼 此 之 间 也 不 平 行 。砂 土 的 icl 与 csl 规 律 相 似 ，被 描 述 为 :2. 1屈 服 面为 了 更 准 确 地 描 述 结 构 性 松 砂 的 屈 服 特 性 ，引 入 姚 仰 平 等11 在 修 正 剑桥模型基础上提出的针对 松 砂 的 屈 服 函 数 :p (p)( 2 )e =ei0a式 中 : ei0 为 p = 0 时 对 应 于 icl 上 孔 隙 比 的 值 ， 是icl 的 斜 率 ，所 有 等 压 固 结 线 在 相 当 大 的 有 效 应 力 范 围 内 最 后 都 交 于 n 点 ，如 图 1 。py+ 1 = 0( ) 2( 7 )f =m p其 中 = q / p式 中 : 为 应 力 比 ; m 为 临 界 状 态 应 力 比 ; py 为 等 压压 缩的回弹点所对应的 平 均 主 应 力 ，如 图 2 所 示 ; 是 基 于 ncl 与 csl 间 塑 性 孔 隙 比 增 量 确 定 的 一 个 结 构 性 参 数 ，用 来 调 整 处 于 不 同 初 始 状 态 的 结 构 性松 砂屈服曲面的形状 。对 式 ( 3 ) 全 微 分 有 :图 1 e ( p / pa ) 平面内正常固结线 、临 界 状 态 线 、等压固结线及最密固结线 fig 1 normally consolidated，critical state，isotropic consolidationpdp y 1 yden= n ( p)( 8 )paa塑 性孔隙比增量为 :and densest consolidation lines in e versus ( p / pa ) spacepy 1 dpypden= ( n e ) ( p)( 9 )从 图 1 可 以 看 出 ，在 所 有 的 等 压 固 结 线 ( icl )中 ，有 两 条 特 殊 线 : 正 常 固 结 线 ( ncl ) 和 最 密 固 结 线 ( dcl) 。而 砂 土 的 ncl 是 icl 中 在 p = 0 时 与 csl 相 交 的 等 压 固 结 线 ，因 为当 砂土的围压接近 0 时 ，它 的 剪 切 体 变 值 也 接 近 于 0 。ncl 表 达 式 为 :paa塑 性 体 应 变 增 量 d p 与 塑 性 孔 隙 比 增 量 dep 间 的 关vn系 为 :depd p = n( 10 )v1 + e将 式 ( 9 ) 代 入 式 ( 10 ) 中 并 进 行 积 分 可 得 : p en = n ( p )en0( 3 )py 1 dpy 1 + e d pa式 中 : en0 为 p = 0 时 对 应 于 ncl 上 孔 隙 比 的 值 ， n是 ncl 的 斜 率 。dcl 是 指 可 能 的 最 密 等 压 固 结 线 ，定 义 为 : v = ()= ppneaa( 11 )pyp0() ()ppaa p 假 定 与 塑性体应变相关的硬 化 参 数 增 量 dh 定 义为 : d ( p)( 4 )ed =ed0a同 样 地 ，ed0 为 p = 0 时 对 应 于 dcl 上 孔 隙 比 的值 ， d 是 dcl 的 斜 率 。 其 他 等 压 固结线的斜率可 通 过 对 ncl 和 dcl 进 行 插 值 得 到 :1 + ed p( 12 )dh =v n e将 硬 化 参 数 式 ( 12 ) 代 入 式 ( 11 ) 可 得 :2 ( e e d ) 1 + e ()py (p0)dh = p d + ( n d )( 5 )d v=( 13 )=2( en ed )papane31基于统一硬化模型的结构性松砂应力 应 变 关 系 姚 仰 平 ，等将 依 赖 于 应 力 路径的塑性体应 变 修 正 系 数r( ) 定 义 为 :m4 m44 fcr( ) =( 21 )m4 m44 cf式 ( 19 ) 所 示 的 硬 化 参 数 h 可 变 形 为 := 1 + e 1 d ph( 22 ) e r( )v结 构 性 松 砂 屈 服 函 数 式 ( 15 ) 中 的 硬 化 参 数 h用 式 ( 19 ) 替 换 。对 于 等 压 固 结 中 过 en0 的 固 结 砂 土 ，图 2 e lg( p / pa ) 平 面 内 砂 土 的等压压缩和回弹曲线 fig 2 isotropic loading and unloading curves of sand in e-lg( p / pa ) space式 中 : p0 为 初 始 平 均 有 效 应 力 。由 式 ( 13 ) 可 得 :塑 性 体 应 变 修 正 系 数 r( )到 塑性体应变的形式 。1 ，硬化参数就退回 =2. 2剪 胀 方 程 和 塑 性 势 面参 照 修 正 剑 桥 模 型 中 的 剪 胀 方 程 ，并 引 入 特 征状 态 应 力 比 m ，建 立 与 松 砂 状 态 相 关 的 剪 胀 表 达c 1 py = ( pa h + p0 )将 式 ( 14 ) 代 入 式 ( 7 ) 可 得 :( 14 )式 为 :d pm2 2vc( 23 )=p2 + 1 1 ( p h + p ) 1d d) 2f = (= 0a0m pp式 中 :d d 为塑性剪应变增量 。( 15 )式 ( 15 ) 是 本 文 模 型 所 采 用 的 松 砂 屈 服 函 数 。在 式 ( 15 ) 基 础 上 结 合 式 ( 10 ) ，对 砂 土 沿 着 等 p 路 径所以塑性势函数可相应 地 表 示 为 :m2c pxg = + m2 2= 0( 24 )cpp = p0 剪 切 至 临 界 状 态 ( m ) 时 ( 图 2 ) :=式 中 : px 为塑性势函数曲线在 p 轴 对 应 的 平 均 主 应力 。决 定 塑 性 应 变 大 小 的 塑 性 因 子 确 定 方 法 如 下 :对 屈 服 函 数 式 ( 15 ) 全 微 分 有 : 1 ep2()+ 1 ( pa )+ 1 1n= 01 n ep0( 16 )在 等 p 路 径 p = p0 上 ，ncl 与 csl 间 的 塑 性 孔 隙 比 增 量 表 示 为 :df = f dp + f d +f d p= 0( 25 )vp pv()p0ep( )( 17 )=gpcn( 25 a)其 中d v=ppa联 立 式 ( 16 ) 和 式 ( 17 ) 可 得砂土的结构性参数 :由 屈 服 函 数 式 ( 15 ) 和 硬 化 参 数 式 ( 22 ) 可 得 : 1 1 = 1 ( c e ) 1 f f h f 1 + e 1 2( 18 )( 26 )=pvh pvh e r( ) n e将 式 ( 26 ) 代 入 式 ( 25 ) ，整 理 可 得 :为 了 更 合 理 地 反 映 松 砂 剪 缩 体 变 ，借 用 姚 仰 平等9 11 提 出 的 硬 化 参 数 h 的 规 律 ，扩展为如下新形 f fdp +dp 式 ，并 且 对 于 松 砂 ，此 时 m c m f : = ( 27 ) 1 + e f g 1 444m c ( m f )d p1 + eh = h p r( )( 19 )e444v e m f ( m c )对 屈 服 函 数 式 ( 15 ) 变 形 后 可 得 :式 中 : m c 为 特 征 状 态 应 力 比 ; m f 为峰值强度应力 比 。 () p( m )2pppf =+p0 h= 0paaali 和 dafalias5 提 出 m 、m 可 由 式 ( 20 ) 表 示 :( 28 )cfm c = mexp( m)( 20 a)( 20 b)对 式 ( 28 ) 求 导 ，可 得 :mexp( n)m f =fp ( m )+ 1 2(pp)( 29 )=式 中 : m、n 为 模 型 参 数 ，可 由三轴压缩试验得到 。在 剪 切 过 程 中 随 着 当 前 孔 隙 比 e 接 近 ec ， 值 也 慢 慢 趋 近 于 0 ，m c 、m f 就会 相 应地趋近于临界状态 值 m。pa f( p ) ( ) 2 + 1 1 2 m =3pa(m )m ( 30 )32工 业 建 筑 2011 年 第 41 卷 第 9 期塑性体应变增量为 : f( 31 )= 1h g p( 36 )d ij=由 塑 性 势 函 数 式 ( 25 ) 可 得 :ij2g= 1 ( 32 )pm23针 对 toyoura 松 砂 模 型预测与试验结果 的 对 比根据本 文 提 出 的 弹 塑 性 本 构 模 型 对 toyoura松砂在不排水 三 轴 压 缩 条 件 下 进 行 预 测，共 需 10 个c2. 3弹 塑 性 应 变 增 量 的 计 算总的应 变 增 量由弹性 增量与塑性增量组成 :ep( 33 )d ij= d ij + d ij参数，表 1 列出了所需参数。其中 e 、 、m、n、m 和c0 c其 中 弹 性 应 变 增 量 可 由 式 ( 34 ) 确 定 : 是直接 参 考 文 献12确 定 的，泊 松 比 是 给 定 值。参数 ed0 、d 是通过对 toyoura 砂进行等压压缩所得最 密实状态的试 验 数 据 拟 合 得 到，n 是 所 有 等 压 压 缩 试验结果中过 csl 与 p = 0 交点的 icl 的斜率。模 型 中 toyoura 松 砂 的 结 构 性 参 数 确 定 方 法= 1 + d dd e( 34 )ijijm m ijee式 中 : 为 泊 松 比 ; e 为 弹 性 模 量 。弹 性 模 量 e 可 表 示 为 :3 ( 1 2 ) ( 1 + e) pa如 下 : = 0. 008 8 ， = 0 . 0 1 9 ，取 砂 土 最 密 状 态e =( 35 )nc 1 e p表 1模 型 所 用 toyoura 砂 土 10 个 参 数table 1 10 parameters of toyoura sand used in the modelec0ed0 c n dmmn0. 9340. 6150. 0190. 008 80. 003 30. 71. 253. 51. 70. 3的 弹 性 卸 载 斜 率 e= 0. 003 3 ，代 入 结 构 性 参 数 式( 18 ) 可 得 ，最 密 状 态 砂 土 d 0. 46 。图 1 中 横 坐 标 对 应 于 icl 交 点 n 纵 坐 标 对 应 于e0 = en0 = 0. 934 的 点 m，被 假定为最松砂土样本点 。最 松 状 态 的 弹 性 卸 载 斜 率 由 式 ( 6 ) 可 得 : e 0. 007 013 ，代 入 式 ( 18 ) 可 得 ，最 松 状 态 砂 土 l 0. 73 。为 方 便 模 型 对 toyoura 松 砂 进 行 预 测 ，基 于 上 述 两 个 最 密 和 最 松 状 态 砂 土 的 结 构 性 参 数 ，对 所 有 初 始 状 态 位 于 ncl 以 上 的 松 砂 在 d 与 l 之 间 进 行 线 性 插 值 ，插 值 公 式 可 表 示 为 :图 3 砂 土 屈 服 曲 线fig 3 yield curves of sandpc l ( l d )( 37 )=p0 1 ()en0 e0其 中pc =pa n式 中 : pc 为 ncl 上 初 始 孔 隙 比 e0 对应的有效平均 应 力 的 值 ，当 初 始 孔 隙 比 e0 en0 = 0. 934 时 ，松 砂 可 发 生 液 化 ，此时结 构 性参数可直接取 0. 73 。 图 3 所 示是松砂结构 性 参 数 分 别 取 0 、0. 46和 0. 73 时 的 屈 服 曲 线 。根 据 所 提 出 的 模 型 ，对 toyoura 松 砂 进 行 三 轴 不 排 水 试 验 ，试 验 分 两 组 ，一 组的初始孔隙比为 0. 907 ，初 始 围 压 分 别 是 1 ，2 mpa，另一组初始围压 为 1 mpa，初 始 孔 隙 比 分 别 是 0. 921 和 0. 933 。模 型 预 测 与 试 验 结 果 如 图 4 、图 5 所 示 。 另 外 ，在 初 始 围 压 0. 5 mpa 下 ，对 初 始 孔 隙 比 0. 960 的 砂 土 进 行 三3 = 1. 0 mpa; 3 = 2. 0 mpa; 试 验 结 果图 4 不排水条件下不同初始围压的模型 预测与试验结果对比 ( e = 0. 907 )fig 4 comparison between test results and modelpredictions for different initial confining pressure under undrained condition( e = 0. 907 )33基于统一硬化模型的结构性松砂应力 应 变 关 系 姚 仰 平 ，等轴 排 水 试 验 ，预测曲 线 与试验结果对比见图 6 。力 应 变 特 性 。本 文 提 出 的 松 砂 模 型 可 采 用 yao 等14 15 提 出 的 变换应力三维 化 方 法 ，将 模 型 预 测 扩 展 到 三 维 应 力 状 态 。4结 语本 文 在 临 界 状 态 理 论 的 框 架 内 ，提 出 了 一 个 用 结 构性参数调整屈服面的结 构 性 松 砂本 构模 型 。1 ) 模 型通过引入结构性参数 来调整松砂的屈 服 曲 面 ，并采取插值函数 ，从 而 更 好 地 描 述 松 砂 在 剪 切 过程中的屈服特性和塑性 体 变 。2 ) 引 入 状 态 参 量 来 调 节 砂 土 的 特 征 状 态 应 力 比 和 破 坏 应 力 比 ，进 一 步 使 硬 化 参 量 h 和 剪 胀 方 程 依赖于状态参 数 而 变 化 ，从 而 更 好 地 反 映 松 砂 的 剪 缩特性和对密度 、围 压 的 双 重 依 赖 性 。 e0 = 0. 921 ;e0 = 0. 933 ; 试 验 结 果3 ) 等 压 固 结线和临界状态线在 e ( p / pa )平图 5 不排水条件下不同初始孔隙比的模型预测 与试验结果对比 ( 3 = 1 mpa)fig 5 comparison between test results and model predictions for different initial void ratio under undrained condition( 3 = 1 mpa)面 内采用线性关系 ，所 有 的 等 压 固 结 线 最 后 都 交 于一 点 ，在 大 范 围 的 有 效 应 力 作 用 下 更 准 确 地 反 映 砂 土 的 应 力 应 变 特 性 。4 ) 模 型 预 测 与 试 验 结 果 的 对 比 分 析 表 明 ，所 提 出 的模型仅用一 套 材 料 参 数 ，即 可 较 好 地 描 述 不 同 密 度不同初始有效应力下松 砂 的 应 力 应 变 特 性 。参 考 文 献been k， jefferies m g a state parameter for sands j geotechnique，1985 ，35 ( 2 ) : 99 112 cai z y，li x s deformation characteristics and critical state ofsandj chinese journal of geotechnical engineering，2004 ，26( 5 ) jefferies m g nor-sand: a simple critical state model for sandj geotechnique，1993 ，43 ( 1 ) : 91 103 manzari m t，dafalias y f a critical state two-surface plasticitymodel for sandj geotechnique，1997 ，47 ( 2 ) : 255 272 li x s， dafalias y f dilatancy for cohesionless soils j geotechnique，2000 ，50 ( 4 ) : 449 460 li x s，dafalias y f a constitutive framework for anisotropicsand including non-proportional loading j geotechnique，2004 ，54 ( 1 ) : 41 55 ling h i， yue d， kafalias v， et al anisotropic elatoplastic bounding surface model for cohesive soils j journal of geotechnical and geoenvironmental engineering，2006 ，128 ( 7 ) :748 758 stamatopoulos c a an experimental study of the liquefaction strength of silty sands in terms of the state parameterj soil dynamics and earthquake engineering，2010，30( 8) : 662 678 yao y p，matsuoka h，sun d a a unified elastoplastic model for clay and sand with the smp criterion c proc ， 8 th australia-new zealand conf on geomechanics 1999 ，2 : 997 1003 ( 下 转 第 55 页 )123456 e0 = 0. 960 ; 试 验 结 果图 6 排水三轴压缩条件下试验结果与 模 型 预 测 对 比 ( 3= 0. 5 mpa)7fig 6 comparison between test results and model predictionsfor 3= 0. 5 mpa under drained condition图 4 、图 5 表 明 : 对 应 于 同 样 的 初 始 孔 隙 比 ，砂土 的 初 始 围 压 越 大 ，越 容 易 发 生 剪 缩 ，而 对 应 于 相 同 的 初 始 围 压 ，初 始 孔隙比越大不排 水剪切时更易剪 缩 。同时也 反 映 了松砂不 排水剪切时的结构性 。图6 也 同 样 表 明 了 松 砂 在 排 水 剪 切 时 的 剪 缩 特 性 。 通 过 以 上 的 模 型 预 测与试验数 据对比可以看出 ，本 文 所 建 立 的 松 砂 本 构 模 型 能 够 较 好地描述松砂的应 8934工 业 建 筑 2011 年 第 41 卷 第 9 期图 6 效 果 图fig 6 the rendering得 益 于 昆 明 独 一 无 二 的 气 候 条 件 ，顺 城 项 目 设计 为完全开敞的现代化大型 购 物 娱 乐中 心 。自 然 和 生 态成为重要的设计元素 。花 ，作 为 昆 明 的 特 色 ，被 有 机 地 浓 缩 为 建 筑环境与人之间的情感交流载体 ， 同 时 也 成 为 昆 明 市 的 新 地 标 。 建 筑 立 面 的 处 理 ，充 分 考 虑 到当地历史的 继 承和延续 ，引 入“花 ”的 符 号 ，通 过 现 代 金 属 材 质 的 演 绎 ，使 之 具 有 时 代 性 ，同 时 利 用 传 统 材 质 和当地材料来表现现代商业立面 ， 进 一 步 打 破 传 统 与 现 代 的 界 限 ，让 项 目 体 现“城 市 的 演 化 ”。新 与 旧 共 生 共 存 ，为 城 市 空 间 做 出 贡 献 。与 imax 影 院 作 为 商 业 形 象 的 中 心 相 衬 托 ，西 南 的双塔住宅作为住宅建筑 的 轴 心 ，会 所 置 于 其 间 ， 从 而在空间上引导住户的活 动 领 域 。影 院 门 厅 建 筑 上 作 出 相 应 处 理 ，为 观 众 的 购 票 、等 待 、休 息 和 消 费 提 供 场 所 。新 与 旧 的 元 素 在 和 谐 中 共 生 共 存 ，形 成 完 整 的 统 一 。图 4 二 层 平 面fig 4 the plan of the 2 nd floor连 廊 、跑 道 和 漫 步 小 径 连 接 新 月 厅 ( 住 宅 活 动 区 或行 政 会 所 ) 至 高 消费 温 泉 浴场和健身中心 ，并 连 接 一 个 室 外 游 泳 池 。 建 筑 风 格 被 定 义 为 纯 朴 、现 代 的 住 宅 楼 稳 健 地 坐 落 于 裙 房 之 上 ，并 享受东南和西面 的 最 大 景 观 优 势 。 办公楼也 以相同的设计理念 ，以 简 单 突 出 的 形 式 ，运 用 质感和发光的透明玻璃屏幕 ， 使 都 市 的 夜 晚 绚 丽 多 彩 ( 图 5 、图 6 ) 。3结 语昆明是一个美丽传统与现代文明相融 合 的 城 市 ，