悬臂梁有限元优化程序.doc

悬臂梁优化分析 班级： 姓名： 学号： 指导老师： 目录 一、 条件分析 1 二、 分析步骤 1 （一）前处理阶段: 1 （二）求解阶段 3 （三）后处理阶段 4 （四）优化阶段 9 三、 优化结果 13 （一）读取优化结果列表 13 （二）选择优化结果 13 （三） 代入结果分析 14 四、整理命令流 14 参考文献 17 条件分析 由题可知：悬臂梁中的平均应力小于，且梁的挠度小于1厘米。而且横截面积约束条件为：，。（考虑学号系数），连杆的材料属性为：杨氏模量,泊松比为0.3。由于梁的长度一定，若要使梁的重量最小，则要求体积最小，进而可知要求横截面积，所以可确定体积是所求目标，因此可确定： 设计变量 状态变量 平均应力 挠度 目标函数 体积V 分析步骤 1. 定义工作文件名和工作标题 (1) 执行[Utility Menu]\File\change Jobname。弹出对话框，输入panjiafeng12，并选择复选框，单击“OK”按钮。 (2) 执行[Utility Menu]\File\Change Title。弹出的对话框，输入panjiafeng12，单击“OK”按钮。 （一）前处理阶段: 1. 初始化设计变量 执行[Utility Menu]\File\Parameters\Scalar Parameter,弹出对话框，输入X1=0.1cm，X2=0.3cm。 2.定义单元类型，面积，转动惯量 执行[Utility Menu]\Preprocessor\Element Type\Add\Edit\Delete 弹出对话框，选择Structural Beam中的2D elastic 3 单击“OK”单击“Close”。 3.执行[Utility Menu]\File\Parameters\Scalar Parameter,弹出对话框，输入AREA=X1*X2,IZZ=(X1*(X2**3))/12,定义梁的横截面积和惯性矩 4.定义实常数和材料属性 (1) 执行[Utility Menu]\Preprocessor\Real Constants\Add\Edit\Delete,弹出如图所示的对话框，分别输入AREA，IZZ，X2。 执行[Utility Menu]\Preprocessor\Material Props\MaterialModels\Structural\Linear\Elastic\Isotropic\弹出如图所示的对话框，输入和0.3，点击 OK和Close。 5.建立几何模型 （1）执行[Utility Menu]\Preprocessor\Modeling\Create\Nodes\In Active CS,定义三个节点如图所示，坐标分别为（0，0，0）、（2.5，0，0）、（5，0，0） （2）执行[Utility Menu]\Preprocessor\Modeling\Create\Elements\Auto Numbered\Thru Nodes\点击1(0,0),2(2.5,0) → Apply →点击2(2.5,0),3(5,0) 最后点击 OK ，定义单元1包含节点1，2。定义单元2包含节点2，3，如图所示 （二）求解阶段 1. 执行[Utility Menu]Solution\Define Loads\Apply\Structural\Displacement\On Nodes\pick the 点1（0，0）点击OK 点击select ?ALL DOF点击OK；和\Structural\Force/Moment\On Nodes选择点2（2.5，0）点击OK →select ?FY 500 → Apply 选择点3（5，0）点击OK →select ?FY 500 点击OK。如图所示，在节点1处施加约束，所有位移均为0。在节点2和节点3处施加Y方向集中载荷500N。 2. 执行[Utility Menu]Solution\Solve\Current LS点击OK，求解 （三）后处理阶段 执行[Utility Menu]General Postproc\List Results\Sorted Li