版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
典型环节传递函数及伯德图控制系统频率响应分析方法与工程应用Contents目录控制系统频域分析的完整知识脉络,从数学建模到工程实践。01控制系统数学模型基础02伯德图基本原理与优势03八大典型环节传递函数与伯德图04系统伯德图绘制与稳定性分析05工程应用与课程总结Chapter01控制系统数学模型基础从物理系统到数学表达的建模方法论MATHEMATICALMODELS控制系统数学模型的三种形式控制系统数学模型是分析设计的基石,时域微分方程、复域传递函数和频域频率特性三种形式各有优势,其中传递函数因能将微分运算转化为代数运算而成为经典控制理论的核心工具。TIMEDOMAIN时域模型01用微分方程描述系统变量随时间的演化关系,给定初始条件可求解输出响应02直观反映物理过程,但参数变化时需重新列方程求解,不便于设计分析COMPLEXDOMAIN复域模型01通过拉普拉斯变换将微分方程转化为传递函数G(s),微分运算简化为代数运算02传递函数可表征系统动态性能,便于研究结构或参数变化对系统行为的影响FREQUENCYDOMAIN频域模型01将传递函数中s替换为jω得到频率特性G(jω),描述幅值放大与相位偏移02频率响应分析不依赖精确数学模型,可通过实验直接测量,工程实用性强TransferFunction传递函数的定义与核心性质传递函数G(s)=C(s)/R(s)是线性定常系统零初始条件下输出与输入拉氏变换之比,其零极点分布完全决定系统动态特性,且串联环节的传递函数可直接相乘,为典型环节分解提供了理论基础。01定义:线性定常系统在零初始条件下,输出拉氏变换C(s)与输入拉氏变换R(s)之比,记为G(s)=C(s)/R(s)02零极点:分子多项式的根为零点,分母多项式的根为极点,零极点在s平面的分布决定系统稳定性与动态性能03固有属性:传递函数仅取决于系统本身的结构和参数,与输入信号的形式和大小无关04串联分解:串联环节的等效传递函数等于各环节传递函数的乘积,复杂系统可分解为典型环节逐一分析核心公式G(s)=C(s)/R(s)输出拉氏变换与输入拉氏变换之比系统结构示意R(s)G(s)C(s)输入→传递函数→输出G₁(s)·G₂(s)串联相乘G(s)/[1+G(s)]反馈形式TypicalLinkDecomposition开环传递函数的典型环节分解任意线性系统开环传递函数经因式分解均可表示为八种典型环节的乘积,这种分解方法使复杂系统的频率响应分析转化为对基本单元的逐一研究再叠加,是伯德图方法的理论基石。01将开环传递函数G(s)H(s)的分子分母进行因式分解,可转化为K、1/s、s、1/(Ts+1)、Ts+1等典型环节的乘积形式02八种典型环节覆盖了一阶和二阶系统的全部动态特征,更高阶系统也可通过部分分式展开归结为这些基本单元03分解后各环节在伯德图上的幅频特性以分贝相加、相频特性以角度相加,将乘法运算转化为加法运算04最小相位系统的幅频特性与相频特性存在一一对应关系,只需分析幅频即可推断相频,进一步简化分析自动控制原理教材·传递函数因式分解CHAPTER02伯德图基本原理与优势对数坐标下的频率响应图示方法及其工程价值FrequencyResponse伯德图的定义与双图结构伯德图由对数幅频特性图和相频特性图两张图组成,横轴采用频率对数刻度,纵轴分别采用分贝和角度线性刻度,这种设计将串联环节的乘法转化为加法,极大简化了手工绘图与系统分析。频率响应分析实验场景01对数幅频特性图:纵轴为20lg|G(jω)|(单位dB),横轴为lgω,展示系统对各频率正弦信号的放大或衰减程度02相频特性图:纵轴为相角φ(ω)(单位:度),横轴为lgω,展示系统对各频率正弦信号造成的相位超前或滞后03对数坐标的核心优势:将传递函数的乘除运算转化为加减运算,串联环节只需将各自曲线纵向叠加即可得到总特性04宽频率覆盖:对数频率轴可在同一张图上覆盖极宽频率范围(如0.01到1000rad/s),兼顾低频稳态精度与高频动态特性ORIGIN&LEGACY伯德图的发明与工程地位伯德图由贝尔实验室科学家H.W.Bode于1940年提出,其核心创新是用对数坐标将复杂频率响应计算简化为直线近似,至今仍是控制工程领域最主流的频率分析工具。贝尔实验室·H.W.Bode工作场景01发明起源:HendrikWadeBode(1905-1982)在贝尔实验室研究反馈放大器稳定性时发明了对数频率响应图示法,后以他名字命名02核心创新:用渐近线近似替代精确曲线,使工程师无需计算机即可快速手绘系统频率响应,极大降低了分析门槛03工程地位:与Nyquist图、Nichols图并列为三大频域分析工具,因能同时直观展示幅值裕度和相位裕度而应用最广04持久生命力:现代虽有MATLAB等软件辅助精确绘制,但渐近线近似法仍是工程师快速估算系统特性的基本功BODEPLOT伯德图的三大工程优势伯德图之所以成为控制工程首选分析工具,源于其三重优势:对数运算简化计算、裕度判读直观可靠、实验测量无需精确模型,三者共同构成了从理论到实践的完整分析闭环。计算简化对数坐标将环节串联的乘除运算转化为加减运算,多环节叠加只需将曲线纵向移动后相加渐近线近似使手工绘图成为可能,转折频率处修正量固定且可查表,工程精度足够log→±稳定性直观判读增益裕度:在相角为-180°的频率处读取幅值距0dB的距离,正值表示系统稳定相位裕度:在增益穿越频率(幅值为0dB)处读取相角距-180°的距离,正值表示稳定GM·PM实验友好频率响应可通过实验直接测量:输入不同频率正弦信号,记录输出幅值比和相位差即可绘图无需事先知道精确传递函数,对难以建模的复杂系统(如化工过程)尤其适用无模型CHAPTER03八大典型环节传递函数与伯德图逐一剖析每种环节的数学表达、频率特性与图形特征TYPICALELEMENT比例环节:纯增益的静态映射比例环节G(s)=K是最基本的典型环节,其伯德图幅频为水平线20lgKdB、相频恒为0°,虽无动态特性但在系统设计中承担增益调节功能,直接影响稳态精度与系统稳定裕度。运算放大器比例放大电路实物01传递函数G(s)=K,频率特性G(jω)=K,幅值|G(jω)|=K为常数,不随频率变化02对数幅频特性L(ω)=20lgK为水平线:K>1时位于0dB上方表示放大,K<1时位于0dB下方表示衰减03相频特性φ(ω)=0°恒为零,比例环节不引入任何相位超前或滞后,对系统稳定性无额外影响04工程意义:增大K可提高稳态精度但降低稳定裕度,设计时需在精度与稳定性之间取得平衡INTEGRALELEMENT积分环节:累积效应与稳态误差消除积分环节G(s)=1/s的幅频特性以-20dB/dec斜率衰减高频信号,相频恒为-90°,其核心价值在于消除系统稳态误差,但代价是引入固定相位滞后,是控制系统精度与稳定性矛盾的关键节点。传递函数G(s)=1/s,频率特性G(jω)=-j/ω,幅值|G(jω)|=1/ω随频率增大而反比衰减对数幅频L(ω)=-20lgω过(1,0)点,斜率-20dB/dec,频率每增大10倍幅值下降20dB相频特性φ(ω)=-90°恒定不变,对所有频率信号均造成90°相位滞后工程意义消除稳态误差但降低相位裕度,过量积分导致系统振荡甚至不稳定RC积分电路实验实物DIFFERENTIATION微分环节:变化率检测与相位超前微分环节G(s)=s的幅频以+20dB/dec斜率放大高频信号,相频恒为+90°,能预判信号变化趋势提供相位超前补偿,但纯微分会无限放大高频噪声,工程中须用近似微分实现。信号微分电路·电子实验台01传递函数G(s)=s,频率特性G(jω)=jω,幅值|G(jω)|=ω随频率正比增大,高频增益趋向无穷G(s)=s02对数幅频特性L(ω)=20lgω为过(1,0)点、斜率+20dB/dec的直线,与积分环节关于0dB轴镜像对称+20dB/dec03相频特性φ(ω)=arg(jω)=+90°恒定不变,对所有频率信号均提供90°相位超前φ=+90°04工程限制:纯微分环节物理不可实现(高频噪声被无限放大),实际使用一阶微分环节Ts+1或带低通滤波的近似微分器Ts+1CONTROLTHEORY·TRANSFERFUNCTION惯性环节:一阶低通与转折频率惯性环节G(s)=1/(Ts+1)以一阶低通特性广泛应用于系统建模,转折频率ω=1/T将伯德图分为低频0dB平台与高频-20dB/dec衰减段,时间常数T决定了系统带宽与响应速度。传递函数:G(s)=1/(Ts+1),T为时间常数,转折频率ω=1/T是幅频特性由平坦转为衰减的临界点幅频渐近线:低频段(ω≪1/T)为0dB水平线,高频段(ω≫1/T)以-20dB/dec下降,转折频率处实际值比渐近线低约3dB相频特性:φ(ω)=-arctan(ωT)从0°单调下降至-90°,转折频率处φ=-45°,变化速率与T值直接相关物理意义:T越大则转折频率越低、系统带宽越窄、阶跃响应越慢,典型实现为RC低通滤波器和热惯性系统RC低通滤波器·惯性环节典型物理实现FREQUENCYRESPONSE一阶微分环节:高通特性与相位超前补偿一阶微分环节G(s)=Ts+1与惯性环节互为对偶,低频段直通、高频段以+20dB/dec放大并提供0°到+90°的相位超前,是PD控制器改善系统动态性能的核心机制。PD控制器工业应用设备01传递函数G(s)=Ts+1,转折频率ω=1/T,T为微分时间常数,与惯性环节1/(Ts+1)构成零极点对偶关系02幅频渐近线:低频段(ω≪1/T)为0dB水平线,高频段以+20dB/dec上升,转折频率处实际值高约3dB03相频特性φ(ω)=arctan(ωT):从0°单调上升至+90°,转折频率处φ=+45°,为系统提供正相位补偿04工程应用:PD控制器的微分作用即体现为一阶微分特性,在穿越频率附近提供相位超前以提高相位裕度和响应速度FREQUENCYRESPONSE振荡环节:二阶谐振与阻尼效应振荡环节G(s)=ωn²/(s²+2ζωns+ωn²)是二阶系统核心,阻尼比ζ决定谐振峰高度和相位变化陡度,ωn决定转折频率位置,二者共同刻画了系统从过阻尼到欠阻尼的全部动态行为。01传递函数—G(s)=1/((s/ωn)²+2ζ(s/ωn)+1),ωn为无阻尼自然频率,ζ为阻尼比,0<ζ<1时为欠阻尼状态02幅频特性—低频0dB,高频以−40dB/dec下降;ζ<0.707时出现谐振峰,ωr=ωn√(1−2ζ²),峰值Mr=1/(2ζ√(1−ζ²))03相频特性—从0°下降至−180°,在ω=ωn处恒为−90°;ζ越小,相位在ωn附近变化越陡峭04物理实现—弹簧-质量-阻尼系统、RLC串联电路、直流电机调速系统均可建模为振荡环节弹簧-质量-阻尼系统振动实验装置BodeDiagramAnalysis阻尼比ζ对振荡环节伯德图的影响阻尼比ζ是振荡环节最关键的参数,决定了谐振峰的存在与否和高度:ζ<0.707出现谐振峰且ζ越小峰越高,ζ=0.707为最佳阻尼比无谐振峰,ζ>1时系统过阻尼退化为两个惯性环节串联。不同阻尼比ζ下的二阶系统阶跃响应波形对比01ζ<0.5谐振峰显著ζ=0.1时峰值Mr≈5倍(14dB),ζ=0.3时Mr≈1.8倍(5dB),谐振频率ωr=ωn√(1−2ζ²)Mr≈5×02ζ=0.707最佳阻尼比谐振峰刚好消失,阶跃响应超调量约4.3%,兼顾速度与平稳性σ≈4.3%03ζ≥1系统过阻尼振荡环节退化为两个惯性环节串联,传递函数分解为1/(T₁s+1)·1/(T₂s+1)T₁·T₂串联04渐近线修正振荡环节渐近线在转折频率附近的误差与ζ密切相关,须根据ζ值查表或计算精确修正ω修正SECOND-ORDERDIFFERENTIATION二阶微分环节:振荡环节的对偶二阶微分环节G(s)=(s/ωn)²+2ζ(s/ωn)+1与振荡环节互为对偶,高频段以+40dB/dec上升、相位从0°超前到+180°,在工程中较少单独出现但对分析系统零点和非最小相位特性不可或缺。传递函数G(s)=(s/ωn)²+2ζ(s/ωn)+1,ωn为转折频率,ζ为阻尼比,与振荡环节互为倒数关系幅频特性:低频段0dB,高频段以+40dB/dec上升;ζ<0.707时在ωn附近可能出现幅值凹陷或突起的反谐振特征相频特性:从0°上升至+180°,在ω=ωn处为+90°,变化趋势与振荡环节关于0°轴对称工程场景:常见于非最小相位系统分析、高阶控制器设计中的零点对效应,以及振动主动控制中的加速度反馈振动控制系统·加速度传感器工业应用FrequencyResponse·DelayElement延迟环节:纯时滞与相位侵蚀延迟环节G(s)=e^(-τs)的幅频恒为0dB但相频随频率线性下降(-τω),高频段的巨大相位滞后严重侵蚀系统稳定裕度,是过程控制中最棘手的环节,催生了Smith预估器等专门补偿方法。01传递函数G(s)=e−τs,τ为纯延迟时间,频率特性G(jω)=e−jωτ,幅值|G(jω)|=1恒为常数02对数幅频L(ω)=20lg1=0dB恒为零,纯延迟不改变任何频率信号的幅值大小03相频特性φ(ω)=−τω(弧度)=−57.3τω(度),随频率线性无限下降,高频段可产生远超−180°的相位滞后04工程影响延迟使Nyquist图螺旋内卷、伯德图相位裕度急剧缩小,Smith预估器通过内部模型补偿延迟改善控制性能化工过程控制中的管道传输——纯延迟环节的典型物理场景,流体从入口到出口需要经过固定时间τFREQUENCYANALYSIS八大典型环节伯德图特征对比八种典型环节构成完整的频率分析基本单元库,其中积分/微分、惯性/一阶微分、振荡/二阶微分三组对偶关系的伯德图分别关于0dB和0°轴对称。环节名称传递函数幅频特征相频范围转折频率比例环节K20lgK(水平线)0°无积分环节1/s−20dB/dec−90°无微分环节s+20dB/dec+90°无惯性环节1/(Ts+1)0dB→−20dB/dec0°→−90°1/T一阶微分Ts+10dB→+20dB/dec0°→+90°1/T振荡环节ωn²/(s²+2ζωns+ωn²)0dB→−40dB/dec0°→−180°ωn二阶微分(s/ωn)²+2ζs/ωn+10dB→+40dB/dec0°→+180°ωn延迟环节e^(−τs)0dB(水平线)−τω(线性下降)无对偶环节关于0dB轴和0°轴对称,掌握特征即可通过叠加法分析任意系统CHAPTER04系统伯德图绘制与稳定性分析从环节叠加到系统级频率响应的完整分析方法绘制流程系统伯德图的标准绘制流程系统伯德图绘制遵循"分解→标定→起始→叠加→修正"五步法,核心是利用对数坐标的加法性质将各环节渐近线逐段叠加,最终在转折频率处做修正即可得到精确曲线。01标准形式分解—将G(s)H(s)化为典型环节乘积的标准形式,提取增益K、各环节时间常数T和自然频率ωn02转折频率标定—计算所有转折频率(1/T、ωn等),在对数频率轴上从低到高标注,确定需要关注的频率范围03低频起始渐近线—绘制低频段起始渐近线,起始高度由20lgK决定,起始斜率由积分/微分环节个数决定(每个积分−20dB/dec)04斜率逐段叠加—从低频向高频扫描,每经过一个转折频率改变斜率——惯性环节−20dB/dec、一阶微分+20dB/dec、振荡环节−40dB/dec05相频合成与幅频修正—各环节相频曲线角度相加得总相频,转折频率附近根据实际公式修正幅频曲线(振荡环节按ζ值修正谐振峰)FREQUENCYRESPONSE·LOW-FREQASYMPTOTE低频段渐近线与系统型别低频渐近线由增益K和积分环节数v共同决定:起始高度20lgK、斜率-20vdB/dec,系统型别v=0,1,2分别对应阶跃、斜坡、加速度输入的零稳态误差。v=00型系统无积分环节,低频段为水平线20lgKdB,斜率为0。Bode图低频段呈现恒定幅值特性,系统对低频信号具有恒定增益放大能力。跟踪阶跃输入有稳态误差ess=1/(1+K),跟踪斜坡输入稳态误差无穷大,无法跟踪加速度输入。ess=1/(1+K)v=1I型系统含1个积分环节,低频段斜率-20dB/dec,过(1,20lgK)点。积分作用使系统具有记忆特性,能够消除阶跃输入的稳态误差。跟踪阶跃输入稳态误差为零,跟踪斜坡输入有稳态误差ess=1/K,开环增益K越大跟踪精度越高。ess=1/Kv=2II型系统含2个积分环节,低频段斜率-40dB/dec,过(1,20lgK)点。双重积分提供更强的低频衰减和更高的系统型别。跟踪阶跃和斜坡输入均无稳态误差,跟踪加速度输入有稳态误差ess=1/K,适用于高精度伺服控制系统。ess=1/KBODESTABILITYANALYSIS增益裕度与相位裕度的定义及判读增益裕度GM和相位裕度PM是伯德图稳定性分析的两大核心指标,分别从幅值和相位维度量化系统距不稳定边界的距离,工程上通常要求GM>6dB且PM>30°~60°以保证足够的稳定储备。增益裕度GMDEFINITION在相角穿越频率ωg处(φ(ωg)=−180°),幅值距0dB的距离,GM=−L(ωg)dBSTABILITYGM>0表示系统稳定(幅值在0dB以下),GM越大系统对增益变化的容忍度越高PHYSICALMEANING系统开环增益可再增大10(GM/20)倍才会到达临界不稳定状态相位裕度PMDEFINITION在增益穿越频率ωc处(L(ωc)=0dB),相位距−180°的距离,PM=180°+φ(ωc)STABILITYPM>0表示系统稳定(相位在−180°以上),PM越大系统瞬态响应的超调越小ENGINEERINGPM≈30°~60°时系统响应平稳且速度适中,PM过小则振荡剧烈、过大则响应迟缓BodePlot·WorkedExample实例演练:三阶系统伯德图绘制以G(s)=100/(s(s+10)(s+100))为例演示完整绘制流程:标准化后识别出K=0.1、1个积分环节和2个惯性环节,转折频率10和100rad/s,渐近线斜率从-20依次变为-40、-60dB/dec。控制工程课堂·伯德图手绘练习Step01·标准化G(s)=0.1/[s(0.1s+1)(0.01s+1)],增益K=0.1,积分环节v=1,转折频率ω₁=10、ω₂=100Step02·低频段ω=1处幅值20lg0.1=-20dB,斜率由积分环节决定为-20dB/decStep03·斜率变化ω<10→-20dB/dec;10~100→-40;ω>100→-60dB/decStep04·相频叠加φ(ω)=-90°-arctan(0.1ω)-arctan(0.01ω),从-90°降至-270°,读出ωc、ωg计算PM和GMBodeDiagramAnalysis从开环伯德图推断闭环性能开环伯德图的低、中、高三个频段分别对应闭环系统的稳态精度、动态性能和抗噪能力,穿越频率ωc近似等于闭环带宽,相位裕度PM直接决定阶跃响应超调量。LOWFREQ低频段→稳态精度低频增益越高、积分环节数越多,闭环稳态误差越小系统型别v决定了可无差跟踪的输入类型(阶跃/斜坡/加速度)稳态误差essMIDFREQ中频段→动态性能穿越频率ωc近似等于闭环带宽,ωc越大响应越快、调节时间越短相位裕度PM与超调量直接相关:PM=45°时超调约25%,PM=60°时超调约10%相位裕度PMHIGHFREQ高频段→抗噪能力高频段衰减速率越快,系统对高频测量噪声的抑制能力越强但过快的衰减意味着带宽窄,需在响应速度与抗噪性之间权衡衰减斜率CHAPTER05工程应用与课程总结从理论分析到控制器设计的工程实践闭环BodePerspectivePID控制器的伯德图视角三个作用分别对应比例、积分和一阶微分环节的叠加,参数整定本质是调节低频增益、中频穿越频率和高频衰减。Kp比例增益系数ProportionalP比例作用对应比例环节Kp,在伯德图上整体抬高或降低幅频曲线增大Kp提高低频增益改善稳态精度,但穿越频率右移可能降低相位裕度IntegralI积分作用对应积分环节Ki/s,在低频段以-20dB/dec提升增益,消除稳态误差引入-90°恒定相位滞后,过量积分会使系统不稳定,需与D作用配合DerivativeD微分作用近似一阶微分特性,在中频段提供+20dB/dec斜率和正相位超前补偿积分环节的相位滞后,提高穿越频率处的相位裕度,改善动态响应FREQUENCYRESPONSE·COMPENSATION基于伯德图的频率法校正设计频率法校正通过改变伯德图形状满足性能指标:超前校正在中频段提供相位超前提高裕度和带宽,滞后校正提高低频增益改善稳态精度同时降低穿越频率,超前-滞后校正兼顾两者优势。超前校正传递函数Gc(s)=(aTs+1)/(Ts+1),a>1,本质为一阶微分与惯性环节的串联,在中高频段产生超前相位特性核心作用在穿越频率附近提供最大相位超前φm=arcsin((a−1)/(a+1)),有效提高系统相位裕度,扩展带宽加快响应适用场景稳态精度已满足但相位裕度不足、响应速度偏慢,或需要改善动态性能的控制系统滞后校正传递函数Gc(s)=(Ts+1)/(bTs+1),b>1,本质为惯性与一阶微分环节的串联,在低频段呈现积分特性核心作用低频增益提高b倍改善稳态误差,高频衰减使穿越频率左移至相位更充裕区域,牺牲带宽换取精度适用场景动态性能可接受但稳态误差不满足要求,或低频增益不足、需要抑制高频噪声的场合ControlToolboxMATLAB辅助伯德图绘制实践MATLAB的bode、margin等函数将伯德图绘制从手工渐近
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2026年城镇燃气安全管理人员试题库及答案
- 医学检验凝血知识测试题与答案
- 2026黑龙江伊春市乌翠区消防救援大队招聘消防文员4人备考题库
- 2026电工技师考试题库及答案
- 常见幼儿数学试题及答案呈现
- 四川省成都市高中化学 溶度积以及沉淀反应的应用教学设计 新人教版选修4
- 辽宁省沈阳市第二十一中学高中英语 Module 4 Fine Arts Western Chinese and Pop Period 1 IntroductionFunctionCultural Corner教案 外研版必修2
- 化水学习:填空试题及参考答案
- 广西玉林小升初考试试题及答案
- 2026年商业综合体建筑轮廓AI夜间勾勒设计
- 湖北邮政公司招聘笔试题目
- LY/T 1575-2023汽车车厢底板用竹胶合板
- 被执行人生活费申请书范文
- DZ-T 0270-2014地下水监测井建设规范
- KWFB无密封自控自吸泵(服务)课件
- 现代综合评价方法和案例-配套教材
- GB/T 4450-1995船用盲板钢法兰
- 新《义务教育法》解读
- 5.5M焦炉(1×55孔)施工方案
- 大断面方圆坯连铸机的设计特点课件
- 探伤报告格式模板
评论
0/150
提交评论