| 查看: 1995 | 回复: 14 | |||
z770428
|
[求助]
MATLAB 定义外部函数变量
|
|
各位高手，我用MATLAB求解热力耦合问题时，把温度场分析得到的热应变striant作为importdata输入到结构分析的workspace中运行结构分析主程序时，老是提示 ??? Undefined function or variable 'straint'. straint是一100*3矩阵， 我想问如何在结构分析主程序中定义straint？ |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有43人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab高手帮忙做个自定义复杂函数拟合，含有复数，
已经有25人回复
- 求matlab画一个函数的图
已经有6人回复
- MATLAB函数速查手册
已经有211人回复
- matlab如何改变函数变量值在调用函数重新计算
已经有12人回复
- matlab的函数调用问题
已经有3人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- matlab做函数图象的问题
已经有6人回复
wgm096350
银虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1229.2
- 散金: 120
- 红花: 2
- 帖子: 565
- 在线: 72.1小时
- 虫号: 2071589
- 注册: 2012-10-19
- 性别: GG
- 专业: 固体力学
2楼2012-10-19 09:14:24
wgm096350
银虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1229.2
- 散金: 120
- 红花: 2
- 帖子: 565
- 在线: 72.1小时
- 虫号: 2071589
- 注册: 2012-10-19
- 性别: GG
- 专业: 固体力学
3楼2012-10-19 09:14:42
z770428 
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1508.1
- 红花: 2
- 帖子: 416
- 在线: 178.3小时
- 虫号: 1360748
- 注册: 2011-08-04
- 性别: GG
- 专业: 固体力学
4楼2012-10-20 11:15:20
zhizi888888
银虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 5607.9
- 帖子: 17
- 在线: 58.6小时
- 虫号: 1774648
- 注册: 2012-04-23
- 性别: GG
- 专业: 半导体电子器件
5楼2012-10-20 12:23:21
z770428
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1508.1
- 红花: 2
- 帖子: 416
- 在线: 178.3小时
- 虫号: 1360748
- 注册: 2011-08-04
- 性别: GG
- 专业: 固体力学
6楼2012-10-20 12:52:18
z770428
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1508.1
- 红花: 2
- 帖子: 416
- 在线: 178.3小时
- 虫号: 1360748
- 注册: 2011-08-04
- 性别: GG
- 专业: 固体力学
7楼2012-10-21 08:57:35
z770428
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1508.1
- 红花: 2
- 帖子: 416
- 在线: 178.3小时
- 虫号: 1360748
- 注册: 2011-08-04
- 性别: GG
- 专业: 固体力学
8楼2012-10-21 08:57:52
wgm096350
银虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1229.2
- 散金: 120
- 红花: 2
- 帖子: 565
- 在线: 72.1小时
- 虫号: 2071589
- 注册: 2012-10-19
- 性别: GG
- 专业: 固体力学
9楼2012-10-21 15:48:30
z770428
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1508.1
- 红花: 2
- 帖子: 416
- 在线: 178.3小时
- 虫号: 1360748
- 注册: 2011-08-04
- 性别: GG
- 专业: 固体力学
xiegangmai: 帖中代码部分可使用“插入代码”功能 2012-10-23 22:03:05
|
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 |
10楼2012-10-21 21:17:37
