C语言编程实现求解应力向量,基于MATLAB与Visual+C%2b%2b实现主应力的编程求解.pdf...

您所在位置：网站首页 > 海量文档

&nbsp>&nbsp计算机&nbsp>&nbspC/C++资料

基于MATLAB与Visual+C%2b%2b实现主应力的编程求解.pdf4页

本文档一共被下载：次,您可全文免费在线阅读后下载本文档。

下载提示

1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。

2.该文档所得收入(下载+内容+预览三)归上传者、原创者。

3.登录后可充值，立即自动返金币，充值渠道很便利

第15卷第3期 湖南工程学院学报 Vbl．15．No．3

20()5年9月 o￡Huna：aInstituteof Sept．2005

Journal Engineering

基于MATLAB与Visual

莫江涛，王静文

(中南大学机电工程学院，湖南长沙410083)

现一点应力状态主应力及应力主方向的计算方法．

关键词：MATLAB；VC++；应力状态；主应力；应力主方向

文献标识码：A

中图分类号：THl23．4

确定弹性体中一点应力状态的主应力及应力主 ，为3阶单位阵，解出主应力后，将3个主应力分别

方向，是应力分析的主要问题之一．许多问题都要放 代入如下的线性方程组：

在主应力空间内进行讨论，如：屈服条件、塑性力学 (S—Af)V=0 (3)

问题的有限元分析等，都需要求出单元体主应 即可求出与各个主应力对应的应力主方向V．

力…．本文研究在单元体上六个应力分量已知的情

况下，通过MATLAB与VC++混合编程实现主应量．

力及应力主方向计算机求解的原理和方法． 综上可知，已知应力张量，求解主应力及应力主

方向只需分别求解方程(2)和方程(3)．

1应力分析的矩阵方法

2程序实现

文献[2]对应力张量在坐标变换下的相似性原

理进行了详细研究，结果表明：处于不同坐标系下的 文献[2]及上文已说明，求解一点应力状态的主

应力张量是相似矩阵；相似矩阵有相同的特征值，即 应力及主方向实质是求解该点应力张量的特征值与

应力张量的特征值不随坐标变换而改变，这些特征 特征向量．关于矩阵特征值与特征向量的计算已有

值是应力张量的主应力，其相应的特征向量是该点

的应力主方向．这一论述的实际意义在于，若已知某 方法、QR方法等，读者可参考文献[4]．运用文献

点的应力分量，求该点的主应力及主方向只需求解

应力张量的特征值及特征向量，这为主应力的计算 都能实现主应力及主方向的计算机求解，但考虑到

机求解提供了理论依据． 编程难度、运算效率、程序可靠性及方便用户交互等

即，已知应力张量的6个分量为：如，Sy，S。，Sxy，

s。，s。，写成矩阵形式为： 实现主应力及主方向的计算机求解．

『如嘶‰] 2．1程序功能设计与实现方法

S=SxySy Syz (1) 本程序的主要设计功能为：①根据用户输入的

b。跏。：J 应力张量，计算主应力及应力主方向；②实现与用户

求解应力张量S的主应力，只需求解如下的一元三 交互；③输出计算结果；等．

次方程： 实现方法：用VC++制作图形界面，实现与用

S—A，l (2) 户的交互功能，将用户输入的参数传递给MATLAB

计算引擎完成主应力及主方向的求解，将结果返回

上式称为应力张量的特征多项式[3]3．式(2)中，

收稿日期：2005—04—22

作者简介：莫江涛(1979一)，男，硕士研究生，研究方向：机械设计及制造

万方数据

24

发表评论

请自觉遵守互联网相关的政策法规，严禁发布色情、暴力、反动的言论。

用户名:

验证码:

匿名?

发表评论