ANSYS活塞温度场分析全过程

合集下载

基于ANSYS的焊接温度场和应力的数值模拟研究

基于ANSYS的焊接温度场和应力的数值模拟研究

基于ANSYS的焊接温度场和应力的数值模拟研究一、本文概述随着现代工业技术的飞速发展,焊接作为一种重要的连接工艺,在航空、汽车、船舶、石油化工等领域的应用日益广泛。

然而,焊接过程中产生的温度场和应力场对焊接结构的性能有着至关重要的影响。

为了深入理解焊接过程中的热-力行为,预测焊接结构的变形和残余应力,进而优化焊接工艺参数和提高产品质量,本文旨在利用ANSYS有限元分析软件,对焊接过程中的温度场和应力场进行数值模拟研究。

本文首先简要介绍了焊接数值模拟的意义和现状,包括焊接数值模拟的重要性、国内外研究现状和存在的问题等。

随后,详细阐述了ANSYS 软件在焊接数值模拟中的应用,包括其基本原理、分析流程、模型建立、参数设置等方面。

在此基础上,本文以某典型焊接结构为例,详细阐述了焊接温度场和应力场的数值模拟过程,包括模型的建立、边界条件的设定、求解参数的选择、结果的后处理等。

对模拟结果进行了详细的分析和讨论,验证了数值模拟方法的准确性和可靠性,为实际工程应用提供了有益的参考。

本文的研究不仅有助于深入理解焊接过程中的热-力行为,为优化焊接工艺参数和提高产品质量提供理论支持,同时也为ANSYS软件在焊接数值模拟领域的应用推广和进一步发展奠定了基础。

二、焊接理论基础焊接是一种通过加热、加压或两者并用,使两块或多块金属在原子层面结合形成永久性连接的工艺过程。

焊接过程涉及复杂的物理和化学变化,包括金属的熔化、凝固、相变以及应力和变形的产生等。

因此,深入了解焊接过程的理论基础对于准确模拟焊接过程中的温度场和应力分布至关重要。

焊接过程中,热源将能量传递给工件,导致工件局部快速升温并熔化。

熔池形成后,随着热源的移动,熔池中的液态金属逐渐凝固形成焊缝。

焊接热源的类型和移动速度、工件的材质和厚度等因素都会影响焊接过程的温度场分布。

为了准确模拟这一过程,需要了解各种热源模型(如移动热源模型、体积热源模型等)及其适用范围,并选择合适的模型进行数值模拟。

ANSYS计算大体积混凝土温度场的关键技术

ANSYS计算大体积混凝土温度场的关键技术

ANSYS计算大体积混凝土温度场的关键技术在大型混凝土结构建设的过程中,温度场分析对保障混凝土结构的安全性至关重要。

ANSYS作为工程领域中常用的数值模拟工具,能够对混凝土温度场进行准确的计算,为混凝土结构的设计和施工提供科学依据。

但是,对于大体积混凝土的温度场计算,存在一些关键技术需要考虑,下面将进行详细介绍。

1. 混凝土的物理性质混凝土的物理性质是温度场模拟中的关键因素之一。

混凝土在浇筑后的初凝期、成型期、硬化期、老化期等各个阶段的物理性质都存在巨大的变化。

因此,在进行混凝土的温度场计算前,需要准确地测量混凝土在不同时间点的物理性质,如热导率、比热容、密度等。

2. 热源的模拟混凝土的温度场计算需要考虑混凝土内部的各种热源对温度场的影响。

建筑中的热源包括太阳辐射、室内外温度差、人体热辐射等,需要对这些热源进行准确的模拟。

3. 初始条件的设置混凝土温度场计算的初始条件设置直接影响计算结果的准确性。

混凝土在浇筑后的初始温度值、初始变形状态的设置等都需要进行准确、合理的处理。

4. 传热模型的选择对于大体积混凝土的温度场计算,需要选择合适的传热模型。

传热模型可以根据混凝土的物理性质和热源的模拟情况,选择适用于不同情况下的传热模型,如动态传热模型、静态传热模型等。

5. 计算方法的选择针对大体积混凝土温度场的计算,需要选择合适的计算方法。

常用的方法有有限元法、有限差分法等,需要根据混凝土内部温度场、变形场等的变化情况,选择合适的计算方法。

6. 数值模拟使用ANSYS进行混凝土温度场计算,需要进行数值模拟。

数值模拟是对真实物理系统的数学模拟,通过建立数学模型,利用计算机运算获得物理系统的各种行为特性,如温度场、应力场、变形场等。

7. 计算结果的验证在进行混凝土温度场计算后,需要对计算结果进行验证。

验证结果通常采用实验测试的方式进行验证,如温度测试、原位应力测试、变形测试等。

,针对大体积混凝土温度场计算,需要考虑混凝土的物理性质、热源模拟、初始条件设置、传热模型选择、计算方法选择、数值模拟和计算结果验证等方面的关键技术,以保证计算结果的准确性和可靠性。

用ANSYS来计算和分析CEFR堆容器温度场和热应力

用ANSYS来计算和分析CEFR堆容器温度场和热应力

用ANSYS来计算和分析CEFR堆容器温度场和热应力要介绍应用ANSYS程序计算CEFR堆容器的温度场、热应力、地震应力等等,完成CEFR堆容器的整个计算和评定,凸显了ANSYS作为大型通用的有限元程序,对我们核工程实际工作的强大支持,尤其是ANSYS在CEFR堆容器温度场和热应力方面的计算,在一定程度上已经超越了专门热工软件所能完成的范围。

中国实验快堆(CEFR)是国家重点工程,是国家863计划能源领域最重要的项目。

在CEFR以前的初步设计阶段,我们虽然购买了ANSYS 程序,但应用方面没有体现出ANSYS作为国际上大型的通用的有限元分析软件的特点,但是在CEFR施工设计阶段,我们运用ANSYS计算了CEFR堆容器的温度场、热应力、地震应力等等,完成了CEFR堆容器的整个计算和评定,凸显了ANSYS作为大型通用的有限元程序,对我们核工程实际工作的强大支持,尤其是ANSYS在CEFR堆容器温度场和热应力方面的计算,在一定程度上已经超越了专门热工软件能完成的范围。

堆容器是CEFR堆的核心组成部件,按照规范要求,对堆容器的应力强度评定需要计算温度场和热应力,本文利用ANSYS程序的热分析功能,建立堆容器各个部分的分析模型,采用变参数输入、等效热导率等方法计算了堆容器在不同工况下的温度场和热应力,为堆容器的应力强度评定提供各种热应力结果。

1 结构和材料CEFR堆容器由主容器和保护容器双层壳以及堆容器支承组成,外层为保护容器,保护容器外部覆盖有相应厚度的保温层,内层为主容器。

在材料部分,遇到的困难是氩气层如何处理的问题。

按规范要求,氩气层应考虑传导、对流、辐射三种传热方式的混合换热,我们运用ANSYS程序的温度场计算功能,仅由热工提供温度边界条件(环境温度、对流换热系数等),完成了堆容器轴对称模型的温度场计算,这项工作替代了热工的计算工作。

同时,为了解决占用内存、计算时间等问题,我们在计算中采用了等效热导率的方法:将氩气层作为一种实体材料,与堆容器同时建模分网。

基于Ansys的永磁同步电机温度场分析

基于Ansys的永磁同步电机温度场分析

近些年来 ,由于永磁同步电机在工农业 的生产生活中应用越来越广泛 ,关于永磁同步电机的各类 问 题也越来越引起人们的重视 .由于电机的应用逐渐广泛 ,电机的负载逐渐增大 ,电机运行时的损耗也逐渐 增加 ,导致 电机 内部温升也 同时增大.而过高的温升容易使得电机部件出现过早老化 的现象 ,所以电机 内 部温升是电机设计 中必须考虑 的参数 .因此 ,电机的发热问题也引起了人们 的重视.但是由于电机内部结 构 复杂 、磁场分 布不 均匀 等 因素 ,使 得人 们对 于 电机 内部 温度 场分 布 的认识 大 多停 留于公 式计 算 出的 电 机 温升 得平 均值 ,而 没有 具体 到 电机 的温 度分 布 规律 .但 利用 Ansys有 限元 分析 软件 ,可 以模 拟 出 电机 内 部温度场的分布情况 ,从而解决这一问题.本文通过Ansys软件中的热分析模块 ,对实际电机温度场进行 二维分析 ,得 出温度分布云图,模拟 电机在运行过程 中温度场分布 隋况.
摘要 :利用Ansys ̄-f限元分析软件对永磁 同步电机的平面二维模型进行温度场分析 ,将 分析得到的结果与理
论计算结果进行 比较 ,寻找 出一种有效的计算 电机温度场的方法 .
关键词 :永磁 同步 电机 ;Ansys软件 ;瞬态热分析 ;温度场
中图分类号 :TM 341
文献标识码 :A
文章 编号 :2095—2481(2013)02.0138.03
第 25卷第 2期 2013年 5月
宁德师范学院学报(自然科 学版) Journal of Ningde Normal University(Natural Science)
Vol_25No.2 M ay2013
基于Ansys的永磁 同步电机温度场分析

《2024年基于ANSYS的焊接温度场和应力的数值模拟研究》范文

《2024年基于ANSYS的焊接温度场和应力的数值模拟研究》范文

《基于ANSYS的焊接温度场和应力的数值模拟研究》篇一一、引言随着科技的发展,焊接技术作为制造行业中的关键工艺之一,其质量和效率直接关系到产品的性能和寿命。

因此,对焊接过程中的温度场和应力分布进行精确的数值模拟显得尤为重要。

ANSYS作为一种功能强大的工程仿真软件,被广泛应用于焊接过程的数值模拟。

本文将基于ANSYS,对焊接温度场和应力进行数值模拟研究,以期为实际生产提供理论依据。

二、焊接温度场的数值模拟1. 模型建立在ANSYS中建立焊接过程的有限元模型,包括焊件、焊缝、热源等部分。

其中,焊件采用实体单元进行建模,焊缝则通过线单元进行描述。

热源模型的选择对于模拟结果的准确性至关重要,应根据具体的焊接工艺选择合适的热源模型。

2. 材料属性及边界条件根据实际材料,设定焊件和焊缝的热导率、比热容、热扩散率等物理参数。

同时,设定初始温度、环境温度等边界条件。

3. 数值模拟过程根据焊接过程的实际情况,设定加载步和时间步长,模拟焊接过程中的温度变化。

通过ANSYS的热分析模块,得到焊接过程中的温度场分布。

三、焊接应力的数值模拟1. 耦合分析焊接过程中,温度场的变化会导致应力的产生。

因此,在ANSYS中,需要将在热分析中得到的温度场结果作为应力分析的输入条件,进行热-结构耦合分析。

2. 本构关系与材料模型根据材料的本构关系和力学性能,设定材料的弹性模量、泊松比、热膨胀系数等参数。

同时,选择合适的材料模型,如各向同性模型或各向异性模型。

3. 应力分析通过ANSYS的结构分析模块,结合耦合后的温度场结果,进行应力分析。

得到焊接过程中的应力分布和变化情况。

四、结果与讨论1. 温度场结果分析通过ANSYS的后处理功能,可以得到焊接过程中的温度场分布图。

分析温度场的分布情况,可以了解焊接过程中的热传导和热扩散情况,为优化焊接工艺提供依据。

2. 应力结果分析同样,通过后处理功能可以得到焊接过程中的应力分布图。

分析应力的分布和变化情况,可以了解焊接过程中产生的残余应力和变形情况。

ANSYS热分析指南

ANSYS热分析指南

第一章简介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量〕等。

热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。

二、ANSYS的热分析∙在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。

∙ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。

∙ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。

此外,还可以分析相变、有内热源、接触热阻等问题。

三、ANSYS 热分析分类∙稳态传热:系统的温度场不随时间变化∙瞬态传热:系统的温度场随时间明显变化四、耦合分析∙热-结构耦合∙热-流体耦合∙热-电耦合∙热-磁耦合∙热-电-磁-结构耦合等第二章 基础知识一、符号与单位W/m 2-℃ 3二、传热学经典理论回顾热分析遵循热力学第一定律,即能量守恒定律:●对于一个封闭的系统(没有质量的流入或流出〕PE KE U W Q ∆+∆+∆=-式中: Q —— 热量;W —— 作功;∆U ——系统内能; ∆KE ——系统动能; ∆PE ——系统势能;● 对于大多数工程传热问题:0==PE KE ∆∆; ● 通常考虑没有做功:0=W , 则:U Q ∆=;●对于稳态热分析:0=∆=U Q ,即流入系统的热量等于流出的热量; ●对于瞬态热分析:dtdUq =,即流入或流出的热传递速率q 等于系统内能的变化。

三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。

热传导遵循付里叶定律:dxdTkq -='',式中''q 为热流密度(W/m 2),k 为导热系数(W/m-℃),“-”表示热量流向温度降低的方向。

基于ANSYS的温度场仿真分析

基于ANSYS的温度场仿真分析
测试实验 , 来 判 定 绝 缘 材 料 的绝 缘 性 能 。 因 此 , 有 必 要 研 究 恒 温
式中 : [ C] 为 比热 矩 阵 , 考 虑 系 统 内能 的增 加 : [ K] 为 传 导 矩
阵, 包含导热系数 、 对流 系数及辐射率 和形 状 系数 ; { T} 为节点温
度向量 ; { T} 为温度对时间的导数 ; { Q( t ) } 为 节 点 热 流 率 向量 , 包 含 热 生成 。 温 度箱 中温 度 场 的热 分 析 属 于 瞬态 热 传 导过 程 。
绝 缘 材 料 在 高 温 条 件 下 长 期 运 行 会 造 成 热 老 化 ,使 绝 缘 性 能下降 , 在 高压 作 用 下 易 击 穿 。 为 进 一 步 研 究 绝 缘 材 料 热 老 化 特 性, 需 将 绝 缘 材 料 放 置 在 温 度 箱 内长 时 间加 热 , 再 进 行 绝 缘 特 性
元 类 型 、定 义 材 料 相
导热 系数 比热 容
( Wl m・ x、 J / k g ‘ K
0 . 1 2 1 0 3 0
生热 率
W/ m3
密度
k g / m
1 . 2 型 通 用 有 限元 分 析 软 件 , 能 够 进 行 机 械 应
1 0 4
基 于 AN S Y S的 温 度 场 仿 真 分 析
基于 A N S Y S的温度场仿真分析
T e mp e r a t u r e F i el d Si mu l a t i o n An a l y s i s B a s e d o n ANS YS
潘从 芳 娄 毅 蔺 红 张起 瑞 杨 一 胡 贺 明
( 新 疆 大学 电 气工程 学 院 , 新疆 鸟 鲁木 齐 8 3 0 0 4 7 )

Ansys Workbench在电机温度场分析的实际运用

Ansys Workbench在电机温度场分析的实际运用

Ansys Workbench在电机温度场分析的实际运用摘要:温升高是电机最为主要的故障原因,而电机的种类很多,不同种类有着多种多样冷却方式,因此,电机的温度分析较为复杂,传统方法是以热负荷作为基准根据试验结果类比电机的设计温升,对于一些特殊结构的电机,热负荷类比法就不能满足设计需要。

采用Ansys Workbench仿真软件通过FEA有限元分析(Finite Element Analysis),可以对特殊结构电机定转子热源分布、以及传导、对流、辐射等要素进行网格化分析。

本文以具体案例的设计分析过程,论述Ansys Workbench稳态温度场在电机设计中的实际运用。

关键词:温升电机温度场有限元Ansys1引言我们以一台低压变频异步电动机YVF400-6-315KW、380V、50HZ为研究对象,对其定转子温度场进行仿真分析,对比求解结果与最后型式试验的偏差,从而验证Ansys Workbench仿真软件在特殊电机设计的实际运用。

电动机主要的设计参数如下:2 2D建模2.1 在Ansys Workbench程序界面下,通过ToolBox,进入稳态温度分析Steady-state-Thermal工作平台;2.2 在Steady-state-Thermal工作平台点取Geometry进行几何形状设置;或者在Analysis Systems树状下右侧窗口,右键选取Geometry--Import导入定子或转子三维部件的stp、sat、step等格式;2.3 在Steady-state-Thermal工作平台右键选取Geometry-第二行Edit Geometry in DesignModeler----可进入DM-右键点取Import1,选取生成-Generate;可获得每个部件的建模信息。

3、材料设置3.1 回到workbench对第三行Model右键Edit---进入Mechancial Enterprise机械单元,在Outline下的Model点取Geometry项下的每个零部件,左下表格中可以查看体积、面积、重量,并且设置材料名称、密度、热导率等;3.2 定转子材料设置,40度左右的热导率按下表:4、划分网格在workbench对第三行Model右键Edit进入Mechancial Enterprise机械单元,菜单栏点取Generate Mesh ,在树状mesh下表格relevance设置网格的相关属性,数字越小,节点和单元数越少。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

目录
1 概述------------------------------------------ 2
2 CATIA建模过程--------------------------------- 3
3 ANSYS分析过程------------------------------- 10
4 结果分析-----------------------------------------14
5 参考文献--------------------------------------- 15
1.概述
1.16125柴油机活塞基本条件:
缸径D=125mm,6缸。

活塞是发动机的重要部件之一,与连杆构成发动机的心脏,活塞通过运动将燃气压力传递给连杆再至曲轴输出,工作时受力非常复杂。

随着发动机向高速度、低能耗方向发展,采用优异的活塞材料尤为重要。

目前车用发动机活塞材料以铝合金为主,其他还有铸铁、铸钢、陶瓷材料等。

铝合金的突出优点是密度小,可降低活塞质量及往复运动惯性,因此铝合金活塞常用于中、小缸径的中、高速发动机上。

与铸铁活塞相比,铝合金活塞导热性好,工作表面温度低,顶部的积碳也较少。

活塞由活塞顶、头部、群部构成。

活塞顶的形状分为平顶、凸顶、凹顶。

平顶活塞结构简单、制造容易、受热面积小、应力分布较均匀、多用在汽油机上;凸顶活塞顶部突起成球状、顶部强度高、起导向作用、有利于改善换气过程。

凹顶活塞可改变可燃混合气的形成和燃烧,还可以调节压缩比。

活塞工作时温度很高,顶部可达600 ~700K,且温度分布很不均匀;活塞顶部承受的气体压力很大,特别是作功行程压力最大,柴油机活塞顶燃烧最高压力5~9Mpa,这就使活塞产生冲击和侧压力的作用;根据活塞实际最大爆发压力工况添加边界条件,选用压力为5MPa
便于做有限元分析,此方案采用w顶活塞,用于六缸发动机
2.设计的初步准备:
1.选好各个值的长度
2.了解6125柴油机发动机
3.学会catia软件建模
4.学会ansys软件进行对模型的温度场分析
5.了解发动机活塞的性能,并能做出正确分析
2.活塞主要结构尺寸计算:
2.1选定各个比例系数:
D=125mm
H: (0.8~1.3)*D 取 150 H1: (0.5~0.8)*D 取 80
H2: (0.4~0.8)*D 取 100 h1: (0.1~0.2)*D 取 15
h3: (0.3~0.4)* H2 取30
h4: (0.6~0.7)* H2 取 70
d: (0.3~0.38)*D 取 40
C1: (0.04~0.08)*D 取5
其他环岸: (0.025~ 0.045)*D 取 5
B: (0.35~0.42)*D 取44
2.3设定其他参数:
顶部厚度:15mm
活塞度: 16mm
3.活塞三维建模:
1.首先打开catia软件,点击机械设计·草图绘制器·进入xy作平面
进行绘制
图3.1
如图3 .1 绘制之后,退出草图工作平面
2.然后点击回转体按钮,绕y方向短轴旋转360 ,,如图3-2,所示。

图3-2
4.选择xy平面绘制如下图形
图3-3
5.退出草图选择凸台对其进行拉伸,拉伸厚度为48mm如下图
图3-4 6.对上述4中的图形进行凹槽,凹槽深度为22mm
图3-5 7.再点击xy平面绘制圆,其半径为20如下图
8.退出草图对其进行凹槽,凹槽深度为70如图
9.
9绘制去重孔并对其进行凹槽深度为50mm
4ANASYS分析.
1. 将CATIA中的模型导入ANSYS:
2.单位换算:
ANSYA Main Menu→ Preprocessor→ Modeling→ Operate → Scale →V olumes选择实体,在RX,RY,RZ Scale factor 分别中输入,0.001,0.001,0.001
3. 单元类型设置
Main Menu→Preprocessor→Element Type→Add/Edit/Delete
→Add→Solid→Tet 10node 87
4.划分网格
Main Menu→Preprocessor→Meshing→Meshtool,打开Meshtool对话框,勾选Smart Size,设定为6,点击“Mesh”按钮,选择活塞模型进行自由网格划分,如图4-1
图4-1
5.材料属性
活塞材料及其属性:铝合金,弹性模量为70e9,泊松比为0.3
Main Menu→Preprocessor→ Meterial Props→Material Models→Thermal→Conductivity→Isotropic→KXX=160
Main Menu→Preprocessor→ Meterial Props→Material Models→Favorate→Isotropic→EX=70e9 PRXY=0.3
6.定义求解类型:
Main Menu→Solution→Analysis Type→New Analysis,设定分析类型为Steady-State
7.施加载荷:
Main Menu→Solution→Define Loads→Apply→Thermal→Convention→On Areas,分别对活塞的各个表面施加的参数如下表4-2
各部分如下图所示. 表4-2
输入各项数据后,得到如下图所示的模型,如图4-3
图4-3
8.变形
将得到的模型进行变形后得到结果如图4-4
图4-4
9.求解
Main Menu→Solution → Solve → Current LS
10.显示结果
General Postproc → Plot Results → Contour Plot → Nodal Solu → DOF Solution →Nodal Temperature点击OK得到结果如图4-5
图4-5
五.结论及结果分析
由温度场可以看出,活塞顶部与燃气接触的表面最高温度约为377℃,位于活塞顶部中心处,因此活塞的最高温度在活塞顶面位置,这要求活塞有较大的导热截面以防止大的热应力和烧蚀发生。

从总体上考虑,活塞的火力岸区域是受热最为严重的地方,所以要充分的考虑活塞的材质,提高活塞的抗热性,以防活塞顶部发生热变形。

活塞环槽上下面及内侧面的温度差为20~30℃,温度梯度相对较大,加上此处的尖角多,又是活塞散热的重要途径,所以应特别考虑此处的结构。

特别是第一环槽的温度差很大,应予以充分的重视
活塞的轴向温度大约在235℃,所以在设计活塞的时候,要考虑活塞在受热时产生的轴向的热变形
在外侧面,沿轴向从上到下,温度由高到低,第一环槽处温度为380℃左右,
活塞销座处最高温度大约为250℃,内腔面最高温度大约为300℃.从分析结果来看,活塞顶部的温度稍高,特别是第一环槽部分的温度过高,活塞顶部形成热积累。

活塞环区的温度对于发动机的可靠性是极为重要的,环区温度过高,将使润滑油变质甚至碳化,造成活塞环粘结,失去活动性,使环槽迅速磨损,变形,严重时将造成发动机气缸套擦伤,甚至拉缸。

经过对活塞的温度场分析,该活塞顶部热负荷偏大,因此有必要对活塞进行优化设计,以改善活塞本体的传热性能并优化活塞的结构
六.参考文献
[1]丁仁亮.CATIA V5 基础教程[M].北京:机械工业出版社,2006.10 [2]张乐乐.ANSYS辅助分析应用基础教程上机指导[M].北京:清华大学出版社;北京交通大学出版社,2007.12
[3]蔡春源.机械零件设计手册[M].第三版.北京:冶金工业出版社,1994 [4]胡海龙.CATIA V5零件设计实例教程. [M].北京:清华大学出版社,2007.9。

相关文档
最新文档