| 查看: 2009 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
z770428
|
[求助]
MATLAB 定义外部函数变量
|
||
|
各位高手，我用MATLAB求解热力耦合问题时，把温度场分析得到的热应变striant作为importdata输入到结构分析的workspace中运行结构分析主程序时，老是提示 ??? Undefined function or variable 'straint'. straint是一100*3矩阵， 我想问如何在结构分析主程序中定义straint？ |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有6人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab高手帮忙做个自定义复杂函数拟合，含有复数，
已经有25人回复
- 求matlab画一个函数的图
已经有6人回复
- MATLAB函数速查手册
已经有211人回复
- matlab如何改变函数变量值在调用函数重新计算
已经有12人回复
- matlab的函数调用问题
已经有3人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- matlab做函数图象的问题
已经有6人回复
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
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
