 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4999 | 回復(fù): 4 | ||
hustboshao銅蟲 (小有名氣)
|
[求助]
很簡(jiǎn)單的一個(gè)matlab程序運(yùn)行時(shí)間過長(zhǎng)，得不到結(jié)果 已有1人參與
|
|
主程序： clc clear all global kesai omega0 m rou D a1 a3 K C %%主系統(tǒng)基本參數(shù) global epsilong_T s_33 C_pzt C_pzt1 K_pzt K_pzt1 k0 d_33 s_0 l_0 theta L R %壓電材料參數(shù) kesai=0.0013;m=0.44;omega0=62.83;rou=1.2;D=0.015; a1=2.3;a3=-18; K=1736.947916; C=0.07187752; s_0=0.005*0.005*pi; %cross section area of the PZT rod l_0=0.01; %length of the PZT rod epsilong_T=3400*8.854e-12; % permittivity under constant stress 常應(yīng)力下的介電常數(shù) s_33=20.7e-12; %compliance of the piezoceramic under constant electric 恒定電場(chǎng)下的彈性柔順系數(shù) d_33=593e-12; %壓電常數(shù) C_pzt=epsilong_T*s_0/l_0; %capacitance between the electrodes of the PZT rod with no external force無外力下的電容 K_pzt=(1/s_33)*s_0/l_0; %stiffness of the short-circuited PZT rod 短路剛度 k0=d_33*(1/(sqrt(epsilong_T*s_33))); % electromechanical coupling cofficient 機(jī)電耦合系數(shù) %C_pzt1=C_pzt*(1-k0^2);%capacitance of the PZT rod under constant strain 常應(yīng)變下的電容 C_pzt1=1.2e-7; %K_pzt1=K_pzt*(1/(1-k0^2)); % stiffness of the PZT rod with open electrodes 回路剛度 %K_pzt1=C_pzt1/(epsilong_T*s_33*(1-k0^2)^2); K_pzt1=K*0.25; %theta=k0/(1-k0^2)*(1/(sqrt(epsilong_T*s_33)));% electromechanical coupling factor theta=1.55e-3; L=0.000001; % 電感 R=10; % 電阻 step_size=0.01; %時(shí)間步長(zhǎng) precision=1e-9; %精度量級(jí)1 options=odeset('reltol',precision); N=20; NN=N/0.01; tic U=5 %[T1,Y1]=ode45(@uncontrol,[0:step_size:N],[0.001 0],options,U); [T2,Y2]=ode45(@control,[0:step_size:N],[0.001 0 0.001 0],options,U); toc %w_1=max(Y1(NN/2:end,1))/0.015 w_2=max(Y2(NN/2:end,1)) w_2_Q=max(Y2(NN/2:end,3)) %figure(1) %plot(T1,Y1(:,1)/0.015,'r') %hold on %figure(2) %plot(T2,Y2(:,1)/0.015,'b') %hold on %figure(3) %plot(T2,Y2(:,3)/0.015,'G') %hold on figure(4) plot(T2,Y2(:,1)/0.015,'-r',T2,Y2(:,3)/0.015,'--b') xlabel('Time','FontName','Times New Roman','FontSize',20) ylabel('w_1,w_2','FontName','Times New Roman','FontSize',20) legend('uncontrol','control') hold on 子程序： function dy=control(t,y,U) global kesai omega0 m rou D a1 a3 K C %%主系統(tǒng)基本參數(shù) global epsilong_T s_33 C_pzt C_pzt1 K_pzt K_pzt1 k0 d_33 s_0 l_0 theta L R %壓電材料參數(shù) dy=[y(2); (rou*D*U*a1*y(2)*0.5+rou*D*a3*y(2)^3./U*0.5-C*y(2)-K*y(1))./m-K_pzt1*y(1)./m+theta*y(3)./m; y(4); -R*y(4)./L-y(3)./(L*C_pzt1)+theta*y(1)./L]; end |
回復(fù)此樓» 猜你喜歡
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有138人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有1人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 大家看看，我該如何投稿到哪個(gè)期刊？《從金庸的武學(xué)視角談?dòng)?jì)算機(jī)編程思維的培養(yǎng)》
已經(jīng)有10人回復(fù)
getengqing
木蟲 (正式寫手)
- 應(yīng)助: 104 (高中生)
- 金幣: 3759.8
- 散金: 30
- 紅花: 5
- 帖子: 592
- 在線: 110.5小時(shí)
- 蟲號(hào): 2986702
- 注冊(cè): 2014-02-23
- 性別: GG
- 專業(yè): 自然語言理解與機(jī)器翻譯
2樓2015-11-09 16:49:24
hustboshao
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1308.9
- 紅花: 1
- 帖子: 137
- 在線: 78小時(shí)
- 蟲號(hào): 1567118
- 注冊(cè): 2012-01-06
- 性別: GG
- 專業(yè): 固體力學(xué)
3樓2015-11-09 16:59:43
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
hustboshao: 金幣+10, ★★★很有幫助 2015-11-12 10:15:42
感謝參與，應(yīng)助指數(shù) +1
hustboshao: 金幣+10, ★★★很有幫助 2015-11-12 10:15:42
|
換成ode23s或者15s，沒仔細(xì)看你的程序，但是發(fā)現(xiàn)很多很小的數(shù)值，感覺你的問題應(yīng)該是剛性的，剛性方程不能用ode45解。 還有一個(gè)建議就是無量綱化，讓你要解的參數(shù)在 0 1 之間，或者至少在一個(gè)數(shù)量級(jí)上。 不過你要解的方程很少，如果只是用ode45時(shí)間長(zhǎng)，那換成ode23s/15s應(yīng)該會(huì)很快給你一個(gè)結(jié)果。 |
» 本帖已獲得的紅花（最新10朵）
hustboshao4樓2015-11-10 14:54:40
hustboshao
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1308.9
- 紅花: 1
- 帖子: 137
- 在線: 78小時(shí)
- 蟲號(hào): 1567118
- 注冊(cè): 2012-01-06
- 性別: GG
- 專業(yè): 固體力學(xué)
5樓2015-11-12 10:16:14
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hustboshao 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +4 | 大饞小子 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +10 | leaflight 2026-03-24 | 10/500 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +5 | 傾____蕭 2026-03-27 | 6/300 |
|
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4 | 12.27 2026-03-27 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
|
[考研] 一志愿華北電力大學(xué)能動(dòng)專碩，293，求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|