 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5001 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
hustboshao銅蟲 (小有名氣)
|
[求助]
很簡單的一個matlab程序運行時間過長，得不到結(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 恒定電場下的彈性柔順系數(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 機電耦合系數(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; %時間步長 precision=1e-9; %精度量級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ù)此樓» 猜你喜歡
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有241人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有29人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有10人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 大家看看，我該如何投稿到哪個期刊？《從金庸的武學(xué)視角談計算機編程思維的培養(yǎng)》
已經(jīng)有10人回復(fù)
【答案】應(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解。 還有一個建議就是無量綱化，讓你要解的參數(shù)在 0 1 之間，或者至少在一個數(shù)量級上。 不過你要解的方程很少，如果只是用ode45時間長，那換成ode23s/15s應(yīng)該會很快給你一個結(jié)果。 |
hustboshao4樓2015-11-10 14:54:40
getengqing
木蟲 (正式寫手)
- 應(yīng)助: 104 (高中生)
- 金幣: 3759.8
- 散金: 30
- 紅花: 5
- 帖子: 592
- 在線: 110.5小時
- 蟲號: 2986702
- 注冊: 2014-02-23
- 性別: GG
- 專業(yè): 自然語言理解與機器翻譯
2樓2015-11-09 16:49:24
hustboshao
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1308.9
- 紅花: 1
- 帖子: 137
- 在線: 78小時
- 蟲號: 1567118
- 注冊: 2012-01-06
- 性別: GG
- 專業(yè): 固體力學(xué)
3樓2015-11-09 16:59:43
hustboshao
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1308.9
- 紅花: 1
- 帖子: 137
- 在線: 78小時
- 蟲號: 1567118
- 注冊: 2012-01-06
- 性別: GG
- 專業(yè): 固體力學(xué)
5樓2015-11-12 10:16:14
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 材料工程專業(yè)碩士申博 +3 | 麟正宇 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑求調(diào)劑 +5 | 121. 2026-04-02 | 5/250 |
|
|
[考研] 318求調(diào)劑，計算材料方向 +8 | 吸喵有害笙命 2026-04-01 | 9/450 |
|
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +20 | Youth_ 2026-03-30 | 20/1000 |
|
|
|
[考研] 085600 295分求調(diào)劑 +16 | W55j 2026-03-30 | 19/950 |
|
|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +9 | 雙馬尾痞老板2 2026-04-01 | 9/450 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +7 | Ong3 2026-03-27 | 7/350 |
|
|
|
[考研] 材料專碩322分 +11 | 哈哈哈吼吼吼哈 2026-04-01 | 11/550 |
|
|
|
[考研] 324求調(diào)劑 +5 | 想上學(xué)求調(diào) 2026-04-01 | 6/300 |
|
|
|
[考研] 英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑 +12 | 12.27 2026-03-30 | 14/700 |
|
|
|
[考研] 求調(diào)劑 +5 | 零八# 2026-03-27 | 5/250 |
|
|
|
[考研] 0855機械初試280求調(diào)劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
|
[考研] 江蘇蘇北高校誠邀調(diào)劑同學(xué) +3 | zzll406 2026-03-31 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿大連理工大學(xué)，機械工程學(xué)碩，341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 材料專碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|