三片式直角应变花主应力及其夹角计算

#三片式直角应变花主应力及其夹角计算
是对于一些受力复杂的部位,很难判断受力主方向的测点,可采用三片式直角应变花进行测试,本文主要对三片式直角应变花的主应力和主应力方向进行计算。
三片式直角应变花的主应变计算公式如下(式1):在这里插入图片描述
对于任意角度的应变计算公式如下(式2):
在这里插入图片描述

第一第二主应力计算公式如下(式3):
在这里插入图片描述
主应变与0°应变片夹角(式4):
在这里插入图片描述
上述公式计算了第一第二主应变和第一第二主应力的大小,但是在夹角计算公式描述为主应变与0°夹角,而非第一或第二主应力与0°的夹角。为了便于判断计算出第一主应力和0°的夹角。本文按照如下思路:

  1. 通过公式4计算主应变的夹角
  2. 通过公式1计算第一主应变的大小;
  3. 将步骤1计算出的夹角带入公式2中计算出主应变大小
  4. 比较步骤2和步骤3计算出的主应变大小是否相等,进而确定第一主应变的夹角。
    其Matlab的计算代码如下:
function [stress_1, stress_2, Angle]=Prin_stress(data,E,u)
%data的三列分别为0°、90°、45°应变数据,E、u分别为弹性模量和泊松比
%stress_1为第一主应力
%stress_2为第二主应力
%angel为第一主应力方向与0°应变片夹角
	c1=data(:,1)+data(:,2);
    c2=data(:,1)-data(:,2);
    c3=2*data(:,3)-data(:,1)-data(:,2);
    L2=(c2.^2+c3.^2).^0.5;
    stress_1=E*(c1/(1-u)+L2/(1+u))/2;
    stress_2=E*(c1/(1-u)-L2/(1+u))/2;
    angle_1=atan(c3./c2);
    angle_2=angle_1+pi;
    strain_A_1=c1/2+(c2.^2+c3.^2).^0.5*0.5;
    strain_B_1=c1/2+c2.*cos(angle_1)*0.5+c3.*sin(angle_1)*0.5;
    strain_B_2=c1/2+c2.*cos(angle_2)*0.5+c3.*sin(angle_2)*0.5;
     Angle=zeros(length(strain_A_1),1);
    for i=1:1:length(strain_A_1)
        if abs(strain_A_1(i)-strain_B_1(i))<=0.15
            Angle(i)=angle_1(i)*180/pi/2;
        elseif abs(strain_A_1(i)-strain_B_2(i))<=0.15
            Angle(i)=angle_2(i)*180/pi/2;
        else
            msgbox('数据出现问题');
             break;
        end
    end
end
  • 2
    点赞
  • 13
    收藏
    觉得还不错? 一键收藏
  • 2
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值