版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
強度計算.數(shù)值計算方法:非線性分析:非線性分析案例研究1非線性分析基礎(chǔ)1.1非線性分析概述非線性分析是結(jié)構(gòu)工程中一種重要的分析方法,用于解決那些在載荷與響應(yīng)之間存在非線性關(guān)系的問題。這種非線性關(guān)系可能源于材料的非線性行為、結(jié)構(gòu)的幾何非線性、接觸非線性,或是這些因素的組合。非線性分析能夠更準確地預測結(jié)構(gòu)在極端條件下的行為,如地震、爆炸或極端溫度變化等。1.2非線性材料特性1.2.1原理材料的非線性特性通常指材料在應(yīng)力超過一定閾值后,其應(yīng)變與應(yīng)力的關(guān)系不再遵循線性比例。常見的非線性材料模型包括彈塑性模型、粘彈性模型和超彈性模型。彈塑性模型描述了材料在塑性階段的應(yīng)力應(yīng)變關(guān)系,而粘彈性模型則考慮了時間對材料行為的影響。1.2.2內(nèi)容在非線性材料分析中,一個關(guān)鍵的步驟是定義材料的應(yīng)力-應(yīng)變曲線。例如,對于彈塑性材料,可以使用雙線性模型,其中材料在彈性階段遵循胡克定律,而在塑性階段,應(yīng)力保持不變,應(yīng)變繼續(xù)增加。示例代碼importnumpyasnp
importmatplotlib.pyplotasplt
#定義材料參數(shù)
E=200e9#彈性模量,單位:Pa
sigma_y=250e6#屈服強度,單位:Pa
epsilon_0=sigma_y/E#彈性應(yīng)變極限
#定義應(yīng)力-應(yīng)變關(guān)系
defstress_strain(epsilon):
ifepsilon<=epsilon_0:
sigma=E*epsilon
else:
sigma=sigma_y
returnsigma
#生成應(yīng)變數(shù)據(jù)
epsilon=np.linspace(0,0.005,100)
#計算應(yīng)力
sigma=[stress_strain(e)foreinepsilon]
#繪制應(yīng)力-應(yīng)變曲線
plt.plot(epsilon,sigma)
plt.xlabel('應(yīng)變')
plt.ylabel('應(yīng)力')
plt.title('彈塑性材料的應(yīng)力-應(yīng)變曲線')
plt.grid(True)
plt.show()1.3幾何非線性1.3.1原理幾何非線性分析考慮了結(jié)構(gòu)變形對自身幾何形狀的影響。當結(jié)構(gòu)的變形較大時,這種影響變得顯著,傳統(tǒng)的線性分析方法將不再適用。幾何非線性分析通常在大位移、大旋轉(zhuǎn)或薄殼結(jié)構(gòu)的分析中使用。1.3.2內(nèi)容在幾何非線性分析中,結(jié)構(gòu)的變形狀態(tài)需要在每一步迭代中更新。這涉及到使用非線性方程組來描述結(jié)構(gòu)的平衡狀態(tài),其中包含了位移、旋轉(zhuǎn)和應(yīng)力的非線性關(guān)系。示例代碼#假設(shè)一個簡單的梁在大位移下的非線性分析
importsympyassp
#定義變量
x,u=sp.symbols('xu')
#定義非線性位移方程
u_x=u+x**2/(1+u)
#定義外力
F=1000
#定義平衡方程
balance_eq=sp.diff(u_x,x)-F
#解平衡方程
u_solution=sp.dsolve(balance_eq,u)
#打印解
print(u_solution)1.4接觸非線性1.4.1原理接觸非線性分析處理的是兩個或多個物體在接觸時的相互作用。這種分析考慮了接觸面的摩擦、間隙、滑移等非線性效應(yīng),對于預測結(jié)構(gòu)在復雜載荷條件下的行為至關(guān)重要。1.4.2內(nèi)容接觸非線性分析通常需要定義接觸面的屬性,如摩擦系數(shù)、接觸剛度等。在分析過程中,需要檢查每個時間步或載荷步中接觸面的狀態(tài),以確定是否發(fā)生了接觸、滑移或分離。示例代碼#使用Python模擬兩個物體接觸的非線性分析
importnumpyasnp
#定義物體參數(shù)
mass1=1.0
mass2=1.0
k=1000#接觸剛度
c=0.1#摩擦系數(shù)
#定義初始條件
u1=0.0
u2=0.0
v1=1.0
v2=0.0
#定義時間步
dt=0.01
t_end=1.0
#定義時間數(shù)組
t=np.arange(0,t_end,dt)
#定義位移數(shù)組
u1_hist=np.zeros_like(t)
u2_hist=np.zeros_like(t)
#初始化位移
u1_hist[0]=u1
u2_hist[0]=u2
#進行時間積分
foriinrange(1,len(t)):
#計算接觸力
delta_u=u1_hist[i-1]-u2_hist[i-1]
ifdelta_u>0:
F_contact=k*delta_u
else:
F_contact=0
#計算摩擦力
ifv1>0andv2>0:
F_friction=c*(v1-v2)
else:
F_friction=0
#更新速度和位移
v1=v1-F_contact*dt/mass1
v2=v2+F_contact*dt/mass2
u1=u1+v1*dt
u2=u2+v2*dt
#存儲位移歷史
u1_hist[i]=u1
u2_hist[i]=u2
#繪制位移歷史
plt.plot(t,u1_hist,label='物體1')
plt.plot(t,u2_hist,label='物體2')
plt.xlabel('時間')
plt.ylabel('位移')
plt.legend()
plt.show()1.5非線性方程求解方法1.5.1原理非線性方程求解是非線性分析的核心。在非線性分析中,結(jié)構(gòu)的平衡狀態(tài)通常由一組非線性方程描述。這些方程可能非常復雜,無法通過解析方法求解,因此需要使用數(shù)值方法。1.5.2內(nèi)容常用的非線性方程求解方法包括牛頓-拉夫遜法、擬牛頓法和弧長法。牛頓-拉夫遜法是一種迭代方法,通過在當前點計算雅可比矩陣來線性化非線性方程,然后求解線性方程組來更新解。擬牛頓法是牛頓-拉夫遜法的一種變體,它避免了在每一步迭代中重新計算雅可比矩陣,從而提高了計算效率?;¢L法則是一種用于追蹤非線性分析中路徑依賴解的方法,特別適用于屈曲分析。示例代碼#使用牛頓-拉夫遜法求解非線性方程
importnumpyasnp
#定義非線性方程
deff(x):
returnx**3-2*x-5
#定義非線性方程的導數(shù)
defdf(x):
return3*x**2-2
#定義初始猜測
x0=2.0
#定義迭代次數(shù)
max_iter=100
#定義收斂準則
tol=1e-6
#進行迭代
foriinrange(max_iter):
x1=x0-f(x0)/df(x0)
ifabs(x1-x0)<tol:
break
x0=x1
#打印解
print("解為:",x1)以上代碼展示了如何使用牛頓-拉夫遜法求解一個非線性方程。通過迭代更新猜測值,直到滿足收斂準則。這種方法在非線性分析中非常常見,用于求解結(jié)構(gòu)的平衡方程。2非線性有限元分析2.1有限元法基礎(chǔ)回顧在非線性有限元分析中,我們首先回顧有限元法(FEM)的基礎(chǔ)。有限元法是一種數(shù)值方法,用于求解復雜的工程問題,如結(jié)構(gòu)力學、熱傳導、流體力學等。它將連續(xù)體離散化為有限數(shù)量的單元,每個單元用一組節(jié)點來表示,通過在這些節(jié)點上求解偏微分方程的近似解,來獲得整個結(jié)構(gòu)的解。2.1.1原理有限元法基于變分原理和加權(quán)殘值法。在結(jié)構(gòu)力學中,它通常用于求解彈性力學方程。對于一個給定的結(jié)構(gòu),有限元法通過將結(jié)構(gòu)分解為多個小的、簡單的單元,然后在每個單元上應(yīng)用局部的平衡方程和變形協(xié)調(diào)條件,來構(gòu)建整個結(jié)構(gòu)的全局方程。2.1.2代碼示例以下是一個使用Python和scipy庫進行線性有限元分析的簡單示例,用于求解一個簡單的梁的彎曲問題。雖然這里展示的是線性分析,但非線性分析的流程與此類似,只是在求解過程中需要考慮非線性效應(yīng)。importnumpyasnp
fromscipy.sparseimportlil_matrix
fromscipy.sparse.linalgimportspsolve
#定義梁的長度和節(jié)點數(shù)
L=1.0
n=10
h=L/(n-1)
#創(chuàng)建剛度矩陣
K=lil_matrix((n,n),dtype=float)
foriinrange(n-1):
K[i,i]+=12
K[i,i+1]-=6
K[i+1,i]-=6
K[i+1,i+1]+=4
ifi<n-2:
K[i+1,i+2]-=3
K[i+2,i+1]-=3
K[i,i]+=6/h**3
K[i,i+1]-=3/h**3
K[i+1,i]-=3/h**3
K[i+1,i+1]+=4/h**3
ifi<n-2:
K[i+1,i+2]-=1/h**3
K[i+2,i+1]-=1/h**3
#應(yīng)用邊界條件
K[0,:]=0
K[-1,:]=0
K[0,0]=1
K[-1,-1]=1
#定義載荷向量
F=np.zeros(n)
F[n//2]=-1.0
#求解位移向量
U=spsolve(K.tocsr(),F)
#輸出位移向量
print(U)2.2非線性有限元模型建立非線性有限元分析考慮了材料非線性、幾何非線性和接觸非線性等因素。在建立非線性有限元模型時,需要定義材料屬性、幾何形狀、載荷和邊界條件,以及非線性效應(yīng)的類型。2.2.1原理非線性效應(yīng)可能包括材料的塑性、大變形、接觸、摩擦等。這些效應(yīng)使得結(jié)構(gòu)的響應(yīng)不再是載荷的線性函數(shù),因此需要使用迭代方法來求解非線性方程組。2.2.2代碼示例使用Python和FEniCS庫建立一個非線性有限元模型的示例。這里我們考慮一個簡單的平面應(yīng)力問題,其中材料表現(xiàn)出塑性行為。fromfenicsimport*
importmatplotlib.pyplotasplt
#創(chuàng)建網(wǎng)格
mesh=UnitSquareMesh(8,8)
#定義函數(shù)空間
V=VectorFunctionSpace(mesh,'Lagrange',2)
#定義邊界條件
defboundary(x,on_boundary):
returnon_boundary
bc=DirichletBC(V,Constant((0,0)),boundary)
#定義材料屬性
E=1e3
nu=0.3
yield_stress=100
#定義本構(gòu)關(guān)系
defsigma(F):
I=Identity(F.shape[0])
J=det(F)
ifJ>1:
returnE/(1+nu)*(F-inv(F.T))*J**(-1/3)+E*nu/((1+nu)*(1-2*nu))*(J**(1/3)-1)*I
else:
returnE/(1+nu)*(F-inv(F.T))*J**(-1/3)
#定義變分問題
u=TrialFunction(V)
v=TestFunction(V)
F=Identity(2)+grad(u)
P=Constant((0,-1))
T=inner(sigma(F),grad(v))*dx-inner(P,v)*ds
#定義非線性問題
problem=NonlinearVariationalProblem(T,Function(V),bc,J=T)
#定義求解器
solver=NonlinearVariationalSolver(problem)
#設(shè)置求解參數(shù)
parameters=solver.parameters
parameters['newton_solver']['relative_tolerance']=1e-6
parameters['newton_solver']['absolute_tolerance']=1e-6
parameters['newton_solver']['maximum_iterations']=20
#求解
solver.solve()
#可視化結(jié)果
u=Function(V)
plot(u)
plt.show()2.3非線性有限元網(wǎng)格劃分網(wǎng)格劃分是有限元分析中的關(guān)鍵步驟,它直接影響到分析的精度和計算效率。在非線性分析中,網(wǎng)格的細化和自適應(yīng)劃分尤為重要,因為非線性效應(yīng)可能在局部區(qū)域產(chǎn)生顯著影響。2.3.1原理網(wǎng)格劃分應(yīng)考慮結(jié)構(gòu)的幾何特征和載荷分布。對于非線性分析,通常在高應(yīng)力或應(yīng)變區(qū)域使用更細的網(wǎng)格,以捕捉非線性行為的細節(jié)。自適應(yīng)網(wǎng)格劃分技術(shù)可以根據(jù)解的局部誤差自動調(diào)整網(wǎng)格密度。2.3.2代碼示例使用FEniCS庫進行自適應(yīng)網(wǎng)格劃分的示例。fromfenicsimport*
importmatplotlib.pyplotasplt
#創(chuàng)建初始網(wǎng)格
mesh=UnitSquareMesh(8,8)
#定義函數(shù)空間
V=FunctionSpace(mesh,'Lagrange',2)
#定義邊界條件
defboundary(x,on_boundary):
returnon_boundary
bc=DirichletBC(V,Constant(0),boundary)
#定義變分問題
u=TrialFunction(V)
v=TestFunction(V)
f=Constant(-6)
a=dot(grad(u),grad(v))*dx
L=f*v*dx
#求解
u=Function(V)
solve(a==L,u,bc)
#自適應(yīng)網(wǎng)格劃分
error_control=AdaptiveMeshRefinementControl()
error_control.mark(mesh,u.vector().get_local(),0.5)
#重復求解和自適應(yīng)劃分
foriinrange(5):
mesh=refine(mesh,error_control)
V=FunctionSpace(mesh,'Lagrange',2)
bc=DirichletBC(V,Constant(0),boundary)
u=Function(V)
solve(a==L,u,bc)
#可視化最終網(wǎng)格
plot(mesh)
plt.show()2.4非線性有限元邊界條件設(shè)置邊界條件在有限元分析中至關(guān)重要,它們定義了結(jié)構(gòu)的約束和載荷。在非線性分析中,邊界條件可能隨載荷的增加而變化,例如接觸邊界條件。2.4.1原理邊界條件包括位移邊界條件和載荷邊界條件。位移邊界條件用于固定結(jié)構(gòu)的某些部分,而載荷邊界條件則用于施加外部力。在非線性分析中,接觸邊界條件需要特別處理,以確保接觸面的正確行為。2.4.2代碼示例使用FEniCS庫設(shè)置接觸邊界條件的示例。fromfenicsimport*
importmatplotlib.pyplotasplt
#創(chuàng)建網(wǎng)格
mesh=UnitSquareMesh(8,8)
#定義函數(shù)空間
V=VectorFunctionSpace(mesh,'Lagrange',2)
#定義邊界條件
defboundary(x,on_boundary):
returnon_boundary
bc=DirichletBC(V,Constant((0,0)),boundary)
#定義接觸邊界條件
defcontact_boundary(x,on_boundary):
returnnear(x[1],0)andon_boundary
contact_bc=ContactCondition(V,Constant(0),contact_boundary)
#定義變分問題
u=TrialFunction(V)
v=TestFunction(V)
F=inner(sigma(u),grad(v))*dx-inner(Constant((0,-1)),v)*ds
#定義非線性問題
problem=NonlinearVariationalProblem(F,Function(V),bc,J=derivative(F,u))
#設(shè)置接觸邊界條件
problem.set_contact_conditions([contact_bc])
#定義求解器
solver=NonlinearVariationalSolver(problem)
#設(shè)置求解參數(shù)
parameters=solver.parameters
parameters['newton_solver']['relative_tolerance']=1e-6
parameters['newton_solver']['absolute_tolerance']=1e-6
parameters['newton_solver']['maximum_iterations']=20
#求解
solver.solve()
#可視化結(jié)果
plot(Function(V))
plt.show()2.5非線性有限元求解控制非線性有限元分析通常需要迭代求解,直到滿足收斂準則。求解控制包括選擇合適的求解器、設(shè)置收斂準則和迭代策略。2.5.1原理求解控制策略應(yīng)考慮非線性問題的性質(zhì)。對于材料非線性,可能需要使用增量加載和增量迭代。對于幾何非線性,可能需要使用弧長法或Riks法。收斂準則通?;跉埐詈臀灰频母淖兞俊?.5.2代碼示例使用FEniCS庫進行增量迭代求解的示例。fromfenicsimport*
importmatplotlib.pyplotasplt
#創(chuàng)建網(wǎng)格
mesh=UnitSquareMesh(8,8)
#定義函數(shù)空間
V=VectorFunctionSpace(mesh,'Lagrange',2)
#定義邊界條件
defboundary(x,on_boundary):
returnon_boundary
bc=DirichletBC(V,Constant((0,0)),boundary)
#定義變分問題
u=TrialFunction(V)
v=TestFunction(V)
F=inner(sigma(u),grad(v))*dx-inner(Constant((0,-1)),v)*ds
#定義非線性問題
problem=NonlinearVariationalProblem(F,Function(V),bc,J=derivative(F,u))
#定義求解器
solver=NonlinearVariationalSolver(problem)
#設(shè)置求解參數(shù)
parameters=solver.parameters
parameters['newton_solver']['relative_tolerance']=1e-6
parameters['newton_solver']['absolute_tolerance']=1e-6
parameters['newton_solver']['maximum_iterations']=20
#增量迭代求解
u=Function(V)
foriinrange(10):
solver.solve()
u.vector().axpy(1,Function(V).vector())
#可視化結(jié)果
plot(u)
plt.show()以上示例和原理概述了非線性有限元分析的基本流程和關(guān)鍵步驟。在實際應(yīng)用中,這些步驟可能需要根據(jù)具體問題進行調(diào)整和優(yōu)化。3非線性分析案例研究3.1橋梁結(jié)構(gòu)非線性分析3.1.1原理橋梁結(jié)構(gòu)的非線性分析主要考慮材料非線性、幾何非線性以及邊界條件非線性。在實際應(yīng)用中,橋梁可能遭受大荷載、地震等極端條件,導致結(jié)構(gòu)發(fā)生非線性變形。非線性分析能夠更準確地預測結(jié)構(gòu)在這些條件下的行為,評估其安全性和穩(wěn)定性。3.1.2內(nèi)容材料非線性:考慮混凝土、鋼材等材料的塑性、彈塑性或損傷模型。幾何非線性:處理大變形和大位移,如懸索橋的纜索松弛。邊界條件非線性:分析支座、連接件等的非線性行為。3.1.3示例假設(shè)我們有一個簡化的橋梁模型,使用Python和FEniCS庫進行非線性分析。fromfenicsimport*
importmatplotlib.pyplotasplt
#創(chuàng)建網(wǎng)格和函數(shù)空間
mesh=UnitSquareMesh(10,10)
V=VectorFunctionSpace(mesh,'Lagrange',2)
#定義邊界條件
defboundary(x,on_boundary):
returnon_boundary
bc=DirichletBC(V,Constant((0,0)),boundary)
#定義材料參數(shù)
E=1e3
nu=0.3
mu=E/(2*(1+nu))
lmbda=E*nu/((1+nu)*(1-2*nu))
#定義非線性材料模型
defsigma(F):
returnlmbda*(tr(F)-3)*Identity(2)+2*mu*(F-Identity(2))
#定義變分問題
u=TrialFunction(V)
v=TestFunction(V)
f=Constant((0,-10))
T=Constant((1,0))
F=inner(sigma(grad(u)),grad(v))*dx-inner(f,v)*dx-inner(T,v)*ds
#解非線性問題
u=Function(V)
solve(F==0,u,bc)
#可視化結(jié)果
plot(u)
plt.show()此代碼示例展示了如何使用FEniCS庫解決一個非線性彈性問題,模擬橋梁結(jié)構(gòu)在荷載作用下的變形。3.2高層建筑非線性地震響應(yīng)3.2.1原理高層建筑在地震作用下可能經(jīng)歷非線性變形,特別是當?shù)卣饛姸瘸^設(shè)計基準時。非線性地震響應(yīng)分析考慮結(jié)構(gòu)的非線性特性,如材料的塑性變形、結(jié)構(gòu)的幾何非線性以及阻尼的非線性。3.2.2內(nèi)容地震荷載模擬:使用時程分析或反應(yīng)譜分析。結(jié)構(gòu)非線性模型:包括材料非線性、幾何非線性以及阻尼非線性。結(jié)果評估:分析結(jié)構(gòu)的位移、加速度、內(nèi)力等,評估結(jié)構(gòu)的損傷程度。3.2.3示例使用Python和PySDOF庫進行高層建筑的非線性地震響應(yīng)分析。importnumpyasnp
frompysdofimportSdofSystem
#定義單自由度系統(tǒng)參數(shù)
mass=10000#質(zhì)量,單位:kg
stiffness=1e6#剛度,單位:N/m
damping=1000#阻尼,單位:Ns/m
yield_force=1.5e6#屈服力,單位:N
#創(chuàng)建非線性單自由度系統(tǒng)
sys=SdofSystem(mass,stiffness,damping,yield_force=yield_force)
#定義地震加速度時程
time=np.linspace(0,10,1000)
acc=np.sin(2*np.pi*time)#簡化示例,實際應(yīng)用中應(yīng)使用地震記錄
#進行時程分析
response=sys.time_history(acc,time)
#輸出結(jié)果
print(response['displacement'])
print(response['velocity'])
print(response['acceleration'])此代碼示例展示了如何使用PySDOF庫模擬一個高層建筑在地震作用下的非線性響應(yīng),包括位移、速度和加速度。3.3復合材料非線性破壞分析3.3.1原理復合材料的非線性破壞分析通常涉及損傷力學模型,考慮材料的多軸應(yīng)力狀態(tài)和損傷累積。分析的目標是預測復合材料在復雜載荷下的破壞模式和壽命。3.3.2內(nèi)容損傷模型:如最大應(yīng)力理論、最大應(yīng)變理論、Tsai-Wu理論等。非線性本構(gòu)關(guān)系:描述復合材料在不同應(yīng)力狀態(tài)下的行為。壽命預測:基于損傷累積理論預測復合材料的剩余壽命。3.3.3示例使用Python和SciPy庫進行復合材料的非線性破壞分析。importnumpyasnp
fromegrateimportodeint
#定義損傷累積
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2024排水工程施工合同范本
- 親子春游主題活動方案模板五篇
- 2024年環(huán)保型木材料購銷協(xié)議
- 二零二五年度基礎(chǔ)設(shè)施建設(shè)擔保合同樣本3篇
- 原發(fā)性閉角型青光眼護理查房
- 2024年電商快遞代發(fā)服務(wù)協(xié)議
- (四川樂山卷)2022年中考物理第二次模擬考試(參考答案)
- 露在外面的面(說課稿)-2023-2024學年五年級下冊數(shù)學北師大版001
- 一年級數(shù)學100以內(nèi)口算題
- 2024版地暖安裝施工合同協(xié)議書范本
- 農(nóng)業(yè)機械控制系統(tǒng)硬件在環(huán)測試規(guī)范
- 翁潭電站大王山輸水隧洞施工控制網(wǎng)設(shè)計說明書
- 隆胸術(shù)培訓課件
- 鋼筋焊接培訓課件
- 行政內(nèi)勤培訓課件
- 酒吧、夜店、KTV員工管理手冊(2022版)
- 化纖企業(yè)(化學纖維紡織企業(yè))安全生產(chǎn)操作規(guī)程
- 重大事故隱患專項排查檢查表
- 中建住宅二次結(jié)構(gòu)專項施工方案
- 紅薯淀粉加工項目可行性研究報告
- 《工程圖學基礎(chǔ)教程(第4版)》 課件 第7章 零件圖
評論
0/150
提交評論