第三章 构造研究中的应力分析基础
一、应力
内力与面力、体力、外力、内力等概念不同，在固体力学中应用广泛。物体的变形是由内力直接引起的（外力仅是引起内力改变的原因），因而在固体力学中更关心物体内部各部分之间的内力的变化情况，因而引入“应力”这个概念。
为了研究物体内部某点的内力分布状态，通常设想过该点作一个微小面积的截面。设这一微小截面的面积为ΔF ，作用在该截面上的内力为ΔP ，则将 称为该点处该截面上的应力。由于内力是矢量，应力也是矢量。
由于ΔF 截面上的内力P 可以分解为垂直于截面的内力分量N 和平行于截面的内力分量T ，相应地，应力也可以分解为垂直于截面的应力（σ）和平行于截面的应力（τ）
垂直于截面的应力（σ）称为该截面上的正应力，平行于截面的应力（τ）称为该截面上的剪应力。
二、主应力、主方向和主平面
在物体内部的某点处总是可以找到这样一个包含该点的微小的正六面体（立方体），它的三对正交截面上没有剪应力而只有正应力作用，这种情况下的三对正应力称为该点的主应力，分别用σ1、σ2、σ3表示，并规定压应力为正，拉应力为负，在代数值上保持σ1＞σ2＞σ3。 主应力的方向称为该点的应力主方向，三对截面则称为该点的三个主平面。
一点的3个主应力决定了该点的应力状态，当3个主应力中有两个为零时称单轴应力状态；有1个为零时称双轴应力状态或平面应力状态；当3个主应力都不为零时称为三轴应力状态。
三、应力莫尔圆
应力莫尔圆是一种重要的图解方法，可以直观地表示一点的应力状态。
以横坐标代表正应力σ，纵坐标代表剪应力τ，根据σ1、σ2和σ3的大小作出的用以反映一点应力状态的一个圆就是应力莫尔圆。
1．单轴应力状态的二维应力莫尔圆
圆上任一点的坐标代表与主应力σ1呈θ夹角的截面上所受到的正应力和剪应力。（其中θ是过该点的半径与横坐标轴所呈夹角的一半）。
2．双轴应力状态的二维应力莫尔圆
与单轴应力状态类似，圆上任一点的坐标代表与主应力σ1呈θ夹角的截面上所受到的正应力和剪应力。（其中θ是过该点的半径与横坐标轴所呈夹角的一半）。
从单轴和双轴应力莫尔圆上可以看出：
（1）剪应力互等定律：在两个相互垂直的截面上剪应力大小相等，方向相反；
（2）正应力之和守恒：在两个相互垂直的截面上正应力之和不变，等于主应力之和；
（3）在与最大主应力σ1呈45?和135?的截面上所受到的剪应力最大。
3．三轴应力状态的三维应力莫尔圆
三个分别包含σ1和σ2轴、σ2和σ3轴、σ1和σ3轴的3个二维应力莫尔圆共同组成的区域内的任一点的横坐标和纵坐标即代表了三维空间中某截面上的正应力和剪应力。
该图上也可以看出最大剪应力位于σ1和σ3构成的应力圆上，位于与σ1呈45?或135?夹角的截面上。
p dF dP
F P
F ==??→?0lim dF dT
dF dN
F F 0
0lim lim →?→?=＝和τσ
四、应力场、应力轨迹（主应力迹线）与应力扰动（或应力集中）
应力场：受力物体内的每一点都存在与之对应的应力状态，物体内部各点的应力状态在物体占据的空间内组成的总体称为应力场。“场”是分布在空间区域内的物理量或数学函数。
物体内各点应力状态相同的应力场是均匀应力场，否则构成非均匀应力场。
由于地壳内部上覆岩石的压力σh=ρgh （式中ρ是岩石密度，g是重力加速度，h是距地表的深度）随深度而变化，以及地壳岩石的非均匀性，地壳中不存在理想的均匀应力场。
应力迹线：主应力迹线又称应力轨迹，或应力网络，它定性地表示了主应力和最大剪应力的作用方位。应力轨迹可通过光弹模拟实验和计算方法求得。
应力扰动（或应力集中）：由于岩块或地块内部的局部不均匀性和不连续性，在岩体内部造成的应力场的局部变化。包括应力迹线的偏移和应力值的局部集中或变异。（在单轴压缩的情况下，在圆孔附近可以产生拉张应力，也可以出现单向均匀压应力三倍的应力）。
第三章 常用传感器 一、知识要点及要求 (1)掌握常用传感器的分类方法； (2)掌握常用传感器的变换原理； (3)了解常用传感器的主要特点及应用。 二、重点内容及难点 (一)传感器的定义、作用与分类 1、定义：工程上通常把直接作用于被测量，能按一定规律将其转换成同种或别种量值输出的器件，称为传感器。 2、作用：传感器的作用就是将被测量转换为与之相对应的、容易检测、传输或处理的信号。 3、分类：传感器的分类方法很多，主要的分类方法有以下几种： (1)按被测量分类，可分为位移传感器、力传感器、温度传感器等； (2)按传感器的工作原理分类，可分为机械式、电气式、光学式、流体式等； (3)按信号变换特征分类，可概括分为物性型和结构型； (4)根据敏感元件与被测对象之间的能量关系，可分为能量转换型与能量控制型； (5)按输出信号分类，可分为模拟型和数字型。 (二)电阻式传感器 1、分类：变阻式传感器和电阻应变式传感器。而电阻应变式传感器可分为金属电阻应变片式与半导体应变片两类。 2、金属电阻应变片式的工作原理：基于应变片发生机械变形时，其电阻值发生变化。金属电阻应变片式的的灵敏度v S g 21+=。 3、半导体电阻应变片式的工作原理：基于半导体材料的电阻率的变化引起的电阻的变化。半导体电阻应变片式的的灵敏度E S g λ=。 (三)电感式传感器 1、分类：按照变换原理的不同电感式传感器可分为自感型与互感型。其中自感型主要包括可变磁阻式和涡电流式。 2、涡电流式传感器的工作原理：是利用金属体在交变磁场中的涡电流效应。 (四)电容式传感器 1、分类：电容式传感器根据电容器变化的参数，可分为极距变化型、面积变化型、介质变化型三类。 2、极距变化型：灵敏度为201δ εεδA d dC S -==，可以看出，灵敏度S 与极距平方成反比，极距越小灵敏度越高。显然，由于灵敏度随极距而变化，这将引起非线性误差。 3、面积变化型：灵敏度为常数，其输出与输入成线性关系。但与极距变化型相比，灵敏度较低，适用于较大直线位移及角速度的测量。 4、介质变化型：可用来测量电介质的液位或某些材料的厚度、湿度和温度等；也可用于测量空气的湿度。 (五)压电式传感器 1、压电传感器的工作原理是压电效应。
精心整理第一章工程概况
根据**院提供的岩土工程勘察报告，该场地工程地质条件如下：
三、检测桩位示意图 四、钢筋应力计在桩身埋设位置示意图 钢筋应力计在各试桩中位置示意图
二、测试设备及钢筋测力计的埋设 1、每桩钢筋应力计设置在各土层交界面处，每一个截面设2只钢筋测力计（基本呈180°对称布置），各钢筋应力计埋设截面的平、剖面图如前图； 2、JTM-V1000振弦式钢筋应力计采用焊接法固定在钢筋笼主筋上，并与桩身纵轴线平行；
3、连接在应力计的电缆线用柔性材料保护，绑扎在钢筋笼内侧并 引至地面； 4、所有应力计均用明显标记编号； 5、仪器设备：检测仪器设备采用JTM-V1000振弦式钢筋应力计、JTM-V10B 型频率读数仪、集线箱等组成。 三、测试原理 1位2ε c1j = εεs1j 3E cj 、E sj —砼弹性模量、钢筋弹性模量[E s 取2.0×108(kPa)] A cj 、A sj —同一截面处砼面积、钢筋总面积。 εcj 、εsj —同一截面处砼与钢筋的应变 4、钢筋应力计受力的计算公式： ) 2()(' 2 02 ----------------??=-?=Si Sij S i ij Sij A E F F k P ε
式中： P Sij —第i 量测截面处在j 级荷载下应力计所受轴向力（kN ） F ij —第i 量测截面处在j 级荷载下应力计的实测频率值(Hz) F i0—i 截面处钢筋应力计的初始频率值（Hz ） K A si ’—56f ij P ij —i A i 12、弦式钢筋应力计宜放在两种不同性质土层的界面处，以测量桩在不同土层中的分层摩阻力。在地面处(或以上)应设置一个测量断面作为钢筋应力计传感器标定断面。钢筋应力计埋设断面距桩顶和桩底的距离不宜小于1倍桩径。在同一断面处对称设置2个钢筋应力计。钢筋计应按主筋直径大小选择。仪器的可测频率范围应大于桩在最大加载时的频率的1.2倍； 3、使用前应对钢筋计逐个标定，得出压力(拉力)与频率之间的关系。带有接长 ) 3()(' -------------------------?= Si S Sij Sij A E P ε
地质构造应力场分析方法与原则 发表时间：2019-01-04T10:34:05.383Z 来源：《基层建设》2018年第34期作者：郭建锐[导读] 摘要：构造应力场是地球动力学重要组成部分，是地壳动力学的主体部分，其研究对于构造分析研究、地震分析预报、工程抗震等领域都有着十分重要的理论和实践意义。 赤峰市利拓矿业有限公司内蒙古赤峰市 024000摘要：构造应力场是地球动力学重要组成部分，是地壳动力学的主体部分，其研究对于构造分析研究、地震分析预报、工程抗震等领域都有着十分重要的理论和实践意义。本次研究针对地质构造应力场的测量方法水力压裂法、井壁崩落法、磁组构法进行分析，并对地质构造应场力分析原则进行阐述，继而进一步丰富构造应力场的理论。 关键词：地质构造；构造应场力；应场力引言：构造应力场就是在一个空间范围内构造应力的分布。构造应力场是作用在地壳某一地区内部的和由于这一地区某种变形的构造单元的发育而出现的应力总和。应力场是一种物理场，它和其他物理场，如重力场、电滋场、位势场等一样，也是物质存在的一种形式。场不是空间，而是在空间范围内某个物理量的按势分布。随着时间的变化，场内各点的强度和方向也将发生变化。构造应力场是地球动力学重要组成部分，是地壳动力学的主体部分，其研究对于构造分析研究、地震分析预报、工程抗震等领域都有着十分重要的理论和实践意义。 1.地质构造应力场概述 构造应力场概念是由我国地质学家李四光率先提出的。1947年李四光提出用构造形迹反推构造应力场，并研究各种不同力学性质的构造形迹与应力方向、应力作用方式之间的相互关系。1940年格佐夫斯基也提出研究构造应力场，并把用赤平投影求主应力轴方向的方法引进构造应力场的研究。1950年一1996年国内外地质工作者结合地震地质的研究工作开展了构造应力测量，经多年努力，通过野外与室内实测证实了构造应力的存在，并探索、研究了行之有效的构造应力测量技术方法，完善了构造应力测量的理论基础，建立了可靠的测量技术方法和数据处理系统。万天丰（1999）、武红岭（1999，2003）等将矿场构造应力场研究的方法延伸到盆地构造研究领域，取得了人量的研究认识和资料，极大地丰富了构造应力场研究理论，也为盆地构造应力场研究积累了丰富的地质认识和方法。1970年构造应力场的研究有长足进展，逐渐深入到地质学的多个领域。1980年以后，构造应力场问题越来越受到国内外地质学界的重视，研究内容多涉及板块、大陆，大洋地区的构造应力场。1990年以来，全球大陆与海洋科学钻探计划开始研究现今构造应力和古应力状态和岩石圈动力学问题。 2.地质构造应力场分析方法 构造应力场研究的主要内容是在确定各地的点应力状态（应力方向和应力大小）的基础上，研究在一定区域范围内各个构造活动时期的构造应力分布特征。古应力测量可通过构造形迹分析法、古地磁法、节理测量法来确定古构造应力作用方向，利用声发射法。晶格位错法等可确定古地应力值的大小（导致地层变形时的最大水平古应力）。现今应力测量可利用震源机制解法、水力压裂法、井壁崩落法等来确定现今构造应力最大主应力方向，利用声发射法、经验公式法可确定现今地应力大小。 2.1.1水力压裂法 水力压裂测量地应力的方法首先在美国发展起来，1977年B.Haimson在井深5.1Km处进行了水力压裂地应力测量。我国学者葛洪魁（1998）、康红普（2014）均在研究中采用水力压裂测量法进行验证。水力压裂（Hydraulic fracturing）地应力测量是通过在井眼周围地层中诱发人工裂缝来获取地应力的一种方法，测试精度受多种因素的影响，如测试层位筛选、施工仪器设备、施工方案的选择以及测试数据的分析等。 2.1.2井壁崩落法 井壁崩落椭圆法的理论依据为崩落椭圆是由地壳内的构造应力场形成的，所以二者之间存在确定的关系。它的基本原理是，由于地壳内存在水平差应力，致使钻井壁形成应力集中，在垂直于最大水平主应力（压应力为正）方向的井壁端切向应力最大，当该处切向应力达到或超过岩石的破裂极限强度时，即发生破裂，从而形成井壁崩落椭圆。1970年加拿大Bell在研究阿尔伯达油田四臂井径测量的地层倾角测井资料后，发现井眼扩大方向与区域内的最小水平主应力方向平行，Gough等也发现了这种现象。1985年，Zoboek使用井下电视观测证实了Boll的发现，并与B.Haimson等人对井眼崩落机制进行研究，说明了井壁崩落法是测量水平主应力方向的可行方法。shulnberger测井公司研究应用测井资料解释地层压力问题，并用于解释石油工程中的地层破裂压力、地层坍塌压力及油层出砂等问题。这种用测井资料解释地应力剖面的方法，己经能够解决石油工程中的许多问题。 2.1.3磁组构测量法 磁组构是指磁性颗粒或晶格的定向排列或组合，其实质是岩石磁化率各向异性。岩石磁化率各向异性是指岩石的磁化强度随方向的变化性质，包括感应磁化率各向异性与剩余磁化率各向异性。GrahamJ.w（1954）提出，儿乎所有岩石都可以观测到磁各向异性。研究表明，岩石的磁化率一般表现为磁化率数量椭球的形状和方向。椭球可以反映岩石内部铁磁性颗粒长轴的主要分布方向，与沉积搬运和充填方式、岩浆岩流动构造、变质岩类型和变质程度、页理、线理、褶皱轴方向等存在一定对应关系，是地史时期定向应力和温度作用的结果，是岩组分析和有限应变测量的重要手段之一。 3.地质构造应力场分析原则 3.1时间局限性原则 一般认为根据不同构造形变的切错和叠加等关系可以确定构造应力场的分期，即相对活动次序。可以根据组成构造形变的最新地层时代和角度不整合面之上的最老上覆地层的时代，来确定构造应力场作用的大致时间。如果有地层或侵入体同位素年代的资料时，构造应力作用的时间可以确定得更准确些。即使如此，构造应力作用的时间还是不可能确定得十分精确。 如果已知组成某一构造形变的最新地层年代和侵蚀了构造形变的不整合面之上的最老上覆地层的年代，构造形变肯定是在不整合形成期间发生的；但两个沉积地层的年代之间，发生了许多变化：老地层沉积之后要下沉、硬结成岩；受构造应力作用后造成构造形变；隆起遭受剥蚀；地壳重新下降，接受新的沉积。可以看出在整个不整合的形成过程中造成构造形变的构造应力作用只局限在一个较短的时间内。如果再考虑到同位素年代的不精确性（由于采样、测试方法等原因），要准确测定构造应力作用的时间实际上目前还难以实现。 3.2空间动态性原则
管道应力分析基础知识 2009-04-09 13:55 1. 进行应力分析的目的是 1) 使管道应力在规范的许用范围内； 2) 使设备管口载荷符合制造商的要求或公认的标准； 3) 计算出作用在管道支吊架上的荷载； 4) 解决管道动力学问题； 5) 帮助配管优化设计。 2. 管道应力分析主要包括哪些内容？各种分析的目的是什么？ 答：管道应力分析分为静力分析和动力分析。 1) 静力分析包括： (l)压力荷载和持续荷载作用下的一次应力计算――防止塑性变形破坏； (2)管道热胀冷缩以及端点附加位移等位移荷载作用下的二次应力计算――防止疲劳破坏； (3)管道对设备作用力的计算――防止作用力太大，保证设备正常运行； (4)管道支吊架的受力计算――为支吊架设计提供依据； (5)管道上法兰的受力计算――防止法兰泄漏； (6)管系位移计算――防止管道碰撞和支吊点位移过大。 2) 动力分析包括： (l)管道自振频率分析――防止管道系统共振； (2)管道强迫振动响应分析――控制管道振动及应力； (3)往复压缩机气柱频率分析――防止气柱共振； (4)往复压缩机压力脉动分析――控制压力脉动值。 3. 管道应力分析的方法 管道应力分析的方法有：目测法、图表法、公式法、和计算机分析方法。选用什
么分析方法，应根据管道输送的介质、管道操作温度、操作压力、公称直径和所连接的设备类型等设计条件确定。 4. 对管系进行分析计算 1) 建立计算模型（编节点号），进行计算机应力分析时，管道轴测图上需要提供给计算机软件数据的部位和需要计算机软件输出数据的部位称作节点： (1)管道端点 (2)管道约束点、支撑点、给定位移点 (3)管道方向改变点、分支点 (4)管径、壁厚改变点 (5)存在条件变化点（温度、压力变化处） (6)定义边界条件（约束和附加位移） (7)管道材料改变处（包括刚度改变处，如刚性元件） (8)定义节点的荷载条件（保温材料重量、附加力、风载、雪载等） (9)需了解分析结果处(如跨距较长的跨中心点) (10) 动力分析需增设点 2) 初步计算（输入数据符合要求即可进行计算） (1) 利用计算机推荐工况（用CASWARII计算，集中荷载、均布荷载特别加入） (2) 弹簧可由程序自动选取 (3) 计算结果分析 (4) 查看一次应力、二次应力的核算结果 (5) 查看冷态、热态位移 (6) 查看机器设备受力 (7) 查看支吊架受力（垂直荷载、水平荷载） (8) 查看弹簧表
压力容器分析报告
目录 1 设计分析依据 0 1.1 设计参数 0 1.2 计算及评定条件 0 1.3 材料性能参数 0 2 结构有限元分析 (1) 2.1 理论基础 (1) 2.2 有限元模型 (1) 2.3 划分网格 (1) 2.4 边界条件 (2) 3 应力分析及评定 (2) 3.1 应力分析 (2) 3.2 应力强度校核 (2) 4 分析结论 (3) 4.1 上封头接头外侧 (4) 4.2 上封头接头内侧 (5) 4.3 上封头壁厚 (7) 4.4 筒体上 (9) 4.5 筒体左 (10) 4.6 下封头接着外侧 (12) 4.7 下封头壁厚 (14)
1 设计分析依据 （1）压力容器安全技术监察规程 （2）JB4732-1995 《钢制压力容器-分析设计标准》-2005确认版 1.1 设计参数 表1 设备基本设计参数 1.2 计算及评定条件 （1）静强度计算条件 表2 设备载荷参数 载荷进行计算，故采用设计载荷进行强度分析结果是偏安全的。 1.3 材料性能参数 材料性能参数见表3，其中弹性模型取自JB4732-95表G-5，泊松比根据JB4732-95的公式（5-1）计算得到，设计应力强度分别根据JB4732-95的表6-2、表6-4、表6-6确定。 表3 材料性能参数性能
2 结构有限元分析 2.1 理论基础 传统的压力容器标准与规范，一般属于“常规设计”，以弹性失效准则为理论基础，由材料力学方法或经验得到较为简单的适合于工程应用的计算公式，求出容器在载荷作用下的最大主应力，将其限制在许用值以内，即可确认容器的壁厚。对容器局部区域的应力、高应力区的应力不做精细计算，以具体的结构形式限制，在计算公式中引入适当的系数或降低许用应力等方法予以控制，这是一种以弹性失效准则为基础，按最大主应力理论，以长期实践经验为依据而建立的一类标准。 塑性理论指出，由于弹性应力分析求得的各类名义应力对结构破坏的危险性是不同的，随着工艺条件的苛刻和容器的大型化，常规设计标准已经不能满足要求，尤其是在应力集中区域。若不考虑应力集中而只按照简化公式进行设计，不是为安全而过分浪费材料就是安全系数不够。基于各方面的考虑，产生了“分析设计”这种理念。采用以极限载荷、安定载荷和疲劳寿命为界限的“塑性失效”与“弹性失效”相结合的“弹塑性失效”准则，要求对容器所需部位的应力做详细的分析，根据产生应力的原因及应力是否有自限性，分为三类共五种，即一次总体薄膜应力( Pm) 、一次局部薄膜应力( Pc) 、一次弯曲应力( Pb) 、二次应力( Q) 和峰值应力( F) 。 对于压力容器的应力分析，重要的是得到应力沿壁厚的分布规律及大小，可采用沿壁厚方向的“校核线”来代替校核截面。而基于弹性力学理论的有限元分析方法，是一种对结构进行离散化后再求解的方法，为了获得所选“校核线”上的应力分布规律及大小，就必须对节点上的应力值进行后处理，即应力分类，根据对所选“校核线”上的应力进行分类，得出各类应力的值，若满足强度要求，则所设计容器是安全的。 按照JB4732-1995进行分析，整个计算采用ANSYS13.0软件，建立有限元模型，对设备进行强度应力分析。 2.2 有限元模型 由于主要关心容器开孔处的应力分布规律及大小，为减少计算量，只取开孔处作为分析对象，且取其中较为关心的大孔进行分析校核。分析设计所用的几何模型如图1所示。在上下封头和筒体之间存在不连续的壁厚，由于差距和影响量较小，此处统一采用上下封头的设计厚度。 图1 压力容器模型 2.3 划分网格 在结构的应力分析中，采用ANSYS13.0中的solid187单元进行六面体划分，如图2所示。图3~图5
固体中的应力波 李清 中国矿业大学（北京）
参考书： 1 王礼立. 《应力波基础》第2版(2005年8月1日)，国防工业出版社 2 李玉龙. 《应力波基础简明教程》第1版 (2007年4月1日)，西北工业大学 3 丁启财(美国). 《固体中的非线性波》，中国友谊出版公司 4 宋守志. 《固体中的应力波》，煤炭工业出版社 5 杨善元. 《岩石爆破动力学基础》，煤炭工业出版社 6 莱茵哈特(杨善元译). 《固体中的应力瞬变》，煤炭工业出版社 7 徐小荷. 《冲击凿岩的理论基础与电算方法》，东工出版社 8 郭自强. 《固体中的波》，地震出版社
目录 第0章绪论 (1) 1 波动现象 (1) 2 应力波的概念 (1) 3 应力波分类 (3) 4 应力波理论与其它力学理论的关系 (3) 5 应力波理论的发展 (3) 6 应力波理论在岩土工程中的应用 (3) 第1章一维应力波基础 (4) §1.1波动方程及其解 (4) 1.1.1 一维纵波的波动方程 (4) 1.1.2 波的传播速度 (4) 1.1.3 波动方程的解 (5) 1.1.4 解的物理意义 (6) §1.2 应力波的几个基本参量 (7) §1.3 应力波的能量 (7) §1.4 波的衰减 (8) 1.4.1 原因 (8) 1.4.2 度量 (8) 1.4.3 衰减率α的测定 (9) §1.5 考虑杆的横向效应的波动方程 (10) §1.6 杆中的扭转波与弯曲波 (12) 1.6.1 扭转波 (12) 1.6.2 弯曲波 (13) 第2章二维和三维弹性波理论基础 (14) §2.1 弹性体的运动微分方程 (14) §2.2 弹性体的无旋波与等容波 (15) 2.2.1 无旋波(纵波、P波) (15) 2.2.2 等容波(横波、S波) (16) §2.3 平面波的传播 (17) 2.3.1 平面纵波(V//c) (17) 2.3.2 平面横波(V⊥c) (18) §2.4 薄板中的应力波 (19) 2.4.1 控制方程 (19) 2.4.2 纵波 (20) 2.4.3 横波 (21) 2.4.4 各种波速关系 (21) §2.5 球面波 (22) 2.5.1 波动方程及其解 (22) §2.6 柱面波 (23) 第3章应力波的相互作用 (24) §3.1 一维应力波在界面的反射和透射 (24) 3.1.1 应力波在不同介质界面的反射和透射 (25) 3.1.2应力波在变截面杆中的反射和透射 (26)
《有限元分析》报告基本要求： 1. 以个人为单位完成有限元分析计算，并将计算结果上交；（不允许出现相同的分析模型，如相 同两人均为不及格） 2. 以个人为单位撰写计算分析报告； 3. 按下列模板格式完成分析报告； 4. 计算结果要求提交电子版，报告要求提交电子版和纸质版。（以上文字在报告中可删除） 《有限元分析》报告 一、问题描述 （要求：应结合图对问题进行详细描述，同时应清楚阐述所研究问题的受力状况和约束情况。图应清楚、明晰，且有必要的尺寸数据。） 一个平面刚架右端固定，在左端施加一个y 方向的-3000N 的力P1，中间施加一个Y 方向的-1000N 的力P2，试以静力来分析，求解各接点的位移。已知组成刚架的各梁除梁长外，其余的几何特性相同。 横截面积：A=0.0072 m2 横截高度：H=0.42m 惯性矩：I=0.0021028m4ｘ 弹性模量： E=2.06ｘ10n/ m2/ 泊松比：u=0.3 二、数学模型 （要求：针对问题描述给出相应的数学模型，应包含示意图，示意图中应有必要的尺寸数据；如进行了简化等处理，此处还应给出文字说明。） （此图仅为例题）
三、有限元建模（具体步骤以自己实际分析过程为主，需截图操作过程） 用ANSYS 分析平面刚架 1．设定分析模块 选择菜单路径：MainMenu—preference 弹出“PRreferences for GUI Filtering”对话框，如图示，在对话框中选取：Structural”,单击[OK]按钮，完成选择。 2．选择单元类型并定义单元的实常数 （1）新建单元类型并定 （2）定义单元的实常数在”Real Constants for BEAM3”对话框的AREA中输入“0。0072”在IZZ 中输入“0。0002108”，在HEIGHT中输入“0.42”。其他的3个常数不定义。单击[OK]按 钮，完成选择 3．定义材料属性 在”Define Material Model Behavier”对话框的”Material Models Available”中，依次双击“Structural→Linear→Elastic→Isotropic”如图
第三章 弹性波的相互作用 3-3 已知两种材质的弹性杆A 和B 的弹性模量、密度和屈服极限分别为： E A ＝60GPa ， ρA =2.4g/cm 3，Y A =120MPa ，E 1A =E A /5； E B ＝180GPa ，ρB =7.2g/cm 3，Y B =240MPa ，E 1B =E B /5。 试对Ⅵ－10所示四种情况分别画出X －t 及σ—v 图，并确定撞击结束时间、两杆脱开时间以及分离之后各自的整体飞行速度。 解：两种材料的参数计算如下： s m E C A A A /500010104.210606 39 0=???= = -ρ s m C C A A /10005/01==，s m C v A yA yA /105000 10004.210120)(6 00-=???-=-=ρσ s m E C B B B /500010102.7101806 39 0=???= = -ρ s m C C B B /10005/01==，s m C v B yB yB /667.65000 10002.710240)(6 00-=???-=-=ρσ A C )(00ρ＝2.4×10-3×106×5000＝12×106kg/（sm ） B C )(00ρ＝7.2×10-3×106×5000＝36×106kg/（sm ） （1）： A B v . 由上图可知：当左杆波从自由端反射至接触面时，速度，为－4m/s ，应力为0，撞击结束。 撞击结束时间：0.02μs 。
两杆脱开时间即接触到脱开时间：0.02μs。 短杆整体飞行速度：－4 m/s（3区）。 长杆整体飞行速度：2m/s（5区速度）。 （2） 撞击结束应在A点。 撞击结束时间：0.04μs。 两杆脱开时间即接触到脱开时间：0.04μs。 短杆整体飞行速度：2 m/s（7区）。 长杆整体飞行速度：9m/s（6，10区）。 （3） v
XXX球罐应力分析报告 设备名称：XXX球罐 设备位号：XXX 应力分析报告
目录 1基本设计参数 (4) 2计算数据 (6) 2.1 计算条件 (6) 2.2材料性能数据 (7) 3主要受压元件计算 (8) 4整体结构分析计算 (9) 4.1 力学模型和有限元模型 (9) 4.2 载荷工况分析 (11) 4.3 载荷边界条件 (12) 4.4 位移边界条件 (15) 4.5 应力强度分布云图及路径选取 (15) 4.6 应力线性化及强度评定 (20) 4.7 整体结构强度评定汇总 (33) 5局部结构分析计算 (34) 5.1 人孔与接管N1/N4局部结构分析 (34) 5.1.1 力学模型和有限元模型 (34) 5.1.2载荷边界条件 (36) 5.1.3位移边界条件 (38) 5.1.4应力分布云图及路径选取 (39) 5.1.5 应力线性化及强度评定 (40) 5.1.6 人孔与接管N1/N4应力线性化及强度评定 (48) 5.2 人孔与接管V1/K3/K4局部结构分析 (48) 5.2.1 力学模型和有限元模型 (48) 5.2.2载荷边界条件 (51) 5.2.3位移边界条件 (53) 5.2.4应力分布云图及路径选取 (54) 5.2.5 应力线性化及强度评定 (55)
5.2.6 人孔与接管V1/K3/K4应力线性化及强度评定 (63) 5.3 人孔与接管K1/K2局部结构分析 (63) 5.3.1 力学模型和有限元模型 (63) 5.3.2载荷边界条件 (66) 5.3.3位移边界条件 (68) 5.3.4应力分布云图及路径选取 (69) 5.3.5 应力线性化及强度评定 (70) 5.3.6 人孔与接管K1/K2应力线性化及强度评定 (78) 5.4 人孔与接管N2局部结构分析 (78) 5.4.1 力学模型和有限元模型 (78) 5.4.2载荷边界条件 (81) 5.4.3位移边界条件 (83) 5.4.4应力分布云图及路径选取 (84) 5.4.5 应力线性化及强度评定 (85) 5.4.6 人孔与接管N2应力线性化及强度评定 (93) 5.5 人孔与接管N5局部结构分析 (93) 5.5.1 力学模型和有限元模型 (93) 5.5.2载荷边界条件 (96) 5.5.3位移边界条件 (99) 5.5.4应力分布云图及路径选取 (100) 5.5.5 应力线性化及强度评定 (101) 5.5.6 人孔与接管N5应力线性化及强度评定 (109) 6结论 (109) 附录 (109) 球罐SW6计算文件
《应力场分析与裂缝预测》教学大纲 （2004年制定，2012年第二次修改） 课程名称：应力场分析与裂缝预测 课程英文名称：Stress Field Analysis and Fracture Prediction 课内学时：32 课程学分：2 课程性质：学位课开课学期：每学年第一学期 教学方式：课堂讲授考核方式（考试/考查）：考试 大纲执笔人：曾联波主讲教师：曾联波 师资队伍：曾联波、童亨茂、陈书平 一、课程内容简介 《应力场分析与裂缝预测》是地质学专业和资源勘探与地质工程专业硕士研究生的一门专门课程。讲授古、现应力场和储层裂缝的研究方法及其在油气勘探与开发中的应用，包括应力与应力场的基础概念、古构造应力场分析方法、现今地应力测量方法、裂缝的基础知识，裂缝定量预测方法、古应力场在油气勘探中的应用、现今地应力和裂缝在低渗透油气田开发中的应用。本门课程为32学时，2学分。 二、课程目的和基本要求 课程的目的是培养学生掌握古、今应力场分析与储层裂缝预测的基本理论和方法分析油田应力场分布及进行储层裂缝预测的基本能力。《应力场分析与裂缝预测》课程涉及构造地质学、地质力学、储层地质学、岩石力学、石油地质学和油气藏工程等多方面的基本知识，要求学生系统学习了大学本科地质类专业的构造地质学、固体力学、石油地质学和储层地质学等课程。 学完本课程后，应达到以下基本要求： 1．了解应力、应力场和裂缝的基本概念及基本分布特征； 2．掌握古应力场研究方法及进展，并能运用这些基本方法分析油田古应力场分布和指导油气勘探； 3．掌握现今地应力测量方法，并能运用这些方法分析低渗透油气田的地应力分布和指导油气田开发。 4. 掌握储层裂缝的研究和预测方法，并能运用这些方法研究和预测低渗透储层裂缝的分布规律。 三、课程主要内容 §1. 应力场分析和裂缝预测的基础知识（4学时） §1.1应力、应力场和裂缝的基本概念。 §1.2应力场和裂缝研究的基本内容与方法。 §1.3应力场分析和裂缝预测的研究现状与发展趋势。 §1.4应力场分析和裂缝预测的研究意义。 §2. 现今地应力测量方法（4学时） §2.1现场地应力测量方法。 §2.2岩心地应力测量方法。 §2.3测井地应力分析方法。 §2.4地应力的分布规律及影响因素 §3. 古构造应力场分析方法（6学时） §3.1古构造应力方向分析方法。
ANSYS基础教程——应力分析 关键字：ANSYS 应力分析 ANSYS教程 信息化调查找茬投稿收藏评论好文推荐打印社区分享 应力分析是用来描述包括应力和应变在的结果量分析的通用术语，也就是结构分析，应力分析包括如下几个类型:静态分析瞬态动力分析、模态分析谱分析、谐响应分析显示动力学，本文主要是以线性静态分析为例来描述分析，主要容有：分析步骤、几何建模、网格划分。 应力分析概述 ·应力分析是用来描述包括应力和应变在的结果量分析的通用术语，也就是结构分析。 ANSYS 的应力分析包括如下几个类型: ●静态分析 ●瞬态动力分析 ●模态分析 ●谱分析 ●谐响应分析 ●显示动力学 本文以一个线性静态分析为例来描述分析步骤,只要掌握了这个分析步骤，很快就会作其他分析。 A. 分析步骤 每个分析包含三个主要步骤:
·前处理 –创建或输入几何模型 –对几何模型划分网格 ·求解 –施加载荷 –求解 ·后处理 –结果评价 –检查结果的正确性 ·注意！ANSYS 的主菜单也是按照前处理、求解、后处理来组织的；
·前处理器(在ANSYS中称为PREP7)提供了对程序的主要输入； ·前处理的主要功能是生成有限元模型，主要包括节点、单元和材料属性等的定义。也可以使用前处理器PREP7 施加载荷。 ·通常先定义分析对象的几何模型。 ·典型方法是用实体模型模拟几何模型。 –以CAD-类型的数学描述定义结构的几何模型。 –可能是实体或表面，这取决于分析对象的模型。 B. 几何模型 ·典型的实体模型是由体、面、线和关键点组成的。 –体由面围成，用来描述实体物体。 –面由线围成，用来描述物体的表面或者块、壳等。 –线由关键点组成，用来描述物体的边。 –关键点是三维空间的位置，用来描述物体的顶点。
第五章 刚性卸载近似 5－2 有弹性――线性硬化材料的两有限长杆（长度为L ），其弹性波速C 0相等，塑性波速C 1分别为C 0/3和C 0/5。它们均以5倍的屈服速度分别撞击刚性靶，如图Ⅳ－34所示。试分别画出x －t 图和v －σ图，并确定这两根杆脱离靶板的时间。 解：图解如下： （1），塑性波速C 1为C 0/3。 σ y 由图中看出，v －σ图上的17点已为拉应力，故应脱靶。由x －t 图可得其脱靶时间为3.5L/C 0。 （2）、 塑性波速C 1为C 0/5。 y X 图Ⅳ－34 长为v
由图中看出，v －σ图上的17已为拉应力，故应脱靶。由x －t 图可得其脱靶时间为38L/（9C 0）。 5－5半无限长杆的材料为弹性－线性硬化材料，其弹性波速C 0和塑性波速C 1均已知，且C 1＝C 0/10。若在杆端作用一如图Ⅳ－35所示的应力载荷)(t σ，试采用刚性卸载近似来确定杆中残余应变段的长度。 解：如图所示。 ∵ ??? ? ??? -==-C d dv d dv x m m t m u m 000ρσρσσ ∴ 0)(0=-+u m m dt d C X σσσ 又： )(m l t C X στ-= ∴ 0)()] ([0=-+-u m m m l dt d t σσσστ 在本题中，残余应变段的长度尽头的应力应等于Y 。 由公式（5－27）： )]([)(00 11 max t t t t t t t t m --- = σσ 在本题中：01max 2,4t t Y ==σ Y t t t t Y t m =--= )](2[4)(00 σ 解得：002)237(4 1 t t t >+= ∴ 应采用公式（5－28）： )()(1020 max t t t t t t m --= σσ 在本题中：01max 2,4t t Y ==σ Y t t t t Y t m =--= )2(4)(2 020σ 解得：08 33t t = 代入公式（5－26）得： 1021)(t t t C t X -==? 的应力载荷 O O 图Ⅳ－35 作用于杆端
管道应力分析概述 CAESARII软件介绍 CAESARII管道应力分析软件是由美国COADE公司研发的压力管道应力分析专业软件。它既可以分析计算静态分析，也可进行动态分析。CAESARII向用户提供完备的国际上的通用管道设计规范，使用方便快捷。交互式数据输入图形输出，使用户可直观查看模型（单线、线框，实体图）强大的3D计算结果图形分析功能，丰富的约束类型，对边界条件提供最广泛的支撑类型选择、膨胀节库和法兰库，并且允许用户扩展自己的库。钢结构建模，并提供多种钢结构数据库.结构模型可以同管道模型合并,统一分析膨胀节可通过标准库选取自动建模、冷紧单元/弯头，三通应力强度因子（SIF）的计算、交互式的列表编辑输入格式用户控制和选择的程序运行方式,用户可定义各种工况。 一、管道应力分析的原则 管道应力分析应保证管道在设计条件下具有足够的柔性，防止管道因热胀冷缩、管道支承或端点附加位移造成应力问题。 二、管道应力分析的主要内容 管道应力分析分为静力分析和动力分析。 静力分析包括： 1）压力荷载和持续荷载作用下的一次应力计算——防止塑性变形破坏； 2）管道热胀冷缩以及端点附加位移等位移荷载作用下的二次应力计算——防止疲劳破坏； 3）管道对设备作用力的计算——防止作用力太大，保证设备正常运行； 4）管道支吊架的受力计算——为支吊架设计提供依据； 5）管道上法兰的受力计算——防止法兰汇漏。 动力分析包括：
l）管道自振频率分析——防止管道系统共振； 2）管道强迫振动响应分析——控制管道振动及应力； 3）往复压缩机（泵）气（液）柱频率分析——防止气柱共振； 4）往复压缩机（泵）压力脉动分析——控制压力脉动值。 三、管道上可能承受的荷载 （1）重力荷载：包括管道自重、保温重、介质重和积雪重等； （2）压力荷载：压力载荷包括内压力和外压力； （3）位移荷载：位移载荷包括管道热胀冷缩位移、端点附加位移、支承沉降等； （4）风荷载； （5）地震荷载； （6）瞬变流冲击荷载：如安全阀启跳或阀门的快速启闭时的压力冲击： （7）两相流脉动荷载； （8）压力脉动荷载：如往复压缩机往复运动所产生的压力脉动； （9）机械振动荷载：如回转设备的振动。 四、管道应力分析的目的 1）为了使管道和管件内的应力不超过许用应力值； 2）为了使与管系相连的设备的管口荷载在制造商或国际规范（如 NEMA SM-23、API-610、API-6 17等）规定的许用范围内； 3）为了使与管系相连的设备管口的局部应力在 ASME Vlll的允许范围内； 4）为了计算管系中支架和约束的设计荷载；
科技信息2011年第27期 SCIENCE &TECHNOLOGY INFORMATION 华北地区构造应力场研究 李富涛1孟昭焕2贾宝刚1 (1.山东省煤田地质局物探测量队山东泰安271021；2.莱芜市国土资源局山东莱芜271100) 【摘要】本文结合相关数据、模型和软件分别利用重力场、重力垂线偏差与构造应力场的内在关系式对华北地区陆地构造应力值进行了计算，通过对相关数据结果进行对比分析，总结并得出了华北地区重力总水平梯度、构造应力场和研究方法本身的一些规律和特征。 【关键词】重力场；重力水平梯度；垂线偏差；构造应力 0引言构造应力场是地球动力学研究领域一个重要的组成部分。由于我 们不能直接测量得到浅层地表以外的岩石圈构造应力场，一些学者于 是另辟蹊径，以可以直接测量得到的相关区域重力数据为参考，通过 研究构造应力场与重力场之间的内在关系的方法而最终获得构造应 力场数据。在这方面，典型的代表人物有游永雄、向文、方剑等。游永雄 曾利用重力场研究了包括华北地区在内的多个地区的构造应力场情 况，本文即利用近似方法专门针对华北地区东经[106°，124°]、北纬[31° 43°]范围的大陆构造应力场进行更加细致地研究[1]，以期使得对该区 域构造应力场及其变化规律和研究手段本身认识得更加详尽。 1 重力和垂线偏差场转换构造应力场公式1.1利用重力场计算华北地区构造应力场 游永雄推导了重力场转换构造应力场的公式即[2]：Δσxx ＝g ρx ，y m g x （1）其中，Δσxx 代表构造应力；g 为正常重力；f 为引力常量；ρx ，y 为均衡 改正的单位均衡柱体密度；ρm 为地幔密度；g x 为重力总水平梯度，其水 平分量Δg x 和Δg y 值可用下面公式计算[2][3]： Δg x ＝－1＋∞－∞乙＋∞－∞乙（x－x＇）Δg z ［（x－x＇）2＋（y－y＇）2＋H 2］32dx＇dy 'Δg y ＝－12π＋∞－∞乙＋∞－∞乙（y －y ＇）Δg z ［（x－x＇）2＋（y－y＇）2＋H 2］ 32dx＇dy 乙乙乙乙乙乙乙乙乙乙乙'（2）Δg z 是得到的重力异常值，x＇和y '是流动坐标，遍及整个测量区 域，H 是空间延拓高度， 积分面积可以有限化和离散化，以适应计算，本文即以离散化后 2度范围为积分区域来计算。 求g x 的值的计算式为：g x =(Δg x )2+(Δg y )2 姨（3）1.2利用重力垂线偏差计算构造应力场公式 利用垂线偏差计算构造应力场公式如下[2]： Δσxx ＝-g 24πf ρx ，y ·u ρm ·ρ （4）式中，u 代表重力垂线偏差；ρ=206265rad ·s 。 u 的值根据下式计算[4]：u =（ξ2+η2） 1/2（5）其中，ξ为南北垂线偏差（垂线偏差子午圈分量）；η为西东垂线偏 差（垂线偏差卯酉分量）。 2 计算华北地区构造应力场2.1利用重力场计算华北地区构造应力场 本文利用华北地区5′×5′分辨率的DTM 数据、360阶重力场模型 EGM96并借助于PALGrav1.0软件[5]求得该区域布格重力异常值Δg z ， 然后计算得到重力总水平梯度g x 。在此基础上，再利用重力延拓知识[6]， 并根据式（1）分别计算得到了华北地区地表、20公里和40公里深度 处的构造应力值。以下分别是该区5′×5′分辨率DTM 图、重力总水平 梯度图和地表、20公里、40公里深度构造应力场图。 2.2利用重力垂线偏差计算华北地区构造应力场 利用上述同样DTM 数据、重力场模型和软件计算南北垂线偏差 ξ和东西垂线偏差η，然后计算重力垂线偏差u 。根据公式（4）进一步 计算得到该区构造应力值。以下分别是利用垂线偏差计算得到的华北 地区重力总水平梯度图和地表构造应力场图。3分析和讨论 图1华北地区DTM 图（单位：m ）图 2 华北地区重力总水平梯度矢量图（单位：E ） 图3华北地区重力总水平梯度等值线图（单位：E ）图4华北地区地表构造应力场矢量图（单位：MPa ） ○科教前沿○
《化工设备机械基础》习题解答 第三章 内压薄壁容器的应力分析 一、名词解释 A 组： ⒈薄壁容器：容器的壁厚与其最大截面圆的内径之比小于的容器。 ⒉回转壳体：壳体的中间面是直线或平面曲线绕其同平面内的固定轴线旋转360°而成的壳体。 ⒊经线：若通过回转轴作一纵截面与壳体曲面相交所得的交线。 ⒋薄膜理论：薄膜应力是只有拉压正应力没有弯曲正应力的一种两向应力状态，也称为无力矩理论。 ⒌第一曲率半径：中间面上任一点M 处经线的曲率半径。 / ⒍小位移假设：壳体受力以后，各点位移都远小于壁厚。 ⒎区域平衡方程式：计算回转壳体在任意纬线上径向应力的公式。 ⒏边缘应力：内压圆筒壁上的弯曲应力及连接边缘区的变形与应力。 ⒐边缘应力的自限性：当边缘处的局部材料发生屈服进入塑性变形阶段时，弹性约束开始缓解，原来不同的薄膜变形便趋于协调，边缘应力就自动限制。 二、判断题（对者画√，错着画╳） A 组： 1. 下列直立薄壁容器，受均匀气体内压力作用，哪些能用薄膜理论求解壁内应力哪些不能 （1） 横截面为正六角形的柱壳。（×） （2） 横截面为圆的轴对称柱壳。（√） （3） 横截面为椭圆的柱壳。 （×） · （4） 横截面为圆的椭球壳。 （√） （5） 横截面为半圆的柱壳。 （×） （6） 横截面为圆的锥形壳。 （√） 2. 在承受内压的圆筒形容器上开椭圆孔，应使椭圆的长轴与筒体轴线平行。（×） 3. 薄壁回转壳体中任一点，只要该点的两个曲率半径R R 2 1 =，则该点的两向应力σσθ=m 。 （√） 4. 因为内压薄壁圆筒的两向应力与壁厚成反比，当材质与介质压力一定时，则壁厚大的容器， 壁内的应力总是小于壁厚小的容器。（×） 5. 按无力矩理论求得的应力称为薄膜应力，薄膜应力是沿壁厚均匀分布的。（√） B 组： 1. 卧式圆筒形容器，其内介质压力，只充满液体，因为圆筒内液体静载荷不是沿轴线对称分布 的，所以不能用薄膜理论应力公式求解。（√） 2. 由于圆锥形容器锥顶部分应力最小，所以开空宜在锥顶部分。（√） - 3. 凡薄壁壳体，只要其几何形状和所受载荷对称于旋转轴，则壳体上任何一点用薄膜理论应力 公式求解的应力都是真实的。（×） 4. 椭球壳的长，短轴之比a/b 越小，其形状越接近球壳，其应力分布也就越趋于均匀。（√） 5. 因为从受力分析角度来说，半球形封头最好，所以不论在任何情况下，都必须首先考虑采用 半球形封头。（×）
固体中的应力波
李清
中国矿业大学（北京）
参考书：
1 王礼立. 《应力波基础》第 2 版(2005 年 8 月 1 日)，国防工业出版社 2 李玉龙. 《应力波基础简明教程》第 1 版 (2007 年 4 月 1 日)，西北工业大学 3 丁启财(美国). 《固体中的非线性波》，中国友谊出版公司 4 宋守志. 《固体中的应力波》，煤炭工业出版社 5 杨善元. 《岩石爆破动力学基础》，煤炭工业出版社 6 莱茵哈特(杨善元译). 《固体中的应力瞬变》，煤炭工业出版社 7 徐小荷. 《冲击凿岩的理论基础与电算方法》，东工出版社 8 郭自强. 《固体中的波》，地震出版社
目录
第 0 章 绪论..................................................................................................................错误!未定义书签。 1 波动现象..............................................................................................................错误!未定义书签。 2 应力波的概念......................................................................................................错误!未定义书签。 3 应力波分类..........................................................................................................错误!未定义书签。 4 应力波理论与其它力学理论的关系..................................................................错误!未定义书签。 5 应力波理论的发展..............................................................................................错误!未定义书签。 6 应力波理论在岩土工程中的应用......................................................................错误!未定义书签。
第 1 章 一维应力波基础..............................................................................................错误!未定义书签。 §波动方程及其解 .................................................................................................... 错误!未定义书签。 一维纵波的波动方程 .................................................................................... 错误!未定义书签。 波的传播速度 ................................................................................................ 错误!未定义书签。 波动方程的解 ................................................................................................ 错误!未定义书签。 解的物理意义 ................................................................................................ 错误!未定义书签。 § 应力波的几个基本参量......................................................................................错误!未定义书签。 § 应力波的能量......................................................................................................错误!未定义书签。 § 波的衰减..............................................................................................................错误!未定义书签。 原因 ................................................................................................................ 错误!未定义书签。 度量 ................................................................................................................ 错误!未定义书签。 衰减率α的测定 ............................................................................................ 错误!未定义书签。 § 考虑杆的横向效应的波动方程..........................................................................错误!未定义书签。 § 杆中的扭转波与弯曲波......................................................................................错误!未定义书签。 扭转波 ............................................................................................................ 错误!未定义书签。 弯曲波 ............................................................................................................ 错误!未定义书签。
第 2 章 二维和三维弹性波理论基础..........................................................................错误!未定义书签。 § 弹性体的运动微分方程....................................................................................错误!未定义书签。 § 弹性体的无旋波与等容波..................................................................................错误!未定义书签。 无旋波(纵波、P 波) ......................................................................................错误!未定义书签。 等容波(横波、S 波).......................................................................................错误!未定义书签。 § 平面波的传播....................................................................................................错误!未定义书签。
平面纵波(V 错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指
定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定
书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书
签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。
错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错
误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!
未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未
指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指
定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定
书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书签。错误!未指定书