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

user: 1. 引言 隐身能力已经成为衡量现代军事飞行器性能的重要指标之一，等离子体被称为物质第四态，由于在电磁方面具有特殊的物理性质，成为了隐身材料研究中热门的课题。目前，等离子体对电磁波的隐身性能的研究主要有三个方面：折射隐身、吸收隐身和反射隐身。 文献 [ 1 ] 详细讨论了电磁波在不同折射率时等离子体的折射隐身规律，利用费马原理和变分法给出了不同折射率时等离子体中的电波轨迹，文献 [ 2 ] 讨论了电磁波在等离子体中吸收隐身规律。当电磁波在等离子体中传播时，相当大部分的能量将被等离子体吸收，其机理是电磁波的电场使等离子体中电子加速获得动量，相互之间产生电磁振荡和粒子的有效碰撞，成而使电磁波的能量发生衰减，这就是所谓的碰撞吸收。由于碰撞吸收可以使电磁波的能量产生较多的衰减，因此吸收隐身成为了研究的重点。文献 [ 3 ] [ 4 ] 从等离子体对电磁波碰撞吸收的角度讨论了当等离子体密度呈线性和指数分布时，等离子体对电磁波衰减的影响。文献 [ 5 ] [ 6 ] [ 7 ] 讨论了等离子体密度呈抛物线分布时，讨论了电磁波在等离子体中衰减与等离子体碰撞频率、电磁波频率、等离子体最大密度、电磁波入射角等的关系。然而在实验室中产生的等离子的电子密度常取Epstein分布，文献 [ 2 ] 虽然讨论了等离子体密度呈Epstein分布时电磁波与等离子体之间的相互作用，但是作者从反射角度讨论了电磁波在等离子体中的传播特性，而且理论基础也比较简单，因此有必要从吸收角度在理论上进行深入探讨其微观作用机理，同时对这三种等离子体密度分布衰减进行分析比较，为未来实际应用提供参考。 WKB是一种普遍采用计算束缚态能量和势垒穿透率的近似计算方法。在研究电磁波在非均匀等离子体中传播碰撞吸收特性时，可将等离子体的介电常数看成随着电磁波的波长变化，同时忽略在传播过程中的耦合效应。本文首先将具有一定厚度的等离子体看成密度分布是不均匀的，然后将其构建为不均匀分布的等离子体层，利用Matlab进行编程和数值模拟，计算了三种不同类型等离子体密度分布对电磁波碰撞吸收衰减的影响，并且对相关结果进行了分析比较。
system: [SYSTEM]

user: 2. 模型建立和波动方程求解 [ 5 ] [ 6 ] [ 7 ] 图1给出了等离子体与电磁波相互作用的几何模型。金属目标物表面被包裹着一层厚度为 r 0 的非磁性不均匀等离子体。将等离子体看作是平面分层各向同性介质。沿着z轴方向，等离子体呈不均匀分布，而在垂直于z轴即在xy平面内分布是均匀的。 当等离子体内部收到外部扰动时，等离子体将产生集体振荡，相应的振荡频率称为等离子体频率 ω p ，它表征等离子体物性的重要参数之一。其表达式为 ω p = n e e 2 / ε 0 m e (1) 式中e为电子的电荷量， ε 0 为真空介电常数， m e 为电子质量。 图1. 电磁波在金属衬底等离子体中传播示意图 一般来说，电子和中性分子的碰撞频率 ν 是电子密度 n e 和温度T的函数，实际应用中通常取经验值，用 ν e 表示。当考虑碰撞时，非磁化等离子体的介电常数 ε 可表示为 ε = ε r + i ε i = ε 0 ( 1 − ω p 2 / ( ω 2 + ν e 2 ) + i ( ν e / ω ) ω p 2 / ( ω 2 + ν e 2 ) ) (2) 在耗散介质中，波数 k = k r + i k i 为复数，实部 k r 表现了波在传播过程中的相位延迟，虚部 k i 表现了波在传播过程中的幅度衰减，结合 k 2 = ω 2 μ ε 以及式(2)得到 k r = ω ( μ / 2 ) 1 / 2 [ ε r + ( ε r 2 + ε i 2 ) 1 / 2 ] 1 / 2 (3) k i = ω ( μ / 2 ) 1 / 2 [ − ε r + ( ε r 2 + ε i 2 ) 1 / 2 ] 1 / 2 (4) 可看出，耗散介质中波数k为入射电磁波频率 ω 和等离子体密度 n e 的函数。 根据模型仅考虑垂直入射面线性极化的平面波，则电场E仅在y轴方向有分量。此时波动方程可表示为： d 2 E y d z 2 + Q 2 ( z ) E y = 0 (5) 式中： Q 2 ( z ) = k 2 ( z ) − ω 2 c 2 sin θ n (6) 根据菲涅尔定律： 1 ⋅ sin θ i = n n sin θ n (7) 其中 n n 、 θ n 和 θ i 表示第n层等离子体的折射率、入射角和透射角。 根据方程(6)、(7)解方程(5)，解得斜入射情况下的WKB近似解为 E y = E 0 exp [ ∓ i ∫ 0 r 0 ( k 2 ( z ) − ω 2 c 2 sin 2 θ n ) d z ] (8) “−”代表正向波，“+”代表反向波。WKB有效条件为 1 k 0 2 | 3 4 ( 1 Q 2 d Q d z ) 2 − 1 2 1 Q 3 d 2 Q d z 2 | ≪ 1 (9) 出射到空气中的电磁波电场强度可以由式(8)得出 E = E 0 exp [ − 2 i ∫ 0 r 0 ( k 2 ( z ) − ω 2 c 2 sin 2 θ n ) d z ] (10) 当电磁波从 z = r 0 处以入射角 θ n 传播到等离子体内部并在等离子体表面即 z = 0 处被反射回来，当电磁波再次传播到等离子体的表面时其能量可表示为 P ( r 0 ) = P 0 exp [ − 4 Im ( ∫ 0 r 0 ( k 2 ( z ) − ω 2 c 2 sin 2 θ n ) d z ) ] (11) 这里 P 0 是电磁波入射时 z = r 0 处的能量。 电磁波的双程衰减为 A t t = | 10 lg P ( r 0 ) P 0 | = | 17.37 Im ( ∫ 0 r 0 ( k 2 ( z ) − ω 2 c 2 sin 2 θ n ) d z ) | (12) c为真空中光速。
system: [SYSTEM]

user: 3. 计算模拟与结果分析 根据以上理论分析可知，等离子体隐身的效果与电磁波的频率具有密切的关系。因此，为了更直观的揭示它们之间的相互作用机理，本文从以下角度进行了分析和研究。首先我们将等离子体密度分别取以下三种分布： 线性分布： n e ( z ) = n e 0 ( z / r 0 ) 。 抛物线分布： n e ( z ) = n e 0 ( z / r 0 ) 1 / 2 。 Epstein分布： n e = n e 0 [ 1 + exp ( − ( z − r 0 / 2 ) / z 0 ) ] − 1 。 计算中，等离子体层厚度 r 0 = 0.1   m ， z 0 = r 0 4.0 × ln 10 ≈ 0.0109 ， n e 0 为等离子体峰值密度，e和 m e 分别为单电子的电量和质量。 3.1. 电磁波入射频率对衰减的影响 图2揭示了垂直入射的电磁波在等离子体密度分别为线性、抛物线和Epstein分布时相互作用的特点。其中 n e 0 取值10 11 m −3 ， ν e 取值10 11 Hz。当入射波频率较小时，等离子体对电磁波的吸收能力都较强。随着入射波的频率不断增大，存在明显的最大吸收峰。随着入射波频率不断变大，等离子体对电磁波的吸收在减弱。以上现象产生的原因可解释为，当入射波的频率接近或等于等离子体振荡截止频率时，等离子体和电磁波产生共振吸收效应，从而对电磁波的吸收达到最大值；当入射波频率小于或大于等离子体振荡截止频率时，等离子体中的电子对电磁波的相应不同步，等离子体内的电子极化现象较弱，从而使等离子体中的电子对电磁波能量的吸收大大减弱。通过图2进一步看出，当等离子体密度为抛物线分布时，入射电磁波的衰减最为明显。 3.2. 等离子体电子碰撞频率对衰减的影响 将等离子体三种分布设定了相同的峰值密度( n e 0 = 10 17     m − 3 )，考虑垂直入射到非磁化等离子体中电磁波的碰撞吸收情况。图3反映了当碰撞频率较小时，等离子体对电磁波的吸收频谱较窄，随着等离子体振荡频率增大，电磁波的吸收衰减系数增加，等离子体吸收峰值向更低的频率移动，电磁波的吸收频谱范围增加，当等离子体电子碰撞频率增大时，等离子体对电磁波的吸收能力不增反减，这表明对于一定的等离子体密度分布模型，等离子体对电磁波的吸收能力存在一个最优化的电子碰撞频率，在该碰撞频率下等离子体层对电磁波的吸收能力最强。通过对图3三种分布的类比可以知道，当等离子体呈抛物线分布时，等离子体对电磁波具有更好的吸收效果。 图2. 电磁波入射频率对衰减的影响 图3. 电磁波的衰减与电磁波频率的关系 3.3. 等离子体峰值密度对衰减的影响 将线性分布、抛物线分布和Epstein分布设定了相同的碰撞频率( ν e = 10 10     Hz )，考虑垂直入射到非磁化等离子体中电磁波的碰撞吸收情况。 图4揭示了当入射到等离子体层中的电磁波频率较低和等离子体峰值密度越大时，电磁波在等离子体层中衰减越明显。同时可以发现，当等离子体峰值密度较高时，对电磁波的吸收衰减越强，吸收频谱也越宽。从图4三种分布数值模拟的结果可知，抛物线分布等离子体层比另外两种分布的吸收效果更好。 图4. 等离子密度不同时，衰减随入射波频率的变化关系 3.4. 入射角度对衰减的影响 将线性分布、抛物线分布和Epstein分布设定了相同的峰值密度 n e 0 = 10 17     m − 3 、碰撞频率 ν e = 10 10     Hz ，考虑斜入射到非磁化等离子体中电磁波的碰撞吸收情况。从图5中可以看出，当电磁波的入射角增大，等离子体对电磁波的衰减能力随电磁波频率的变化趋势是一致的。当电磁波的入射角度较小时，电磁波在等离子体中传播的距离也相对较短，因此其衰减程度也越缓慢。当等离子体密度呈抛物线分布时，等离子体对电磁波具有更好的吸收效果。 图5. 入射角度不同时，电磁波的衰减与电磁波频率的关系
system: [SYSTEM]

user: 4. 结论 在构建等离子体与电磁波相互作用的几何模型基础上，利用WKB近似方法进行了理论推导，得出了电磁波在非磁性不均匀等离子体中传播时的双程吸收衰减和其他参数之间的联系，利用Matlab软件进行了数值模拟。研究发现，如果调整等离子体的峰值密度，当入射电磁波的频率较小时，等离子体峰值密度越大，等离子体对电磁波的衰减能力越强；当等离子体的碰撞频率取较低的值时，对电磁波的吸收频谱将变窄，当等离子体电子碰撞频率数值增加时，对电磁波的吸收衰减也增加；当电磁波的入射角度不同时，电磁波入射角度越小，衰减越慢。通过以上分析可知，无论是改变等离子体峰值密度、等离子体碰撞频率，还是改变电磁波的入射角，三种分布中，等离子体呈抛物线分布时，等离子体对电磁波具有更好的吸收效果。这主要是因为等离子体对电磁波的吸收衰减能力取决于边缘电子密度大小和电子密度梯度之间的关系，当等离子体的边缘密度缓慢增加而不是陡增时，使更多电磁波容易进入等离子层参与了衰减，因此等离子体对电磁波的吸收衰减能力增强。
system: [SYSTEM]

user: 摘要 为探明不同肥料品种对耕地的培肥效果，选用羊粪、商品有机肥作为有机改良剂，尿素、包膜尿素、磷酸二铵、过磷酸钙、硫酸钾、氯化钾，硅肥、硼肥等化肥作为无机改良剂，以单施尿素 + 磷酸二铵 + 硫酸钾为对照，通过水稻–油菜盆栽试验，研究其对土壤理化性质及水稻、油菜产量的影响，得出以下结论：1) 在2年水稻–油菜复种中，随着种植年限的增加，土壤有效磷、速效钾含量均呈增加趋势；2) 施用有机肥处理土壤有机质含量高于单施化肥处理，施用羊粪处理与施用商品有机肥处理土壤有机质含量差异不显著。3) 施用尿素处理土壤硝态氮与施用包膜尿素处理差异不显著；4) 施用磷酸二铵处理土壤有效磷含量高于施用过磷酸钙处理，但二者之间无显著差异；5) 综合土壤养分含量及作物需求评价各肥料组合，包膜尿素 + 磷酸二铵 + 氯化钾 + 羊粪为最佳的肥料组合。
system: [SYSTEM]

user: 1. 引言 耕地是关乎国家粮食安全的重要资源，虽然我国拥有960万平方公里的广袤国土，但是在14亿人口的基数下，人均土地面积严重不足，而适用于农业种植的耕地更显稀少。随着我国科技、经济的发展，城市的扩张进一步侵占农业用地面积，耕地呈现逐年减少的趋势 [ 1 ]。在坚守18亿亩耕地红线不动摇的政策下，如何合理地开发与保护耕地，提升耕地土壤质量，是维护我国粮食安全的重大问题。 在我国耕地无法大幅增长的前提下，提高耕地质量是保障农作物产量的有效方法。经过连年耕作或者是新增耕地的条件下，耕地往往存在肥力低下、养分不平衡等问题，如土壤有机质含量低，氮磷钾等养分缺乏或者不平衡，严重影响着耕地质量，需要通过适宜的培肥措施调节土壤养分，以保障作物的高产、稳产。肥料种类是影响培肥效果的关键因素之一，肥料种类不同，养分释放的速率、在土壤中的运移速度、作物吸收利用的难易程度不一，直接影响着土壤肥力的累积速度及作物生长发育。研究表明，施用有机肥可以为作物提供必要的氮磷钾等养分 [ 2 ]，同时有机肥的施用可以改良土壤空隙结构，提高土壤通气性和透水性，有利于农作物根系在土壤中的延伸，从而提高作物对水分的养分的吸收能力，进而提高作物产量 [ 3 ]。氮素是植物生长所必需的重要元素，不同的氮肥形态会对植物的生长发育产生影响 [ 4 ]，以铵态氮为主要氮源时，植物叶片生长一般会受到抑制，而硝态氮作为氮源的情况下，植物的总生物量和根冠比较高 [ 5 ]，根系中糖类和代谢酶活性较强 [ 6 ]，氮素的供应水平也会极大的影响作物产量。研究表明，磷钾的供应也是影响植物生产和果实产量的重要因素，磷素供应不足会对水稻根系产生负面影响，显著增加水稻的根冠比 [ 7 ]，抑制地上部发育。而钾素作为产量要素，供应不足往往会导致作物减产 [ 8 ]。微量元素硅和硼在水稻和油菜的生长发育过程中也起着重要作用，如促进茎秆的发育和提高抗逆性能等。 本研究拟在重构土体的基础上，通过添加不同类型的肥料提高土体肥力，以保证水稻的正常生长需求。研究通过盆栽试验探明不同肥料类型对土壤理化性质、水稻、油菜产量和养分含量的影响，从而确定重构土体上适宜添加的肥料组合，为黄河滩地重构土体快速培肥提供理论依据。
system: [SYSTEM]

user: 2. 材料与方法 1) 试验设置 盆栽试验采用高50 cm、直径30 cm的塑料桶，下部填充10 cm沙土，中部为5 cm容重为1.7 g/cm 3 黄土构建的人造犁底层，上部为25 cm容重为1.4 g/cm 3 的黄土耕作层，供试土壤理化性质见表1。每盆种植水稻3穴，每穴种4株，根据秧苗大小保持1~5 cm水面(以不淹苗为准)。油菜季每盆种植油菜3株。种植期间各处理保持灌溉量一致。 试验以施用常规肥料尿素、磷酸二铵、硫酸钾为对照(CK)，选取羊粪、商品有机肥2种有机改良剂和包膜尿素、过磷酸钙、氯化钾、硅、硼5种肥料，组成12个处理，每个处理重复3次(处理设置见表2)，共计36盆，随机区组排列。硅、硼两因素分别设置不施用硅硼处理(CK)、单施硅肥(8 mg/kg)处理、单施硼肥处理(4 mg/kg)和硅肥(8 mg/kg)、硼肥(4 mg/kg)均施处理。各处理有机改良剂施用量保持氮、磷、钾养分含量一致(相当于4 g/kg羊粪)，化学改良剂施用量保持氮、磷、钾养分一致(0.3 g/kg N， 0.06 g/k P 2 O 5 ， 0.06 g/kg K 2 O)。水稻种植时，磷肥、钾肥与有机肥作为基肥在水稻种植前施入土壤。插秧后2~3天施入20%的氮肥，返青期(插秧后10~12天)施入50%氮肥，抽穗期将剩余氮肥施入土壤，50%的硅肥、硼肥作为基肥施入土壤，剩下的作为追肥在拔节期施入土壤。油菜种植时，氮、钾肥总量的60%、全部的磷肥作基肥施用，20%的氮、钾肥和全部的硼肥作腊肥(移栽后50天)，20%的氮、钾肥作薹肥(移栽后80天)施用。 Table 1 Physical and chemical properties of tested soi 土层 有机质(g/kg) 硝态氮(mg/kg) 铵态氮(mg/kg) 有效磷(mg/kg) 速效钾(mg/kg) 0~20 cm 7.91 13.23 6.74 9.89 126.72 表1. 供试土壤理化性质 Table 2 Types of fertilizers for treatmen 处理 肥料种类 硅酸钠(以Si计) 硼酸钠(以B计) 有机肥 CK 尿素 磷酸二铵 硫酸钾 — — — 1 尿素 磷酸二铵 硫酸钾 — — 羊粪 2 包膜尿素 磷酸二铵 硫酸钾 — — 羊粪 3 尿素 过磷酸钙 硫酸钾 — — 羊粪 4 尿素 磷酸二铵 氯化钾 — — 羊粪 5 尿素 磷酸二铵 硫酸钾 — — 商品有机肥 6 包膜尿素 磷酸二铵 硫酸钾 — — 商品有机肥 7 尿素 过磷酸钙 硫酸钾 — — 商品有机肥 8 尿素 磷酸二铵 氯化钾 — — 商品有机肥 9 尿素 磷酸二铵 硫酸钾 8 mg/kg — — 10 尿素 磷酸二铵 硫酸钾 — 4 mg/kg — 11 尿素 磷酸二铵 硫酸钾 8 mg/kg 4 mg/kg — 表2. 化肥种类处理设置 2) 测定项目 土壤指标：水稻/油菜收获后采集表层20 cm土壤样品，测定有机质、硝态氮、铵态氮、有效磷、速效钾。有机质采用重铬酸钾–外加热法测定，硝铵态氮采用流动分析仪测定，有效磷采用钼锑抗比色法测定，速效钾采用火焰光度计测定。 水稻和油菜成熟后，整株挖出，根系用清水冲洗。样品晾干后在40℃下烘干，将样品分为根、茎、叶和稻米。测定水稻根系和地上部分生物量及水稻产量。
system: [SYSTEM]

user: 4. 结论 1) 在2年水稻–油菜轮作中，随着种植年限的增加，土壤有效磷、速效钾、含量均呈增加趋势。经过2年的水稻–油菜轮作，土壤有效磷的增幅最大，相比2017年油菜收获后，2017年水稻收获后、2018年油菜收获后、2018年水稻收获后土壤有效磷含量分别增加108.9%、47.7%、17.9%。 2) 施用有机肥处理土壤有机质含量高于单施化肥处理，施用羊粪与施用商品有机肥处理土壤有机质含量差异不显著。 3) 施用尿素处理土壤硝铵态氮含量高于施用包膜尿素处理，但二者无显著差异。种植水稻时，土壤长期处于淹水状态，土壤中的硝态氮易发生反硝化作用或随水淋失，因此建议施用包膜尿素，使土壤氮素缓慢释放，降低氮素转化与损失。 4) 施用硫酸钾处理土壤速效钾含量略高于氯化钾，但差异不显著，考虑到种植水稻时土壤长期处于淹水状态，硫酸根易转化为硫化氢，对作物产生毒害，建议选择氯化钾作为钾肥。
system: [SYSTEM]

user: 摘要 随着我国人民收入和生活水平的提高，以及全面二胎政策的出台和老龄化问题的加重，对家政服务业发展提出更高要求。同时，家政服务业作为第三产业，对解决就业和消除贫困具有积极意义。本文根据我国家政服务业在人员构成、培训、法律等方面存在的问题，指出政府在行业发展中作用的发挥。
system: [SYSTEM]

user: 1. 引言 改革开放以来，我国各项事业均取得伟大成功，社会经济得到长足发展，人民收入日益增长。根据《2018年国民经济和社会发展统计公报》显示 [ 1 ]，2018年我国居民人均可支配收入达到28，228元，比上年增长8.7%，人们的物质生活极大丰富的同时，人们对生活的质量提出更高要求，为我国家政服务业，尤其是高端家政服务提供了广大市场。
system: [SYSTEM]

user: 2. 我国家政服务行业发展趋势与必要性 2016年，二胎生育政策全面实施，虽然生育率在当年出现高峰后连续两年下降，但在政策鼓励下，新生儿数量势必会攀升。人口老龄化问题日趋突出，老年人口数量逐年增长，截至2018年年底，我国60岁及以上老年人口约2.49亿，占总人口的17.9%。加之伴随我国家庭小型化及知识女性逐渐走出家庭主妇定位等现实原因，使得月嫂、育儿嫂、养老护理员、家庭教师、家庭财务规划师、家庭管家等专业化高端家政服务需求日益增加。 家政服务业作为市场需求量巨大的第三产业，具有强大的就业吸纳能力。根据《2017年中国家政服务行业发展现状分析及未来发展趋势预测》数据测算，2016年，全国家政服务业从业人数为2542万人，同比增长9.3%。劳动与社会保障部的调查显示，2017年家政行业市场从业人数近3000万人，但家政服务业供应量仍存在巨大缺口。面对家政服务的大量需求，家政服务业将继续为吸纳社会剩余劳动力贡献力量。 此外，工信部指出我国农村家庭收入的三分之一来自进城务工，而家政服务业是农村进城务工人员职业的主要选择之一。2020年是脱贫攻坚的收官之年，家政服务业在解决就业问题的同时，也是政府消除贫困的渠道之一。而家政服务业在我国发展总体处于起步阶段，在发展过程中仍存在诸多问题，面对市场对家政服务业更加专业化和职业化的要求，政府应采取多方面积极措施，促进家政服务业的发展。
system: [SYSTEM]

user: 3. 家政服务业在发展中存在的问题 3.1. 人员构成复杂，综合素质普遍不高 我国家政服务人员主要来源于城镇下岗职工、退休人员、农村外来务工人员以及一部分失学的未成年人 [ 2 ]。这些从业人员绝大部分文化水平不高，在从事工作时大多也未经过相关技能培训，完成工作仅依靠个人生活经验。从事家政服务行业人员年龄普遍偏高，在从事专业性要求较高的岗位并对其进行培训时，对知识的接受能力差，导致从事对专业性要求较高的岗位如月嫂、育儿嫂、护理等人员紧缺。 3.2. 培训内容不规范，考核制度不统一 家政服务业的规范专业培训教材匮乏，我国各地对家政服务行业的从业人员进行的相关培训内容各异，甚至同一城市不同公司培训内容也不尽相同。培训机构为获得政府补贴，在培训过程中往往流于形式，大多只进行简单的岗前培训，时间短，内容单一，只重视理论，忽视实际操作，在课时分配上没有参照国家规定。而且在培训过程中没有重视家政服务人员的道德培养、人际沟通能力、安全防范能力 [ 3 ]。 3.3. 法律法规不完善，行业发展不规范 现阶段我国家政服务业没有明确的能够完全适用的法律，主要依据的法律仅为《劳动法》、《劳动合同法》，而这也仅限于是会员制和员工制的家政公司 [ 4 ]。我国现阶段家政服务业发展形式主要还是中介制，而中介在发生纠纷时不承担任何法律责任，家政服务人员与客户之间的劳动关系往往只依靠口头约定，而不签订劳动合同，二者一旦发生纠纷，劳动合同法也不在适用范围内，且家政服务人员的各项社会保障也难以得到保障。 3.4. 社会认可度低，人员流动性大 在传统观念中，家政服务往往被人们简单的认为等同于保姆，保姆的工作包括洗衣、做饭、收拾屋子等一切室内家务，有时甚至还涵盖照顾儿童或老人及外出采买食材。而受封建思想的影响，保姆工作就是旧社会中的“老妈子”，所以其在社会中的地位比较低下，在工作过程中常常不被人们尊重，社会认可度低。虽然现在家政服务逐渐趋于专业化，技术含量不断增加，但社会现状仍然很难改变。正因为职业的社会认可度低，加之社会保障和薪酬难以达到与付出相对等，导致人员流动性极大。 3.5. 政府扶持政策落实不详，行业监管措施不到位 近年来，国家针对家政服务业的发展下发了许多相关文件，各省市也根据自身发展现状及特点制定了相应的发展措施和扶持政策，这一定程度地促进了家政市场的活跃程度。但是，政府在对家政服务业进行政策扶持时，对相关优惠政策是否落到实处缺乏了解，导致优惠政策“有名无实”。此外，家政服务行业在发展过程中缺乏完善的监督体系，监管措施不到位，作为具有直接监督责任的行业协会往往不作为，导致家政服务行业乱象丛生，发展缓慢。
system: [SYSTEM]

user: 4. 政府在家政服务业发展中作用的发挥 4.1. 严格把关市场和行业准入，明确服务人员准入条件 为促进家政服务市场发展，政府降低了家政服务企业的市场准入门槛，在财政及企业审批方面出台许多优惠政策，促使大量企业涌入该行业。因此，政府在各类家政企业申请进入该行业时，应严格地对家政市场和行业准入进行把关，对注册公司的相关资质细致审查，在源头上保证市场和行业有序合规。 为适应家政服务业发展现状及未来趋势，针对当前家政服务人员低文化、低技能、低素质等现状，政府应明确家政服务人员进入该行业的条件，规定只有经过岗前培训和相关岗位技能及素质等培训并取得证书的人员才能进入该服务岗位，保证人员持证上岗。 4.2. 指定培训机构，规范培训内容及考核标准 为提升家政服务人员综合素质，提高家政服务质量，政府对家政服务的培训机构出台诸多优惠政策，在财政方面采取财政补贴或减免税收，使得许多培训机构为获得政府财政补贴或减税进入家政服务培训市场。面对层出不穷的培训机构，地方政府应指定一家或几家权威培训机构，统一培训内容及考核标准，注重理论与实践相结合。培训和考核结束后，由权威培训机构发放统一的从业及等级证书，提高证书含金量，规范家政培训市场。 4.3. 完善相关法律法规，切实保障各方权益 我国缺少完全适用于家政行业的法律，致使发生纠纷时，在处理案件时无法可依，各方权益无法得到保障，从而不利于社会的稳定。故政府应积极推动家政服务业专门立法，建立符合我国国情的家政服务合同，将家政服务人员纳入到法律保护范围内的同时，也同样保护客户合法权益。 4.4. 推动消除社会偏见，树立正确家政观念 我国当前社会对家政服务普遍存在偏见，甚至家政服务人员自身也对自己所从事的行业不认可。面对社会刻板印象，政府应利用电视、广播、报纸以及新媒体加强对家政服务业的宣传，引导大众摒除落后的思想观念，以进步和平等的眼光看待家政服务业和家政服务人员。将家政教育引进中小学课堂，使人们从小开始接触和了解家政，既可以培养学生独立的生活技能，又使之树立正确的家政观念，并且在高考填报志愿时可以以平常心报考家政专业，从而扩大家政服务群体，使家政服务人员整体素质提高。大中专及高等院校开设高端的家政专业和家政管理专业，为家政服务业培养高素质人才。 4.5. 加强行业监管，确保政策落实 确立家政服务行业管理机构，明确行业责任主体，完善家政服务行业协会发展，加强行业的监管，不定时对家政公司、家政中介进行抽检，加大对市场上违法违规经营的家政公司和中介的处罚力度，促进家政服务行业健康有序发展。对培训机构进行评估考核，定时核验其培训资质。此外，政府在对培训机构财政补贴后，应对资金流向进行监督，保证资金落在实处。 4.6. 鼓励家政方面科学研究，实行家政品牌战略 我国对家政服务方面的研究起步较晚，科研成果和建树相对较少，在研究深度和广度上均落后于美国、日本被等发达国家及以菲佣为主要产业的菲律宾。政府应为高校中家政服务专业相关的师生、家政领域的学者、从事家政服务业相关的工作人员等创造良好的科研环境，加大在财政上给予的资金支持，为有突出贡献者给予科研奖励，鼓励各方人士对家政方面进行科学研究，促进家政服务业在理论方面的丰富与进步。 同时，各地政府应根据城市发展的现状和特点，人民的生活水平和需求，打造符合本城市的家政品牌。随着人民生活水平的不断提高，人们对生活水准和服务质量提出更高要求，在家政服务方面也更注重服务水平和口碑。此外，近年我国在科学技术方面取得长足进步，各类生活均向智能化发展，O2O线上家政服务已逐渐走进人们生活，各地政府应根据人们需求，顺应未来发展趋势，实行家政品牌战略，促进家政服务业在实践方面的发展。
system: [SYSTEM]

user: 5. 结语 家政服务业作为第三产业的重要组成部分，具有广阔的发展前景，属于未来社会发展的朝阳产业，家政服务业也是满足人民日益增长的美好生活需要的重要载体。政府应在家政服务业发展中发挥积极作用，从而优化产业结构，扩大就业渠道，进一步拉动内需，改善民生，提升百姓生活福祉。
system: [SYSTEM]

user: 摘要 本文利用2017年川渝地区71个地面站点的实测降水数据，在年、月和日尺度上验证了TRMM 3B42V7降水数据的精度，分析TRMM卫星降水数据在该地区的适用性。结果表明，在年、月和日尺度上，TRMM 3B42V7与实测降水数据之间均呈显著线性关系，但回归方程斜率都较小，分别为0.249、0.290和0.198，且降水量相对误差较高，分别为−68.60%、−65.88%和−65.77%。这说明TRMM 3B42V7降水数据显著低估实际降水量，但通过线性回归方程可估算出实际降水量。就单个站点而言，月尺度上绝大部分站点降水实测数据与TRMM 3B42V7的相关系数较高，平均值高达0.915；而在日尺度上气象站点的降水实测数据与TRMM 3B42V7降水量之间的相关系数平均值下降为0.622。无论是月尺度还是日尺度上，峨眉山站的相关系数最小，分别为0.628和0.191。
system: [SYSTEM]

user: 1. 引言 降水是研究区域和全球气候变化的基本要素之一，也是地表与大气间能量与水循环过程的重要组成部分 [ 1 ]，因此准确掌握降水信息对研究不同时空尺度下的大气循环和水循环过程都起着关键作用 [ 2 ]，特别是对于区域防灾减灾具有重要指导作用。 目前降水数据的获取方式主要包括地面气象站点观测、地面测雨雷达和卫星遥感。尽管地面气象台站和地面测雨雷达精度都比较高，但均受制于观测点稀少、空间代表性差，且其空间插值精度非常有限，严重限制了其应用范围。相比之下，卫星遥感获取的降水数据具有高时空分辨率、大覆盖范围以及不受地形地貌条件限制的特点，其提供的降水时空分布信息具有明显的优势 [ 3 ]，其中典型的是由美国和日本联合研制的热带降雨测量卫星计划TRMM (Tropical Rainfall Measuring Mission)卫星在日本发射，为全球降水的高时空分辨率探测提供了新的途径和方法 [ 4 ]，目前已广泛应用于水文过程模拟 [ 5 ]、降雨分布特征 [ 6 ]、天气过程分析 [ 7 ] 等研究中。与此同时，TRMM卫星降水数据的可靠性与精度问题也受到了广泛关注。 川渝地区地貌地形复杂且受多个气候系统的影响，降水空间分布极不均匀，是我国降水局部区域差异最大、变化最复杂的地方之一。有研究表明在年、季和月时间尺度上，TRMM 3B43降水数据在川渝地区有较高精度，具有良好的适用性 [ 8 ]。该结论是基于2000~2011年川渝地区72个气象站点观测数据，对卫星降水数据进行检验得到的，是对多年降水评估得到的平均结果。但是川渝地区年际间、季节内降水差异均较大，TRMM卫星数据在某一年内精度如何，关于这方面的研究比较少见。本文通过比较2017年TRMM卫星数据与实测降水数据，以期进一步评估TRMM卫星数据在川渝地区的可靠性与适用性，从而为应用该卫星数据对该区域水文、气象和防洪减灾等模拟及预报工作奠定基础。
system: [SYSTEM]

user: 2. 资料与方法 本文所用数据包括地面气象站点实测降水数据和TRMM 3B42V7数据。实测降水观测数据来自中国气象数据共享服务中心，选取2017年川渝地区71个气象站点观测数据，川渝地区地势概况以及71个气象站点的分布如图1所示。对卫星降水数据进行检验。降水数据可以通过气象站观测、卫星遥感、大气模式模拟等获取，但地面气象站点观测的降水数据是最准确的 [ 9 ]，因此在本文中将该实测数据作为“真值”。TRMM 3B42V7数据从http://trmm.gsfc.nasa.gov/免费下载，空间分辨率为0.25˚ × 0.25˚，其时间分辨率为3 h，时间范围为2017年1月1日至2017年12月31日。TRMM 3B42V7数据为网格数据，站点实测值为单点数据，本文通过最临近插值法得到卫星降水数据在各个地面气象站点的降水估计值，同时将TRMM卫星3 h时间尺度的降水数据进行累加得到逐日降水数据，进一步与气象站日降水数据进行匹配。在年、月和日尺度上将TRMM 3B42V7降水产品与之进行线性回归，利用相关系数和相对误差来检验TRMM卫星降水数据的精度，从而明确其在川渝地区的适用性。 图1. 川渝地区概况以及71个气象站点的分布图
system: [SYSTEM]

user: 3. TRMM卫星降水量精确检验 3.1. 年降水量精度检验 如图2所示，年尺度上TRMM 3B42V7降水数据与实测降水数据之间具有很好的相关性，相关系数R达到0.66，但线性回归方程斜率仅为0.14。研究区71个站点实测的年降水量平均值为1068.24 (±306.27) mm，其对应的TRMM 3B42V7降水平均值为335.38 (±65.77) mm，年降水量相对误差为−68.60%。这些 图2. 2017年TRMM 3B42V7与地面观测站年降水量散点图 都表明TRMM 3B42V7降水数据明显小于实测降水数据，这可能与川渝地区地形复杂、气候系统多样有关。TRMM卫星数据不能够直接替代实际降水数据，但是二者之间的线性回归方程通过了0.01的显著性水平检验，说明通过回归方程可以估算出可信度比较高、与实际降水量接近的降水数据。 为进一步分析TRMM 3B42V7降水空间分布特征，利用克里金插值法将站点降水数据进行空间插值，分别绘制基于TRMM 3B42V7和实测降水数据的降水空间分布图。如图3所示，实测和TRMM卫星来源的降水量在空间分布上尽管数量大小存在一定差异，但其空间分布特征基本相同，从西部高原区向东部盆地区均呈现为增加趋势，且中部地区的四川盆地降水量明显大于四周山地降水量。 图3. 2017年川渝地区不同降水数据来源的年降水量空间分布图：(a) TRMM 3B42V7和(b) 实测降水量 3.2. 月降水量精度检验 如图4所示，月尺度上TRMM 3B42V7降水数据与实测降水数据相关性更好，相关系数R达到0.91，线性回归方程斜率增加到了0.31。研究区71个站点实测的月降水量平均值为97.88 (±88.23) mm，其对应的TRMM 3B42V7降水平均值为33.40 (±29.83) mm，月降水量相对误差仍然达到−65.88%，这表明月尺度上TRMM卫星估算的降水量同样明显小于实际降水量。因此，月尺度上TRMM卫星数据与实际数据相差很大，但由于二者之间的线性回归方程通过了0.01的显著性水平检验，且决定系数R 2 高达0.83，说明利用该方程根据TRMM 3B42V7卫星数据可估算出实际月降水量。 图4. 2017年TRMM 3B42V7与地面观测站月降水量散点图 3.3. 日降水量精度检验 如图5所示，日尺度上TRMM 3B42V7降水数据与实测降水数据散点分布明显比较分散，二者的相关系数与年尺度上的基本相同，为0.65。研究区71个站点实测的日降水量平均值为3.17 (±8.70) mm，其对应的TRMM 3B42V7降水平均值为1.08 (±3.14) mm，日降水量相对误差为−65.77%，与月降水量相对误差基本相同。说明TRMM 3B42V7降水数据在日尺度上同样显著低估了实际降水量。日尺度上TRMM卫星数据与实际数据之间的线性回归方程也通过了0.01的显著性水平检验，决定系数R 2 达到了0.42，利用回归方程也可以根据TRMM 3B42V7卫星数据来推算实际日降水量，但是显然估计值整体来看变异系数较大，原因是实测的日降水量和TRMM 3B42V7降水量的变异分别达到了274.4%和290.7%。 3.4. 最大降水量精度检验 由于本文研究是基于一年(2017年)的降水量数据，所以只对比了各站点月最大和日最大降水量与其相应的TRMM卫星降水量，结果如图6所示。实测的月和日最大降水量与其对应的TRMM 3B42V7卫星数据均具有很好的相关性，但月尺度上其相关系数达到了0.83，明显高于日尺度上的相关系数0.45。月和日尺度上的实测最大降水量与其对应的TRMM 3B42V7卫星数据的线性回归方程均通过了0.01的显 图5. 2017年TRMM 3B42V7与地面观测站日降水量散点图 图6. 2017年地面观测站实测降水量最大值与相应的TRMM 3B42V7卫星数据散点图：(a) 月最大降水量，(b) 日最大降水量 著性水平检验，但回归方程的斜率均明显小于理想值1，分别为0.20和0.12。此外，各站点实测月降水量最大值的平均值为255.79 (±98.33) mm，其对应的TRMM 3B42V7降水量平均值为80.54 (±23.12) mm，月降水量最大值相对误差为−68.51%；各站点实测日降水量最大值的平均值为82.57 (±63.85) mm，其对应的TRMM 3B42V7降水量平均值为27.24 (±12.60) mm，月降水量最大值相对误差为−67.02%。这进一步说明TRMM 3B42V7降水数据明显小于实测降水数据。
system: [SYSTEM]

user: 4. 单个站点数据精度检验 对年、月、日尺度降水数据的检验结果可知，TRMM3B42V7卫星降水数据明显小于站点实测数据，但是二者间相关性较好，但整体检验难免会掩盖少数站点数据与对应TRMM 3B42V7降水数据之间的差异，原因在于降水受地形、经纬度、海拔、坡度、坡向、大气环流、海陆位置等要素的综合影响，空间分布上具有不连续性和高度非线性特征。因此仅对降水数据进行总体精度评估是不够的，还需要对其各个站点进行验证。 以71个气象站点2017年实测日降水量和由此累计得到月降水量为自变量，以气象站点对应的TRMM 3B42V7卫星降水数据为因变量，作一元线性回归分析，分别得到TRMM3B42V7数据与气象站点之间在月尺度和日尺度上的相关系数R，如表1所示。大部分气象站点的月降水实测数据与TRMM 3B42V7的月降水数据之间具有较强的相关性，相关系数平均值高达0.915 (标准差0.074)，相关系数R > 0.9的台站数占总数的71.8%，相关系数 > 0.8的台站数占总数的90.1%。气象站点实测数据与TRMM月降水数据之间的相关系数在空间分布上表现出不均匀性(变异系数为8.0%)，其中峨眉山站的相关系数最小(仅为0.628)，这可能与峨眉山站点地形起伏明显、海拔落差大有关。相比之下，日尺度上实测降水量与TRMM 3B42V7降水量之间的相关系数平均值下降为0.622 (标准差0.121)，说明二者之间相关性较弱。日尺度上峨眉山站点的实测降水量与TRMM 3B42V7降水量之间的相关系数亦为最小，仅为0.191，说明二者之间没有相关性。日尺度上相关系数的变异系数高达19.4%，说明相关系数空间分布的不均匀性非常显著。 Table 1 Correlation coefficient of average monthly and daily precipitations between TRMM 3B42V7 and corresponding meteorological station data in Sichuan-Chongqing region in 201 站点R 月尺度 /R 日尺度 站点R 月尺度 /R 日尺度 站点R 月尺度 /R 日尺度 站点R 月尺度 /R 日尺度 站点R 月尺度 /R 日尺度 站点R 月尺度 /R 日尺度 石渠0.973/0.660 巴塘0.922/0.780 雷波0.964/0.653 开县0.938/0.674 大足0.793/0.488 南川0.876/0.570 若尔盖0.938/0.672 新龙0.977/0.561 宜宾0.943/0.730 云阳0.935/0.800 荣昌0.906/0.660 长寿0.861/0.592 德格0.900/0.615 理塘0.985/0.716 盐源0.940/0.567 巫溪0.978/0.756 永川0.940/0.602 涪陵0.809/0.441 甘孜0.964/0.532 雅安0.937/0.520 西昌0.970/0.632 奉节0.962/0.771 万盛0.848/0.621 丰都0.853/0.518 色达0.968/0.616 稻城0.988/0.615 攀枝花0.909/0.603 巫山0.975/0.751 铜梁0.923/0.558 武隆0.863/0.580 道孚0.782/0.254 康定0.953/0.747 会理0.993/0.788 遂宁0.919/0.506 北培0.928/0.541 黔江0.978/0.725 马尔康0.998/0.720 峨眉山0.628/0.191 广元0.928/0.603 潼南0.780/0.633 合川0.792/0.615 彭水0.896/0.602 红原0.961/0.734 乐山0.967/0.508 万源0.983/0.743 垫江0.757/0.661 渝北0.923/0.438 叙永0.683/0.508 小金0.917/0.621 木里0.891/0.679 阆中0.860/0.498 梁平0.976/0.777 璧山0.886/0.409 綦江0.930/0.545 松潘0.993/0.636 九龙0.975/0.688 巴中0.929/0.614 万州0.944/0.829 沙坪坝0.937/0.485 酉阳0.955/0.797 都江堰0.976/0.602 越西0.947/0.688 达川0.917/0.516 忠县0.879/0.715 江津0.802/0.674 秀山0.964/0.838 绵阳0.961/0.527 昭觉0.957/0.666 城口0.972/0.727 石柱0.950/0.610 巴南0.858/0.621 表1. 2017年川渝地区基于月平均和日平均降水量的TRMM 3B42V7数据与站点实测值的相关系数R
system: [SYSTEM]

user: 5. 结论 (1) 年尺度上，TRMM 3B42V7与实测降水数据之间具有很好的线性关系，但回归方程斜率仅为0.14，TRMM 3B42V7卫星数据与实测降水数据相比，年降水量相对误差达到了−68.60%，这表明TRMM 3B42V7数据明显小于实测降水数据。二者的空间分布特征基本相同：西部地区降水偏少，东部地区降水偏多，中部四川盆地降水明显大于四周山地降水。 (2) 月尺度上，TRMM 3B42V7与实测降水数据相关性最高，相关系数达到0.91，回归方程斜率增加到0.31，但月降水量相对误差仍高达−65.88%；日尺度上，TRMM 3B42V7与实测降水数据相关性降低到0.65，而且日降水量相对误差为−65.77%。这表明在月和日尺度上TRMM卫星估算的降水量明显小于实际降水量，对月和日实测最大降水量与对应的TRMM 3B42V7数据的分析也得到了相同的结论。 (3) 对71个气象站点的3B42V7降水数据进行逐一精度验证，发现在月尺度上，绝大大部分气象站点的降水实测数据与TRMM 3B42V7的相关性较好，相关系数平均值高达0.915，变异系数为8.0%，相关系数 > 0.8的台站数占总台站数的90.1%；而在日尺度上降水实测数据与TRMM 3B42V7降水量之间的相关系数平均值下降为0.622，变异系数增大到19.4%，说明相关系数空间分布显著不均匀。无论是月尺度还是日尺度上，峨眉山站的相关系数最小(分别为0.628和0.191)，这可能与峨眉山站点地形地势特征有关。
system: [SYSTEM]

user: 摘要 目的：对双螺钉和三螺钉固定治疗胫骨平台Schatzker I型骨折的生物力学进行有限元分析，为临床手术设计以及选择提供思路。方法：利用1位健康男性志愿者下肢CT数据构建胫骨平台Schatzker I型骨折的三维有限元模型。采用不同螺钉固定方式对骨折模型进行固定后，胫骨平台上方设置600 N的力量载荷同时对胫骨模型下方坚强固定来模拟成人膝关节的受力。比较模型在不同固定方式下的等效应力以及剪切力情况。结果：双螺钉固定系统的最大von Mises应力值为23.014 MPa，三螺钉固定系统的最大von Mises应力值为6.252 MPa。双螺钉固定系统与三螺钉固定系统中骨折块的von Mises应力最大值分别是0.736 MPa，2.158 MPa；胫骨的von Mises应力最大值为1.933 MPa 1.991 MPa。结论：使用三螺钉对骨折处进行固定优于双螺钉固定治疗胫骨平台Schatzker I型骨折，可降低固定物固定失败的风险，有利于愈合。
system: [SYSTEM]

user: 1. 背景 胫骨平台骨折属于临床上常见的下肢骨折之一，约占所有胫骨骨折的9.2%，常见于高能量损伤如高处坠落以及交通意外，尤以外侧平台骨折多见 [ 1 ] [ 2 ]。由于胫骨平台骨折通常累及平台关节面，因此往往需要通过手术来获得解剖复位并坚强固定来提供足够的生物力学支撑，从而避免严重并发症的发生 [ 3 ] [ 4 ] [ 5 ]。因此固定区域的生物力学强度至关重要，由于生物力学强度不足而引起的骨块间相互滑移是导致骨折愈合不良甚至不愈合的主要原因，经皮螺钉内固定是治疗胫骨平台骨折的手段之一，由于术后往往能获得较为满意的疗效，在临床上较为常用 [ 6 ] [ 7 ] [ 8 ]。双螺钉固定作为临床治疗胫骨外侧平台骨折的标准固定方式，近几年为了追求更加坚固的手术固定，防止骨折块滑移，三螺钉固定方式也被广大术者所使用 [ 6 ]。然而，关于采用不同螺钉固定方式治疗胫骨平台Schatzker I型骨折的数字生物力学优劣尚未得到证实。 三维有限元分析作为近些年出现的计算机数字分析方法，在医学领域中得到了广泛应用，主要集中于骨骼等的静力力学分析以及血流方面的流体力学分析 [ 9 ] [ 10 ] [ 11 ]。基于有限元分析的生物力学研究可以为临床提供了一个更深入的视角来研究生物力学分布及位移分布，为临床的术式设计以及器械改进等提供理论依据 [ 12 ]。因此在本研究中，笔者采用有限元方法对不同螺钉固定方式治疗Schatzker I型骨折的应力与位移进行比较与分析。本研究假设，三螺钉相较于双螺钉固定为Schatzker I型胫骨平台骨折提供更加坚强的固定，更适合临床的骨折治疗。
system: [SYSTEM]

user: 2. 方法 排除如骨性关节炎，骨肿瘤以及骨质疏松等骨性疾病后，选取1名30岁、身高175 cm、体重65 kg的健康男性志愿者，对其下肢进行CT扫描，获得健康成年男性胫骨CT数据。 2.1. 实验模型建立 将层厚1 mm的CT原始数据以DICOM格式导入到有限元建模软件Mimics 15.0 (比利时Materialize公司)中，利用软件中指定阈值及区域三维重建功能，获得健康成年男性胫骨三维有限元模型。将获得的胫骨模型导入Geomagic Studio软件(美国3D System公司)进行进一步的预处理以减少模型的重叠，利用裁剪功能对胫骨平台内侧约5 cm处进行无缝隙分割，制造出Schatzker I型胫骨平台骨折有限元模型。根据螺钉制造商提供的数据，利用Creo 3.0软件(美国Parametric技术公司)绘制直径6.5cm螺钉的有限元模型。如图1所示完成装搭之后，导入ANSYS软件(美国ANSYS技术公司)进行下一步的分析运算。在本研究中，单元类型选择为四节点四面体单元，以获得更准确的几何非线性分析结果。双螺钉模型与三螺钉模型的节点数分别为5222个、5763个；单元数分别为26，725个、29，561个。 图1. 模拟固定胫骨平台Schatzker I型骨折及两种螺钉固定方式的有限元模型。A.重建的完整胫骨平台；B.重建的胫骨平台Schatzker I型骨折有限元模型；C.双螺钉固定模型；D.三螺钉固定模型 2.2. 材料属性赋值 骨骼与螺钉模型的单元均定义为线性弹性材料属性，皮质骨，松质骨以及固定物的杨氏模量分别设置为17 Gpa，5 Gpa与110 Gpa；皮质骨，松质骨以及固定物的松柏比分别0.33，0.33与0.3 [ 13 ]。本研究中的所有材料，均被简化为均质且各向性质相同。 2.3. 边界条件及载荷 胫骨模型下方进行稳定固定，在上方处加载载荷，如图2所示进行设置。胫骨平台上方加载600 N方向垂直向下的力来模拟正常成人双足站立时膝关节处生理载荷，其中根据生理分布，载荷分为60%与40%均匀分布于胫骨内侧平台与外侧平台以最大程度模拟生理状态。为合理简化模型及减少不必要的运算量，螺钉与骨骼的固定简化为摩擦系数为0.3的摩擦接触 [ 9 ]。 图2. 胫骨有限元模型受力载荷分布及边界示意图 2.4. 实验观察指标 对各项参数完成设置之后，利用ANSYS软件对模型进行分析运算，结果中的等效应力以及各个方向的剪切力导出后作为实验观察指标作进一步的分析；对模型中每个部件的最大应力点的数值以及位置也进行收集整理，以分析模型载荷后的应力集中情况；对模型中每个部件的各方向最大剪切力进行整理，分析模型在和后的局部受力情况；并将模型整体力学以及各向剪切力分布绘制云图。 2.5. 统计学分析 利用SPSS17.0统计学软件，对实验结果进行分析整理，结果采用两独立样本t检验，P值小于0.05被认为差异具有统计学意义。
system: [SYSTEM]

user: 3. 结果 3.1. 两种模型von Mises应力的比较 双螺钉固定系统的螺钉最大应力为23.014 MPa，出现在靠近前方螺钉的骨折块间隙处位置，见图3；骨折块最大应力为0.736 MPa，出现在靠近前方螺钉孔上方与螺钉接触的位置，见图4；胫骨最大应力为1.933 MPa，出现在下方螺钉孔下方与螺钉接触的位置，见图5。三螺钉固定系统的螺钉最大应力为6.252 MPa，出现在下方螺钉的骨折块间隙处位置。骨折块最大应力为2.158 MPa，出现在靠近前方螺钉孔螺帽与骨块接触的位置；胫骨最大应力为1.991 MPa，出现在下方螺钉孔下方的位置。 3.2. 两种模型XY轴、YZ轴及XZ轴上剪切力的比较 两个模型XY轴、YZ轴以及XZ轴上的最大应力(见表1)。三螺钉系统的螺钉部分以及胫骨部分的各向剪切力均小于双螺钉固定系统的各向剪切力；而骨折块部分的XY轴与YZ轴方向剪切力则是三螺钉系统稍大于双螺钉固定系统。上述所有最大剪切力均集中在骨折块间隙周围位置。 Table 1 Maximum stress of different components in two fixation method model 双螺钉系统(MPa) 三螺钉系统(MPa) 螺钉系统 von Mises应力 23.014 6.252 XY轴剪切力 3.594 1.145 YZ轴剪切力 2.223 0.771 XZ轴剪切力 2.644 0.705 骨折块 von Mises应力 0.736 2.158 XY轴剪切力 0.053 0.562 YZ轴剪切力 0.256 0.508 XZ轴剪切力 0.903 0.124 胫骨部分 von Mises应力 1.933 1.991 XY轴剪切力 0.175 0.163 YZ轴剪切力 0.77 0.739 XZ轴剪切力 0.962 0.636 表1. 两种固定方式模型各部件的最大应力值 图3. 两种固定模型中螺钉系统的应力分布云图。双螺钉固定模型中螺钉系统的(A) von Mises力分布云图，(B) XY轴剪切力分布云图，(C) YZ轴剪切力分布云图，(D) XZ轴剪切力分布云图；三螺钉固定模型中螺钉系统的(E) von Mises力分布云图，(F) XY轴剪切力分布云图，(G) YZ轴剪切力分布云图，(H) XZ轴剪切力分布云图 图4. 两种固定模型中骨折块的应力分布云图。双螺钉固定模型中骨折块的(A) von Mises力分布云图，(B) XY轴剪切力分布云图，(C) YZ轴剪切力分布云图，(D) XZ轴剪切力分布云图；三螺钉固定模型中骨折块的(E) von Mises力分布云图，(F) XY轴剪切力分布云图，(G) YZ轴剪切力分布云图，(H) XZ轴剪切力分布云图 图5. 两种固定模型中胫骨的应力分布云图。双螺钉固定模型中胫骨的(A) von Mises力分布云图，(B) XY轴剪切力分布云图，(C) YZ轴剪切力分布云图，(D) XZ轴剪切力分布云图；三螺钉固定模型中胫骨的(E) von Mises力分布云图，(F) XY轴剪切力分布云图，(G) YZ轴剪切力分布云图，(H) XZ轴剪切力分布云图 3.3. 两种模型平均应力的比较 两种模型平均应力的比较中，三螺钉固定系统的螺钉平均应力明显小于双螺钉系统，三螺钉固定系统的总平均应力也小于双螺钉系统的总平均应力，而三螺钉固定系统的骨块平均应力与双螺钉系统的骨块平均应力差距没有前述显著，见图6。 图6. 两种固定模型中平均应力分布。A.螺钉系统；B.骨骼系统；C.总并平均应力
system: [SYSTEM]

user: 4. 讨论 胫骨平台骨折是一种常见的关节内骨折，占全身骨折发生率的1%~2%，对于老年人占到更高的全身骨折发生率 [ 14 ]。由于是关节内骨折，往往需要进行手术以达到满意的解剖复位，避免由于骨折所引起的延迟愈合甚至不愈合等并发症而影响膝关节的功能。使用拉力或皮质骨螺钉对骨折处进行固定是临床上治疗胫骨平台Schatzker I型骨折常见的治疗手段之一。然而，螺钉使用的数量增加往往与手术时间，术中出血量已经手术并发症风险密切相关，然而螺钉使用熟料不足则会导致固定不理想，螺钉松动以及骨折不愈合等并发症 [ 15 ] [ 16 ]。目前尚未有关于比较不同数量螺钉固定治疗胫骨平台Schatzker I型骨折的生物力学研究。有限元分析目前已广泛应用于医学领域，研究者可通过构建不同研究对象的三维有限元模型，在计算机中模拟受力环境，进而分析比较不同力学环境对于对象模型的力学影响 [ 17 ]。本研究中，笔者利用三维有限元分析方法，重建胫骨平台Schatzker I型骨折三维有限元模型，对采用双螺钉以及三螺钉固定的两种固定方式进行了比较。 从两个螺钉固定模型的应力分布分析结果可以看出，对于螺钉系统，主要应力集中于螺钉的骨折间隙处，表明螺钉提供了对抗骨块滑动以支撑骨块的作用。三螺钉系统中螺钉最高应力是双螺钉系统的三分之一，且最高应力出现在增加的下方螺钉中段，说明增加的下方螺钉分担了较多的应力，这种现象在其他位置骨折的固定物中同样出现 [ 18 ]。过度的应力集中不仅可能会引起固定物的断裂，另外也会由于应力的分布不均匀以及固定物的物理属性有别于骨骼从而导致应力遮挡效应的出现，进一步导致骨折处附近部位疼痛，骨折延迟愈合甚至不愈合 [ 19 ] [ 20 ]。然而目前尚未有针对研究引起骨折延迟愈合及不愈合的应力耐受极限，故该方面还有待进一步的研究。 同时，在剪切力与平均应力方面，三螺钉系统也体现出来了较为明显的优势。三螺钉固定方式模型在螺钉系统以及胫骨上的剪切力最大值均小于双螺钉固定模型，可能是由于三螺钉固定具有更大的受力面积，可以减缓分担剪切力从而使得最大剪切力的数值下降。较大的剪切力同样不利于骨骼的愈合以及固定物的稳定，从这一方面而言也说明三螺钉固定系统更安全稳定。这个研究结果与Moran等的研究结果相符，他们发现在实体标本的力学试验中，采用三螺钉固定的模型相较于双螺钉固定的模型具有更加强的支撑能力以及抗滑移能力 [ 6 ]。 当然，采用双螺钉进行固定也有一定的价值。一方面可以降低手术时间以及术中出血，同时可以避免置入下方螺钉时对周围神经的损伤，另外也可以减轻患者的经济负担。虽然双螺钉固定其受力系统中最大应力值相比较大，但对于恢复能力强且骨骼强度较高的青少年患者，若能适当推迟术后的下地负重时间，亦可作为备选的手术方式。
system: [SYSTEM]

user: 5. 结论 综上所述，对于胫骨平台Schatzker I型骨折的患者而言，使用双螺钉进行固定更加简单，创伤相比较小，感染风险更低，但是相比于三螺钉固定存在更大的固定物断裂可能。三螺钉固定系统中应力分布更加均匀，有利于患者骨折的愈合，尤其对于可能存在骨质疏松等基础疾病的老年患者而言，选用三螺钉甚至联合钢板固定更有利于患者。 本研究的局限：一方面，模型的材料均设置为了各向同性的材料属性，与实际情况中的骨骼性质稍稍有所不同；同时本研究没有将肌肉，筋膜等软组织对于骨骼力学环境的影响考虑在内；另外，本研究中只是单纯了研究了模型的静态载荷，而与实际环境中的动态力学环境有所出入。
system: [SYSTEM]

user: 摘要 为了有效利用沿海土地资源，我国广大沿海地区开展了吹填造陆。由此产生的人工吹填土作为建设工程地基的一部分，其性质也是近年来工程研究热点之一。本文分析了影响人工吹填土工程性质的四方面因素，对比分析了我国现行规范中采用标准贯入试验锤击数和静力触探测试数据的液化判别方法间的区别和适用性，得到了人工吹填土场地主要液化土层和大面积地基处理的分析方法。本文还根据曹妃甸滦南工业园地区工程实例中相关勘察、测试和地基处理数据计算得到的液化判别结果，分析了不同地基土体的液化性质，提供了针对人工吹填土地基处理的建议。
system: [SYSTEM]

user: 1. 引言 吹填造陆作为解决土地稀缺和有效利用沿海丰富滩涂资源的有效途径，在我院国沿海经济发达地区广泛开展。由此产生的吹填土由于沉积物和沉积环境，以及距离吹填口远近的不同，吹填土的成分、含砂量和工程性质也各不相同。 人工吹填土没有经过地基处理，大都为可液化土层。而地震液化是导致地基失稳，建筑物上部结构破坏的直接原因之一，因此，为了保证场地稳定，人工吹填土场地需要根据地基液化程度和设计条件，选择适合的地基处理方法消除地基液化。如何准确的评价人工吹填土的液化势和液化程度，为地基处理提供可靠数据，是岩土工程中重要的研究课题。 本文通过分析人工吹填土的成分构成和性质，讨论现行规范中主要地基液化判别方法的差异，采用原位测试成果细化地基土层液化程度，结合工程实例提出人工吹填土场地消除地基液化的建议。