将弧度化为秒的换算系数

将弧度化为秒的换算系数 ρ=206262

参考:《仪器精度理论》 马宏主编 北京航空航天大学出版社




clc;clear;close all; gamma=17.65; D=6; Q=0; cc1=0; fai=deg2rad(34); delta=fai; Ka=(1-sin(fai))./(1+sin(fai)); figure; hold on; theta11=pi/2-1/2.*(asin(sin(delta)./sin(fai)))+delta./2; Kanw1=3.*(cos(theta11).^2+Ka.*sin(theta11).^2)./(3-(1-Ka).*cos(theta11).^2);%小主应力侧压力系数 Kanw2=0; A1=Kanw1.*tan(delta); A2=3*(1-Ka)*cos(fai)/(2*(3-(1-Ka)*cos(theta11)^2)); A3=3*(1-Ka)*cos(fai)^2./(2*sin(fai)*(3-(1-Ka)*cos(theta11)^2)); B1=cc1.*(1+(1-Ka).*(((Kanw1*cos(theta11).^2)/3)-sin(theta11).^2))./tan(fai); B2=cc1*(1-Ka)*cos(fai).^2./2*sin(fai); B2=B2*(1+(1-Ka)*cos(theta11)^2./3*(3-(1-Ka)*cos(theta11)^2)); B3=cc1*(1-Ka)*cos(fai).^2./(2*tan(fai)*sin(fai)); B3=B3*(1+(1-Ka)*cos(theta11)^2./3*(3-(1-Ka)*cos(theta11)^2))-cc1./tan(fai); B4=cc1*(1-Ka).*(((Kanw1*cos(theta11).^2)/3)-sin(theta11).^2)./tan(fai); % beita_value=deg2rad(50):deg2rad(1):deg2rad(60); beita_value=deg2rad(89); e=deg2rad(0.1); m=cc1*(1-Ka).*(Kanw1*cos(theta11).^2/3-sin(theta11).^2)./tan(fai); for j=1:length(beita_value) %循环beita从fai—90°,β从φ+1°开始避免除零 beita=beita_value(j); R0=D./(sin(beita).*exp((beita-fai).*tan(fai))-sin(fai)); z(j)=R0.*sin(fai);%筒仓高度 La=cos(beita).*R0.*exp((beita-fai).*tan(fai)); %Lb=cos(fai).*R0-La; %J=2*(B+Lb)./(B*Lb); q2=0; %上部筒仓带来的力 Fai_value=fai:e:beita-e; %Fai_value=fai:e:deg2rad(50); mm=length(Fai_value); for ii=1:mm Fai=Fai_value(ii); omiga=@(Fai) pi/2+fai-Fai; B=@(Fai) D./tan(omiga(Fai)); R=@(Fai) R0.*(exp((Fai-fai).*tan(fai))); %z0=@(Fai) R(Fai).*sin(Fai)-z(j); L=@(Fai) R(Fai).*cos(Fai)-La; %土条的长度 z0(ii)=R(Fai)*sin(Fai)-z(j); rou_1=@(Fai) A1-tan(fai-Fai)+2.*Kanw2.*tan(fai).*sin(omiga(Fai)).*L(Fai)./B(Fai)+(2.*Kanw2.*tan(fai).*cos(omiga(Fai)).*L(Fai)+B(Fai).*Kanw1).*(tan(fai)+cot(omiga(Fai)))./(B(Fai).*(1-tan(fai).*cot(omiga(Fai)))); rou_2=@(Fai) gamma.*L(Fai)+2.*cc1.*sin(omiga(Fai)).*L(Fai)./B(Fai)+(tan(fai)+cot(omiga(Fai))).*((m+cc1.*cot(omiga(Fai))).*B(Fai)+2.*cc1.*cos(omiga(Fai)).*L(Fai))./(B(Fai).*(1-tan(fai).*cot(omiga(Fai)))); rou_3=@(Fai) R(Fai).*(tan(fai).*sin(Fai)+cos(Fai))./L(Fai); end tspan=[fai beita-e]; f=@(Fai,q) [rou_3(Fai).*(rou_1(Fai).*q+rou_2(Fai))]; %q为平均竖向应力 q0=q2; [z,q]=ode45(f, tspan, q0) seigema_h=Kanw1*q+cc1*(1-Ka).*(Kanw1*cos(theta11).^2/3-sin(theta11).^2)./tan(fai);%水平土压力 plot(seigema_h./(gamma.*D), z0(ii)./D, '-k', 'Linewidth', 1.5); end hold on; grid off box off set(gca,'LineWidth',1.5,'FontSize',10,'FontName','Times New Roman'); xlabel('$\sigma_{\mathrm{h}} / \gamma D$', 'Interpreter', 'latex','FontSize', 14,'FontName', 'Times New Roman'); ylabel('D/(m)','FontSize',14,'FontName','Times New Roman'); axis ij;为什么最后算出来的值好多数量级,太大了,不符合常理,帮我检查原因。正确的规律应该是从0先增大后减小至0
03-19
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

东山一角

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值