《金属塑性加工原》课件

金属塑性成形原理材料科学与工程专业核心基础课程Contents课程目录金属塑性成形原理——从基础理论到工程应用的系统性知识框架01绪论与金属塑性加工基础02塑性变形的应力分析03塑性变形的应变分析04屈服准则与塑性条件05塑性变形本构关系06塑性成形工程解法CHAPTER01绪论与金属塑性加工基础从工业文明的高度认识塑性成形的地位与课程知识版图METALPLASTICFORMING金属塑性成形的基本概念与发展历程从公元前5000年天然金锻打到现代万吨级液压机，人类利用塑性加工的历史跨越七千年金属锻压车间实景基本概念与加工分类金属在外力作用下产生塑性变形，获得所需形状、尺寸和组织性能的基本加工技术涵盖锻造、轧制、挤压、拉拔等七大类，按温度分为冷、温、热加工相比切削加工具有材料利用率高、力学性能改善、适合大批量生产等优势古代铁器锻造工艺历史发展与工业演进公元前5000年已有天然金锻打制品，商代出土铁刃铜钺证明早期锻造技术已存在1766年出现二辊轧机，1836年蒸汽机驱动轧机标志进入近代大工业生产阶段中国1871年建立最早轧钢厂，1997年钢产量达1.08亿吨居世界第一CourseFramework课程知识体系与学习目标《金属塑性成形原理》至20世纪40年代已发展为完整的应用基础科学，由五大板块构成，是连接材料微观机理与工程宏观应用的桥梁课程。五大理论板块塑性加工力学：研究变形体中应力、应变分布规律及力学平衡方程，是成形力能计算的数学基础塑性变形材料学：关注变形过程中材料组织结构与力学性能的演变规律，揭示加工硬化与再结晶机理变形体物理学与摩擦学：从微观层面解释塑性变形机理，并分析工件与模具接触面的摩擦行为5Modules核心学习目标理解金属塑性加工的基本概念，体会塑性变形的微观机理与宏观规律之间的内在联系掌握变形材料与工模具的组织结构和性能变化规律，为工艺设计提供理论依据获知工件尺寸精度、表面质量及性能控制的基本原理和技术途径3Objectives预备知识要求材料学基础（晶体结构、位错理论）、力学基础（静力学、动力学）和流变学基础（应力应变关系）机械加工基础（冷热加工、压力加工、锻造工艺）以及微积分、线性代数等数学工具FoundationMECHANICSFUNDAMENTALS弹性变形、塑性变形与应力基本概念弹性变形与塑性变形是金属受力后的两种本质不同的响应模式。弹性变形可逆、遵循胡克定律，塑性变形不可逆、伴随晶格滑移与位错运动。在塑性加工过程中，工件同时经历弹塑性变形，而应力作为描述内力分布强度的基本物理量，以张量形式完整表征了材料内部各点的受力状态。01弹性变形外力去除后完全恢复，应力与应变呈线性关系（胡克定律），是塑性成形中回弹计算的物理基础胡克定律02塑性变形不可恢复的永久变形，微观上通过晶格滑移和孪生实现，宏观上表现为工件几何形状的永久改变晶格滑移03外力与内力外力分为体积力（重力、惯性力等）和表面力（接触力、分布力），内力是材料抵抗外力产生的内部作用力体积力·表面力04应力张量应力是单位面积上的内力，通过9个分量（3个正应力+6个切应力）构成应力张量σ_ij，完整描述点的应力状态σ_ij·9分量05正应力与切应力正应力σ垂直于截面，引起材料拉伸或压缩；切应力τ平行于截面，是驱动塑性滑移变形的主要力学因素σ·τChapter02塑性变形的应力分析建立描述变形体内力分布的完整数学框架StressAnalysis主应力、主切应力与应力张量不变量主应力是切应力为零的特殊截面上的正应力，代表材料内部受力最纯粹的方向。应力张量的三个不变量I₁、I₂、I₃与坐标系无关，是后续屈服判据和塑性条件的数学基石。主应力与特征方程主应力σ₁≥σ₂≥σ₃作用于切应力为零的主平面上，通过求解应力张量的特征方程σ³−I₁σ²+I₂σ−I₃=0获得。σ³−I₁σ²+I₂σ−I₃=0三个应力不变量第一不变量I₁=σ₁+σ₂+σ₃反映静水压力大小，第二不变量I₂与形状改变相关，第三不变量I₃影响屈服面形状。I₁·I₂·I₃主切应力与最大切应力主切应力τ₁₂、τ₂₃、τ₁₃分别作用于两主应力方向的45°截面上，最大切应力τmax=(σ₁−σ₃)/2是Tresca准则的核心。τ_max=(σ₁−σ₃)/2坐标无关性与客观性应力张量不变量的坐标无关性确保了屈服判据和塑性条件的客观性，不随观察者坐标系旋转而改变。ObjectivityContinuumMechanics应力偏张量、球张量与等效应力应力张量可唯一分解为球张量与偏张量。球张量仅引起体积弹性变化；偏张量才是驱动塑性流动的根本力学因素。01应力球张量σ_m=(σ₁+σ₂+σ₃)/3代表各向等压的静水压力，仅引起弹性体积变化，金属塑性变形中体积不变，故球张量不贡献塑性功02应力偏张量s_ij=σ_ij−σ_m·δ_ij是应力张量减去球张量的剩余部分，其第二不变量J₂直接决定Mises等效应力的大小03八面体切应力τ₈=√(2J₂/3)作用于与三主轴等倾的八面体面上，物理意义是各方向剪切强度的统计平均04Mises等效应力σ̄=√(3J₂)=√[(σ₁−σ₂)²+(σ₂−σ₃)²+(σ₃−σ₁)²]/√2，将三向应力等效为单轴拉伸应力值SolidMechanics·StressAnalysis应力莫尔圆与应力平衡微分方程应力莫尔圆以几何图形直观表达了任意截面上的正应力与切应力关系；应力平衡微分方程从微元体力平衡出发建立应力分量的空间分布约束，三个方程含六个未知量构成静不定问题。01莫尔圆以正应力σ为横轴、切应力τ为纵轴，圆心坐标为((σ₁+σ₃)/2,0)，半径为(σ₁−σ₃)/2，直观表达应力状态。σ,τAxes02三向应力状态下形成三个嵌套莫尔圆，任意截面的应力状态点落在最大圆内、小圆外的阴影区域。3NestedCircles03直角坐标系下应力平衡微分方程∂σ_x/∂x+∂τ_yx/∂y+∂τ_zx/∂z+X=0（及y、z方向类似）描述微元体力平衡。∂σ/∂xEquilibrium04平衡方程3个、未知应力分量6个，问题为一次静不定，需引入屈服准则（提供1个方程）和本构关系补充方程求解。3eq./6UnknownsCHAPTER03塑性变形的应变分析量化描述金属塑性变形中的几何变化与运动学关系CONTINUUMMECHANICS位移、应变张量与小应变几何方程应变是位移场空间梯度的对称化度量，反映了材料局部几何形状的变化程度。在小变形条件下，应变张量由位移梯度的线性对称部分组成，其数学结构与应力张量完全对称——同样存在主应变、应变不变量和莫尔圆等分析工具，形成应力-应变的完整对偶体系。01位移场变量：位移向量u(x,y,z)描述各点从初始位置到变形后位置的移动，是变形运动学的基本场变量02几何方程：小应变几何方程εx=∂u/∂x、γxy=∂u/∂y+∂v/∂x等六个方程将应变分量与位移梯度联系起来03物理含义：正应变εi描述沿坐标轴方向的线段伸缩率，切应变γij描述两坐标轴之间直角的改变量04对偶结构：应变张量与应力张量数学结构对称，同样存在三个主应变ε₁、ε₂、ε₃和三个应变不变量，可绘制应变莫尔圆KINEMATICS·CONSTITUTIVE塑性变形体积不变条件与应力-应变状态对比金属塑性变形的体积不变条件（ε₁+ε₂+ε₃=0）将三个主应变的独立变量从3个降为2个，显著简化了工程问题的求解。应力状态与应变状态需通过本构方程建立定量联系。01体积不变条件ε₁+ε₂+ε₃=0（应变第一不变量为零），物理本质是金属塑性滑移不改变原子间距和晶胞体积ε₁+ε₂+ε₃=002独立变量降维将三个主应变的独立变量从3个减少到2个，是平面应变和轴对称问题分析的关键约束方程3→2独立变量03Levy-Mises本构方程应力偏张量sij与应变增量dεij成正比，而非应力与应变直接对应sij∝dεij04偏张量驱动塑性流动三向等压下偏应力为零、不产生塑性应变，球张量不参与塑性变形，仅偏张量驱动流动偏张量驱动CHAPTER04屈服准则与塑性条件判定材料从弹性状态进入塑性变形的力学临界条件YieldCriteriaTresca与Mises屈服准则Tresca准则基于最大切应力判据，物理图像直观但需预知主应力排序；Mises准则基于畸变能判据，数学形式光滑且无需主应力排序，与多数金属材料的实验数据吻合更优。Tresca最大切应力准则当最大切应力τmax=(σ₁−σ₃)/2达到材料临界值k时发生屈服，即σ₁−σ₃=σs。使用前需预知主应力的大小顺序。σ₁−σ₃=σsMises畸变能准则当等效应力σ̄=√[(σ₁−σ₂)²+(σ₂−σ₃)²+(σ₃−σ₁)²]/√2达到σs时屈服，数学形式光滑，无需预知主应力排序。σ̄=σsπ-Planeπ平面几何对比在π平面上，Tresca准则呈正六边形，Mises准则为其外接圆。两者在纯剪切状态下差异最大，单轴拉伸时完全重合。≈15.5%Engineering工程应用首选大量金属材料（低碳钢、铜、铝等）的实验数据更接近Mises圆，工程计算中优先采用Mises准则作为屈服判据。MisesPreferredYIELDCRITERIA屈服准则的几何表达与工程比较在主应力空间的π平面上，Tresca准则呈正六边形、Mises准则为外接圆，两者在单轴拉伸时重合，在中间主应力等于两端均值时差异最大（约15.5%）。01主应力空间中屈服面是以静水压力轴（σ₁=σ₂=σ₃线）为中心线的无限长柱面，π平面截面即为屈服轨迹π平面·屈服轨迹02单轴拉伸时σ₁=σs、σ₂=σ₃=0，两准则均给出σ₁=σs；纯剪切时Tresca给出τ=k=σs/2，Mises给出τ=σs/√3，差异约15.5%≈15.5%03加工硬化材料的屈服面会随塑性变形累积而扩张（各向同性硬化）或平移（随动硬化），屈服准则需配合硬化模型使用各向同性·随动硬化04工程设计中Mises准则因数学光滑（处处可微）更适合数值计算，Tresca准则因形式简单在手算分析中仍有优势Mises光滑·Tresca简洁CHAPTER05塑性变形本构关系建立应力与塑性应变之间的定量数学模型ConstitutiveTheory弹性本构关系与塑性增量理论弹性阶段应力-应变遵循广义胡克定律的线性关系，由杨氏模量和泊松比完全确定。进入塑性阶段后，增量理论（Levy-Mises方程dε_ij=dλ·s_ij）揭示塑性应变增量方向与应力偏张量方向一致，建立了塑性流动方向与当前应力状态的定量联系，是塑性成形数值模拟的理论核心。01[(1+ν)/E]σij−(ν/E)σkk·δij广义胡克定律：εij=[(1+ν)/E]σij−(ν/E)σkk·δij描述各向同性线弹性行为，由E和ν两个常数完全确定E·ν02Levy-Mises增量方程：dεij=dλ·sij表明塑性应变增量张量与应力偏张量成比例，dλ为非负标量乘子由加载条件确定dλ≥003核心假设：塑性流动方向仅取决于当前应力偏张量方向，与变形路径历史无关（适用于比例加载）比例加载04Prandtl-Reuss方程：在Levy-Mises基础上增加弹性应变增量，形成弹塑性增量本构，适用于需要考虑回弹的精密成形分析弹塑性PlasticityTheory全量理论与真实应力-应变曲线全量理论在比例加载下与增量理论等价，适合简单加载工程问题；幂律硬化模型σ=K·εⁿ是成形力计算的关键。01Hencky全量方程εᵖᵢⱼ=(3ε̄ₚ/2σ̄)·sᵢⱼ直接联系总塑性应变与应力偏张量，比例加载下与增量理论等价02简单加载适用各应力分量按比例增长时精度高；非比例加载（如先拉后扭）下精度降低，工程慎用03真实应力σ=P/A（A为瞬时截面积）消除颈缩面积变化影响，均匀变形阶段始终单调递增04幂律硬化σ=K·εⁿ中K为强度系数，n为硬化指数（0<n<1），n越大成形性越好金属材料拉伸试验·应力-应变曲线实测CHAPTER06塑性成形工程解法主应力法、上限法与滑移线法的理论原理与工程应用PLASTICFORMING·SLIP-LINEMETHOD主应力法的基本原理与镦粗力计算主应力法通过假设某一方向为主应力方向，沿该方向取微元体建立力平衡方程并配合屈服准则求解接触压力分布。以平砧镦粗为例，可解析求解径向压力分布并积分得到总成形力。01核心假设变形体某一方向为主应力方向，沿该方向取微元体建立平衡微分方程，将三维问题简化为一维求解。3D→1D02平衡方程取径向微元环，建立dσ_r/dr+(σ_r−σ_θ)/r=0的平衡方程，配合屈服准则σ_r−σ_z=σ_s求解接触面压力分布。σ_r−σ_z=σ_s03摩擦山效应无摩擦时压力均匀为σ_s；有摩擦时呈"摩擦山"形分布，中心压力最高、边缘降至σ_s。摩擦山04总镦粗力P=∫p·dA通过对接触面压力积分获得，是选择锻压设备吨位和设计模具强度的核心依据。P=∫p·dAPLASTICFORMINGMECHANICS主应力法在挤压与拉拔工艺中的应用主应力法在挤压和拉拔工艺中通过沿轴向建立微元体平衡方程，综合考虑模壁摩擦、锥面正压力和材料屈服强度，可解析计算成形力与工艺参数的定量关系。01正向挤压分析中在锥形变形区取微元体，沿轴向建立平衡方程，考虑模壁摩擦和锥面正压力的轴向分量求解挤压力P=f(σs,λ,α,μ)02挤压力与挤压比、模具半角和摩擦系数密切相关，主应力法可给出解析表达式λ=A₀/A₁03拉拔工艺通过主应力法计算拉拔力，确定不产生断裂的最大道次压缩率，为多道次拉拔工艺规程制定提供依据Multi-passDrawing04主应力法计算结果偏于保守，工程上作为设备吨位选择的安全下限，实际生产中以实验校准修正SafetyLowerBoundPLASTICITYTHEORY虚功原理与上下限定理上下限定理以虚功原理建立成形力边界：静可容应力场给出力下限，动可容速度场给出力上限。上限法构造速度场计算各功率之和估计成形力上限，适用广且灵活。虚功原理外力在虚速度场上做的虚功率等于内力在虚应变率场上做的虚功率，是能量守恒在连续介质中的表达能量守恒上限定理任意动可容速度场算出的外力功率≥真实外力功率，给出成形力的上限估计P*≥P下限定理任意静可容应力场算出的外力功率≤真实外力功率，给出成形力的下限估计P*≤P上限法计算步骤构造动可容速度场→计算变形功率+剪切功率+摩擦功率→总功率J*→上限力P*=J*/v₀J*=ΣẆUPPERBOUNDMETHOD上限法的速度场构造与典型应用上限法的核心技巧在于构造合理的动可容速度场。连续速度场用解析函数描述变形区速度分布，适合简单几何；三角形刚性块速度场将变形区离散为若干匀速运动刚性块，通过块间速度不连续面的剪切耗散逼近真实解。后者因简单灵活而成为平面应变问题中最常用的上限法工具。连续速度场用数学函数（如多项式、三角函数）描述变形区速度分布，需满足∇·v=0（体积不变）和边界速度条件解析函数三角形刚性块模型将变形区分成若干刚性三角形，块内速度均匀、块间存在速度不连续面，通过速度矢量图确定各面速度跳跃量速度矢量图收敛精度平面应变压缩问题中，三角形速度场的上限解随块数增加逐步逼近真实解（滑移线解），4块以上即可获得满意精度4块以上工程应用在锻造、挤压、轧制等工艺中广泛使用，通过优化速度场参数（如三角形角度、分块数）获得最紧的上限估计锻造·挤压·轧制SLIP-LINEFIELDTHEORY滑移线场理论的基本概念与Hencky方程滑移线是塑性变形区内最大切应力方向的轨迹，形成α线和β线两族正交曲线网络。沿滑移线的Hencky方程（p±2kθ=常数）建立了静水压力与滑移线转角之间的定量关系，使得从已知边界条件出发逐步递推求解整个变形区应力分布成为可能，是平面应变问题的精确解析工具。正交滑移线网络滑移线是塑性变形区内最大切应力方向的轨迹，α线和β线互相正交，在Tresca准则下k=σs/2为最大切应力Hencky方程沿α线p+2kθ=Cα、沿β线p−2kθ=Cβ，其中p为静水压力、θ为滑移线与参考轴的夹角Hencky第一定理α线沿β线移动时两端的转角相等，β线沿α线移动时亦然，保证滑移线场几何的一致性逐点递推求解滑移线法从已知应力边界条件出发，利用Hencky方程逐点递推构建完整滑移线场，进而求得全场的应力和速度分布SECTION04·PLASTICITY滑移线法典型应用与塑性成形中的摩擦滑移线法在平砧压缩、楔形模挤压等经典平面应变问题中可给出精确解析解，是验证其他近似方法精度的基准。滑移线法典型应用平砧压缩薄件（h/b<1）：滑移线场含中心扇形区和均匀应力区，精确解p=2k(1+b/4h)，是主应力法和上限法的精度校验基准p=2k(1+b/4h)楔形模正挤压：滑移线场由死区边界和变形区扇形场组成，可精确求解挤压力与模具半角、缩减率的关系精确解析解摩擦的影响与控制摩擦增加成形力（可达无摩擦的2–3倍）、加速模具磨损、导致变形不均匀和工件表面质量下降轧制和拉拔需适度摩擦驱动金属流动（咬入条件μ>tanα），过大摩擦导致能耗增加和表面缺陷润滑按温度选择：冷加工用矿物油/乳化液，温加工用MoS₂/石墨乳，热加工用玻璃润滑剂或涂层冷·温·热Chapter07塑性变形材料学变形过程中金属微观组织演变与宏观性能变化规律DislocationTheory加工硬化机理与位错理论加工硬化是金属塑性变形中流变应力随变形量增加而持续升高的现象，微观根源在于位错密度从退火态10⁶/cm²增至变形态10¹²/cm²，位错间的缠结、塞积和反应阻碍了进一步滑移。加工硬化既增加成形的力能需求，又是冷加工强化的核心机制，需通过合理的中间退火制度平衡成形性与工件强度。01Schmid定律与位错滑移金属塑性变形以位错滑移为主，位错沿密排面在密排方向上运动，临界分切应力τcr=σ·cosφ·cosλ。02位错密度与流变应力变形量增加使位错密度从10⁶–10⁸/cm²（退火态）增至10¹¹–10¹²/cm²（大变形），流变应力τ=τ₀+αGb√ρ。03三阶段硬化曲线易滑移阶段（单系滑移）、线性硬化阶段（多系交滑）与抛物线硬化阶段（交滑移启动）依次演进。04中间退火制度冷成形累积变形量过大时需中间退火消除加工硬化，退火温度和时间需精确控制以恢复塑性而不粗化晶粒。RECOVERY&RECRYSTALLIZATION回复与再结晶及冷热加工界限加工硬化金属加热时依次经历回复和再结晶两个软化过程，再结晶温度TR≈0.4Tm是冷热加工的分界线。0.2–0.3Tm回复阶段点缺陷湮灭和位错重排形成亚晶结构，内应力显著降低但强度和硬度变化不大0.4–0.5Tm再结晶阶段无畸变新晶粒在变形组织中形核长大，完全消除加工硬化，恢复退火态力学性能TR≈0.4Tm再结晶温度界限铁约450°C、铜约200°C、铝约100°C，是区分冷加工与热加工的物理分界线DYNAMICRECRYSTALLIZATION热加工动态再结晶动态再结晶与变形同步进行，稳态流变应力远低于冷加工；温度过高导致晶粒粗化TEXTURE&ANISOTROPY变形织构、各向异性与工程影响塑性变形使晶粒取向趋于集中形成织构，导致材料力学性能呈现方向依赖性（各向异性）。织构在板材深冲中引发"制耳"缺陷，但在硅钢片等产品中可被主动利用以优化磁学性能。理解和控制织构是高品质塑性成形件生产中的关键技术挑战。MECHANISM变形织构是晶粒在滑移过程中取向逐渐集中的结果。冷轧板材常见两种典型织构：{110}<112>黄铜织构和{123}<634>S织构，它们共同决定了板材的各向异性特征。{110}<112>黄铜织构EFFECT织构导致材料各向异性：不同方向的屈服强度、塑性和弹性模量各不相同。工程上用塑性应变比r值（Lankford系数）量化各向异性程度，r值越高深冲性能越好。r值量化各向异性DEFECT深冲制耳是典型缺陷：板材各向异性导致杯形件边缘高度不均匀。当Δr>0时出现4个制耳，严重影响材料利用率和零件精度，需通过退火工艺优化消除。Δr>0产生4制耳APPLICATION硅钢片Goss织构{110}<001>使易磁化方向沿轧向排列，变压器铁芯磁导率提升30%以上，铁损显著降低。这是织构工程最成功的工业应用案例。磁导率提升+30%QUALITYDEFECTS·质量工程塑性成形件质量缺陷分析与控制塑性成形件的典型质量缺陷包括折叠、裂纹和残余应力，其根源在于金属流动不均匀和应力分布不合理。质量控制需从工艺参数优化、模具几何改进和润滑条件调控三方面综合施策，现代有限元数值模拟（ABAQUS/DEFORM）已成为预测缺陷形成和实现工艺虚拟优化的核心工具。锻造件表面质量缺陷实拍·裂纹与折叠缺陷FOLD折叠是金属流动不均匀导致表面金属被压入工件内部形成的折缝缺陷，常见于模锻件飞边桥部和挤压件前端CRACK裂纹由局部拉应力超过材料断裂强度引起，预防方法包括控制变形量、提高变形温度、优化模具圆角半径以减小应力集中STRESS残余应力由变形不均匀导致，卸载后形成自平衡内应力场，可通过退火消除或振动时效降低，防止工件变形和应力腐蚀开裂FEA有限元模拟（ABAQUS、DEFORM、Forge）可预测金属流动轨迹、温度场和损伤分布，实现模具和工艺参数的虚拟优化FutureTrends现代塑性成形技术发展趋势沿数值模拟深化、智能化与绿色化三大方向演进，将传统经验驱动转变为数据驱动和模型驱动的精密制造。有限元仿真·金属成形模拟场景智能锻造生产线·自动化工业场景SIMULATION数值模拟与多尺度仿真热力耦合有限元模拟可预测锻造过程中的温度场、应变场和微观组织演变，实现模具与工艺参数的虚拟优化多尺度模拟将宏观有限元与微观晶体塑性模型耦合，可在介观层面预测织构演化和各向异性行为SMART&GREEN智能制造与绿色成形基于传感器和机器学习的在线质量监控系统可实时预测成形件质量并自动调整工艺参数，实现自适应控制近净成形技术（精密锻造、楔横轧等）将材料利用率从传统工艺的60-70%提升至90%以上，显著减少加工余量和废料MATERIALSCIENCE常用金属材料塑性成形性能参数不同金属材料因晶体结构和位错行为差异，在屈服强度、硬化指数、延伸率等塑性成形关键参数上表现出显著差异。低碳钢综合成形性优良、纯铝深冲性能突出、钛合金需热成形，工艺设计必须依据材料特性参数选择合适的变形温度、速率和润滑方案。常用金属材料室温塑性成形关键参数对比材料