z770428

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1508.1
紅花: 2
帖子: 416
在線: 178.3小時
蟲號: 1360748
注冊: 2011-08-04
性別: GG
專業(yè): 固體力學(xué)

[求助] MATLAB 定義外部函數(shù)變量

各位高手，我用MATLAB求解熱力耦合問題時，把溫度場分析得到的熱應(yīng)變striant作為importdata輸入到結(jié)構(gòu)分析的workspace中運(yùn)行結(jié)構(gòu)分析主程序時，老是提示
??? Undefined function or variable 'straint'.
straint是一100*3矩陣，

我想問如何在結(jié)構(gòu)分析主程序中定義straint？

回復(fù)此樓

» 猜你喜歡

283分材料與化工求調(diào)劑已經(jīng)有12人回復(fù)
349求調(diào)劑已經(jīng)有7人回復(fù)
296材料專碩求調(diào)劑已經(jīng)有11人回復(fù)
一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè) 已經(jīng)有7人回復(fù)
材料340分調(diào)劑已經(jīng)有8人回復(fù)
316求調(diào)劑已經(jīng)有13人回復(fù)
321求調(diào)劑已經(jīng)有11人回復(fù)
315求調(diào)劑已經(jīng)有6人回復(fù)
330分求調(diào)劑已經(jīng)有13人回復(fù)
化工調(diào)劑303分，過四級已經(jīng)有9人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab高手幫忙做個自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，已經(jīng)有25人回復(fù)
求matlab畫一個函數(shù)的圖已經(jīng)有6人回復(fù)
MATLAB函數(shù)速查手冊已經(jīng)有211人回復(fù)
matlab如何改變函數(shù)變量值在調(diào)用函數(shù)重新計算已經(jīng)有12人回復(fù)
matlab的函數(shù)調(diào)用問題已經(jīng)有3人回復(fù)
已知一隨機(jī)信號的功率譜密度函數(shù)，如何利用Matlab產(chǎn)生其單次時域樣本已經(jīng)有6人回復(fù)
問個關(guān)于matlab里find函數(shù)的問題！已經(jīng)有4人回復(fù)
matlab 計算過程中工作空間變量保存及讀取求解決內(nèi)存不足的方法已經(jīng)有5人回復(fù)
matlab做函數(shù)圖象的問題已經(jīng)有6人回復(fù)

1樓 2012-10-18 17:51:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wgm096350

銀蟲 (正式寫手)

應(yīng)助: 6 (幼兒園)
金幣: 1228.4
散金: 120
紅花: 2
帖子: 565
在線: 72.1小時
蟲號: 2071589
注冊: 2012-10-19
性別: GG
專業(yè): 固體力學(xué)

發(fā)圖

回復(fù)此樓

愿有前程可奔赴，亦有歲月可回首

2樓2012-10-19 09:14:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wgm096350

銀蟲 (正式寫手)

應(yīng)助: 6 (幼兒園)
金幣: 1228.4
散金: 120
紅花: 2
帖子: 565
在線: 72.1小時
蟲號: 2071589
注冊: 2012-10-19
性別: GG
專業(yè): 固體力學(xué)

請發(fā)圖

回復(fù)此樓

愿有前程可奔赴，亦有歲月可回首

3樓2012-10-19 09:14:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

z770428

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1508.1
紅花: 2
帖子: 416
在線: 178.3小時
蟲號: 1360748
注冊: 2011-08-04
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

3樓: Originally posted by wgm096350 at 2012-10-19 09:14:42
請發(fā)圖

您好，能幫幫我么

回復(fù)此樓

4樓2012-10-20 11:15:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhizi888888

銀蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 5607.9
帖子: 17
在線: 58.6小時
蟲號: 1774648
注冊: 2012-04-23
性別: GG
專業(yè): 半導(dǎo)體電子器件

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝參與 2012-10-21 20:06:49

不知道對不對啊，
Global
straint

贊一下

回復(fù)此樓

踏踏實(shí)實(shí)做人，認(rèn)認(rèn)真真做事

5樓2012-10-20 12:23:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

z770428

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1508.1
紅花: 2
帖子: 416
在線: 178.3小時
蟲號: 1360748
注冊: 2011-08-04
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

5樓: Originally posted by zhizi888888 at 2012-10-20 12:23:21
不知道對不對啊，
Global
straint

Global
straint

這個我試了，不行，在WORKSPACE沒結(jié)果

贊一下

回復(fù)此樓

6樓2012-10-20 12:52:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

z770428

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1508.1
紅花: 2
帖子: 416
在線: 178.3小時
蟲號: 1360748
注冊: 2011-08-04
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

6樓: Originally posted by z770428 at 2012-10-20 12:52:18
Global
straint

這個我試了，不行，在WORKSPACE沒結(jié)果...

??? Attempted to access straint(1); index out of bounds because numel(straint)=0.

贊一下

回復(fù)此樓

7樓2012-10-21 08:57:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

z770428

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1508.1
紅花: 2
帖子: 416
在線: 178.3小時
蟲號: 1360748
注冊: 2011-08-04
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

5樓: Originally posted by zhizi888888 at 2012-10-20 12:23:21
不知道對不對啊，
Global
straint

??? Attempted to access straint(1); index out of bounds because numel(straint)=0.

贊一下

回復(fù)此樓

8樓2012-10-21 08:57:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wgm096350

銀蟲 (正式寫手)

應(yīng)助: 6 (幼兒園)
金幣: 1228.4
散金: 120
紅花: 2
帖子: 565
在線: 72.1小時
蟲號: 2071589
注冊: 2012-10-19
性別: GG
專業(yè): 固體力學(xué)

能把程序發(fā)一下不？？？這樣子也不知道到底錯在哪里啊

贊一下

回復(fù)此樓

愿有前程可奔赴，亦有歲月可回首

9樓2012-10-21 15:48:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

z770428

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1508.1
紅花: 2
帖子: 416
在線: 178.3小時
蟲號: 1360748
注冊: 2011-08-04
性別: GG
專業(yè): 固體力學(xué)

xiegangmai: 帖中代碼部分可使用“插入代碼”功能 2012-10-23 22:03:05

引用回帖:

9樓: Originally posted by wgm096350 at 2012-10-21 15:48:30
能把程序發(fā)一下不？？？這樣子也不知道到底錯在哪里啊

clear all; close all; clc; format compact; tic; global CRACK NODES straint ;

%%%% Pre-Processing %%%%%
inputQuasiStatic;                                                          % Define the geometry, materials, discontinuities

iter = numIterations;
for i = 1:iter

%%%%% Processing %%%%%
if i == 1, connectivity; pHDOF = []; else NODES(:,4:29) = 0; end       % Define connectivity
omega          = levelSet(i);                                        % Create phi and psi, define enriched elements

straint=?
[DOF,DISP] = calcDOF;                                                 % Total degrees of freedom
[updElem] = enrElem(i,pHDOF);                                        % Find enriched elements, inclusion elements

if i == 1,
      globalK = displacementstiffnessMatrix(omega,DOF,iter,updElem);
                                                                           % Construct global stiffness matrix  %

else globalK = updatedisplacementStiffness(globalK,omega,DOF,updElem,pHDOF);
end                                                                   % Update the global stiffness matrix
freeDOF       = boundaryCond(DOF);
globalF       = forceVector(DOF,i,straint);                         % Construct global force vector
                                                                              % Solve for the degrees of freedom

DISP(freeDOF,

= globalK(freeDOF,freeDOF)\globalF(freeDOF,

;       % Find the nodal displacement values

%%% Post-Processing %%%%%
if i == iter, plotMain(omega,DISP); end                               % Make plots
if isempty(CRACK) == 0
      pHDOF = 2*max(NODES(:,2));                                     % Maximum constant DOF at current iteration
      [KI,KII] = JIntegral(omega,DISP);                                  % Calculate the stress intensity factors
      exit    = growCrack(KI,KII,omega);                               % Advance crack for quasi-static growth
      globalF       = forceVector(DOF,i,straint);
      if strcmp(exit,'YES') == 1
         disp('WARNING: No crack growth, iterations exited early.')
         plotMain(omega,DISP); break
      end
end

disp(['Iteration ',num2str(i),' completed. Elapsed time is ',num2str(toc,'%0.4f'),'.'])
end

贊一下

回復(fù)此樓

10樓2012-10-21 21:17:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 318求調(diào)劑 +3	篤行致遠(yuǎn). 2026-03-31	4/200	2026-04-02 15:56 by Jaylen.
[考研] 071000生物學(xué)調(diào)劑 +6	知昭蔓 2026-04-02	6/300	2026-04-02 15:28 by sophie2180
[考研] 348求調(diào)劑 +6	吳彥祖24k 2026-04-02	6/300	2026-04-02 14:07 by 給你你注意休息
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +6	1段時間 2026-03-31	7/350	2026-04-02 13:57 by 3041
[考研] 314求調(diào)劑 +11	1xiaojun23 2026-03-31	12/600	2026-04-02 12:31 by 1xiaojun23
[考研] 311求調(diào)劑 +14	藍(lán)月亮亮 2026-03-30	14/700	2026-04-02 12:18 by 1753564080
[考研] 377求調(diào)劑 +3	RASKIN 2026-04-02	3/150	2026-04-02 09:45 by zzchen2000
[考研] 367求調(diào)劑 +8	芋泥啵�！� 2026-03-28	8/400	2026-04-02 09:20 by olim
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學(xué) +8	圣日耳曼條 2026-04-01	11/550	2026-04-02 07:40 by chemdavid
[考研] 求調(diào)劑0703 +5	周嘉堯 2026-03-31	8/400	2026-04-01 20:32 by ltltkkk
[考研] 材料科學(xué)與工程339求調(diào)劑 +11	hyz0119 2026-03-31	12/600	2026-04-01 18:40 by 偉大河北
[考研] 070300求調(diào)劑306分 +5	26要上岸 2026-03-27	5/250	2026-04-01 11:09 by oooqiao
[考研] 322求調(diào)劑 +8	三水sss 2026-04-01	8/400	2026-04-01 10:19 by 唐沐兒
[考研] 375求調(diào)劑 +7	雨夏整夜 2026-03-29	7/350	2026-03-31 18:52 by xhai2011
[考研] 生物考研337分求調(diào)劑 +4	cgxin 2026-03-30	6/300	2026-03-31 14:18 by 記事本2026
[考研] 266分，求材料冶金能源化工等調(diào)劑 +8	哇呼哼呼哼 2026-03-27	10/500	2026-03-31 13:35 by Huaxue_Wang
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3	小圓圈圈ooo 2026-03-30	3/150	2026-03-30 21:05 by 余震yz
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無際的草原
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790