悬臂梁有限元分析作业matlab,

　　【物理应用】Matlab实现两端固支梁热力耦合的有限元分析_MATLAB仿真博客的技术博客_博客

　　一简介矩阵实验室实现两端固支梁热力耦合的有限元分析

　　2部分代码% - %基本参数（桥状结构)% -清晰；水平总计=100%水平总计=70%垂直总计=0.5%初始体积分数vol frac=0.25%允许的体积分数var imin=1e-3；%设计变量的最小值obj scale=14e-3；rmin=2；%过滤器半径_ uni=35%分布式负载loadqnVari=0.3% STM应用于设计变量=28刚度的%斜率=16热应力系数的%斜率tala fa 0=12.1 e-6；%热膨胀系数两个设备=[110]；% - %离散化特征%-%-单元刚度-a=0.5；b=0.5h=1；%元素厚度v 0=(2 * a)*(2 * b)* h；固体元素的体积百分比=2.1e5%弹性模量min=0 * E0NU=0.3 si=-1；ti=1；SJ=1；TJ=1；sm=1；TM=-1；sp=-1；TP=-1；ID=1；% ID=1:平面应力问题，ID=2:平面应变问题dce=零(Syselem，1)；for loopi=1:sys elem c=c 0.5 *(xPhys(loopi)/(1 RE *(1-xPhys(loopi)))* nodes _ d(loopi，)* E0 * KE * nodes_d(loopi，)；%结构符合性dce(loopi)=nodes_d(loopi，)(((1 rbt)/(1 rbt*(1-xphys(loopi)))^2 * E0 *塔拉法0 ) * Feth0(:loopi)).-0.5 *(1 re)/(1 re *(1-xphys(loopi)))^2 * nodes _ d(loopi，)* E0 * KE * nodes_d(loopi，)；结构柔度和柔度的灵敏度百分比=c；dce(:)=H * (dce(:)./Hs)；dve=ones(Syselem，1).* v 0/(sys elem * v 0 * vol frac)；%音量约束dve(:)=H * (dve(:)./Hs)；%-移动渐近线法-m=1；约束函数的%数量a1=零(m，1)；项a _ i * z. c _ MMA=10000 *个(m，1)中常数为a_i的%列向量；项c_i*y_i. d=零(m，1)中常数为c_i的%列向量；%列向量带有0.5*d_i*(y_i)^2.项中的常数d _ IX val=xdesf 0 val=c；%顺应性最小化df0dx=DCE(:)；fval=sum(xPhys .* v 0)/(sys elem * v 0 * vol frac)-1；%音量约束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)；%更新MMA变量xnew=shape(xmma，Syselem，1)；xold 2=xold 1(:)；xold 1=xDes(:)；xDes=xnewxDes=xold 1 qn vari *(xDes-xold 1)；change vari=max(ABS(xDes-xold 1))；%-更新模型数据- xPhys=(H * xDes(:)./Hs；%物理密度vol frac=100 * mean(x phys(:))；%体积分数v=sum(x phys)* v 0 * obj scale；%-打印结果和打印密度- disp([ It . sprintf(%-4 I ，loop) c . sprintf(%-7.2f ，Compli). 沃尔弗拉克 sprintf(%-7.2f ，volfrac)音量 sprintf(%-7.2f ，v). 奇瓦里 sprintf(%-7.3f ，changeVari ) MaxTdif . sprintf(%-7.1f)，max(TEdiffe) ) Penal . sprintf(%-7.0f ，RE)])xPhys fig=shape(xPhys，nely，nelx)；彩色地图（灰色);imagesc(1-xphyspig)；轴相等；轴紧；轴关闭；暂停(1e-6)；end3仿真结果

　　编辑

　　四参考文献[1]黄雅洁。圆筒热力耦合冲击下热应力场的有限元分析及矩阵实验室编程d .昆明理工大学， 2015.

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

　　编辑

郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。