【物理应用】Matlab实现两端固支梁热力耦合的有限元分析

1 简介

Matlab实现两端固支梁热力耦合的有限元分析

2 部分代码

% -------------------------------------------------------------------------%                  BASIC PARAMETERS (BRIDGE-LIKE STRUCTURE)    % -------------------------------------------------------------------------clear ; clc ;HoriTotal  = 100 ;          % horizontalVeriTotal  = 70;            % verticalIniVolfrac = 0.5  ;         % initial volume fractionVolfrac    = 0.25 ;         % allowable volume fractionvarimin    = 1e-3 ;         % smallest value of design variableObjScale   = 14e-3 ;        rmin       = 2 ;            % filter radiusF_uni      = 35 ;           % distributed loadqnVari     = 0.3 ;          % STM applied to design variablesRE         = 28 ;           % RAMP for stiffnessRbt        = 16 ;           % RAMP for thermal stress coefficientTalafa0   = 12.1e-6 ;       % thermal expansion coefficientTwoDeVec  = [ 1 1 0 ];    % -------------------------------------------------------------------------%                         DISCRETIZATION FEATURES % -------------------------------------------------------------------------%-ELEMENT STIFFNESS-a = 0.5 ;                                              b = 0.5;                                              h = 1 ;                    % element thicknessv0 = (2*a) * (2*b) * h ;   % volume of solid elementE0 = 2.1e5 ;               % elastic modulusEmin = 0 * E0 ;                                        NU = 0.3 ;                                             si = -1 ; ti = 1 ;                                     sj = 1 ; tj = 1 ;sm = 1 ; tm = -1 ;sp = -1 ; tp = -1 ;ID = 1 ;                   % ID = 1: plane stress problem,   ID = 2: plane strain problem    dce    = zeros( Syselem ,1 ) ;           for loopi = 1 : Syselem                      c = c + 0.5 * (xPhys(loopi)/(1 + RE * (1-xPhys(loopi)))) * nodes_d(loopi,:) * E0 * KE * nodes_d(loopi,:)' ;  % structural compliance        dce(loopi) = nodes_d(loopi,:) * (((1 + Rbt)/(1+Rbt*(1-xPhys(loopi)))^2 * E0 * Talafa0 ) * Feth0(:,loopi))...            - 0.5 * (1+RE)/(1 + RE * (1-xPhys(loopi)))^2 * nodes_d(loopi,:) * E0 * KE * nodes_d(loopi,:)' ;  % sensitivity of structural compliance         end    Compli = c ;    dce(:) = H * (dce(:)./Hs);          dve = ones(Syselem,1).* v0/(Syselem * v0 * Volfrac);  % volume constraint    dve(:) = H * (dve(:)./Hs);                                %-METHOD OF MOVING ASYMPTOTES-    m     = 1;                % number of constraint functions     a1    = zeros(m,1);       % Column vector with the constants a_i in the terms a_i*z.    c_MMA = 10000*ones(m,1);  % Column vector with the constants c_i in the terms c_i*y_i.    d     = zeros(m,1);       % Columns vector with the constants d_i in the terms 0.5*d_i*(y_i)^2.    xval  = xDes;        f0val = c;          % compliance minimization    df0dx = dce(:);               fval = sum(xPhys.*v0)/(Syselem * v0 * Volfrac)-1 ;    % volume constraint    dfdx = dve' ;                                            [xmma, ~, ~, ~, ~, ~, ~, ~, ~, low,upp] = ...        mmasub(m, n, loop, xval, xmin, xmax, xold1, xold2, ...        f0val,df0dx,fval,dfdx,low,upp,a0,a1,c_MMA,d);      % Update MMA Variables        xnew     = reshape(xmma , Syselem , 1 );    xold2    = xold1(:);    xold1    = xDes(:);    xDes      = xnew ;    xDes      = xold1 + qnVari * (xDes-xold1) ;                                     changeVari = max(abs(xDes-xold1)) ;    %-UPDATE MODEL DATA-    xPhys = (H * xDes(:))./Hs;                    % physical density    volfrac = 100 * mean(xPhys(:)) ;              % volume fraction    v      = sum(xPhys) * v0 * ObjScale ;    %-PRINT RESULTS AND PLOT DENSITIES-    disp([' It.: ' sprintf('%-4i',loop) ' c.: ' sprintf('%-7.2f',Compli)...        ' Volfrac.: ' sprintf('%-7.2f',volfrac) 'Volume.: ' sprintf('%-7.2f',v ) ...        ' chVari.: ' sprintf('%-7.3f',changeVari ) 'MaxTdif.: ' sprintf('%-7.1f',max(TEdiffe) ) 'Penal.: ' sprintf('%-7.0f', RE )])    xPhysFig = reshape(xPhys ,nely,nelx) ;    colormap(gray); imagesc(1-xPhysFig); axis equal; axis tight; axis off;pause(1e-6);end

3 仿真结果

4 参考文献

[1]黄雅洁. 圆筒热力耦合冲击下热应力场的有限元分析及MATLAB编程[D]. 昆明理工大学, 2015.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

5 代码下载

【物理应用】Matlab实现两端固支梁热力耦合的有限元分析相关推荐

matlab梁应力分析视频,有限元分析及应用清华大学曾攀主讲视频教程
0.0 有限的单元无限的能力 1.1 力学的分类:质点.刚体.变形体的力学 1.2 变形体力学的要点 1.3 微分方程求解的方法 1.4 关于函数逼近的方式 1.5 针对复杂几何域上的函数表征及逼近 ...
两端固定弦的自由振动 | 分离变量法（一）| 偏微分方程（十三）
分离变量法是一种可用来求某些典型区域上定解问题精确解的经典方法,本章将通过典型例子,介绍分离变量法的基本思想和具体步骤,并提出方法的理论核心--固有值问题,进而从Fourier展开角度来认识和应用分离 ...
matlab模拟n维谐振子,n维耦合谐振子的能量谱条件数理论研究
n维耦合谐振子的能量谱条件数理论研究 Theoretical Study of Energy Spectrum Condition Number of n-Dimension Coupled Harm ...
三次固支样条matlab,matlab连续梁程序的编制与使用
<matlab连续梁程序的编制与使用>由会员分享,可在线阅读,更多相关<matlab连续梁程序的编制与使用(21页珍藏版)>请在人人文库网上搜索. 1.第三章连续梁程序的编制 ...
matlab 索力迭代,索梁组合体系桥梁施工索力确定的割线迭代法
摘要索梁组合体系桥梁属于高次超静定结构,由于各种非线性因素的影响,该类桥梁施工索力的确定有着很高的技术要求,针对索梁组合体系桥梁中主梁采用满堂支架施工的情况,提出了割线迭代算法,该方法思路清晰,计算 ...
固高倒立摆 matlab,MATLAB-Control-Software 针对固高直线一级倒立摆 - 下载 - 搜珍网...
MATLAB版实时控制软件/ MATLAB版实时控制软件/MATLAB Patch/ MATLAB版实时控制软件/MATLAB Patch/Readme-Important.txt MATLAB版实时 ...
matlab物理运动,MATLAB软件在中学物理运动学教学中的应用
卞救软件研制软件在中学物理运动学教学中的应用魏青 (北京市海淀区教师进修学校附属实验学校 ,北京 ) 摘要 :计算机技术被广泛应用于 ...
用matlab求解物理方程,MATLAB计算四类数学物理方程的举例求解题库.ppt
数学物理建模与计算机辅助设计第5章四类数学物理方程的求解举例本章内容 §5.1 求解本征值型数学物理方程 §5.2 求解稳定型数学物理方程 §5.3 求解热传导型数学物理方程 §5.4 求解波动 ...
matlab运行时风扇,TCFD和CAESES耦合优化案例-轴流风扇
TCFD是CFD Support团队为我们带来的新一代叶轮机械专用CFD模拟工具.TCFD不受用户人数和核数的限制,具有完全自动化的流程,极大地提升了CFD模拟的效率:同时,它保持求解器的开源,可以由 ...

【物理应用】Matlab实现两端固支梁热力耦合的有限元分析

1 简介

2 部分代码

3 仿真结果

4 参考文献

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

5 代码下载

【物理应用】Matlab实现两端固支梁热力耦合的有限元分析相关推荐

最新文章

热门文章