DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 2. 理论阐述 随机森林 [ 7 ] 指的是利用多棵决策树对样本训练集进行训练并预测的一种分类器，在机器学习中，随机森林是一个包含多个决策树的分类器，并且其输出的类别是由个别树输出的类别的众数决定。在预测的时候，输入一个样本，每个决策树都会给出一个预测结果，根据所有决策树的预测结果的众数来确定样本最终的类别值。在生成随机森林的过程中，每一棵决策树都是通过一个独立抽取的样本建立的，决策树之间的相关性会影响最终误差，因此要尽量保证决策树之间不相关。 随机森林构建的具体步骤如下 [ 8 ] ： 1) 输入参数：原始训练集样本的个数 N ，变量的数目 M ，节点分类变量 m ； 2) 用自助抽样法有放回地随机抽取 k 个自助样本集，根据自助样本集，构建 k 棵决策树，每次未被抽到的样本组成了 k 个袋外数据，即简称OOB； 3) 在训练的过程中，自助样本集逐渐生成一棵棵决策树，在每个节点处，随机选取 M 个特征值中的 m 个( m 小于 M )，随后从 m 个特征中选取一个特征进行分支，此过程要保证结点不纯度最小，直至每个节点不纯度达到最小； 4) 生成多个决策树分类器，输入需要预测的元组，每个决策树都给出了一个分类结果，根据所有决策树的分类结果得出最终的类别。 随机森林的优点如下 [ 9 ] ： 可以处理高维度的数据，并且不用进行指标的选择；模型训练完毕之后，也可以得出数据每个指标的重要程度；模型的训练速度很快，并且可以进行并行化操作；最小化了决策树之间的相关性，所以模型的分类精度较高，与决策树相比，可以产生高精度的分类器；在训练过程中可以很好地估计缺失值，因此缺失值对分类模型的准确率不会产生很大影响。 随机森林的缺点如下 [ 10 ] ： 对于噪声较大的数据进行分类时，随机森林可能会产生过拟合现象从而影响分类的结果；对于低维度数据，可能不能产生很好的分类。
system: [SYSTEM]

user: 4. 政策建议 基于本文的结论，并结合目前房地产上市公司的盈利能力现状，对房地产行业发展提出如下建议： 1) 坚持宏观调控方向 目前，在国家宏观政策调控的大背景下，房地产行业发展下行，房地产投资放缓，但是房地产作为刚性需求，政策因素稍有变动，房价仍然面临高速增长的趋势。所以，国家政府有关部门对于房地产行业的调控政策必须坚持实行，同时还要保证经济的持续平稳发展，从而有效抑制房地产投资过热，控制房价，杜绝炒房。 2) 加快实施房地产税收制度改革 房地产的税收制度改革对于房地产市场的调节非常重要。首先，要以降低政府干预为原则，以市场化为导向，在此基础上建立起完善的税收制度。推进房产税的实施，从而使房地产市场达到一种供需平衡的状态，通过税收的调节控制房价，使房地产市场能够健康发展。 3) 上市公司进驻三四线城市，开发有潜力的地产 目前，京沪深三地房地产上市公司开发房地产以一线城市以及发达的二线城市为主，这些发达地区的房产已经被过度开发，发展空间小，要注重三四线城市的房地产建设，这些地区多数处于高速发展阶段，发展基数小，要将这些地区的房地产潜力开发出来，不但能在很大程度上提高公司的盈利能力，促进公司发展，而且最大化房产公司的盈利潜力，从而很大程度上促进整个行业的发展。
system: [SYSTEM]

user: 摘 要 本文介绍了黎曼猜想等数学伟大猜想对数学方法论发展的巨大推动作用。同时，简述了黎曼猜想与H值138的巨大关联作用，供广大读者和研究者参考。 关键词 :数学猜想，黎曼猜想，费尔马猜想，H值138，胡氏约等式，思想桥梁 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 给任何科学下一个十分完备的定义，都存在一定的困难。所谓数学，简洁地说，就是关于数量关系及其规律的科学，是一切科学基础的基础。现代物理学的前沿已与数学前沿几乎重叠在一起了！数学领域中许多奇妙猜想及其探索，有力推动了数学的发展。 数学的许多猜想，包括哥德巴赫猜想、费尔马猜想、庞加莱猜想和黎曼猜想等，十分引人入胜！证明的过程，发展了许多新方法，构成了数学王国中的瑰宝！ 1900年，伟大数学家希尔伯特(Hilbert)在巴黎举办的第二届国际数学家大会上提出了23个影响深远数学的问题，它为整个20世纪的数学发展指明了方向。大名鼎鼎的黎曼猜想作为第8个问题的一部分而被世人所知，走过了沧桑百年的历程。百年轮回，时至今日，23个问题中已经有19个确定解决，还有3个部分解决。值千禧年之际，美国克雷研究所提出了7个世纪性的数学难题包括黎曼猜想在内，并慷慨地为每个问题设置了100万美元的奖金，多年来，这些难题吸引了许多数学家，但还没有人获奖。 黎曼猜想依然如巍峨的奇山，矗立在人类的智力“喜马拉雅山”巅峰之上。 关于哥德巴赫猜想，国内介绍文献较多，在此不再赘述了。 庞加莱猜想：所谓庞加莱猜想，是指法国数学家庞加莱在1904年提出的一个猜想，即如果一个封闭空间中所有的封闭曲线都可以收缩成一点，那么这个空间一定是三维圆球。当我们在拜读庞加莱猜想已得到证明的消息时，我们想推广庞加莱猜想到高维的情况，结果又看到了推广的高维庞加莱猜想已于20世纪60和80年代得到了证明。
system: [SYSTEM]

user: 2. 由费尔马大定理引出的猜想 20世纪90年代之前的几个世纪，数学领域存在9大世界难题，包括：庞加莱猜想，黎曼假设，霍奇猜想，杨–米尔理论，P与NP问题，波奇·斯温纳顿–戴雅猜想，纳威厄–斯托克斯方程，歌德巴赫猜想和费尔马大定理。20世纪60年代，我国著名数学家陈景润将歌德巴赫猜想推进了一大步，证明了1 + 2。20世纪90年代，英国天才数学家安德鲁·维尔斯证明了费尔马大定理。2006年，我国数学家朱熹平和曹怀东教授在哈佛大学丘成桐教授的指导下，证明了庞加莱猜想。我们相信其他难题在本世纪内将被一个个攻克。 1637年，爱好数学的大法官费尔马在一本书的页边写下了他对一个问题的看法：他发现了一个简洁的证明，但是由于纸张太小无法写下来。这就是被后世称为费尔马猜想的问题，其完整的证明直到358年后的1995年才由英国数学家怀尔斯借助最艰深的现代工具所完成。 费尔马大定理起源于两千多年前，挑战人类三个多世纪，多次震惊全世界，耗尽人类众多最杰出大脑的精力，也让千千万万业余者痴迷，终于在1995年被英国天才数学家安德鲁·维尔斯(Andrew J. Wiles)攻克，成为惊世传奇。1975年安德鲁·维尔斯在剑桥大学攻读研究生，导师为他幸运选择了椭圆曲线的研究方向，掌握了他实现梦想的工具。通过数位天才数学家的努力和他自己千百次的峰回路转，终于利用“科利瓦金–弗莱切”方法证明了日本数学家谷山丰与志村五郎的猜想：每一个椭圆方程的E–序列，都对应一个模形式的M–序列。依据德国数学家格德·费赖关于谷山丰与志村五郎猜想和费尔马大定理之间的联系及肯·里贝特的证明，安德鲁·维尔斯在证明谷山丰与志村五郎猜想的基础上，用反证法显然证明了费尔马大定理。 古希腊的丢番图写过一本著名的“算术”一书，经中世纪的愚昧黑暗时代到文艺复兴的时候，“算术”的残本重新被发现研究，其中有关于勾股定理的描述： X 1 2 + X 2 2 = Y 2 ，西方称为毕达哥拉斯定理，我国古代科学家比西方更早就发现了这一定理。1637年法国业余大数学家费尔马(Pierre de Fermat)在阅读“算术”时，在关于勾股数问题的页边上，写下了著名的猜想： 是不能成立的(这里n > 2，X 1 ，X 2 ，Y，n都是非零整数)，此猜想后来称为费尔马大定理。费尔马在书中写到：当时想到了此猜想的绝妙证明，但因为书页空白太小，无法写下。从此，费尔马大定理成为困扰人类几个世纪的数学难题。据有的数学家猜测，费尔马当时并未想到理想的证明。 费尔马大定理的最终证明是数学发展历史上的一座丰碑。如果能证明下列费尔马大定理的推广的猜想，数论的发展可能会达到一个更新的境界。 X 1 n + X 2 n = Y n (1) 当n = 2时，这就是著名的勾股定理，或西方人称为毕达哥拉斯定理。 但当 n > 2 ， X 1 ， X 2 ， Y ， n 都是正整数，这个方程没有解，这就是著名的费尔马大定理。 X 1 n + X 2 n + X 3 n = Y n (2) 这里 n > 2 ， X 1 ， X 2 ， X 3 ， Y ， n 都是正整数，这个方程没有解，但有解的可能性不能完全排除。 X 1 n + X 2 n + X 3 n + ⋯ + X m n = Y n (3) 这里 m < n ， n > 2 ， X 1 ， X 2 ⋯ X m ， Y ， n ， m 都是正整数。 这个方程没有解，但存在有解的可能性。 上述猜想可能会较费尔马大定理容易得多，这与庞加莱猜想的情况有些类似。不论上述关于费尔马大定理推广猜想是否正确，其有关证明会发展许多数学方法，可望对数论相关领域的发展起重要的推动作用。
system: [SYSTEM]

user: 4. 黎曼猜想与H值138之关系 据英国《新科学家》周刊网站2018年9月24日报道，尽管有人声称可以证明，但数学中最著名的未解难题之一“黎曼猜想”可能仍未被解决。在9月24日的海德堡国际数学与计算机科学获奖者论坛上，退休数学家迈克尔·阿提亚发表了其所谓的证明“黎曼猜想”的惊人思路，他的同行们在近160年的时间里都未能证明“黎曼猜想”。 阿提亚说：“证明‘黎曼猜想’会让你成名。如果你已经成名，那就会臭名远扬。没人相信有人能证明‘黎曼猜想’，因为它太难了。还没有人能证明它，所以现在为什么你能呢？当然，除非你有了一个全新思路。” 阿提亚所谓的“简单证明”是建立在两位20世纪著名数学家约翰·冯诺伊曼和弗里德里希·希策布鲁赫的研究之上的。通过结合他们的研究成果，并假设“黎曼猜想”不成立，阿提亚声称得出了一个逻辑矛盾，这表明该假说事实上肯定是正确的。阿提亚说：“这看起来很神奇。但我要说所有的辛苦工作都是70年前完成的。” 阿提亚在发言中回顾了冯诺伊曼、希策布鲁赫以及数学史上其他著名人物的成果。他仅用几页幻灯片讲解了证明“黎曼猜想”的过程，并称这可与精细结构常数(≈1/H)建立联系。 阿提亚出生于1929年，是英国最著名的数学家之一，他曾获得通常被称为“数学界诺贝尔奖”的两个奖项——菲尔茨奖和阿贝尔奖。他还曾担任伦敦数学学会和英国皇家学会的会长。阿提亚近年来发表了多篇论文，得出了一些令人瞩目的结论，但这些结论尚未令同行信服。尽管他的最新结论尚未经过严格的同行评审以检验其正确性，但外界最初的反应是持谨慎怀疑态度。 英国华威大学的尼古拉斯·杰克逊说：“黎曼猜想是一道出了名的难题。多年来，很多一流数学家接近却未能成功证明这一猜想。”“希望能从阿提亚的研究中得出一些有用的见解，即使证明黎曼猜想的思路站不住脚。” 然而，阿提亚希望，如果得到证实，他的证明思路将激励年轻一代将研究成果的应用范围延伸至“黎曼猜想”的更普遍情况，以及似乎与之不相关的数学领域，甚至包括物理领域。”
system: [SYSTEM]

user: 5. H值138是通向微观、介观、宏观和宇观世界的桥梁 胡氏约等式(胡文祥约等式)可以表述为： 现代交往朋友人数 ≈ 猿人洞里的人数 ≈ 原始部落的人数 ≈ 母系氏族人数 ≈ 早期村落的人数 ≈ 动物单群里的个数 ≈ 现代学术交流会议最佳人数 ≈ 现代军事单位连队的人数 ≈ 现代村组(过去的生产队)的平均人数 ≈ 社会基层组织平均人数 ≈ 人类理想平均寿命 ≈ 精细结构常数的倒数 ≈ 宇宙中原子序数的上限 ≈ 哈勃时间 ≈ 宇宙年龄(亿年) ≈ 黎曼猜想特征值 ≈ 138。 上述胡氏约等式中138称为胡氏数(胡文祥数，或称为H值)，胡氏数是笔者推导出的宇宙中原子序数的上限，接近精细结构常数1/137的倒数，精细结构常数是物理学中一个非常重要的无量纲数，表示电子在第一玻尔轨道上的运动速度和真空中光速的比值，是微观世界的一个常数，却在数学世界、无机世界、有机世界、生物和人类社会中扮演了重要角色 [ 3 ] 。英国杰出数学家阿提亚在上文中指出：黎曼猜想可与精细结构常数(≈1/137)建立联系。 上述胡氏约等式中16个约等号(还可以更多)蕴含了丰富的内容，充分表明：H值138是宇宙中的一个特征数值，是通向微观、介观、宏观和宇观世界的伟大桥梁!数学世界、原子世界、无机世界、有机世界和生物社会界及人类社会界乃至整个宇宙等都遵守共同的宇宙基本规律 [ 4 ] [ 5 ] [ 6 ] 。 这些研究成果再次印证了伟大的物理学家伽利略曾说过的一句至理名言：数学是上帝用来书写宇宙的文字。
system: [SYSTEM]

user: 摘 要 在由单一供应商和单一零售商组成的供应链中讨论零售商单位生产成本以及不完全信息对供应链成员决策的影响，包括供应商批发价格决策以及零售商的生产和采购策略决策。构建了零售商在完全信息和不完全信息两种情形下选择生产和采购两种产品策略时的数学模型，分别对比两种产品策略下的供应链成员收益得出结论。研究发现，当供应商提供的批发价格高于零售商的生产成本时，零售商也会选择以更高的价格向供应商进行采购。同时不完全信息会削弱“低类型零售商”采购的意愿(不完全信息下低类型零售商仅选择生产策略)，提高“高类型零售商”选择采购策略的意愿。 关键词 :生产和采购，不完全信息，供应链管理 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 生产和采购策略在现在企业管理与供应链管理中是重要的研究主题之一，几乎所有企业都会对该策略定期地进行重新评估。传统观点认为，企业采购决策可能仅仅是将内部生产成本与外部供应商收取的批发价格进行比较，进而选择更便宜的策略。然而，已有的研究成果表明企业的生产或采购策略选择取决于多种因素，会更加复杂。 国外早期的经典文献对影响供应链管理中生产和采购策略的因素进行了分析研究。Demski (1994) [ 1 ] 研究了供需双方的长期动态交互对供应链决策的影响，Anderson和Parker (2002) [ 2 ] 考虑了边做边学(learning-by-doing)会影响供应链生产和采购决策的可能性。Demski和Sappington (1993) [ 3 ] 还强调了可能削弱企业采购策略的实际因素，包括需要充分地保证投入质量以及避免在采购策略中披露私有信息的考虑。此外，Van Long (2005) [ 4 ] 的研究结果已经指出，企业的采购策略会造成技术溢出，这种技术溢出对竞争对手更有利。与此同时，Shy和Stenbacka (2003) [ 5 ] 还得出向供应商采购可以避免企业过多固定成本的结论。Arya等人(2008) [ 6 ] 的研究强调产出市场的竞争与投入市场定价之间的相互作用也可能使采购策略复杂化。特别是，如果企业的主要竞争对手依赖于供应商投入，那么企业可能更加依赖向供应商采购，即使这样做比企业自己生产更昂贵。 由此可见，早期研究表明企业是选择生产还是采购取决于多种因素。国内也有许多学者对这一问题进行了研究，陈月梅(2006) [ 7 ] 认为企业采购是降低成本，提高核心竞争力的重要措施，并且分析得出企业战略、资产专用性程度和供应商选择等是影响企业采购策略的主要因素。尹小勇和全裕吉(2009) [ 8 ] 分析了跨国公司业务外包的动因及其必要的条件，认为企业采购可以降低成本和风险，增加企业灵活性和响应市场的速度。张汉江，张媛和李俊萍等(2011) [ 9 ] 人的研究发现供应链下游企业选择采购时的上游企业利润总比下游企业选择生产时更高，对于下游企业来说，其自制获得的利润总是高于采购获得的利润。 现有文献大多是在信息完全状态下进行讨论，在供应链背景下，市场的复杂多变使得企业之间往往具备信息非对称，因此本文将不完全信息纳入到研究模型中。本文中供应商自有的销售渠道可以看做零售商的依赖于供应商投入的竞争对手，研究结果证明即使供应商提供的批发价格高于自己生产的成本，零售商也会选择采购策略。零售商更愿意选择采购策略是因为它可以降低供应商自有销售渠道的既得利益，从而减少供应商对自有销售渠道的产品投入，提高零售商自身的竞争力。早期文献也有对于采购战略优势的研究，Schrader和Martin (1998) [ 10 ] 在上游和下游市场的Cournot竞争模型中，表明了采购策略的价值，用以减少下游竞争对手可获得的投入供应，从而阻碍竞争对手的零售业务。 现有文献对零售商生产和采购策略的研究大部分是在完全信息条件下进行的，本文以完全信息为基准，讨论了不完全信息下零售商生产和采购决策的变化。本文下面的章节安排如下，第2章对供应链模型和模型符号进行说明，第3章分析了作为讨论基准的完全信息下零售商的生产和采购策略选择；第4章分析了不完全信息下零售商的生产和采购策略选择，同时对比研究了不完全信息对零售商决策带来的影响；第5章总结了本文研究结论以及未来研究方向。
system: [SYSTEM]

user: 2. 模型说明 首先对论文使用的数学模型模型及模型的变量和参数进行说明。假设市场中存在两家公司：供应商S (Supplier)和零售商R (Retailer)。其中供应商有两种产品销售渠道，一是以批发价格w出售给零售商R，二是自己销售，销售数量为Q，将供应商的单位生产成本和运营成本归一化为0，零售商R以单位成本c自己生产产品或者从供应商处以批发价格w采购产品。供应商和零售商面临的市场需求满足： p = a − Q − q 。 其中p表示产品的销售价格，Q与q分别表示供应商和零售商的销售数量。a表示市场需求。本文 假设产品的市场需求 a i 有两种类型：高需求和低需求类型，即 i ∈ { H ， L } ， a H 代表产品市场需求为高， a L 代表产品市场需求为低，且假设 a H > a L 。产品市场需求类型的概率分布满足 ， P r ( a L ) = 1 − β ，且 0 < β < 1 ，则市场需求的平均值为 μ = β a H + ( 1 − β ) a L 。 在不完全信息下，由于零售商R更接近消费者市场，因此可以确切地了解产品实际市场需求 ，而 对于供应商则是未知的，只能了解 的概率分布和市场需求均值 μ 。本文用 表示供应商对市场需求的信念，其中 j ∈ { H ， L } 。即当 时，表示供应商对产品市场需求的判断与实际情况相符，反之亦然。因此产品市场需求 a i 是零售商R的私有信息(为了表述方便，下文将市场需求为低类型时的零售商R简称“低类型R”，市场需求为高类型时的零售商R简称“高类型R”)。
system: [SYSTEM]

user: 5. 结论 本文在完全信息和不完全信息下，分别研究零售商的生产和采购决策以及供应商的批发价格决策，本文所取得的研究成果主要有以下三点。首先，由于供应商想要保证其直销渠道在与零售商的竞争中处于有利地位，所以当零售商选择生产策略时会提高供应商直销渠道的产品销售量，而零售商的采购策略会降低竞争对手供应商的销售量。其次，与传统观点不符的是，零售商为了在与供应商直销渠道的序贯博弈中获得leader优势，更愿意支付超过其内部生产成本的批发价格进行采购。最后，在不完全信息下，低类型零售商只选择自己生产产品，也就是说不完全信息削弱了“低类型零售商”进行采购的意愿，同时提高“高类型零售商”采购产品的意愿。 本研究为供应链中的生产或采购决策补充了理论依据和实践启示，并且将不完全信息等因素纳入考虑，使得研究内容更加丰富。同时由于对市场需求和供应链模型的刻画进行了一定程度的简化，对其进行更详实和现实的描述从而完善研究结论是今后的研究方向。
system: [SYSTEM]

user: 摘 要 供应链当中由于各个节点企业间存在信息传递缓慢、信息不对称导致的需求预测不准、订货时间不稳定、订货量不稳定等多种因素导致库存和安全库存高居不下，库存周转率相对较低、库存成本相对较高等常见问题。本文从系统动力学的角度，分析导致库存过高的原因。文中运用vensim建立供应商、配送中心、门店这一常见的供应链仿真模型，最后通过某一品类常见的快消品进行实例分析。文章最后通过对库存调整周期、配送延迟时间、配送中心库存、各个节点的安全库存进行敏感性分析，结果表明：通过改变库存策略，可以在满足超市的需求的情况下，使超市和配送中心保持较低的库存水平，大大降低了库存成本。该系统能够有效地反映库存及其成本的变化，可以为供应链库存管理提供参考。 关键词 :SD，Vensim，供应链库存控制 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着现在科技的发展，网购因为其“物美价廉”、“24小时营业”等特性越来越受大众的普遍欢迎，然而，由于商品的种类相对较多，同时人们对不同的商品的消费习惯的不同，社区便利店/小型的连锁超市并没有受到太大的影响，反而在这场角逐之战中，体现出其独特的生命力。在整个连锁超市的供应链当中，库存占据着重要的地位。所以库存管理是整个连锁零售企业运营管理中的不可或缺的一个重要的环节。我们都知道，库存存在的主要目的是为了防止短缺，快速满足用户的需求，库存实际上是人们在无法准确预测未来需求变化的一种无奈之举。 然而，由于过多的库存导致资金的占用、仓库管理成本的增加、由于预测不准导致产品滞销，因为产品滞销带来的资金、和库房占用等等造成费用在整个零售企业运行费用占有主导地位 [ 1 ] [ 2 ] [ 3 ] ，所以，如何在能够满足用户需求的前提之下，尽量的降低库存至关重要。在连锁超市这条供应链中库存存在于供应商、配送中心和门店等多方。连锁零售企业多级库存具有明显的非线性、多环路、动态性、延迟性等复杂动力学特征，系统动力学是Forrest教授于1958年提出的一种分析企业生产管理、库存管理等问题的系统仿真方法，系统动力学擅长处理高阶复杂、多环路、非线性及时间延迟的动态问题 [ 4 ] ，因此系统动力学具有应用于库存控制领域的基础条件 [ 5 ] [ 6 ] 。应用系统动力学建模的步骤： ① 问题的识别，明确研究目的； ② 确定系统的边界，即系统分析的对象和范围； ③ 构建模型，绘制系统流程图，并建立方程； ④ 模型模拟； ⑤ 结果分析，对模拟结果进行分析，预测，设计，测试选择最优方案。
system: [SYSTEM]

user: 2. 文献及理论分析 在供应链库存控制领域的系统动力学应用中，不同学者从不同角度进行了大量的研究。福瑞斯特(Forrest) [ 7 ] 根据系统动力学的基本原理，分析出了消费需求波动沿供应链向上游企业逐级放大的特性，并通过灵敏性分析了因决策延迟、需求波动、随机噪声、生产能力限制、广告等因素对库存的影响。多维尔(Towill) [ 8 ] 在系统动力学基础上建立了著名的库存——订货控制系统(IOBPCS)，解决了由需求波动引起的库存量增加等一系列问题。波勒斯(Poles) [ 9 ] 从制造供应链角度建立了生产–库存的系统动力学模型，分析库存量、交货时间等因素对系统的影响，通过灵敏度分析了生产提前期、再制造提前期整个系统的影响系数。巴拉吉和詹姆斯(Balaji & James) [ 10 ] 研究建立供应商和制造商协同的供应链库存管理系统动力学模型，简称CPFR，通过比较不同随机需求下该模型的效果，验证了协同式供应链库存管理可以很好的降低库存成本。 在国内学者中，张力菠、韩玉启等 [ 11 ] 通过总结系统动力学的理论研究与应用发展的基础上，对供应链管理相关问题进行分类汇总，指出了在供应链管理领域中系统动力学具有独特的应用方向。廖诺、张毕西、吴小结 [ 12 ] 建立了一个供应链节点企业的系统动力学模型，通过比较不同外部需求条件下模型的动态性能，提出了供应链优化的方案。汪小京、刘志学、徐娟 [ 13 ] 将第三方物流(TPL)引入VMI中，建立了第三方物流管理库存(TMI)的系统动力学模型，通过仿真分析发现，在需求相同的情况下，TMI比VMI能够进一步的减低库存水平。成琼文和周璐 [ 14 ] 则从绿色供应链角度，建立了绿色供应链管理系统，其中主要包括采购、库存、设计和管理等四个方面，最后通过实践分析发现，各个因素之间相互影响，互为联系。以上文献主要是从传统供应链的角度研究采购、生产、配送、销售各个环节之间的联系，但在连锁零售供应链多级库存这个方向上的研究还比较缺乏。本文从系统动力学的角度，考虑到连锁零售企业多级库存网络化、连锁性、订货批量小、交易频繁和需求变化大等特点，优先考虑了连锁零售企业多级库存管理系统中供应商、配送中心和门店库存之间的内在信息流因果关系，然后，根据因果关系构造了构造一个动态库存管理模型，最后，通过敏感性实验对模型进行了优化，力争为连锁零售企业多级库存管理提供有利的决策依据 [ 15 ] 。
system: [SYSTEM]

user: 3. 确定模型边界和相关假设 3.1. 建模目的 对于一个涉及供应商、配送中心、超市的供应链系统，如何实现供应链各个节点之间的协调发展，对整个供应链的效率有着重要的影响。通过建立系统动力学仿真模型，研究了在消费者随机需求下各节点企业的供应商、配送中心和超市的库存水平，试图找出影响库存动态变化的因素，降低整条供应链的库存水平，进而缓解库存费用。 3.2. 模型边界 文中研究的是超市供应链系统，商品从供应商到配送中心，之后有配送中心配送到门店，并最终到消费者手中。在整条供应链的上游，考虑配送中心向向供应商递交订单进行采购，只考虑供应商的延迟交货，下游考虑的便是便利店向用户提供产品。 3.3. 模型假设 为了模拟研究问题，在建立模型之前必须对模型做如下假设： · 在模型中，只考虑某一品类的产品。 · 研究对象是配送中心和超市商店组成的库存管理系统。 · 配送中心的上游供应商可以随时满足其产品采购要求，不考虑供应商的生产和供货能力，假设供应商的生产和运输能力不受限制。 · 供应商与超市配送中心、超市配送中心和超市商店之间的交货延迟的时间是固定的。 · 假设在整个销售的过程中不存在退货，即没有逆向物流。 · 根据现实情况，在连锁零售的商品，一般都会有一定的保质期，为了保证商品能够在保质期内售卖完，文中模拟超市通过降价促销的手段来增加销售额。
system: [SYSTEM]

user: 4. 库存管理系统模型流程图构建 多级库存管理流程图 在分析因果图的基础上，构建了超市库存系统的系统动态流程图，如图1所示。 图1. 配送中心和超市库存流程图
system: [SYSTEM]

user: 5. 案例分析 5.1. 模型模拟和结果分析 本文以啤酒这一类快消品为例。在Vensim进行仿真并得到以下结果：在现实生活中，顾客的需求是波动的，为了尽可能的模拟现实的波动情况，所以我们将消费者的波动需求设置为阶跃函数 [ 16 ] ：ECD = 5 + RAMP(0.3，0，60)。具体的变量赋值情况请如表1所示。 Table 1 Model variable assignment tabl 变量名称 赋值 SI 60 SSCT 2 ECD INPUT CDST 1 SPI 0 DT1 1 SIAC 3 SOC 1 SOB 0 DI 300 DPI 0 DOB 0 DT2 2 DIAC 5 DOC 2 DSCT 5 表1. 模型变量赋值表 图2. SI随延迟时间的变化 图3. SPI随延迟时间的变化 图4. DI的延迟时间 图5. DPI随延迟时间的变化 延迟时间1是影响商店到达率的时间。延迟时间2是影响配送中心到达率的时间。从图2~5中，我们可以看到随着延迟时间2的增加，配送中心的库存和配送中心的库存正在减少。随着延迟时间1的缩短，商店中商店的库存减少，商店库存增加。门店库存处于短缺状态，当客户需求稳定时，实现稳定状态。之后，它直接影响商店的正常销售和客户的满意度。 如图6和图7所示，随着商店的库存调整周期缩短，商店和配送中心供应商的每个库存调整的振荡幅度和频率都有所降低，但是库存水平上升。因此，在确定商店的库存调整周期时，应注意确保库存保持在低水平并且幅度波动小。当配送中心库存调整周期延长时，配送中心的库存波动增加，商店库存不受影响。 图6. IAC随SI的变化情况 图7. IAC随DI的变化情况 5.2. 模型改进和效果分析 1) 流程优化 从库存管理流程和模拟模型结果的模型分析，门店通过往期消费者的消费情况，来预测消费者的需求量，再根据门店现有库存和门店的在途库存，来制定下一期的订货方案；当订货需求传递到配送中心之后，配送中心根据超市的订单量、配送中心在途库存、安全库存等因素综合制下一期的配送中心的订货量；供应商根据配送中心的订单量和所了解的市场波动情况以及安全库存等因素制定生产表，进行批量生产。在供应链的整个链条之中，存在信息不对称等问题，通过上述系统动力学的模型很好的解释了由于信息孤岛导致的“牛鞭效应”。 牛鞭效应产生的恶果是延长了供货期，导致供应商很容易出现当前的生产力、库存和订单需求匹配度不高的问题；增加了供应链上各级的库存成本等。 可以通过一下方式来缓解“牛鞭效应”的危害： a) 调整供应链上的信息结构。 建立有效的信息共享机制，其中最典型的就是VMI (供应链库存管理)。以各方的最低成本为目标，在一个共同的协议之下，通过供应商来管理客户的库存，通过不断的修改协议的内容，从而到达库存成本最小，进而形成战略上合作伙伴关系。 b) 建立现代化管理信息系统 充分利用供应链各个节点之间的企业订单自动处理系统、电子商务系统、扫码技术等进行集成，通过互联网，能够很好的进行采购、运输、仓库等各个环节更好的衔接起来，进而提升效率，降低成本。 2) 动力学方程的优化 SRR = MIN ( DRR ， CAD ) (1) DDI = CAD ∗ DSCT (2) DOC = IF _ THEN _ ELSE ( DIG = O ， CAD ， DIG / DIAC + CAD ) (3) SRR = MIN ( CAD ， DI ) (4) 3) 参数优化 优化现有的可控决策变量，并重置决策变量的参数，如表2所示。 Table 2 Parameter optimizatio 变量名称 赋值 DIAC 2 DT1 2 DT2 1 表2. 参数优化 结合表2的参数优化，通过仿真软件对优化策略进行仿真和仿真，仿真结果如下，从图8和图9中可以看出： 文中的延迟时间1和延迟时间2，主要是值物品运输时间和入库时间，随着交通的发展，再加之连锁超市的实际情况，不考虑天气的影响，目前路上的运输时间基本上是固定的，没有办法进行缩短，所以在实际的管理过程中调整库存调整周期来合理的观测库存的变化，最终确定最优的各级库存调整订货周期的时间。
system: [SYSTEM]

user: 6. 结论和不足 6.1. 结论 本文通过系统动力学的方法，从复杂性系统的角度进行分析，首先定义连锁超市库存管理系统边界，建立了连锁超市库存的仿真模型；然后，分析系统中各因素之间的关系，经过测试，得出仿真模型与实际系统更加一致；最后，对超市库存管理系统的仿真模型进行了优化和改进。 本文最后发现可以从以下两个方面进行连锁零售企业库存的管理：一是提高供应商的供应效率，减少配送中心的到货时间，降低配送中心的库存水平和波动，降低配送中心的库存管理效益；二是减少中间环节，提高信息共享度。终端客户需要直接通过信息系统到配送中心，配送中心分析需求信息以产生对上游供应商的购买请求。一方面，它可以有效提高库存响应速度，提高客户服务水平和满意度。另一方面，配送中心直接维护所有店铺的库存，店铺库存仅满足市场日的需求，可以节省店铺管理成本，发挥配送中心的综合管理作用。 图8. 优化后超市库存的模拟结果 图9. 优化配送中心库存的模拟结果 6.2. 不足 由于连锁超市库存管理系统因素众多，关系复杂，需要进一步完善。例如，系统外部供应商的供应能力和人力资源等因素。
system: [SYSTEM]

user: 摘 要 基于《综合幸福感问卷》 (MHQ)量表和SF-12生存质量量表进行实证研究，分析比较流动人口与城市居民幸福感的差异，探索流动人口幸福感与生存质量的相关性，为提高流动人口幸福感提供参考依据。结果显示：城市居民主观幸福感(SWB)、友好关系得分高于流动人口，负性情感得分流动人口高于城市居民，差异有统计学意义(p < 0.05)。流动人口幸福感与躯体健康和心理健康方面以及多个维度得分存在相关性(p < 0.05)。幸福感与生存质量存在正相关，可通过提高流动人口生存质量提升其幸福感。 关键词 :流动人口，幸福感，生存质量 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着我国经济和工业的迅速发展，我国社会的主要矛盾转化为人民日益增长的美好生活需要和不平衡不充分的发展之间的矛盾，这充分表现出当代人民对生活质量的追求，以及我们对人民生活的幸福感提升的欠缺。在城市快速化的发展中我国大规模的流动人口做出了巨大贡献，他们是城市化建设中的重要力量，据《中国流动人口发展报告2016》显示，未来一段时间我国人口流动迁移仍将持续活跃，且新生代流动人口占比持续提升，因此对流动人口的生存质量和幸福感问题的关注和解决必须提上日程。幸福感的研究主要分为两种取向：主观幸福感(SWB)和心理幸福感(PWB)，主观幸福感是衡量个人生活质量的重要综合性心理指标，具有多层次结构模型，与心理幸福感有着相互依赖的关系，将两者整合能更好更全面地反应流动人口的幸福感(严标宾，郑雪，& 邱林，2004)。生存质量是个体在所处社会文化背景下，依照自己的目标、期望和关注，对自身生活状态的感受与评价，包括心理、生理、社会功能和物质状态四个方面，是对流动人口在异地生活情况的综合反映，对流动人口的幸福感、生活满意度及社会支持的建立和发展都具有极其重要和深远的意义。 本研究比较流动人口和城市居民综合幸福感和生存质量的差异，探索流动人口幸福感与生存质量的相互关系，为缩小流动人口与城市居民差距，提高流动人口整体生存质量和幸福感等提供参考依据。
system: [SYSTEM]

user: 2. 对象与方法 2.1. 研究对象 采用分层抽样的方法，将浙江省按照经济水平分为高、中、低三类，每类各随机抽取1个市，对流动人口集中的企业进行现场面对面问卷调查。最终获得样本量流动人口768人，平均年龄为33.95 ± 9.95岁，其中男性占62.4%、女性占37.6%，同时选择相应调查单位工作3个月及以上的本地居民为城市对照组，共计350名，平均年龄为36.69 ± 10.18岁，其中男性占59.9%、女性占40.1%。 2.2. 研究方法 问卷主要包括苗元江的《综合幸福问卷》 (MHQ)、SF-12生存质量量表和一般情况问卷三部分。 苗元江的《综合幸福问卷》(MHQ)：该量表包含50个条目，9个维度(生活满意、正性情感、负性情感、生命活力、健康关注、利他行为、自我价值、友好关系、人格成长)，其中前3个维度构成主观幸福感，其余6个维度构成心理幸福感(苗元江，2003)。 SF-12生存质量量表：由躯体健康(PCS)和心理健康(MCS)两部分组成，共有8个维度：生理功能(PF)、生理职能(RP)、躯体疼痛(BP)、总体健康(GH)、活力(VT)、社会功能(SF)、情感职能(RE)、心里健康(MH) (张莎，田晋，刘巧兰，等，2011)。 一般情况问卷：内容包括出生年月、性别、原居住地、学历、婚姻状况、所在行业、每天工作时间、平均月收入、自觉社会地位等。 2.3. 统计分析 采用SPSS18.0进行统计分析，统计方法包括统计描述、t检验、相关分析等。
system: [SYSTEM]

user: 3. 结果与分析 3.1. 流动人口与城市对照幸福感的比较 流动人口主观幸福感(SWB)、友好关系的得分低于城市居民，负性情感得分则高于城市居民，差异有统计学意义(p < 0.05)；流动人口总幸福感、心理幸福感、生活满意、正性情感、健康关注、生命活力、自我价值、人格成长、利他关系与城市居民无统计学意义差异(p > 0.05)。见表1。 Table 1 Comparison of the happiness of floating population and urban contro 流动人口 城市对照 t值 p值 总幸福感 4.81 ± 0.74 4.88 ± 0.78 −1.402 0.161 主观幸福感(SWB) 4.82 ± 0.83 4.94 ± 0.89 −2.132 0.033 心理幸福感(PWB) 4.80 ± 0.83 4.84 ± 0.87 −0.816 0.415 生活满意 4.34 ± 1.21 4.45 ± 1.23 −1.369 0.171 正性情感 4.64 ± 1.21 4.60 ± 1.28 0.416 0.677 负性情感 2.56 ± 1.28 2.23 ± 1.13 4.329 <0.001 健康关注 5.21 ± 1.12 5.19 ± 1.16 0.242 0.809 生命活力 4.60 ± 1.16 4.67 ± 1.09 −1.000 0.318 自我价值 5.04 ± 1.04 5.04 ± 1.07 −0.017 0.986 人格成长 4.56 ± 0.81 4.56 ± 0.86 −0.063 0.949 友好关系 4.92 ± 1.20 5.09 ± 1.13 −2.224 0.026 利他行为 4.74 ± 1.12 4.85 ± 1.14 −1.568 0.117 表1. 流动人口与城市对照幸福感的比较 3.2. 流动人口幸福感与生存质量的相关性 流动人口幸福感与躯体健康和心理健康方面以及多个维度得分存在相关性(p < 0.05)。负性情感与生活质量及各个维度呈负相关，总幸福感、主观幸福感、心里幸福感、正性情感、健康关注、生命活力、自我价值、人格成长、友好关系与生活质量及各维度呈正相关(详见表2)。
system: [SYSTEM]

user: 4. 讨论 随着我国城市化进程的进一步加快，流动人口数量不断增加。本研究结果显示，流动人口主观幸福感和友好关系低于城市居民，负性情感则高于城市居民(p < 0.05)。由于户籍制度的限制，教育程度和劳动能力的差异，流动人口的生存质量往往低于城市居民。职业的可选择性、社会医疗保障等与城市居民的差异也使流动人口缺乏安全感，难以融入城市(袁爱清，2014)，影响了流动人口的主观幸福感。在自我群体归属上，流动人口存在难以准确将自己定义为“农”或“工”，群体归属困难，缺乏归属感的状况，又难以融入城市生活，加之房租、抚养小孩、对家乡的牵挂等一系列的生活压力，使其无法安逸快乐地享受生活，影响了友好关系，负性情感得分较高。据《中国流动人口发展报告2016》指出基本公共卫生服务、计划生育服务和社会医疗保险的流动人口覆盖面不断扩大(国家卫生和计划生育委员会流动人口司2016)，这些制度措施的实施，将有利于加强流动人口的安全感，提升幸福指数，但政策的落实和辐射效应以及流动人口幸福感的提升尚需一段时间。 Table 2 Correlation between the happiness of floating population and quality of lif PCS MCS PF RP BP GH VT RE SF MH 总幸福感 0.293 ** 0.248 ** 0.246 ** 0.204 ** 0.205 ** 0.282 ** 0.191 ** 0.208 ** 0.200 ** 0.245 ** 主观幸福感 0.258 ** 0.303 ** 0.200 ** 0.226 ** 0.247 ** 0.262 ** 0.150 ** 0.267 ** 0.199 ** 0.307 ** 心理幸福感 0.259 ** 0.183 ** 0.226 ** 0.157 ** 0.147 ** 0.245 ** 0.183 ** 0.144 ** 0.162 ** 0.176 ** 生活满意 0.107 * 0.142 ** 0.090 * 0.060 0.071 0.192 ** 0.043 0.094 ** 0.093 * 0.158 ** 正性情感 0.262 ** 0.207 ** 0.185 ** 0.178 ** 0.130 ** 0.295 ** 0.179 ** 0.134 ** 0.067 0.233 ** 负性情感 −0.164 ** −0.243 ** −0.121 −0.222 ** −0.300 ** −0.042 −0.074 * −0.300 ** −0.241 ** −0.198 ** 健康关注 0.310 ** 0.138 ** 0.250 ** 0.185 ** 0.176 ** 0.256 ** 0.162 ** 0.167 ** 0.176 ** 0.116 ** 生命活力 0.262 ** 0.077 * 0.253 ** 0.085 * 0.112 ** 0.230 ** 0.122 ** 0.078 * 0.095 ** 0.082 * 自我价值 0.181 ** 0.189 ** 0.174 ** 0.139 ** 0.125 ** 0.155 ** 0.169 ** 0.125 ** 0.161 ** 0.178 ** 人格成长 0.196 ** 0.264 ** 0.143 ** 0.167 ** 0.152 ** 0.232 ** 0.231 ** 0.160 ** 0.166 ** 0.255 ** 友好关系 0.120 ** 0.097 * 0.078 * 0.097 ** 0.105 ** 0.107 ** 0.079 * 0.068 0.036 0.131 ** 利他行为 0.150 ** 0.086 * 0.156 ** 0.079 * 0.040 0.144 ** 0.079 * 0.084 * 0.108 ** 0.066 表2. 流动人口幸福感与生存质量的相关性 *P < 0.05，**P < 0.01。 研究表明，流动人口幸福感与躯体健康和心理健康方面以及多个维度得分均存在相关性，MCS和PCS得分越高，幸福感越高。心理和躯体越健康，情绪越积极向上，越有机会接触新事物，越有能力和活力去创造生活，能更好地融入新环境，增加归属感、幸福感。
system: [SYSTEM]

user: 5. 结论 综上所述，为进一步提高流动人口幸福感，改善流动人口生存质量，应加强对流动人口健康的关注，尽所能去满足健康卫生方面的需求，完善社会保障制度、完善社会基础医疗设施；切实保障农民工的基本权益，更好地贯彻落实相关法律；应加强对流动人口的人文关怀，号召公司、社区居委等积极组织集体活动，提升流动人口归属感与积极的情绪，帮助流动人口消除“外来”的陌生感，温暖积极地与他们相处，营造人人平等、和谐互助的社会氛围，促进社会的融合。
system: [SYSTEM]

user: 摘 要 为研究钢筋混凝土柱墩的抗震性能，进行了箱型截面混凝土柱墩水平双向往复拟静力试验，利用ABAQUS有限元分析软件在PQ-Fiber子程序的基础上进行了pushover数值模拟，模拟结果与实验基本吻合。在实验基础上截断部分钢筋进行双塑性铰的模型分析，表明可以在墩底形成双塑性铰，通过对比单双塑性铰柱墩的骨架曲线、滞回曲线、应变曲线等，表明双塑性铰柱墩承载力没有明显变化，而极限变形能力明显提高，桥墩的抗震能力得到增强。 关键词 :钢筋混凝土柱墩，双塑性铰，有限元分析，滞回曲线，骨架曲线，抗震性能 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 钢筋混凝土柱墩作为桥梁结构中的主要承载地震力的构件，震损破坏通常出现在柱墩的底部或顶部，国内外学者研究发现通过加固或使用特殊材料可以使塑性铰发生转移。邓江东 [ 1 ] 研究了一种新型复合立柱，通过在钢筋混凝土柱外表面粘贴一定的高度和厚度的钢板并将其焊接在基座的嵌入式工字钢上，发现柱底和加强钢板的上边缘可以形成双塑性铰，双塑性铰桥墩的最大承载力、极限变形、耗能能力都显著提高。美国加州大学Panagiotou [ 2 ] 对高层建筑中剪力墙结构采用双塑性铰DPH设计方法，即在结构的底部和1/2处分别设置塑性铰，表明双塑性铰链可消耗大部分的地震能量，延长结构周期，减小地震力的作用，而且合理的布置还可降低经济成本。梁兴文 [ 3 ] 建立了钢筋混凝土剪力墙模型，对6种不同的塑性铰数量和布置位置进行了非线性时程分析，表明多塑性铰区的塑性变形可减小高阶振型的影响，提高抗震性能。蒋永生 [ 4 ] 采用2种人工塑性铰转移铰方案，将塑性铰设置在离柱边一定距离的梁腹中下部，实验结果表明塑性铰转移后节点抗震性能优良。闻洋 [ 5 ] 通过对三层框架矩形截面钢管混凝土进行低周反复试验，试验发现塑性铰出现的顺序从中间层梁端到底层再到顶层梁端，最后出现在柱底部，体现了“强柱弱梁”的破坏机制，采用ABAQUS模拟得到的滞回曲线相吻合，抗震性能良好。 但总体来说对混凝土柱墩产生双铰的研究很少，在模型试验的基础上，本文采用原有钢筋截断的方法在截面的1/4处及柱底形成双塑性铰，并用ABAQUS进行模拟计算分析双塑性铰的滞回性能、骨架曲线，表明形成双塑性铰对结构的抗震更加有益。
system: [SYSTEM]

user: 2. 实验模型 2.1. 截面设计 本次桥墩实验共包括两个试件，A1试件加载高度为3.2米，竖向轴力为420 KN，轴压比0.2。B1试件加载高度为2.2米，轴压比0.1。为了在柱的顶部进行两个方向的往复推拉，设计了540 × 700 × 720 mm的柱帽，为了锚固墩底，设计了尺寸为250 × 800 × 1300 mm的基座。采用C50细石混凝土，模型柱纵筋配28ϕ8，箍筋ϕ6，间距通长100 mm，具体截面尺寸和配筋如图1，力学参数见表1。 2.2. 加载方式 模型竖直和水平方向采用不同的加载方式，竖直方向实施力控制加载，采用千斤顶预加载并保持恒定的力。水平方向位移加载，采用双向受力同步加载制度，加载曲线如图2 (单位：mm)。 图1. 截面配筋 图2. 水平加载幅值曲线 Table 1 Mechanical parameters of concrete (mpa 强度等级 抗压强度 抗拉强度 弹性模量 C50 37 3.3 3.54 × 10 4 钢筋的力学参数 Mechanical parameters of steel bars 直径 屈服强度 极限强度 弹性模量 6 375.27 503.81 2.11 × 10 5 8 326.51 456.6 2.12 × 10 5 表1. 混凝土的力学参数(mpa)
system: [SYSTEM]

user: 3. 基于ABAQUS的数值模拟 3.1. 混凝土材料本构模型 3.1.1. 受压本构 ABAQUS提供了三种自带的混凝土模型，本文采用能模拟往复加载的塑性损伤模型(Concrete Damage Plasticity)，混凝土的应力应变关系曲线多种多样，选用我国规范(GB 50010-2010)本构 [ 6 ] ，由80年代早期清华大学过镇海等人提出的计算式能较好的用于ABAQUS模拟。其受压曲线表达式如下： 上升段x升段其受压曲线表达式如下： y = a a x + ( 3 − 2 a a ) x 2 + ( a a − 2 ) x 3 下降段x降段受压曲线表达式如下： y = x a d ( x − 1 ) 2 + x ε d = ( 700 + 172 f c k 0.5 ) * 10 − 6 α a = 2.4 − 0.0125 f c k α d = 0.157 f c k 0.785 − 0.905 其中， x = ε / ε d ， y = σ / f c k ， f c k 为混凝土轴心抗压强度，若有实测值则取实测。 α a 为上升段参数， α d 为下降段参数， ε d 为混凝土峰值应力所对应的应变。 E c 为弹性模量， E 0 为0.7 f c 对应的割线模量，根据计算得出表2数据。曲线应变截断值取3倍峰值应力对应的应变 [ 7 ] [ 8 ] 。 Table 2 Compressive constitutive parameters of concrete (unit: Mpa f c ε d α a α d E 0 37 0.00175 1.9375 1.7676 3.22765 × 10 4 表2. 混凝土的受压本构参数(单位：Mpa) 3.1.2. 受压损伤因子d c 采用由sidiroff根据能量等价原理，应力作用在受损材料产生的弹性余能与作用在无损材料的余能在形式上相同，将弹性模量改为损伤时的等效弹性模量，得出的无损伤方程： W 0 e = σ 2 2 E 0 等效有损伤弹性余能： W d e = σ ′ 2 2 E d 化简得到： d = 1 − ( σ / E 0 ε ) 0.5 。 3.1.3. 受拉本构 ABAQUS中提供了3种定义混凝土单轴受拉应力–应变关系的接口，包括应力–开裂应变关系；应力–裂缝宽度关系；直接输入受拉断裂能。由混凝土断裂能准则可以看出混凝土的应力–开裂应变关系与单元尺寸有关，网格大小不同时，应力应变关系也不同，本文通过建模分析时发现钢筋应变达到0.066作为柱的破坏准则。因此，为减少建模工作量及结果分析的准确性，采用直接输入断裂能，按欧洲规范CEB-FIP MC90 [ 9 ] 建议计算： G f = α ( 0.1 f c ) 0.7 式中： f c 单位为MPa； G f 单位为N/mm；系数α与混凝土最大骨料直径D max 有关，CEB-FIP MC90中建议，D max = 8 mm时，α = 0.025；D max = 16 mm时，α = 0.03；D max = 32 mm时，α = 0.058。 3.1.4. 其他参数 混凝土粘滞参数影响着ABAQUS计算收敛，太大结果不精细，太小不易收敛，本次模型取值0.05，具体参数如表3所示。 Table 3 Concrete plastic properties parameter 膨胀角 偏心率 f b 0 /f c 0 k c 粘性参数 30 0.1 1.16 0.667 0.05 表3. 混凝土塑性性能参数 3.2.钢筋本构 对于ABAQUS自带的钢筋本构，属于纯钢筋滞回规则，不能很好的模拟出“捏拢”效应，因此本文通过子程序PQ-Fiber [ 10 ] [ 11 ] [ 12 ] 引入钢筋本构通过削弱钢筋卸载刚度来模拟这种效果。PQ-Fiber是清华大学土木工程系结构工程研究所潘鹏副教授和曲哲博士基于大型通用有限元程序ABAQUS开发的一组材料单轴滞回本构模型的集合。选用再加载刚度按Clough本构退化的随动硬化单轴本构模型USteel02，如图3所示。 图3. 往复加载时的单轴应力应变关系 3.3. 模型的建立 有限元ABAQUS混凝土单元选用C3D8R，钢筋选用桁架单元T3D2，钢筋与混凝土之间的约束采用“Embeded Region”。为防止RF1，RF2两个参考点直接加载在混凝土柱墩顶部出现应力集中及过约束，建立了正六面体钢板，通过参考点同时施加X、Y方向荷载，耦合到钢板的侧面底边线中点。模型初始分析步约束基座侧面及底面的各方向位移，模拟固定端。第一个分析步，同时施加柱顶的轴压及重力，使关系平稳的建立起来，第二个分析步施加两方向的水平推力。柱墩网格划分先分割部件，再按120 mm设置全局种子，钢筋按100 mm分网，其他辅助加载部分适当加大网格尺寸，A1共3651个单元，模型建立如图4所示。
system: [SYSTEM]

user: 4. 双塑性铰位置分析 塑性铰是指在弯矩基本保持不变的情况下，截面曲率激增，形成了一个能转动的“铰”。塑性铰是在破坏开始时形成的，故将其作为试件破坏的标志。双塑性铰即在一个设计截面形成了两个能转动的铰，通过实验与模拟对比分析，本文判定双塑性铰产生的条件是钢筋应变在破坏处同时达到0.0667。A2，B2组试件在原实验的基础上进行了部分钢筋截断 [ 13 ] [ 14 ] 如图5，使其在截断处与柱底产生双塑性铰，截断钢筋用红色标记，轴压比等其他参数与实验的参数一致，具体信息见表4。 图4. 有限元模型 图5. 双塑性铰设计截面 Table 4 Column design parameter 组 编号 轴压比 加载高度(m) 截断高度(m) 绞的个数 配筋率 实验与模拟 C C1 0.2 3.2 0 1 1.437% 实验及模拟 C2 3.2 0.8 2 柱底1.437%， 截断处0.617% 模拟 E E1 0.1 2.2 0 1 1.437% 实验及模拟 E2 2.2 0.55 2 柱底1.437%， 截断处1.026% 模拟 表4. 柱的设计参数 4.1. 应变曲线 纵筋的应变显示A2、B2在墩底和1/4柱高处，即分别距墩底800 mm、550 mm处，形成了两个塑性铰(图6)。两个塑性铰位置钢筋最大应变均达到0.066。 图6. 双塑性铰与柱底塑性铰应变值对比 4.2. 混凝土应力云图 计算所得混凝土桥墩A1和B1的应力云图如图7，可以看出当出现塑性铰时混凝土应力约达到抗压强度f c 的50%，此时钢筋已经屈服破坏出现在柱底及截断处的对角边，与实验结果相符。A2和B2桥墩的应力云图显示，在底部区域形成了两个塑性铰，和图6中的应变数据对应。 图7. 双塑性铰与柱底塑性铰应力云图
system: [SYSTEM]

user: 5. 模拟结果与实验结果对比分析 5.1. 滞回曲线 计算所得曲线如图8所示。为了和模型实验对比，模拟的结果加了试验与千斤顶之间滑动摩擦力，摩擦系数为0.03。 图8. 双塑性铰与柱底滞回曲线对比 由上图可以看出，实验与计算曲线整体吻合程度较好。模拟曲线初始刚度偏大，可能是混凝土材料本构曲线较为复杂，初始弹模取值较大造成的。滞回曲线的捏拢效应没有试验的明显，可能原因是有限元没有考虑钢筋与混凝土之间的粘结滑移，后续有待研究。对比单双铰模型结果，可以看出双塑性铰模型滞回圈更饱满，抗震性能更好，但最大抗剪承载力无明显变化。 5.2. 骨架曲线 骨架曲线是通过连接滞回曲线中各级位移加载峰值时对应的荷载值，平滑连接，形成的一条滞回曲线的包络线。如图9所示，A1的极限位移为60 mm，A2的极限位移提高了75%，达到105 mm。B1的极限位移为40 mm，B2为60 mm，B2的极限位移比B1提高了50%。可以看出双塑性铰相比传统柱底出现的单铰增大了柱墩的极限位移。 图9. 双塑性铰与柱底塑性铰骨架对比
system: [SYSTEM]

user: 6. 结论 在模型试验的基础上，本文利用ABAQUS建立了双向往复荷载作用下混凝土桥墩柱的非线性有限元模型，进而设计了双塑性铰柱墩，并与单塑性铰柱墩进行了对比，得到如下结论： 1) 通过有限元计算与实验结果比较，ABAQUS计算所得的滞回曲线、骨架曲线与实验大致吻合，表明通过本文建立的非线性有限元模型能较好地反映钢筋混凝土墩柱的受力性能。 2) 通过设计截断钢筋截面的不同的位置、数量和配筋率，表明可以在混凝土柱墩底部形成两个塑性铰。 3) 通过对比单双塑性铰柱墩，结果表明在墩底形成双塑性铰的情况下桥墩具有更好的极限变形能力，混凝土柱墩的抗震性能得到增强，但极限承载力无明显变化。