下载文档原格式
传热过程中稳态与非稳态的传热性能研究热传导是物体内部或物体之间通过分子间碰撞来传递热量的过程。
在传热过程中，往往会涉及到稳态和非稳态传热性能的研究。
稳态传热是指系统内部的温度分布和热流密度保持不变的传热过程，而非稳态传热则是指系统内部的温度分布和热流密度随时间变化的传热过程。
本文将探讨传热过程中稳态与非稳态的传热性能研究。
一、传热过程中的稳态传热性能研究稳态传热是热传导的一种重要情况，其特点是温度分布和热流密度在系统内部保持不变。
在稳态传热过程中，可以通过一系列的实验和数学模型来研究传热性能。
1. 实验方法稳态传热性能的研究通常需要通过实验来获得数据。
实验中需要测量物体的温度和热流密度等参数，并利用传热方程来计算传热速率。
同时，还可以进行不同条件下的对比实验，以研究传热性能随温度差、材料特性等因素的变化。
2. 数学模型在稳态传热性能研究中，数学模型起到了重要作用。
通过分析热传导方程、输运方程以及边界条件等，可以建立物体的温度分布模型，从而计算热流密度和传热速率。
常见的数学方法包括有限元法、有限差分法等，这些方法可以通过计算机模拟和数值计算来得到稳态传热性能的结果。
二、传热过程中的非稳态传热性能研究非稳态传热是热传导的另一种情况，其特点是温度分布和热流密度随时间变化。
非稳态传热性能的研究对于理解瞬态过程、优化传热设备等具有重要意义。
1. 实验方法对于非稳态传热性能的研究，实验方法同样起到了关键作用。
实验中需要测量物体的温度随时间的变化，并利用实验得到的数据建立传热模型。
此外，还可以通过改变传热边界条件、加热方式等来研究传热性能的变化。
2. 数学模型非稳态传热性能的研究同样离不开数学模型的建立。
根据热传导方程和边界条件，可以建立非稳态传热模型。
此外，还可以利用转换方法将非稳态问题转化为稳态问题进行求解。
例如，可以采用拉普拉斯变换等数学方法来分析非稳态传热性能。
三、稳态与非稳态传热性能的比较与应用稳态传热性能和非稳态传热性能在实际应用中都具有重要意义。
热传导方程反问题热传导方程反问题是指在已知温度分布的情况下，通过测量边界上的温度来确定材料的热传导系数。
这个问题可以用数学模型来描述，即热传导方程。
热传导方程是描述物质内部温度分布随时间和空间变化的偏微分方程。
它可以用以下形式表示：∂u/∂t = α∇^2u其中，u表示温度分布，t表示时间，α表示热传导系数，∇^2表示拉普拉斯算子。
在反问题中，我们已知边界上的温度分布和时间变化情况，需要求解未知的热传导系数α。
为了解决这个问题，可以采用逆问题方法。
逆问题方法是一种数学处理方法，在已知输出数据和输入模型之间寻找最优解。
在热传导方程反问题中，逆问题方法可以通过以下步骤进行：1. 建立正问题模型：根据已知条件建立热传导方程，并求解出温度分布。
2. 确定目标函数：目标函数是一个衡量模型输出与实际观测值之间差异的指标。
在本例中，目标函数可以定义为测量值与模拟值之间的平均误差。
3. 选择逆问题方法：逆问题方法有很多种，包括正则化方法、贝叶斯方法、遗传算法等。
在本例中，可以采用最小二乘法。
4. 求解逆问题：根据正问题模型和目标函数，使用最小二乘法求解未知的热传导系数α。
热传导方程反问题的求解过程中需要注意以下几点：1. 数据收集：在进行反问题求解前需要收集足够的数据，包括边界上的温度分布和时间变化情况。
2. 正确建立模型：建立正问题模型时需要考虑材料的物理特性和实际情况，并进行合理简化。
3. 选择合适的逆问题方法：不同的逆问题方法适用于不同类型的反问题，需要根据具体情况选择合适的方法。
4. 对结果进行验证：求解出热传导系数后需要对结果进行验证，比较模拟值与实际观测值之间的差异，以评估求解结果的可靠性和精度。
总之，热传导方程反问题是一种重要的数学处理方法，在工程领域中具有广泛应用。
通过正确建立模型、选择合适的逆问题方法和对结果进行验证，可以求解出未知的热传导系数，为工程设计和优化提供有力支持。
导热系数单位的物理意义-概述说明以及解释1.引言1.1 概述导热系数是一个重要的物理量，用来描述物质导热性能的好坏。
在热传导过程中，物质内部的热能会通过导热介质的传递而传递到外部环境。
导热系数是衡量这种热传导效果的一个指标，它能够告诉我们物质导热的快慢程度。
导热系数的单位是热导率除以厚度，通常用瓦特/米·开尔文（W/m·K）来表示。
这个单位的物理意义可以解释为：单位时间内单位面积上热流通过的热量，对于单位温度差的变化。
这说明导热系数与物质的热导率、厚度以及温度差有关。
导热系数的物理意义可以从热传导的角度来理解。
在一个温度差存在的物质中，高温区域的分子具有较高的热能，而低温区域的分子则具有较低的热能。
通过导热介质，高温区域的热能会传递给低温区域的分子，使得温度逐渐均匀。
导热系数就是描述这种传递过程的特性之一，它越大则说明物质的导热性能越好，能够更有效地将热量传递出去。
在实际应用中，导热系数的物理意义也是非常重要的。
例如，在建筑材料的选择和设计中，需要考虑到材料的导热性能，以确定其在不同温度环境下的导热效果。
此外，在电子器件的散热设计中，导热系数的大小也直接关系到器件的温度分布和热量的传递效率。
因此，理解导热系数的物理意义对于研究和应用领域具有重要的指导意义。
综上所述，导热系数是一个描述物质导热性能的重要物理量，其单位瓦特/米·开尔文表示了单位时间内单位面积上热流通过的热量对于单位温度差的变化。
通过理解导热系数的物理意义，我们能够更好地理解热传导的过程，并在应用中更好地利用和设计物质的导热性能。
1.2 文章结构本篇长文旨在探讨导热系数单位的物理意义。
文章主要分为引言、正文和结论三个部分。
引言部分将概述导热系数的背景和重要性，介绍读者对于导热系数有一个初步的了解。
在引言的结尾，将明确文章的目的，即探讨导热系数单位的物理意义。
正文部分将分为三个子部分。
首先，我们将对导热系数进行定义，明确导热系数的意义和计算方法。
《Nb3Sn超导材料力-电磁-热多场耦合本构关系的多尺度计算方法》篇一一、引言随着科技的不断进步，超导材料在能源、医疗、交通等众多领域的应用日益广泛。
其中，Nb3Sn超导材料因其卓越的超导性能和良好的机械性能，在核磁共振成像仪、超导磁体以及超导电力电缆等设备中扮演着重要角色。
本文旨在探讨Nb3Sn超导材料在力-电磁-热多场耦合下的本构关系，并介绍一种多尺度计算方法。
二、Nb3Sn超导材料的特性Nb3Sn超导材料具有较高的临界电流密度和临界磁场强度，使其在超导应用中具有显著优势。
然而，由于力-电磁-热多场耦合的存在，其本构关系变得复杂。
在力场作用下，Nb3Sn超导材料的机械性能会发生变化；在电磁场作用下，其超导性能会受到影响；而在热场作用下，其相变行为和热稳定性将发生变化。
因此，研究其多场耦合下的本构关系具有重要意义。
三、多场耦合本构关系在力-电磁-热多场耦合下，Nb3Sn超导材料的本构关系涉及多个物理过程的相互影响和作用。
为准确描述这一过程，需要建立多尺度计算方法。
该方法包括微观尺度的电子结构和原子结构分析，中观尺度的材料力学性能和电磁性能分析，以及宏观尺度的多场耦合效应分析。
通过综合分析这些尺度的数据，可以得出Nb3Sn超导材料在多场耦合下的本构关系。
四、多尺度计算方法1. 微观尺度计算：利用第一性原理计算和分子动力学模拟等方法，研究Nb3Sn超导材料的电子结构和原子结构。
通过分析电子的能级结构、态密度以及原子间的相互作用等，揭示其在力、电、磁、热等作用下的微观行为。
2. 中观尺度计算：基于微观尺度的分析结果，结合连续介质力学和电磁学理论，建立Nb3Sn超导材料的中观模型。
通过分析材料的力学性能、电磁性能等，得出其在中观尺度下的本构关系。
3. 宏观尺度计算：在考虑多场耦合效应的基础上，建立Nb3Sn超导材料的宏观模型。
通过综合分析力、电、磁、热等多场作用下的宏观行为，得出其在宏观尺度下的本构关系。
张洪武等：非经典热传导问题多尺度分析方法研究非经典热传导问题多尺度分析方法研究古张洪武张盛郭旭大连理工大学工业装备结构分析国家重点实验室，工程力学系．大连１１６０２４
摘要根据时空间尺度的高阶均匀化理论．建立分析尉期性结构中非傅立叶热传导问题的时间一空训多尺度分析方法，通过引入放大空间尺度和缩小时间尺度，研究了由空间非均匀性引起的非傅立叶热传导的波动效应和非局部效应，得到具有非局部效应的四阶微分方程，对高阶非局部热传导方程进行修正．使问题的求解避免了对有限元离散的Ｃ—ｌ连续性要求。给出三种不同材料参数情况的计算结果．验证方法正确性的同时，对存在的问题进行了讨沦。关键词非傅立叶热传导；多尺度方法；均匀化方法：非局部模型
１引言传统的傅立叶（Ｆｏｕｒｉｅｒ）导热定律是导热现象规律性的经验总结，它是建立在大量常规传热实验的基础上的。傅立叶导热定律不涉及热传导时间项，定律本身隐含了热传播速度为无限大的假设。但对极端热传导条件下的非稳态传热过程，如激光表面热处理、脉冲干燥及微时间或微空间尺度条件下的传热问题等，热传播速度的有限性却必须考虑，此时会出现一些不同于常规传热过程的物理现象．这种热传导效应称为非傅立叶（ｎｏｎ．Ｆｏｕｒｉｅｒ）热传导效应。传统傅立叶导热定律的本构方程描述了热流量和温度空间梯度分布之间的关系，其数学表达式为抛物线型偏微分方程。而非傅立叶导热定律还考虑了热波的时间迟滞，其数学表达式为双曲线型偏微分方程。热传导问题多尺度分析方法的研究具有极大的学术探讨价值和广泛的工程应用领域“’…。本文主要目的是根据非傅立叶导热定律本构方程，研究不同材料组成的多相结构中热传导问题的时间．空间多尺度分析方法ｈ５１ｏ在非均匀介质的分界面上存在的反射和折射作用影响了脉冲激励的传播，在宏观上出现了勃、散，衰减等现象。为了解决这一勃、散效应，本文采用了多时间一空间尺度的高阶均匀化理论对问题进行分析。通过引入放大空间尺度和缩小时间尺度，从数学上获得不同阶次的时间一空间问题的均匀化方程，对这些具有不同阶次的均匀化方程进行合并整理，最终得到用于结构宏观多尺度分析的高阶均匀化方程。
２时间一空间多尺度渐进分析的基本方法如图ｌ所示，假设宏观的特征尺寸￡远大于非均匀性尺寸，。在空间上引入两个尺度：一个是宏观或整体空间尺度ｘ，另一个是微观或局部放大空间尺度ｙ，且Ｙ＝Ｉ／ｓ，其中ｓ＜＜１。在时间上引入
一个一般时间尺度，即ｔｏ＝ｒ，同时还引入一个缩小时间尺度：ｔｌ＝８２ｆ，以进行时间域的多尺度分析。因为瞬态温度场≯与ｘ、Ｙ、ｔｏ、ｔ１相关，对≯采取近似多尺度渐进展开，得ｄ（ｘ，Ｙ，ｆ）＝如（ｘ，Ｙ，ｔｏ，＾）＋印【（ｘ，Ｙ，ｔｏ，ｒ１）＋占２≯２（ｘ，ｙ，ｆｏ，‘１）＋－…・・（１）
Ｉ鼻＾ｌｅ：＾黝材料＝
图ｌ一维杆和单元结构所研究的结构右侧施加热源加），其余表面绝热，其特征长度为ｆ（在ｘ尺度上）和盎（在ｙ尺度上）
＋幽家自然科学基金５０１７８０１６、杰出青年科学基金１０２２５２１２资助项目堡兰查兰！！！丝苎垫竺兰望兰兰墨垦坌堑三堡翌墨——————————————型且盎：，／ｓ。并认为热流量的温度空间梯度项ｑ’＝々（ｘ）屯在壶上有周期性。瞬杰非傅立叶热传导微分方程＾（√ｓ）扛∽ｓ礁。＋九｝一诲∽ｓ谚，Ｉ，＝０（２）
宏观边界上的边界条件，舳）－ｏ，织（“）＝音鼯
≯－，ｏ）＝，０），畦，ｘ，ｏ）＝ｇｂ）其中：≯（，，：／。，ｆ）表示瞬态温度场；＾Ｇ／ｓ）表示单位体积物质的比热：ｒｂ肛）表示驰豫时间：女ｂ肛）表示热传导系数（々∽ｓ）、ｒ（∥ｓ）和ｚ∽ｓ）具有局部周期性）；Ｆ表示横截面积；Ｑ０）表示热源密度；（１，
和（１，分别表示空间和时间的导数。为建立均匀化模型．定义平均算子（．）：１卅－１ｆ．ｄＹ。利用链式法则，空间和时间导数表示为
畦：＝破，＋ｓ一‘九．≯＝谚ｂ＋ｓ２哆＾，≯＝谚协＋２ｅ２谚＾ｂ＋ｓ４声ｖ。热流量温度空间梯度Ｉ贝为ｑ’＝ｔ（九十ｓ。幻）（３）式（２）热传导微分方程变为州“杀＋ｓ２毒ｘ誉“２静＋（翌ｃ３ｔｏｗ２丝Ｏｔ。）Ｊ＿瞰＋ｓ－ｌｑｊＪ
Ａ・【（哦０ｆ０＋虬）＋ｓ２（２嘛＾＋移＾）＋ｓ４嗍＾】－眈＋ｓ。ｑ：ｙ】（４）把≯（ｘ，Ｙ，ｆ）的渐进展开式（１）代入式（４）左边，得（４）左边２Ａ・【（砒，。０ｆｏ＋九山）＋ｄ确，，。～＋Ａ，‰）＋ｓ‘（晚，。ｍ＋如，ｆ０）＋……＋
ｏｅ２（２ｒ耐ｏ＾ｆＩ＋如，ｆ．）＋ｓ３（２硝Ｊ以＋氟＾）＋占４（２确山ｆ．＋晚．ｆ．）＋’・…’＋４确＾ｆＩ＋Ｆ５吮ｖ１＋ｓ６确，¨＋…‘＿】把≯（ｘ，Ｙ，ｔ）的渐进展开式（Ｉ）代入式（３）ｔ得ｑ’＝七｛【九一＋ｓ一１九∥】＋【渐．，＋卉∥】＋【占２如，，＋锄∥１＋‘‘‘‘。｝
＝七ｐ－１≯¨＋【九。＋蛾，，１＋４＃ｈＦ＋如，，】＋¨…）＝Ｆ一１ｑ：ｌ＋蕊＋叼：＋ｃＺｑ；＋……（５）其中ｑ！ｌ＝ｋ‰∥ｑ：＝ｋ【丸．，＋≯¨，ｙ１，ｓ＝ｏ，ｌ，２……
（６）
把ｑ’的渐进展开（５）代入式（４）右边，得（４）右边＝ｓ＿ｑ！ｔ．ｊ＋酶，，＋唧ｉｄ＋……＋声－２口！ｔ，，＋ｓ叫ｑ；，ｙ＋ｑ０＋……
＝ｓ一２ｑ！ｌ，，＋ｓ一１（ｑ；∥＋ｇ！１．Ｊ）＋（ｑｉ，，＋口；，，）＋‘・・・・＿根据式（４）两边比较（相同的∥项对应系数相等），得到以下各阶传热方程“ｓ“）：目！１．。＝０
０（ｓ一‘）．ｑ如＋ｇ：ｌＪ＝ｏＯ（ｅｏＪ．２－（砒Ｊｍ＋氏．‰）＝叱＋叱
ｏ（ｅ。）．２（硝ｍ＋Ａ，ｆ０）＝吐，＋ｇｉ．，ｏ（ｓ２）：＾（砸２如ｂ＋≯２，～＋２哦山＾＋‰）２ｑ；，ｙ＋ｑ！．，
（７）（８）（９）（１０）（Ｌ１１４８２张洪武等：非经典热传导问题多尺度分析方法研究３对不同阶问题的解
１）ｏ（ｓ一２）问题由传热方程（７），ｇ：¨＝０，两边同乘以九，在单胞域ｒ上分步积分，得ｆ４）ｏｑ＇－ｔ，ｙｄＹ＝Ｉ≯ｏｑ＇－－ｒｉｄｓ一卜（粕，，）２ｄｙ＝ｏｒｄｖｒ由于以Ｙ为周期的单胞边界ａｒ上边界积分项为零，那么式（１２）第一项Ｉ≯ｏｑ＇－。ｒｉｄｓ‘且ｌｔ（九）２ｄ，，＝ｏ，推出ｒ九，，＝０ｊ‰＝中ｏ（ｘ，ｔｏ，ｔ１）
ｑ＇ｔ＝ｔ如。＝０
（ｔ２）０。因为ｋ是止数
ｎ３）ｆ１４１２１ｏ（ｅ１）ｆ司题由传热方程（８），孤，＋ｑ：Ｌ，。＝０，由式（１４）ｑ！Ｉ＝０得ｑ：Ｌ．。＝０，并把式（６）和（１３）代／２（８），得吐ｙ＋ｑ＇－Ｉ．，＝畹ｙ＝【ｋ（Ｏｏ，。＋Ａ，ｙ）】．，＝０（１５）
可以观察到Ａ是关于中¨的线性函数，那么可假设办（ｘ，Ｙ，ｔｏ，ｔｔ）＝Ｏｔ（ｘ，ｔｏ，ｔｔ）＋Ａ（ｙ）（ｂｏ．，（１６）
把式（１６）代入（１５）、（６）得［ｔ（１＋“，ｙ）】．，＝０（１７）
“＝ｍＯ，ｘｋ０＋Ａ，ｙ）（１８）考虑图中的单胞结构。单胞域是由Ａ１和Ａ２微结构子域组成，如Ａ‘”＝圳０＜Ｙ＜蕊】；∥＝叫蕊‘Ｙ＜囱
其中：０≤口≤ｌ是单胞的体积分率方程（１７）和０８）在整个单胞可表示为ｋｊ（１＋Ａｊ，ｙ）＝口Ｊ；甜”＝ｏｏ，，ｋｊ（１＋Ａｊ∥）（Ｊ＝ｔ，２）（１９）
其中：ａ；为常数
由解在单胞上的周期性、界面处的连续性以及均匀化条件可解出Ａｊ，（，）圳２丽（１－ａ）（ｋ２－ｋ０旷，争圳＝百轰瓮护坐竽】
（２０）ｋ＝（州坞））２瓦鬟意（２１）
３）０（ｓｏ）问题由传热方程（９），２－（哦山ｆｏ＋ｑ）Ｏ，ｔｏ）＝乱ｙ＋酰，＝ｔ【哦，，＋如，ｙ】，，＋ｔ【南，，＋Ａ，，】，，根据周期函数的导数在一个周期域上积分值为零和平均算子定义可知（ｑｉ．，）＝０。把式（１３）代入传热方程（９），并对传热方程两边取平均算子。为了进一步推导，取驰豫时间ｆ＝“＾一（面ｏ．ｒ幽＋ｍｏ山）一（ｑｏ．，）＝０（２２）
其中厶＝《五）＝瑾＾＋（１一口）如，“＝（ｆ）＝口１＋（１一ａ）ｒ２把式（１８）代２ｋ（２２）得到宏观传热方程Ｚ一（种。山ｆ０＋Ｏｏ．ｆ０）一ｋＨ＜ｂｏ，。＝０（２３）
ｍ式（２３）得
ｔｏｏ，ｆ。，。＋ｍｏ，ｔ。＝石１ｋ月巾ｏ，肼张洪武等：非经典热传导问题多尺度分析方法研究代入传热方程（９）ｑｌ。＝ｐ（ｙ）ｋ…Ｏｏ—ｑ；，，ｆｌ（Ｙ）＝２（ｙ）／２．
ｆ２４１
把式（１６）代／ｋ（２４）左边，得左边＝ｑｉ．ｙ＝【七（Ａ．，＋如，ｙ）】，ｙ＝｛七［≯２，ｒ＋ｍｌ，，＋４００．“】｝，ｙ把式（１６）４ｔＡ．（２４）右边，得右边２ｆｌ（ｙ）ｋ。ｍｏ，。一（中０，ｘｔ（１＋Ａ．，））．，＝西ｏ．。［ｆｌ（ｙ）ｋ。一ｋ（１＋Ａ．，）】方程变为｛七【疵．，＋中Ｉ．，＋４中ｏ．“】｝．ｐ＝ｏｏ．“【卢（ｙ）一１）】七。（２５）可以观察到也是关于西１．，和ｍｏ。的线性函数，那么可假设≯２（ｘ，ｙ，ｔ０，ｔ１）＝ｑ）２（ｘ，ｔＯ，ｔｉ）＋４（ｙ）ｏＩ，＋Ｂ（ｙ）Ｏｏ“（２６）把（２６）代入（２５）、（６）可知｛ｋ［Ａ＋只，】｝．，＝【ｆｌ（Ｙ）一Ｉｌｋ。（２７）
“＝ｋ（Ｌ＋Ａ．ｙ）ｍｌ一＋ｋ（Ａ＋量Ｐ）ｍｏ．“（２８）同理可解口（，），得州加ｔ鼍ｃ扣一赛辫等，ｙ２＋ｔ一鲁ｃ扣＋篆鬻警岩加
｛一些毪害盟（扣。ａ（１－ⅥａＸ旷ｌ－…２ａ）．ｆｉ２唰（ｋ２－ｋ１）｝２ｋｋ（２９ａ）、ｌ，，、＾。’Ｌ川ｌＩ一口ＩＥ．＋触，ｌ
、’
圳＝ｔ象ｃ砉叫＋而ｏｔ（而ｋ２－ｋ面ｔ）旷斗与警ｃ砉叫＋桨舞笔者¨
卜畦堑型１蒙２ｋｋ生幽监（和＋坐１业２［（１－ａ丝）ｋ警≯｝（２９ｎ）、ｌ２、九
７
ｌ＋础２】’
……
由（２９）可以得到（２，４）＝０（３０）
（ｔ（Ａ＋彤））２０（３１）
４）ｏ（ｅ‘）问题由传热方程（１０），Ａ・（硝州。＋Ａ，～）＝ｑ；．ｙ＋ｑｉ，＝【女（戎．，＋如．，）】，，＋【々（Ａ，，＋如，，）】，，，根据周期函数的导数在一个周期域上积分值为零和平均算子定义可知（吐，）＝０。把式（１６）代入传热方程（１０），并对传
热方程两边取平均算子兄一（ｄｌｌｔＪ０ｆ０＋中ＩＪｏ）＋（＾彳Ｘｒ（ｏｏ，，），，讹＋（ｍｏ．，）ｔｆｏｌ＝（ｇ：，，）（３２）
把式（２８）代ｆｌ，（３２）得＾（椰’¨山＋ｍ‰）＋（＾一Ｘｆ（巾吣）．ｒ山＋（ｏ吣）．ｂｌ＝（七（１＋Ａ，ｙ））ｏＩ．“＋（七（Ａ＋置，））巾ｏ，埘（３３）
把式（２１）、（３０）和（３１）代入（３３）得到宏观传热方程＾。（ａＩｈ，‘＾＋中ｌ，～）一‘ｍｌ，。＝０（３４）把式（１６）、（２６）和（２３）、（３４）代入传热方程（１０）得（ｔ（丸，，＋ｔ１）２，，＋爿中ｌ，。＋脚Ｉｏ．一）｝，ｙ＝ｔ。（卢一ｌ冲ｔ，。＋［ｔ。肛－ｋ（Ａ＋钆）归ｏ．一（３５）可以观察到九是中２∥母１．。和ｍｏ．一的线性函数，那么可假设鸡（ｘ，Ｙ，ｔＯ，ｔｔ）＝０３（ｘ，ｔｏ，ｔｉ）＋爿（ｙ）ｍ２，，＋Ｂ（ｙ）ｍ１．“＋Ｃ（ｙ）ｍｏ．脚（３６）
把式（３６）代入（３５）、（６）得【ｋ（Ｂ＋Ｃ，）ｂ＝ｔ。肛一ｋ（Ａ＋只ｙ）（３７）
醴＝七。西２．，十七（彳＋量ｙ）巾Ｉ．盯＋七（口＋ｃｙ）中ｏ．一（３８）同理可解得ｃ（ｙ），由ｃ（，）可以得到
