哈工大材料力学上机实验资料报告材料
材料力学实验报告
材料力学实验报告材料力学实验报告引言：材料力学是一门研究材料在外力作用下的力学性能和变形规律的学科。
通过实验研究，我们可以深入了解材料的力学性质，为工程设计和材料选择提供依据。
本报告将介绍我们在材料力学实验中的观察和结果，并对实验数据进行分析和讨论。
实验一：拉伸试验拉伸试验是材料力学实验中最常见的一种试验方法，用于研究材料在拉伸载荷下的力学性能。
我们选择了一根标准的金属试样，将其固定在拉伸试验机上，并逐渐施加拉伸力。
通过测量试样的应变和应力，我们得到了应力-应变曲线。
实验结果显示，随着拉伸力的增加，试样开始发生塑性变形。
在这个阶段，应力与应变呈线性关系，即应力随着应变的增加而线性增加。
然而，当拉伸力达到一定程度时，试样出现断裂。
通过观察断裂面的形态，我们可以判断材料的断裂模式，如韧性断裂、脆性断裂等。
进一步分析应力-应变曲线，我们可以得到一些重要的力学参数，如屈服强度、抗拉强度和延伸率。
屈服强度是材料开始发生塑性变形时的应力值，抗拉强度是试样抵抗拉伸力的最大极限，而延伸率则表示试样在断裂前的延展能力。
这些参数对于材料的工程应用和性能评估至关重要。
实验二：硬度测试硬度是材料力学中另一个重要的性能指标，它反映了材料抵抗外力的能力。
我们采用了维氏硬度计进行硬度测试，将金属球压入试样表面并测量压痕的直径。
根据硬度计的原理，我们可以计算出试样的硬度值。
硬度测试的结果显示，不同材料的硬度值存在明显差异。
硬度值高的材料通常具有较好的抗压性能，适用于承载大压力的工程应用。
而硬度值低的材料则更容易受到外力的破坏，适用于需要易变形的应用场景。
实验三：弯曲试验弯曲试验用于研究材料在弯曲载荷下的力学性能。
我们选择了一根长条状的试样，通过在试样两端施加力矩，使试样发生弯曲变形。
通过测量试样的挠度和应力分布，我们可以得到弯曲试验的结果。
实验结果表明，试样的挠度与施加的力矩呈线性关系。
在试样的底部，应力最大，而在试样的顶部，应力最小。
材料力学上机大作业(哈工大)
H a r b i n I n s t i t u t e o f T e c h n o l o g y材料力学上机报告课程名称：材料力学设计题目：二向应力状态分析院系：XXXXXX班级：XXXXXX设计者：XXXXXX学号：XXXXXX设计时间：2013.06.18哈尔滨工业大学二向应力状态分析一：课题要求1.输入：任意一点的应力状态：（σx、σy、τxy）；某截面方位角α2.输出：输入点的主应力（σ1、σ2、σ3），方位角α斜截面上的应力σα、τα。
及主方向角α3.画出应力圆示意图。
4.程序运行时为界面显示形式。
二：程序框图三：所编程序x=str2double(get(handles.edit1,'string'));y=str2double(get(handles.edit2,'string'));xy=str2double(get(handles.edit3,'string'));M=str2double(get(handles.edit4,'string')); %将窗口输入值分别赋给x,y,xy,M b=sqrt((x/2-y/2)^2+xy^2);x1=(x+y)/2+b;x3=(x+y)/2-b;x2=0;if x1<0x2=x1;x1=0;endt=(x1-x3)/2;M=M*pi/180;b1=(x+y)/2+(x-y)*cos(2*M)/2-xy*sin(2*M);b2=(x-y)*sin(2*M)/2+xy*cos(2*M);b3=90*atan((-2*xy)/(x+y))/pi;%计算输出的主切应力大小、方向和截面上的应力并赋值set(handles.edit5,'string',x1);set(handles.edit6,'string',x2);set(handles.edit7,'string',x3);set(handles.edit9,'string',t);set(handles.edit10,'string',b3);set(handles.edit11,'string',b1);set(handles.edit12,'string',b2);%在输出窗口显示主切应力大小、方向和截面上应力b4=sqrt(b.^2+t.^2);v1=(x+y)/2-b4:0.001:(x+y)/2+b4;b11=sqrt(b4.^2-(v1-(x+y)/2).^2);b12=-sqrt(b4.^2-(v1-(x+y)/2).^2);%绘制应力圆上的点axes(handles.axes1); %选择应力圆的输出地址plot(v1,b11,v1,b12);grid on%绘制应力圆以上程序为在matlab中使用GUI编程时的主代码，界面代码请见m文件。
材料力学实验报告
材料力学实验报告引言：材料力学是研究物质在外力作用下的变形和破坏行为的科学。
在工程领域，材料力学实验是非常重要的，它能提供关于材料性能的定量数据，用于设计和优化结构。
本篇实验报告将介绍一项材料力学实验，包括实验目的、实验装置和实验过程，重点关注实验结果的分析和讨论。
实验目的：本次实验旨在研究一种金属材料的拉伸性能，通过对材料在不同载荷下的应力-应变关系曲线的测定，获得材料的力学性能参数，如屈服强度、抗拉强度和延伸率等。
同时，通过断口分析，了解材料的破坏行为和断裂机制。
实验装置：本次实验采用的材料力学实验装置包括拉伸试验机、计算机数据采集系统和金属试样。
拉伸试验机主要包括上夹具和下夹具，通过电机驱动实现上下夹具之间的拉伸和压缩运动。
计算机数据采集系统用于实时记录试验过程中的应变和载荷数据。
金属试样采用标准的矩形横截面形状，制备精细，确保试样的几何尺寸以及表面质量。
实验过程：1. 调整试验机，确保试样正确安装在上下夹具之间，并进行预应力调校。
2. 设置拉伸速率和采样频率，开始实验。
3. 开始加载并进行拉伸实验，直至试样断裂。
4. 实时记录应变和载荷数据，生成应力-应变曲线。
5. 对断口进行分析，观察破坏模式和断裂特征。
实验结果分析：基于实验数据，通过应力-应变曲线的绘制和分析，可以得到材料的力学性能参数。
应力-应变曲线的特点是：一开始，材料的应变随载荷的增加近似线性增加，这是材料的弹性区域。
当应变逐渐超过一定程度时，材料的应变开始迅速增加，即材料进入了屈服区。
进一步增加载荷，材料的应变仍呈线性增加，但增加的速率较之前小，这是材料的塑性区。
除了绘制应力-应变曲线，我们还可以计算出材料的屈服强度、抗拉强度和延伸率等力学性能参数。
屈服强度是指试样开始进入塑性阶段时的应力值，抗拉强度是试样发生破裂时的最大应力值，而延伸率则反映了试样在拉伸过程中的延伸能力。
断口分析是评价材料破坏行为和断裂机制的重要手段。
通过观察断口的形貌特征和变异，可以判断材料的韧性和脆性。
哈尔滨工程大学力学实验报告
Harbin Engineering University材料力学实验报告实验题目：姓名：班级：学号：设备号：实验时间：年月日时分—时分力学实验教学中心材料力学实验上课要求一、实验项目分两部分：第一部分为“必做实验”，是必选的实验；第二部分为“选做实验”，要在提供的实验题目中选做要求的实验个数。
所有实验要在开课周内完成，逾期不补。
二、平时要多留意网上通知。
三、必须在选定的时间前10分钟到实验室进行实验，不许迟到。
如果不按时到课，按缺课处理。
四、进入实验室时，必须首先用学生卡进行刷卡签到。
并带计算器和直尺。
五、课前要认真预习，写好预习报告，注意要把上课时间和设备号写在实验报告的封面上，无预习报告或预习不合格的不能做实验。
六、实验报告要字迹工整，图表规范，书面整洁，在做完实验的一周内，把实验报告按班级统一交给上课教师，逾期按没交处理。
七、进入实验室必须遵守实验室的各项规章制度，不得大声喧哗，保持室内卫生，正确使用仪器设备，注意安全。
八、实验中要详实记录实验数据，不得抄袭和携带它组实验数据。
九、实验结束后应将所用的仪器设备整理好，摆放整齐。
预习报告须经上课老师检查签字后方可离开实验室。
预习报告一、实验题目：二、实验目的：三、实验仪器：四、实验原理（原理图、公式推导和文字说明）：五、实验内容及主要步骤：六、实验数据记录表格（自己设计）：教师签章月日实验报告一、实验题目：二、实验目的：三、实验仪器：四、实验原理（原理图、公式推导和文字说明）：五、实验数据处理（整理表格、计算过程、计算结果）：六、总结及应用（结论、误差分析、讨论问题、本实验的应用）：。
大学材料力学实验报告
大学材料力学实验报告大学材料力学实验报告引言材料力学实验是大学材料科学与工程专业中的一门重要课程。
通过实验，我们可以深入了解材料的力学性质和行为，为材料设计和应用提供基础数据和理论依据。
本次实验旨在通过拉伸试验和硬度测试，探究不同材料的力学性能和硬度特点。
实验一：拉伸试验拉伸试验是一种常用的力学实验方法，用于评估材料的强度、延展性和塑性等性能。
在实验中，我们选择了三种常见的材料进行拉伸试验：钢材、铝材和塑料。
1. 实验步骤首先，我们准备了三个不同材料的试样，分别是圆柱形的钢材、铝材和塑料样品。
然后，将试样固定在拉伸试验机上，并施加逐渐增大的拉力，直到试样断裂为止。
在拉伸过程中，我们记录下拉力和试样的伸长量，以绘制应力-应变曲线。
2. 实验结果通过拉伸试验得到的应力-应变曲线可以反映材料的力学性能。
钢材的应力-应变曲线呈现出明显的弹性区和塑性区，具有较高的屈服强度和延展性。
铝材的应力-应变曲线也呈现出弹性和塑性的特点，但相对于钢材来说，其屈服强度和延展性较低。
而塑料的应力-应变曲线则主要表现为塑性变形，没有明显的弹性区。
实验二：硬度测试硬度是材料力学性能的重要指标之一，用于评估材料的抗压能力和耐磨性。
在实验中，我们选择了三种不同硬度的材料进行硬度测试：钢材、铝材和陶瓷。
1. 实验步骤我们使用了维氏硬度计和洛氏硬度计对试样进行硬度测试。
首先，将试样固定在硬度计上，然后施加一定的压力，观察压头对试样的印痕情况。
根据印痕的大小和形状，我们可以得出试样的硬度数值。
2. 实验结果通过硬度测试，我们发现钢材具有较高的硬度数值，表明其具有较高的抗压能力和耐磨性。
铝材的硬度数值相对较低，说明其相对较软。
而陶瓷的硬度数值最高，表明其具有极高的抗压能力和耐磨性。
结论通过本次实验，我们深入了解了材料的力学性能和硬度特点。
拉伸试验结果表明，钢材具有较高的屈服强度和延展性，铝材次之，而塑料则主要表现为塑性变形。
硬度测试结果显示，钢材具有较高的硬度数值，铝材较低，而陶瓷的硬度最高。
哈工大材力上机实验报告
材料力学电算实验压杆的临界力计算院系：机电工程学院班级：设计者：学号：指导教师：张桂莲软件要求：设计时间：一.概述：本程序使用Microsoft Visual Basic编写，可以对不同材料、不同约束类型、不同截面类型的压杆进行临界力的计算。
杆件的参数可以输入，得出结果之后也可以清零。
二、问题分析及相关公式：1、压杆稳定当短粗杆受压时(图1)，在压力F由小逐渐增大的过程中，杆件始终保持原有的直线平衡形式，直到压力F达到屈服强度载荷F s(或抗压强度载荷F b)，杆件发生强度破坏时为止。
但是，如果用相同的材料，做一根与图1a所示的同样粗细而比较长的杆件(图1b)，当压力F比较小时，这一较长的杆件尚能保持直线的平衡形式，而当压力F逐渐增大至某—数值F1时，杆件将突然变弯，不再保持原有的直线平衡形式，因而丧失了承载能力。
我们把受压直杆突然变弯的现象，称为丧失稳定或失稳。
此时，F1可能远小于F s (或F b)。
可见，细长杆在尚未产生强度破坏时，就因失稳而破坏。
图1在研究压杆稳定时，我们用一微小横向干扰力使处于直线平衡状态的压杆偏离原有的位置，如图1所示。
当轴向压力F由小变大的过程中，可以观察到：1)当压力值F1较小时，给其一横向干扰力，杆件偏离原来的平衡位置。
若去掉横向干扰力后，压杆将在直线平衡位置左右摆动，最终将恢复到原来的直线平衡位置。
2)当压力值F2超过其一限度F cr时，平衡状态的性质发生了质变。
这时，只要有一轻微的横向干扰，压杆就会继续弯曲，不再恢复原状，。
3)界于前二者之间，存在着一种临界状态。
当压力值正好等于F cr时，一旦去掉横向干扰力，压杆将在微弯状态下达到新的平衡，既不恢复原状，也不再继续弯曲，。
临界状态是杆件从稳定平衡向不稳定平衡转化的极限状态。
压杆处于临界状态时的轴向压力称为临界力或临界载荷，用F cr表示。
2、两端铰支细长压杆的临界力图2为一两端为球形铰支的细长压杆，其临界力公式为：图222lEIF cr π=(1)式(1)又称为欧拉公式。
材料力学实验报告模板
材料力学实验报告
报告标题：_________________________
一、实验目的
（简要说明实验的主要目的和预期达到的学习效果）
二、实验原理
（描述实验的理论基础，包括相关的材料力学理论和公式）
三、实验设备和材料
（列出进行实验所需的主要设备、工具和材料）
四、实验步骤
（详细描述实验的操作步骤，包括准备工作和具体的实验流程）
五、实验数据和结果
5.1 实验数据
（记录实验过程中收集的所有数据，可使用表格形式呈现）
5.2 实验结果
（根据实验数据计算出的结果，包括必要的图表和计算过程）
六、结果分析
（分析实验结果，对比理论值和实际值的差异，解释可能的原因）
七、实验结论
（总结实验结果，得出结论，评价实验的成功与否及其科学意义）
八、实验心得和建议
（个人对实验的感想，包括实验过程中的体会、遇到的问题及建议）
九、参考文献
（列出实验报告中引用的所有参考文献）
报告人：_________________________
学号：_________________________
班级：_________________________
日期：__________年__________月__________日。
材料力学实验报告报告
材料力学实验报告报告一、实验目的本实验旨在通过测量不同材料的力学性能参数，了解材料的力学性质，以及分析不同材料的力学性能差异。
二、实验原理1.弹性模量：弹性模量是评价材料抗弯刚性的一个重要指标，可以通过测量材料的拉伸和压缩位移来确定。
拉伸试验时，通过加载材料，测量应力和应变的关系，然后通过斜率求出弹性模量。
2.屈服强度：材料的屈服强度是指材料在拉伸过程中开始出现塑性变形时的抗拉强度，也是一个重要的力学性能参数，通过拉伸试验中的负荷-变形曲线求得。
3.断裂强度：材料的断裂强度是指在材料断裂前能承受的最大负荷，通过拉伸试验中的负荷-变形曲线求得。
三、实验设备与试样准备1.实验设备：拉伸试验机、压缩试验机、材料硬度测试仪等。
2.试样准备：选取不同的材料（如钢材、铝材、铜材等）制作成相同形状、尺寸的试样。
四、实验步骤1.弹性模量测定：(1)将试样固定在拉伸试验机上，设定初始载荷并开始加载。
(2)根据试验机上的位移计和负荷计，测量不同应力水平下的应变，并记录数据。
(3)通过绘制应力-应变曲线，根据直线部分的斜率求得材料的弹性模量。
2.屈服强度测定：(1)将试样固定在拉伸试验机上，设定初始载荷并开始加载。
(2)根据试验机上的压力计和位移计，测量不同载荷下的变形，并记录数据。
(3)通过绘制负荷-变形曲线，找到试样开始出现塑性变形的点，根据载荷计的读数求得材料的屈服强度。
3.断裂强度测定：(1)将试样固定在拉伸试验机上，设定初始载荷并开始加载。
(2)根据试验机上的压力计和位移计，测量试样在拉伸过程中的载荷和位移，并记录数据。
(3)通过绘制负荷-变形曲线，找到试样断裂前的最大负荷，并记录。
五、实验结果与讨论根据实验测量的数据，可以得到不同材料的力学性能参数，如弹性模量、屈服强度和断裂强度。
通过对比不同材料的实验结果，可以得出以下结论：1.钢材的弹性模量较大，机械性能优异。
2.铝材的屈服强度较低，耐腐蚀性能较好。
3.铜材的断裂强度较高，适用于承受较大载荷的工程应用。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
材料力学I上机实验设计报告院系：机电学院班级： 1308***姓名： ***学号： 11308*****指导教师：桂莲时间： 2015年6月一、问题描述1、应力状态分析对于空间或者是平面应力状态的相关计算，如果采用人工计算的方式比较繁琐而且容易出错，对于这种简单的重复计算，编制相应的程序则可以大大提高计算准确度和人工计算强度。
对于平面应力状态，输入量应为（,,x y xy σστ），以及某截面的方位角α，其输出数据应为该单元体所受主应力（123,,σσσ），所受最大剪应力（13max 132σσττ-==），以及方位角为α的斜截面上的应力（,ααστ）以及主方向角σα，同时还要画出其应力圆示意图，以直观的显示其应力状态。
对于空间应力状态，输入量则应该为各应力（,,,,,x y z xy yz xz σσστττ），其输出数据应该为该单元体所受主应力（123,,σσσ），所受最大剪应力（13max 132σσττ-==），同时还要画出其应力圆示意图，以直观的显示其应力状态。
这样，应力状态分析的基本任务就可以完成。
2、常用截面图形几何性质的分析在生活中，有各种各样的几何形状，但是对于工程实际中经常用到的构件，其截面的几何形状则非常有限。
对于不同的截面，其形心位置、对于形心轴的惯性矩也就有所不同，这样在进行如弯曲、扭转等的应力分析时就会到来不便，因此编制相应的程序来计算相关截面的几何性质也就具有了实际应用价值和可行性。
在这部分程序中，截面几何形状分为三角形、矩形、椭圆形、梯形、圆形、扇形等多种形式，对于不同的截面形状，输入量也就不同。
例如，对于扇形应输入直径和圆心角（,d α）；对于梯形则应输入上底、下底和高（,,a b h ）；对于椭圆形，则要输入长轴长和短轴长（,a b ）等等，在此不一一列举，具体输入数据请参看程序运行。
不过对于不同的截面，其输出的量都是相同的，即截面形心的位置、面积、对于形心轴的惯性矩（,,,,C C C C y z y z S I I ），这些输出量就是这些截面的基本几何性质参数，有了这些参数之后则可以对其进行进一步的计算和接下来的分析等问题。
这样，对于工程实际中常见的几种截面形状的几何性质的分析任务就可以用编写程序的方式来完成。
以上两项容即为本次上机实验程序的基本构成，两种功能相互配合、协调运行，以完成程序的稳定运行和人机交互，同时也可以解决工程实际中的一些比较简单的问题和任务，程序编制的任务就可以完成。
二、程序流程本程序分为三个窗体：主界面、应力状态分析界面、截面几何性质分析界面，下面将对各个界面进行程序流程的分析：1、主界面： 开始显示基本信息隐藏计算窗口选择计算任务显示相应的窗口继续运行结束程序运行截图：选择相应的计算程序，然后点击确定即可进入对应的计算程序窗口之中，如果点击退出则退出程序，菜单栏也包含有一些辅助的功能和相关的基本信息。
2、应力状态分析窗口开始定义变量及所用子函数选择应力状态，显示相应的输入提示图输入数据，赋值计算相关的输出量，并将其显示绘制应力圆继续下一步计算结束程序运行截图：点击某应力状态后，相应的输入框可用，输入数据之后，点击确定，输出计算结果，点击绘制应力圆，绘制对应的应力圆及坐标。
点击退出，退出程序；点击返回，返回主界面；点击清除，清除所有数据。
3、截面几何性质分析窗口开始定义变量及所用子函数选择几何形状，显示相应的输入提示图输入数据，赋值计算相关的输出量，并将其显示绘制截面构型图继续下一步计算结束程序运行截图：选择相应的形状类型，点击确定，显示提示图，输入框可用；输入数据，点击计算，输出计算结果和截面构型图；点击清除，清除所有数据和图形显示；点击返回，返回主界面；点击退出，结束程序运行。
按照此流程编制和应用程序，即可实现相应的设计任务。
三、具体某个问题和涉及到的计算公式以及相关理论本程序涉及多条公式及相关理论，现分部分阐述：1、平面应力状态分析平面应力状态任一斜截面上既有正应力又有切应力，则其公式为：平面应力状态主应力大小及方向公式为：（1）主应力计算公式（2）主平面方位计算公式（3）最大切应力：最大切应力发生在与主平面夹45度角的平面方位。
其公式为：应力圆：222222x y x y xy αασσσσσττ+-⎛⎫⎛⎫-+=+ ⎪ ⎪⎝⎭⎝⎭圆心： 2x y a σσ+=半径： 222x y xy R σστ-⎛⎫=+ ⎪⎝⎭根据数学建模结果，编制相应的程序，即可计算出结果。
2、空间应力状态分析三向应力计算公式：特征方程：系数分别为：1x y z I σσσ=++2x yx y zy z xz xy y yz z zx x I στστσττστστσ=++3xyx zx xyy zy xz yz z I στττστττσ=解出此特征方程即可求出单元体的三个主应力，求解该方程可以采用盛金公式，容如下：一元三次方程aX 3＋bX 2＋cX ＋d=0，（a ，b ，c ，d ∈R ，且a≠0）。
重根判别式：A=b 2－3ac ；B=bc －9ad ；C=c 2－3bd ，总判别式：Δ=B 2－4AC 。
当A=B=0时，盛金公式①：X 1=X 2=X 3=－b/(3a)=－c/b=－3d/c 。
当Δ=B 2－4AC>0时，盛金公式②：X 1=(－b －(Y 11/3＋Y 21/3))/(3a)；X 2，3=(－2b ＋Y 11/3＋Y 21/3±31/2 (Y 11/3－Y 21/3)i)/(6a)； 其中Y 1，2=Ab ＋3a (－B±(B 2－4AC)1/2)/2，i 2=－1。
当Δ=B 2－4AC=0时，盛金公式③：X 1=－b/a ＋K ；X 2=X 3=－K/2，其中K=B/A ，(A≠0)。
当Δ=B 2－4AC<0时，盛金公式④：X 1= (－b －2A 1/2cos(θ/3) )/(3a)； X 2，3= (－b ＋A 1/2(cos(θ/3)±31/2sin(θ/3)))/(3a)；其中θ=arccosT ，T= (2Ab －3aB)/(2A 3/2)，(A>0，－1<T<1)。
盛金判别法①：当A=B=0时，方程有一个三重实根；②：当Δ=B 2－4AC>0时，方程有一个实根和一对共轭虚根； ③：当Δ=B 2－4AC=0时，方程有三个实根，其中有一个两重根； ④：当Δ=B 2－4AC<0时，方程有三个不相等的实根。
据此即可解除一元三次方程，得到单元体的主应力。
然后即可根据主应力做出三向应力状态应力圆：121223231313121223231313222222C C C R R R σσσσσσσσσσσσ+=+=+=-=-=-=据此即可做出应力圆及其坐标系。
3、截面几何性质分析本计算程序所用公式较多，如下表所示：按照此表所提供的公式编制程序，即可计算出相应的截面几何性质。
由此，该程序的数学建模过程完成，可以照此编制出相应的程序。
四、程序源代码代码分为三部分：主界面、应力分析窗口、截面几何性质分析窗口：1、主界面'1308***班，11308***10,***'2015年6月21日Private Sub Command2_Click() EndEnd SubPrivate Sub Form_Load()main.Showyingli.Hidejihe.HideEnd SubPrivate Sub jiaoshi_Click()MsgBox "桂莲"End SubPrivate Sub shijian_Click()MsgBox "2015年6月2日" End SubPrivate Sub tuichu_Click()EndEnd SubPrivate Sub zuozhe_Click()MsgBox "11308***10—***"End SubPrivate Sub Command1_Click()If Option1.Value = True Thenyingli.ShowElseIf Option2.Value = True Thenjihe.ShowElseMsgBox "请选择项目"End IfEnd Sub2、应力状态分析窗口Option ExplicitConst pi = 3.1415926Dim σx, σy, σz, σ1, σ2, σ3, τxy, τyz, τxz, τmax, σα, τα, α, ασ As Double Dim temp As DoubleDim I1, I2, I3 As DoubleDim R, C, C12, C23, C31, R12, R23, R31 As DoubleDim j, i As DoublePrivate Sub Command3_Click()EndEnd SubPrivate Sub Command4_Click() yingli.HideEnd SubPrivate Sub Command2_Click() Text1.Text = 0Text2.Text = 0Text3.Text = 0Text4.Text = 0Text5.Text = 0Text6.Text = 0Text7.Text = 0Text8.Text = 0Text9.Text = 0Text10.Text = 0Text11.Text = 0Text12.Text = 0Text13.Text = 0Text14.Text = 0End SubPrivate Sub Command5_Click() Set Picture1.Picture = Nothingσx = Val(Text1.Text)σy = Val(Text2.Text)σz = Val(Text3.Text)τxy = Val(Text4.Text)τyz = Val(Text5.Text)τxz = Val(Text6.Text)α = Val(Text7.Text) / 180 * piIf Option1.Value = True ThenR = Sqr(((σx - σy) / 2) ^ 2 + τxy ^ 2)C = (σx + σy) / 2Picture1.Scale (C - 2 * R - 1, 2 * R + 1)-(C + 2 * R + 1, -2 * R - 1) Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (C - 2 * R - 2, 0)-(C + 2 * R + 2, 0)Picture1.Line (0, 2 * R + 2)-(0, -2 * R - 2)Picture1.DrawWidth = 2Picture1.Circle (C, 0), RPicture1.Line (σx, -τxy)-(σy, τxy)Picture1.DrawWidth = 5Picture1.PSet (σx, -τxy), vbRedPicture1.PSet (σy, τxy), vbRedPicture1.PSet (C + R, 0), vbRedPicture1.PSet (C - R, 0), vbRedPicture1.CurrentX = C - R: Picture1.CurrentY = 0Picture1.Print "σmin"Picture1.CurrentX = C + R: Picture1.CurrentY = 0Picture1.Print "σmax"Pic ture1.CurrentX = σy: Picture1.CurrentY = τxyPicture1.Print "Y(σy,τxy)"Picture1.CurrentX = σx: Picture1.CurrentY = -τxyPicture1.Print "X(σx,-τxy)"Picture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.DrawStyle = 2Picture1.DrawWidth = 1For i = 0 To 200 Step 1Picture1.Line (-20 * R + i * 0.2 * R, -2 * R - 2)-(-20 * R + i * 0.2 * R, 2 * R + 2)Picture1.Line (C - 2 * R - 2, -20 * R + i * 0.2 * R)-(C + 2 * R + 2, -20 * R + i * 0.2 * R)Next iElseIf Option2.Value = True ThenC12 = (σ1 + σ2) / 2C23 = (σ2 + σ3) / 2C31 = (σ3 + σ1) / 2R12 = (σ1 - σ2) / 2R23 = (σ2 - σ3) / 2R31 = (σ1 - σ3) / 2Picture1.Scale (C31 - 2 * R31 - 1, 2 * R31 + 1)-(C31 + 2 * R31 + 1, -2 * R31 - 1)Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (C31 - 2 * R31 - 2, 0)-(C31 + 2 * R31 + 2, 0)Picture1.Line (0, 2 * R31 + 2)-(0, -2 * R31 - 2)Picture1.DrawStyle = 2For i = 0 To 200 Step 1Picture1.Line (-20 * R31 + i * 0.2 * R31, -2 * R31 - 2)-(-20 * R31 + i * 0.2 * R31, 2 * R31 + 2)Picture1.Line (C31 - 2 * R31 - 2, -20 * R31 + i * 0.2 * R31)-(C31 + 2 * R31 + 2, -20 * R31 + i * 0.2 * R31)Next iPicture1.DrawWidth = 2Picture1.Circle (C12, 0), R12Picture1.Circle (C23, 0), R23Picture1.Circle (C31, 0), R31Picture1.CurrentX = C12: Picture1.CurrentY = 0Picture1.Print "C12"Picture1.CurrentX = C23: Picture1.CurrentY = 0Picture1.Print "C23"Picture1.CurrentX = C31: Picture1.CurrentY = 0Picture1.Print "C31"Picture1.CurrentX = σ1: Picture1.CurrentY = 0Picture1.Print "σ1"Picture1.CurrentX = σ2: Picture1.CurrentY = 0Picture1.Print "σ2"Picture1.CurrentX = σ3: Picture1.CurrentY = 0Picture1.Print "σ3"Picture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.DrawWidth = 4Picture1.PSet (C12, 0)Picture1.PSet (C23, 0)Picture1.PSet (C31, 0)Picture1.DrawWidth = 5Picture1.PSet (σ1, 0), vbRedPicture1.PSet (σ2, 0), vbRedPicture1.PSet (σ3, 0), vbRedElseMsgBox "请选择应力状态"End IfEnd SubPrivate Sub Command6_Click()Picture1.Picture = NothingEnd SubPrivate Sub Form_Load()Option1.Value = FalseOption2.Value = False End SubPrivate Sub Option1_Click() Image1.Visible = TrueImage2.Visible = FalseText1.Enabled = TrueText2.Enabled = TrueText3.Enabled = FalseText4.Enabled = TrueText5.Enabled = FalseText6.Enabled = FalseText7.Enabled = TrueText8.Enabled = TrueText9.Enabled = TrueText10.Enabled = TrueText11.Enabled = TrueText12.Enabled = TrueText13.Enabled = TrueText14.Enabled = TrueLabel3.Enabled = FalseLabel5.Enabled = FalseLabel10.Enabled = FalseLabel12.Enabled = FalseLabel13.Enabled = FalseLabel19.Enabled = TrueLabel20.Enabled = TrueLabel21.Enabled = TrueLabel26.Enabled = TrueLabel27.Enabled = TrueLabel28.Enabled = TrueLabel7.Enabled = TrueLabel14.Enabled = True End SubPrivate Sub Option2_Click() Image2.Visible = TrueImage1.Visible = FalseText1.Enabled = TrueText2.Enabled = TrueText3.Enabled = TrueText4.Enabled = TrueText5.Enabled = TrueText6.Enabled = TrueText7.Enabled = FalseText8.Enabled = TrueText9.Enabled = TrueText11.Enabled = TrueText12.Enabled = FalseText13.Enabled = FalseText14.Enabled = FalseLabel3.Enabled = TrueLabel5.Enabled = TrueLabel6.Enabled = TrueLabel7.Enabled = FalseLabel14.Enabled = FalseLabel10.Enabled = TrueLabel12.Enabled = TrueLabel13.Enabled = TrueLabel19.Enabled = FalseLabel20.Enabled = FalseLabel21.Enabled = FalseLabel26.Enabled = FalseLabel27.Enabled = FalseLabel28.Enabled = FalseEnd SubFunction Arccos(X) As Double '反余弦Arccos = Atn(-X / Sqr(-X * X + 1)) + 2 * Atn(1) End FunctionPrivate Sub Command1_Click()If Option1.Value = True Thenσx = Val(Text1.Text)σy = Val(Text2.Text)τxy = Val(Text4.Text)α = Val(Text7.Text) / 180 * piσα = (σx + σy) / 2 + (σx - σy) / 2 * Cos(2 * α) + τxy * Sin(2 * α)τα = -(σx - σy) / 2 * Sin(2 * α) + τxy * Cos(2 * α)If σx <> σy Thenασ = (Atn(2 * τxy / (σx - σy)) / 2) / pi * 180Elseασ = 90End Ifσ1 = (σx + σy) / 2 + Sqr(((σx - σy) / 2) ^ 2 + τxy ^ 2)σ2 = (σx + σy) / 2 - Sqr(((σx - σy) / 2) ^ 2 + τxy ^ 2)σ3 = 0If σ2 > σ1 Thentemp = σ1σ1 = σ2σ2 = tempEnd IfIf σ3 > σ1 Thentemp = σ1σ1 = σ3σ3 = tempEnd IfIf σ3 > σ2 Thentemp = σ2σ2 = σ3σ3 = tempEnd Ifτmax = (σ1 - σ3) / 2Text8.Text = σ1Text9.Text = σ2Text10.Text = σ3Text11.Text = τmaxText12.Text = σαText13.Text = ταText14.Text = ασElseIf Option2.Value = True Thenσx = Val(Text1.Text)σy = Val(Text2.Text)σz = Val(Text3.Text)τxy = Val(Text4.Te xt)τyz = Val(Text5.Text)τxz = Val(Text6.Text)I1 = σx + σy + σzI2 = (σx * σy - τxy ^ 2) + (σy * σz - τyz ^ 2) + (σz * σx - τxz ^ 2)I3 = σx * σy * σz + τxy * τyz * τxz + τxy * τyz * τxz - σx * τyz * τyz - σz * τxy* τxy - σy * τxz * τxzDim aa, bb, cc, dd As DoubleDim a1, b1, c1, derta As Doubleaa = 1bb = -I1cc = I2dd = -I3a1 = bb ^ 2 - 3 * aa * ccb1 = bb * cc - 9 * aa * ddc1 = cc * cc - 3 * bb * ddderta = b1 * b1 - 4 * a1 * c1If (0 = a1) And (0 = b1) Thenσ1 = -bb / (3 * aa)σ2 = -bb / (3 * aa)σ3 = -bb / (3 * aa)'ElseIf derta > 0 Then'Dim y1, y2 As Double'y1 = a1 * b + 3 * a * (-b1 + Sqr(b1 ^ 2 - 4 * a1 * c1)) / 2'y2 = a1 * b + 3 * a * (-b1 - Sqr(b1 ^ 2 - 4 * a1 * c1)) / 2'σ1 = (-b - (y1 ^ (1 / 3) + y2 ^ (1 / 3))) / (3 * a)'σ2 = (-b + (y1 ^ (1 / 3) + y2 ^ (1 / 3)) / 2 + Sqr(3) / 2 * (y1 ^ (1 / 3) - y2 ^ (1 / 3)) * i) / (3 * a)'σ3 = (-b + (y1 ^ (1 / 3) + y2 ^ (1 / 3)) / 2 - Sqr(3) / 2 * (y1 ^ (1 / 3) - y2 ^ (1 / 3)) * i) / (3 * a)ElseIf (0 = derta) Thenσ1 = -bb / aa + b1 / a1σ2 = -b1 / (a1 * 2)σ3 = -b1 / (a1 * 2)ElseIf (derta < 0) ThenDim θ As Doubleθ = Arccos((2 * a1 * bb - 3 * aa * b1) / (2 * Sqr(a1 * a1 * a1))) σ1 = (-bb - 2 * Sqr(a1) * Cos(θ / 3)) / (3 * aa)σ2 = (-bb + Sqr(a1) * (Cos(θ / 3) + Sqr(3) * Sin(θ / 3))) / (3 * aa)σ3 = (-bb + Sqr(a1) * (Cos(θ / 3) - Sqr(3) * Sin(θ / 3))) / (3 * aa) End IfIf σ2 > σ1 Thentemp = σ1σ1 = σ2σ2 = tempEnd IfIf σ3 > σ1 Thentemp = σ1σ1 = σ3σ3 = tempEnd IfIf σ3 > σ2 Thentemp = σ2σ2 = σ3σ3 = tempEnd Ifτmax = (σ1 - σ3) / 2Text8.Text = σ1Text9.Text = σ2Text10.Text = σ3Text11.Text = τmaxElseMsgBox "请选择应力状态"End IfEnd Sub3、截面几何性质分析窗口Option ExplicitDim a, b, h As DoubleDim d, al As DoubleDim iyc, izc, s, yc, zc As Double Const pi = 3.Private Sub Command1_Click()Command5.SetFocusIf Combo1.Text = "三角形" ThenFrame1.Visible = TrueImage1.Visible = TrueImage2.Visible = FalseFrame3.Visible = FalseImage3.Visible = FalseFrame4.Visible = FalseImage5.Visible = FalseFrame5.Visible = FalseImage6.Visible = FalseFrame6.Visible = FalseElseIf Combo1.Text = "矩形" Then Frame1.Visible = TrueImage2.Visible = TrueImage1.Visible = FalseFrame3.Visible = FalseImage3.Visible = FalseFrame4.Visible = FalseImage4.Visible = FalseImage5.Visible = FalseFrame5.Visible = FalseImage6.Visible = FalseFrame6.Visible = FalseElseIf Combo1.Text = "圆形" Then Frame3.Visible = TrueImage3.Visible = TrueImage2.Visible = FalseImage1.Visible = FalseFrame1.Visible = FalseImage4.Visible = FalseImage5.Visible = FalseFrame5.Visible = FalseImage6.Visible = FalseFrame6.Visible = FalseElseIf Combo1.Text = "梯形" Then Image5.Visible = TrueFrame5.Visible = TrueFrame3.Visible = FalseImage3.Visible = FalseImage2.Visible = FalseImage1.Visible = FalseFrame1.Visible = FalseFrame4.Visible = FalseImage4.Visible = FalseImage6.Visible = FalseFrame6.Visible = FalseElseIf Combo1.Text = "椭圆形" Then Image6.Visible = TrueFrame6.Visible = TrueImage1.Visible = FalseFrame1.Visible = FalseImage3.Visible = FalseFrame3.Visible = FalseImage4.Visible = FalseFrame4.Visible = FalseImage5.Visible = FalseFrame5.Visible = FalseElseIf Combo1.Text = "扇形" Then Frame4.Visible = TrueImage4.Visible = TrueImage5.Visible = FalseFrame3.Visible = FalseImage3.Visible = FalseImage2.Visible = FalseImage1.Visible = FalseFrame1.Visible = FalseFrame5.Visible = FalseImage6.Visible = FalseFrame6.Visible = FalseElseIf Combo1.Text = "组合图形" ThenElseMsgBox "请选择截面形状"End SubPrivate Sub Command2_Click() Text1.Text = "0"Text2.Text = "0"Text3.Text = "0"Text4.Text = "0"Text5.Text = "0"Text6.Text = "0"Text7.Text = "0"Text8.Text = "0"Text9.Text = "0"Text10.Text = "0"Text11.Text = "0"Text12.Text = "0"Text13.Text = "0"Text14.Text = "0"Text15.Text = "0"Set Picture1.Picture = Nothing End SubPrivate Sub Command3_Click() jihe.Hidemain.ShowPrivate Sub Command4_Click()EndEnd SubPrivate Sub Command5_Click()Set Picture1.Picture = NothingIf Combo1.Text = "三角形" Thenb = Val(Text1.Text)h = Val(Text2.Text)If b <= 0 Or h <= 0 ThenMsgBox "输入错误"Elses = 0.5 * b * hyc = b / 3zc = h / 3iyc = b * h ^ 3 / 36izc = h * b ^ 3 / 36Text3.Text = ycText4.Text = zcText5.Text = iycText6.Text = iycText7.Text = sPicture1.Scale (-(b + h) / 2, 3 * (b + h) / 2)-(3 * (b + h) / 2, -(b + h) / 2)Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (-(b + h) / 2, 0)-(3 * (b + h) / 2, 0)Picture1.Line (0, -(b + h) / 2)-(0, 3 * (b + h) / 2)Picture1.DrawWidth = 2Picture1.Line (0, 0)-(0, h)Picture1.Line (0, h)-(b, 0)Picture1.Line (b, 0)-(0, 0)Picture1.DrawWidth = 5Picture1.PSet (yc, zc), vbRedPicture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.CurrentX = yc: Picture1.CurrentY = zcPicture1.Print "C"Picture1.CurrentX = yc: Picture1.CurrentY = 3 * h / 2 Picture1.Print "zc"Picture1.CurrentX = 3 * b / 2: Picture1.CurrentY = zc Picture1.Print "yc"Picture1.CurrentX = 0: Picture1.CurrentY = 3 * (h + b) / 2 Picture1.Print "z"Picture1.CurrentX = 17 * (h + b) / 12: Picture1.CurrentY = 0 Picture1.Print "y"Picture1.DrawWidth = 1Picture1.Line (yc, zc)-(yc, 3 * h / 2)Picture1.Line (yc, zc)-(3 * b / 2, zc)End IfElseIf Combo1.Text = "矩形" Thenb = Val(Text1.Text)h = Val(Text2.Text)If b <= 0 Or h <= 0 ThenMsgBox "输入错误"Elses = b * hyc = b / 2zc = h / 2iyc = b * h ^ 3 / 12izc = h * b ^ 3 / 12Text3.Text = ycText4.Text = zcText5.Text = iycText6.Text = iycText7.Text = sPicture1.Scale (-(b + h) / 2, 3 * (b + h) / 2)-(3 * (b + h) / 2, -(b + h) / 2)Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (-(b + h) / 2, 0)-(3 * (b + h) / 2, 0)Picture1.Line (0, -(b + h) / 2)-(0, 3 * (b + h) / 2)Picture1.DrawWidth = 2Picture1.Line (0, 0)-(0, h)Picture1.Line (0, h)-(b, h)Picture1.Line (b, h)-(b, 0)Picture1.Line (b, 0)-(0, 0)Picture1.DrawWidth = 5Picture1.PSet (yc, zc), vbRedPicture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.CurrentX = yc: Picture1.CurrentY = zcPicture1.Print "C"Picture1.CurrentX = yc: Picture1.CurrentY = 3 * h / 2Picture1.Print "zc"Picture1.CurrentX = 3 * b / 2: Picture1.CurrentY = zcPicture1.Print "yc"Picture1.CurrentX = 0: Picture1.CurrentY = 3 * (h + b) / 2Picture1.Print "z"Picture1.CurrentX = 17 * (h + b) / 12: Picture1.CurrentY = 0Picture1.Print "y"Picture1.DrawWidth = 1Picture1.Line (yc, zc)-(yc, 3 * h / 2)Picture1.Line (yc, zc)-(3 * b / 2, zc)End IfElseIf Combo1.Text = "圆形" Thend = Val(Text8.Text)If d <= 0 ThenMsgBox "输入错误"Elses = pi * d ^ 2 / 4yc = 0zc = d / 2iyc = pi * d ^ 4 / 64izc = pi * d ^ 4 / 64Text3.Text = ycText4.Text = zcText5.Text = iycText6.Text = iycText7.Text = sPicture1.Scale (-d, 3 * d / 2)-(d, -d / 2)Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (-d, 0)-(d, 0)Picture1.Line (0, 3 * d / 2)-(0, -d / 2)Picture1.DrawWidth = 2Picture1.Circle (0, d / 2), d / 2Picture1.DrawWidth = 5Picture1.PSet (yc, zc), vbRedPicture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.CurrentX = yc: Picture1.CurrentY = zcPicture1.Print "C"Picture1.CurrentX = yc: Picture1.CurrentY = 17 * d / 12Picture1.Print "zc"Picture1.CurrentX = 0.9 * d: Picture1.CurrentY = zcPicture1.Print "yc"Picture1.CurrentX = 0: Picture1.CurrentY = 3 * d / 2Picture1.Print "z"Picture1.CurrentX = 0.9 * d: Picture1.CurrentY = 0Picture1.Print "y"Picture1.DrawWidth = 1Picture1.Line (yc, zc)-(yc, 17 * d / 12)Picture1.Line (yc, zc)-(0.9 * d, zc)End IfElseIf Combo1.Text = "梯形" Thena = Val(Text11.Text)b = Val(Text12.Text)h = Val(Text13.Text)If a <= 0 Or b <= 0 Or h <= 0 ThenMsgBox "输入错误"Elses = (a + b) * h / 2yc = (a ^ 2 + b ^ 2 + a * b) / (3 * (a + b))zc = (b + 2 * a) * h / (3 * (a + b))iyc = h ^ 3 * (b ^ 2 + 4 * a * b - a ^ 2) / (36 * (a + b))Text3.Text = ycText4.Text = zcText5.Text = iycText6.Text = ""Text7.Text = sPicture1.Scale (-(b + h) / 2, 3 * (b + h) / 2)-(3 * (b + h) / 2, -(b + h) / 2) Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (-(b + h) / 2, 0)-(3 * (b + h) / 2, 0)Picture1.Line (0, -(b + h) / 2)-(0, 3 * (b + h) / 2)Picture1.DrawWidth = 2Picture1.Line (0, 0)-(0, h)Picture1.Line (0, h)-(a, h)Picture1.Line (a, h)-(b, 0)Picture1.Line (b, 0)-(0, 0)Picture1.DrawWidth = 5Picture1.PSet (yc, zc), vbRedPicture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.CurrentX = yc: Picture1.CurrentY = zcPicture1.Print "C"Picture1.CurrentX = yc: Picture1.CurrentY = 3 * h / 2Picture1.Print "zc"Picture1.CurrentX = 3 * b / 2: Picture1.CurrentY = zcPicture1.Print "yc"Picture1.CurrentX = 0: Picture1.CurrentY = 3 * (h + b) / 2Picture1.Print "z"Picture1.CurrentX = 17 * (h + b) / 12: Picture1.CurrentY = 0Picture1.Print "y"Picture1.DrawWidth = 1Picture1.Line (yc, zc)-(yc, 3 * h / 2)Picture1.Line (yc, zc)-(3 * b / 2, zc)End IfElseIf Combo1.Text = "椭圆形" Thena = Val(Text14.Text)b = Val(Text15.Text)If a <= 0 Or b <= 0 ThenMsgBox "输入错误"Elses = pi * a * b / 4yc = 0zc = b / 2iyc = pi * a * b ^ 3 / 64Text3.Text = ycText4.Text = zcText5.Text = iycText6.Text = ""Text7.Text = sPicture1.Scale (-(b + a) / 2, 3 * (b + a) / 2)-(3 * (b + a) / 2, -(b + a) / 2) Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (-(b + a) / 2, 0)-(3 * (b + a) / 2, 0)Picture1.Line (0, -(b + a) / 2)-(0, 3 * (b + a) / 2)Picture1.DrawWidth = 2If a >= b ThenPicture1.Circle (0, b / 2), a / 2, , , , b / aElsePicture1.Circle (0, b / 2), b / 2, , , , b / aEnd IfPicture1.DrawWidth = 5Picture1.PSet (yc, zc), vbRedPicture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.CurrentX = yc: Picture1.CurrentY = zcPicture1.Print "C"Picture1.CurrentX = yc: Picture1.CurrentY = 3 * b / 2Picture1.Print "zc"Picture1.CurrentX = 3 * a / 2: Picture1.CurrentY = zcPicture1.Print "yc"Picture1.CurrentX = 0: Picture1.CurrentY = 3 * (a + b) / 2Picture1.Print "z"Picture1.CurrentX = 17 * (a + b) / 12: Picture1.CurrentY = 0Picture1.Print "y"Picture1.DrawWidth = 1Picture1.Line (yc, zc)-(yc, 3 * b / 2)Picture1.Line (yc, zc)-(3 * a / 2, zc)End IfElseIf Combo1.Text = "扇形" Thend = Val(Text9.Text)al = Val(Text10.Text)If al > pi / 2 Or al <= 0 Or d <= 0 ThenMsgBox "输入错误"Elses = al * d ^ 2 / 4yc = 0zc = d * Sin(al) / 3 / aliyc = d ^ 4 / 64 * (al + Sin(al) * Cos(al) - 16 * Sin(al) * Sin(al) / 9 / al)Text3.Text = ycText4.Text = zcText5.Text = iycText6.Text = ""Text7.Text = sPicture1.Scale (-d, 3 * d / 2)-(d, -d / 2)Picture1.DrawStyle = 0Picture1.DrawWidth = 1Picture1.Line (-d, 0)-(d, 0)Picture1.Line (0, 3 * d / 2)-(0, -d / 2)Picture1.DrawWidth = 2Picture1.Circle (0, 0), d / 2, , (pi / 2 - al), (pi / 2 + al) Picture1.Line (0, 0)-(d / 2 * Sin(al), d / 2 * Cos(al)) Picture1.Line (0, 0)-(-d / 2 * Sin(al), d / 2 * Cos(al)) Picture1.DrawWidth = 5Picture1.PSet (yc, zc), vbRedPicture1.CurrentX = 0: Picture1.CurrentY = 0Picture1.Print "O"Picture1.CurrentX = yc: Picture1.CurrentY = zc Picture1.Print "C"Picture1.CurrentX = yc: Picture1.CurrentY = 17 * d / 12 Picture1.Print "zc"Picture1.CurrentX = 0.9 * d: Picture1.CurrentY = zc Picture1.Print "yc"Picture1.CurrentX = 0: Picture1.CurrentY = 3 * d / 2 Picture1.Print "z"Picture1.CurrentX = 0.9 * d: Picture1.CurrentY = 0 Picture1.Print "y"Picture1.DrawWidth = 1Picture1.Line (yc, zc)-(yc, 17 * d / 12)Picture1.Line (yc, zc)-(0.9 * d, zc)End IfElseIf Combo1.Text = "组合图形" ThenElseMsgBox "请选择截面形状"End IfEnd Sub五、具体问题的计算结果和结果分析（对或者不对）下面对各计算程序进行应用、验证和运行结果的分析：1、平面应力状态分析已知应力状态：σx=56.92MPa,σy=34.63MPa,τxy=66.64MPa，α=52.99°求指定截面的应力，并求出主应力和最大切应力。