| 查看: 2030 | 回复: 14 | ||
z770428
|
[求助]
MATLAB 定义外部函数变量
|
|
各位高手，我用MATLAB求解热力耦合问题时，把温度场分析得到的热应变striant作为importdata输入到结构分析的workspace中运行结构分析主程序时，老是提示 ??? Undefined function or variable 'straint'. straint是一100*3矩阵， 我想问如何在结构分析主程序中定义straint？ |
回复此楼» 猜你喜欢
基金系统什么内容也没有
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab高手帮忙做个自定义复杂函数拟合，含有复数，
已经有25人回复
- 求matlab画一个函数的图
已经有6人回复
- MATLAB函数速查手册
已经有211人回复
- matlab如何改变函数变量值在调用函数重新计算
已经有12人回复
- matlab的函数调用问题
已经有3人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- matlab做函数图象的问题
已经有6人回复
2楼2012-10-19 09:14:24
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
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
30
投票:
