Mathematica入门教程 第1篇 第1章MATHEMATICA概述 (3) 1.1 M ATHEMATICA的启动与运行 (3) 1.2 表达式的输入 (4) 1.3 M ATHEMATICA的联机帮助系统 (6) 第2章MATHEMATICA的基本量 (8) 2.1 数据类型和常数 (8) 2.2 变量 (10) 2.3 函数 (11) 2.4 表 (14) 2.5 表达式 (17) 2.6 常用的符号 (19) 2.7 练习题 (19) 第2篇 第3章微积分的基本操作 (20) 3.1 极限 (20) 3.2 微分 (20) 3.3 计算积分 (22) 3.4 无穷级数 (24) 3.5 练习题 (24) 第4章微分方程的求解 (26) 4.1 微分方程解 (26) 4.2 微分方程的数值解 (26) 4.3 练习题 (27) 第3篇 第5章MATHEMATICA的基本运算 (28) 5.1 多项式的表示形式 (28) 5.2 方程及其根的表示 (29) 5.3 求和与求积 (32) 5.4 练习题 (33) 第6章函数作图 (35) 6.1 基本的二维图形 (35) 6.2 二维图形元素 (40) 6.3 基本三维图形 (42) 6.4 练习题 (46)
第4篇 第7章MATHEMATICA函数大全 (48) 7.1 运算符和一些特殊符号，系统常数 (48) 7.2 代数计算 (49) 7.3 解方程 (50) 7.4 微积分 (50) 7.5 多项式函数 (51) 7.6 随机函数 (52) 7.7 数值函数 (52) 7.8 表相关函数 (53) 7.9 绘图函数 (54) 7.10 流程控制 (57) 第8章MATHEMATICA程序设计 (59) 8.1 模块和块中的变量 (59) 8.2 条件结构 (61) 8.3 循环结构 (63) 8.4 流程控制 (65) 8.5 练习题 (67) --------------习题与答案在68页-------------------
弹性力学简明教程（第四版）课后习题解答 徐芝纶 第一章绪论 1、试举例说明什么是均匀的各向异性体，什么是非均匀的各向同性体？ 【分析】均匀的各项异形体就是满足均匀性假定，但不满足各向同性假定；非均匀的各 向异性体，就是不满足均匀性假定，但满足各向同性假定。 【解答】均匀的各项异形体如：竹材，木材。 非均匀的各向同性体如：混凝土。 1.2 一般的混凝土构件和钢筋混凝土构件能否作为理想弹性体？一般的岩质地基和土质地基能否作为理想弹性体？ 【分析】能否作为理想弹性体，要判定能否满足四个假定：连续性，完全弹性，均匀性, 各向同性假定。 【解答】一般的混凝土构件和土质地基可以作为理想弹性体；一般的钢筋混凝土构件和 岩质地基不可以作为理想弹性体。 1.3五个基本假定在建立弹性力学基本方程时有什么作用？ 【解答】（1）连续性假定：假定物体是连续的，也就是假定整个物体的体积都被组成这个物体的介质所填满，不留下任何空隙。引用这一假定后，物体的应力、形变和位移等物理 量就可以看成是连续的。因此，建立弹性力学的基本方程时就可以用坐标的连续函数来表示他们的变化规律。 完全弹性假定：假定物体是完全弹性的，即物体在对应形变的外力被去除后，能够完全恢复原型而无任何形变。这一假定，还包含形变与引起形变的应力成正比的涵义，亦即两者之间是成线性关系的，即引用这一假定后，应力与形变服从胡克定律，从而使物理方程成为线性的方程，其弹性常数不随应力或形变的大小而变。 均匀性假定：假定物体是均匀的，即整个物体是由同一材料组成的，引用这一假定后整个物体的所有各部分才具有相同的弹性，所研究物体的内部各质点的物理性质都是相同的，因而物体的弹性常数不随位置坐标而变化。 各向同性假定：假定物体是各向同性的，即物体的弹性在所有各个方向都相同，引用此假定后，物体的弹性常数不随方向而变。 小变形假定：假定位移和变形是微小的。亦即，假定物体受力以后整个物体所有各点的 位移都远远小于物体原来的尺寸，而且应变和转角都远小于1。这样在建立物体变形以后的 平衡方程时，就可以方便的用变形以前的尺寸来代替变形以后的尺寸。在考察物体的位移与 形变的关系时，它们的二次幕或乘积相对于其本身都可以略去不计，使得弹性力学中的微分
【2-9】试列出图2-17，图2-18所示问题的全部边界条件。在其端部小边界上，应用圣维南原理列出三个积分的应力边界条件。 x y 2 h 1h b g ρo () 2h b >> h x y l /2/2 h M N F S F 1 q q 图2-17 图2-18 【分析】有约束的边界上可考虑采用位移边界条件，若为小边界也可写成圣维南原理的三个积分形式，大边界上应精确满足公式（2-15）。 【解答】图2-17： 上（y =0） 左(x =0) 右（x =b ） l 0 -1 1 m -1 () x f s () 1g y h ρ+ () 1g y h ρ-+ () y f s 1gh ρ 代入公式（2-15）得 ①在主要边界上x=0，x=b 上精确满足应力边界条件： ()()100(),0;===-+=x xy x x g y h σρτ ()()1b b (),0; ===-+=x xy x x g y h σρτ ②在小边界0y =上，能精确满足下列应力边界条件： () () ,0y xy y y gh σρτ===-= ③在小边界2y h =上，能精确满足下列位移边界条件： ()()2 2 0,0 ====y h y h u v 这两个位移边界条件可以应用圣维南原理，改用三个积分的应力边界条件来代替，当板厚=1δ时，可求得固定端约束反力分别为： 10,,0s N F F gh b M ρ==-=
由于2y h =为正面，故应力分量与面力分量同号，则有： ()()()222 10000 0b y y h b y y h b xy y h dx gh b xdx dx σρστ===?=-???=???=????? ⑵图2-18 ①上下主要边界y=-h/2，y=h/2上，应精确满足公式（2-15） l m x f (s) y f (s) 2h y =- 0 -1 0 q 2 h y = 1 -1q -/2()y y h q σ==-，-/2()0yx y h τ==，/2()0y y h σ==，/21()yx y h q τ==- ②在x =0的小边界上，应用圣维南原理，列出三个积分的应力边界条件：负面上应力与面力符号相反，有 /20/2/2 0/2/20 /2()()()h xy x S h h x x N h h x x h dx F dx F ydx M τσσ=-=-=-?=-??=-???=-???? ③在x=l 的小边界上，可应用位移边界条件0,0====l x l x v u 这两个位移边界条件也可改用三个积分的应力边界条件来代替。 首先，求固定端约束反力，按面力正方向假设画反力，如图所示，列平衡方程求反力： 110,x N N N N F F F q l F q l F ''=+=?=-∑ 0,0y S S S S F F F ql F ql F ''=++=?=--∑ 2 211110,'02222 A S S q lh ql M M M F l ql q lh M M F l =+++-=?=---∑ 由于x=l 为正面，应力分量与面力分量同号，故 M ' N F 'S F '
第一章Mathematica的启动的运行 Mathematica是美国Wolfram公司生产的一种数学分析型的软件，以符号计算见长，也具有高精度的数值计算功能和强大的图形功能。目前最新版本是Mathematica4.0，本附录仅介绍Mathematica4.0的一些常用功能，须深入掌握Mathematica的读者可查阅相关书籍。 在Windows环境下安装好Mathematica4.0，用鼠标双击Mathematica图标（刺球状），在显示器上显示如图1-1的工作窗口，这时可以键入你想计算的东西，比如键入1+1，然后同时按下Shift键和Enter键(数字键盘上只要按Enter键)，这时Mathematica开始工作，计算出结果后，窗口变为图1-2。 图1-1 Mathematica的工作窗口 Mathematica第一次计算时因为要启动核(kernel)，所需时间要长一些，也可以在Mathematica 启动后第一次计算之前，手工启动核，方法是用鼠标点击：Kernel->Start Kernel->Local.这样第一次计算就很快了。
图1-2 完成运算后的Mathematica的窗口 图1-2中的“In[1]:=”表示第一个输入；“Out[1]=”表示第一个输出结果。接下来可键入第二个输入，按这样的方式可利用Mathematica进行“会话式”计算。要注意的是：“In[1]:= ”和“Out[1]=”是系统自动添加的，不需用户键入。Mathematica还提供“批处理”运行方式，即可以将Mathematica作为一种算法语言，编写程序，让计算机执行，这在第七章将会作简要介绍。 第二章 Mathematica的基本运算功能 2.1 算术运算 Mathematica最基本的功能是进行算术运算，包括加（+），减（-），乘（*），除（/），乘方（^），阶乘（！）等。 注意： 1 在Mathematica中，也可用空格代表乘号；数字和字母相乘，乘号可以省去，例如：3*2可写成3 2，2*x可写成2x，但字母和字母相乘，乘号不能省去。建议大家尽可能不要省去乘号，以免引起混乱。 2 在Mathematica中，表达式中用来表示运算的结合次序的括号只允许是圆括号（无论多少层）。例如：4*（2+3/（2-5）） 3 当输入式子中不含小数点，输出结果是完全精确的。例如：输入2/3，输出仍然为2/3。
《弹性力学简明教程》是教育部“十五”国家规划教材。是在第二版的基础上，保持原有的体系和特点，根据教学改革的需要和国家的有关新标准，进行了修订。全书按照由浅入深的原则，安排了平面问题的理论及解答、空间问题的理论及解答和薄板弯曲理论。并着重介绍了弹性力学的数值解法，即差分法、变分法和有限单元法。《弹性力学简明教程》作为弹性力学的入门教材，注重基本理论（基本概念、基本方程和基本解法）的阐述及其应用，以使学生在掌握基本理论的基础上能阅读和应用弹性力学文献，并能初步应用弹性力学的数值解法解决工程实际问题。 主要符号表 第一章绪论1-1 弹性力学的内容1-2 弹性力学中的几个基本概念1-3 弹性力学中的基本假定习题 第二章平面问题的基本理论2-1 平面应力问题与平面应变问题2-2 平衡微分方程2-3 平面问题中一点的应力状态2-4 几何方程刚体位移2-5 物理方程2-6 边界条件2-7 圣维南原理及其应用2-8 按位移求解平面问题2-9 按应力求解平面问题相容方程 2-10 常体力情况下的简化应力函数习题 第三章平面问题的直角坐标解答3-1 逆解法与半逆解法多项式解答 .3-2 矩形梁的纯弯曲3-3 位移分量的求出3-4 简支梁受均布荷载3-5 楔形体受重力和液体压力习题
第四章平面问题的极坐标解答4-1 极坐标中的平衡微分方程4-2 极坐标中的几何方程及物理方程4-3 极坐标中的应力函数与相容方程4-4 应力分量的坐标变换式4-5 轴对称应力和相应的位移4-6 圆环或圆筒受均布压力4-7 压力隧洞4-8 圆孔的孔口应力集中4-9 半平面体在边界上受集中力 4-10 半平面体在边界上受分布力习题 第五章用差分法和变分法解平面问题5-1 差分公式的推导 5-2 应力函数的差分解5-3 应力函数差分解的实例5-4 弹性体的形变势能和外力势能5-5 位移变分方程5-6 位移变分法5-7 位移变分法的例题习题.. 第六章用有限单元法解平面问题6-1 基本量及基本方程的矩阵表示6-2 有限单元法的概念6-3 单元的位移模式与解答的收敛性6-4 单元的应变列阵和应力列阵6-5 单元的结点力列阵与劲度矩阵6-6 荷载向结点移置单元的结点荷载列阵6-7 结构的整体分析结点平衡方程组6-8 解题的具体步骤单元的划分6-9 计算成果的整理6-10 计算实例6-11 应用变分原理导出有限单元法基本方程习题 第七章空间问题的基本理论7-1 平衡微分方程7-2 物体内任一点的应力状态7-3 主应力最大与最小的应力7-4 几何方程及物理方程7-5 轴对称问题的基本方程习题
Mathematica 教程 【Mathematica 简介】 Mathematica 软件是由沃尔夫勒姆研究公司 (Wolfram Research Inc.)研发的。Mathematica 1.0 版发布于1988年6月23日。发布之后，在科学、技术、媒体等领域引起了一片轰动，被认为是一个革命性的进步。几个月后，Mathematica就在世界各地拥有了成千上万的用户。今天，Mathematica 已经在世界各地拥有了数以百万计的忠实用户。 Mathematica已经被工业和教育领域被广泛地采用。实际上，Mathematica负责将高级的数 学和计算引入了传统上非技术的领域，极大的增加了科技软件的市场。一个包含应用、咨询、 书籍、和课程软件的行业支持着国际化的Mathematica用户群，这个行业还在不断地膨胀。 随着沃尔夫勒姆研究公司不断地扩大和Mathematica的使用被不断地扩展到不同的领域， 将会看到Mathematica在全世界范围内对未来产品、重要研究发现、和教学的巨大影响。 数学软件是现在科研工作者的必备的工具，个人比较喜欢用Mathematica，因为它是最接近数学语言的。Mathematica在15日发布，其最显著的变化是允许自由形式的英文输入，而不再需要严格按照Mathematica语法，这类似于Wolfram|Alpha搜索引擎。Mathematica 8 允许用户按照自己习惯的思考过程输入方程式或问题，最令人激动的部分是软件不是逐行执 行命令，而是能理解上下文背景。 1. En ter your queries in pla in En glish using new free-form lin guistic in put 2. Access more tha n 10 trilli on sets of curated, up-to-date, and ready-to-use data 3. Import all your data using a wider array of import/export formats 4. Use the broadest statistics and data visualizati on capabilities on the market 5. Choose from a full suite of engin eeri ng tools, such as wavelets and con trol systems 6. Use more powerful image process ing and an alysis capabilities 7. Create in teractive tools for rapid explorati on of your ideas 8. Develop faster and more powerful applicati ons
弹性力学简明教程（第四版） 习题解答 第一章绪论 【1-1】试举例说明什么是均匀的各向异性体，什么是非均匀的各向同性体？ 【分析】均匀的各项异形体就是满足均匀性假定，但不满足各向同性假定；非均匀的各向异性体，就是不满足均匀性假定，但满足各向同性假定。 【解答】均匀的各项异形体如：竹材，木材。 非均匀的各向同性体如：混凝土。 【1-2】一般的混凝土构件和钢筋混凝土构件能否作为理想弹性体？一般的岩质地基和土质地基能否作为理想弹性体？ 【分析】能否作为理想弹性体，要判定能否满足四个假定：连续性，完全弹性，均匀性，各向同性假定。 【解答】一般的混凝土构件和土质地基可以作为理想弹性体；一般的钢筋混凝土构件和岩质地基不可以作为理想弹性体。 【1-3】五个基本假定在建立弹性力学基本方程时有什么作用？ 【解答】（1）连续性假定：假定物体是连续的，也就是假定整个物体的体积都被组成这个物体的介质所填满，不留下任何空隙。引用这一假定后，物体的应力、形变和位移等物理量就可以看成是连续的。因此，建立弹性力学的基本方程时就可以用坐标的连续函数来表示他们的变化规律。 完全弹性假定：假定物体是完全弹性的，即物体在对应形变的外力被去除后，能够完全恢复原型而无任何形变。这一假定，还包含形变与引起形变的应力成正比的涵义，亦即两者之间是成线性关系的，即引用这一假定后，应力与形变服从胡克定律，从而使物理方程成为线性的方程，其弹性常数不随应力或形变的大小而变。 均匀性假定：假定物体是均匀的，即整个物体是由同一材料组成的，引用这一假定后整个物体的所有各部分才具有相同的弹性，所研究物体的内部各质点的物理性质都是相同的，因而物体的弹性常数不随位置坐标而变化。 各向同性假定：假定物体是各向同性的，即物体的弹性在所有各个方向都相同，引用此假定后，物体的弹性常数不随方向而变。 小变形假定：假定位移和变形是微小的。亦即，假定物体受力以后整个物体所有各点的位移都远远小于物体原来的尺寸，而且应变和转角都远小于1。这样在建立物体变形以后的平衡方程时，就可以方便的用变形以前的尺寸来代替变形
【Mathematica 简介】 Mathematica 软件是由沃尔夫勒姆研究公司（Wolfram Research Inc.）研发的。Mathematica 版发布于1988年6月23日。发布之后，在科学、技术、媒体等领域引起了一片轰动，被认为是一个革命性的进步。几个月后，Mathematica 就在世界各地拥有了成千上万的用户。今天，Mathematica 已经在世界各地拥有了数以百万计的忠实用户。 Mathematica 已经被工业和教育领域被广泛地采用。实际上，Mathematica 负责将高级的数学和计算引入了传统上非技术的领域，极大的增加了科技软件的市场。一个包含应用、咨询、书籍、和课程软件的行业支持着国际化的 Mathematica 用户群，这个行业还在不断地膨胀。随着沃尔夫勒姆研究公司不断地扩大和 Mathematica 的使用被不断地扩展到不同的领域，将会看到 Mathematica 在全世界范围内对未来产品、重要研究发现、和教学的巨大影响。 数学软件是现在科研工作者的必备的工具，个人比较喜欢用Mathematica，因为它是最接近数学语言的。Mathematica 在15日发布，其最显著的变化是允许自由形式的英文输入，而不再需要严格按照Mathematica语法，这类似于Wolfram|Alpha搜索引擎。Mathematica 8允许用户按照自己习惯的思考过程输入方程式或问题，最令人激动的部分是软件不是逐行执行命令，而是能理解上下文背景。 1. Enter your queries in plain English using new free-form linguistic input 2. Access more than 10 trillion sets of curated, up-to-date, and ready-to-use data 3. Import all your data using a wider array of import/export formats 4. Use the broadest statistics and data visualization capabilities on the market 5. Choose from a full suite of engineering tools, such as wavelets and control systems 6. Use more powerful image processing and analysis capabilities 7. Create interactive tools for rapid exploration of your ideas 8. Develop faster and more powerful applications Wolfram Research 的 CEO 和创立者斯蒂芬·沃尔夫勒姆表示：“传统上，让计算机执行任务必须使用计算机语言或者使用点击式界面：前者要求用户掌握它的语法；而后者则限制了可访问函数的范围。”“自由格式语言学能够理解人类的语言，并将其转化为具有特定语法结构的语言。这是产品适用性上的一个突破。 Mathematica 8 是这种创新思想下的第一个产品，但是它已经能够大幅度提高用户的工作效率。” Mathematica简明教程 第1章Mathematica概述 运行和启动：介绍如何启动Mathematica软件，如何输入并运行命令
【2-9】【解答】图2-17： 上（y =0） 左(x =0) 右（x =b ） l -1 1 m -1 () x f s () 1g y h ρ+ () 1g y h ρ-+ () y f s 1gh ρ 代入公式（2-15）得 ①在主要边界上x=0，x=b 上精确满足应力边界条件： ()()100(),0;===-+=x xy x x g y h σρτ()()1b b (),0; ===-+=x xy x x g y h σρτ ②在小边界0y =上，能精确满足下列应力边界条件：() () ,0y xy y y gh σρτ===-= ③在小边界2y h =上，能精确满足下列位移边界条件：()()2 2 0,0 ====y h y h u v 这两个位移边界条件可以应用圣维南原理，改用三个积分的应力边界条件来代替，当板 厚=1δ时，可求得固定端约束反力分别为： 10,,0s N F F gh b M ρ==-= 由于2y h =为正面，故应力分量与面力分量同号，则有： ()()()22210000 0b y y h b y y h b xy y h dx gh b xdx dx σρστ===?=-???=???=?? ??? ⑵图2-18 ①上下主要边界y=-h/2，y=h/2上，应精确满足公式（2-15） l m x f (s) y f (s) 2h y =- 0 -1 0 q 2 h y = 1 -1q -/2()y y h q σ==-，-/2()0yx y h τ==，/2()0y y h σ==，/21()yx y h q τ==- ②在x =0的小边界上，应用圣维南原理，列出三个积分的应力边界条件：负面上应力
Mathematica入门教程 Mathematica的基本语法特征 如果你是第一次使用Mathematica，那么以下几点请你一定牢牢记住： Mathematica中大写小写是有区别的，如Name、name、NAME等是不同的变量名或函数名。 系统所提供的功能大部分以系统函数的形式给出，内部函数一般写全称，而且一定是以大写英文字母开头，如Sin[x],Conjugate[z]等。 乘法即可以用*，又可以用空格表示，如2 3＝2*3＝6 ,x y,2 Sin[x]等；乘幂可以用“^”表示，如x^0.5,Tan[x]^y。 自定义的变量可以取几乎任意的名称，长度不限，但不可以数字开头。 当你赋予变量任何一个值，除非你明显地改变该值或使用Clear[变量名]或“变量名=.”取消该值为止，它将始终保持原值不变。 一定要注意四种括号的用法：()圆括号表示项的结合顺序，如(x+(y^x+1/(2x)));[]方括号表示函数，如Log[x],BesselJ[x,1]；{}大括号表示一个“表”(一组数字、任意表达式、函数等的集合)，如{2x,Sin[12 Pi],{1+A,y*x}}；[[]]双方括号表示“表”或“表达式”的下标，如a[[2,3]]、{1,2,3}[[1]]=1。 Mathematica的语句书写十分方便，一个语句可以分为多行写，同一行可以写多个语句（但要以分号间隔）。当语句以分号结束时，语句计算后不做输出（输出语句除外），否则将输出计算的结果。 一.数的表示及计算 1.在Mathematica中你不必考虑数的精确度，因为除非你指定输出精度，Mathematica总会以绝对精确的形式输出结果。例如：你输入 In[1]:=378/123，系统会输出Out[1]:=126/41，如果想得到近似解，则应输入 In[2]:=N[378/123,5],即求其5位有效数字的数值解，系统会输出Out[2]:=3.073 2，另外Mathematica还可以根据你前面使用的数字的精度自动地设定精度。 Mathematica与众不同之处还在于它可以处理任意大、任意小及任意位精度的数值，如100^7000,2^(-2000)等数值可以很快地求出，但在其他语言或系统中这是不可想象的，你不妨试一试N[Pi,1000]。 Mathematica还定义了一些系统常数，如上面提到的Pi(圆周率的精确值)，还有E(自然对数的底数)、I(复数单位)，Degree(角度一度，Pi/180)，Infinity(无穷大)等，不要小看这些简单的符号，它们包含的信息远远大于我们所熟知的它们的近似值，它们的精度也是无限的。 二.“表”及其用法 “表”是Mathematica中一个相当有用的数据类型，它即可以作为数组，又可以作为矩阵；除此以外，你可以把任意一组表达式用一个或一组{}括起来，进行运算、存储。可以说表是任意对象的一个集合。它可以动态地分配内存，
【3-1】为什么在主要边界（大边界）上必须满足精确的应力边界条件式（2-15），而在小边界上可以应用圣维南原理，用三个积分的应力边界条件（即主矢量、主矩的条件）来代替？如果在主要边界上用三个积分的应力边界条件代替式（2-15），将会发生什么问题？ 【解答】弹性力学问题属于数学物理方程中的边值问题，而要使边界条件完全得到满足，往往比较困难。这时，圣维南原理可为简化局部边界上的应力边界条件提供很大的方便。将物体一小部分边界上的面力换成分布不同，但静力等效的面力（主矢、主矩均相同），只影响近处的应力分布，对远处的应力影响可以忽略不计。如果在占边界绝大部分的主要边界上用三个积分的应力边界条件来代替精确的应力边界条件（公式2-15），就会影响大部分区域的应力分布，会使问题的解答精度不足。 【3-2】如果在某一应力边界问题中，除了一个小边界条件，平衡微分方程和其它的应力边界条件都已满足，试证：在最后的这个小边界上，三个积分的应力边界条件必然是自然满足的，固而可以不必校核。 【解答】区域内的每一微小单元均满足平衡条件，应力边界条件实质上是边界上微分体的平衡条件，即外力（面力）与内力（应力）的平衡条件。研究对象整体的外力是满足平衡条件的，其它应力边界条件也都满足，那么在最后的这个次要边界上，三个积分的应力边界条件是自然满足的,因而可以不必校核。 【3-3】如果某一应力边界问题中有m 个主要边界和n 个小边界，试问在主要边界和小边界上各应满足什么类型的应力边界条件，各有几个条件？ 【解答】在m 个主要边界上，每个边界应有2个精确的应力边界条件，公式（2-15），共2m 个；在n 个次要边界上，如果能满足精确应力边界条件，则有2n 个；如果不能满足公式（2-15）的精确应力边界条件，则可以用三个静力等效的积分边界条件来代替2个精确应力边界条件，共3n 个。 【3-4】试考察应力函数3 ay Φ=在图3-8所示的矩形板和坐标系中能解决什么问题(体力不计)? 【解答】⑴相容条件: 不论系数a 取何值，应力函数3 ay Φ=总能满足应力函数表示的相容方程，式(2-25). ⑵求应力分量 当体力不计时，将应力函数Φ代入公式(2-24)，得 6,0,0x y xy yx ay σσττ==== ⑶考察边界条件 上下边界上应力分量均为零，故上下边界上无面力.
第四章 平面问题的极坐标解答 【4-8】 实心圆盘在r ρ=的周界上受有均布压力q 的作用，试导出其解答。 【解答】实心圆盘是轴对称的，可引用轴对称应力解答，教材中的式（4-11），即 2 2(12ln )2(32ln )20A B C A B C ρ?ρ? σρρσρρτ? =+++? ???=-+++?? ?? =?? (a) 首先，在圆盘的周界（r ρ=）上，有边界条件()=r q ρρσ=-，由此得 -q 2 (12ln )2A B C ρσρρ = +++= (b) 其次，在圆盘的圆心，当0ρ→时，式（a ）中ρσ，?σ的第一、第二项均趋于无限大，这是不可能的。按照有限值条件（即，除了应力集中点以外，弹性体上的应力应为有限值。），当=0ρ时，必须有0A B ==。 把上述条件代入式（b ）中，得 /2C q =-。 所以，得应力的解答为 -q 0ρ?ρ?σστ===。 【4-9】 半平面体表面受有均布水平力q ，试用应力函数 2(sin 2)ΦρB φC φ=+求解应力分量（图4-15）。 【解答】（1）相容条件： 将应力函数Φ代入相容方程40?Φ=，显然满足。 （2）由Φ求应力分量表达式 =-2sin 222sin 222cos 2B C B C B C ρ?ρ?σ?? σ??τ??+?? =+??=--??
（3）考察边界条件：注意本题有两个?面，即2 π ?=± ，分别为?±面。在?±面 上，应力符号以正面正向、负面负向为正。因此，有 2()0,??πσ=±= 得0C =； -q 2 (),ρ??πτ=±= 得2 q B =-。 将各系数代入应力分量表达式，得 sin 2sin 2cos 2q q q ρ?ρ?σ?σ?τ? ?=?? =-??=?? 【4-14】 设有内半径为r 而外半径为R 的圆筒受内压力q ，试求内半径和外半径的改 变量，并求圆筒厚度的改变量。 【解答】本题为轴对称问题，只有径向位移而无环向位移。当圆筒只受内压力q 的情况下，取应力分量表达式，教材中式（4-11），注意到B =0。 内外的应力边界条件要求 r r ()0,()0;(), ()0 R R q ρ?ρρ?ρρρρρττσσ=======-= 由表达式可见，前两个关于ρ?τ的条件是满足的，而后两个条件要求 r 2 22,20A C q A C R ?+=-??? ?+=??。 由上式解得 22 2 ,C () 2() 22 22 qr R qr A R -r R -r =-=。 (a) 把A ，B ，C 值代入轴对称应力状态下对应的位移分离，教材中式（4-12）。 ()()222211cos sin ,(R r )qr R u I K E ρμρμ??ρ?? =-++++??-? ? (b) sin cos 0u H I K ?ρ??=-+=。 (c) 式（c ）中的ρ，?取任何值等式都成立，所以各自由项的系数为零
《弹性力学简明教程》 习题提示和参考答案 第二章习题的提示与答案 2－1 是 2－2 是 2－3 按习题2－1分析。 2－4 按习题2－2分析。 2－5 在的条件中，将出现2、3阶微量。当略去3阶微量后，得出的切应力互等定理完全相同。 2－6 同上题。在平面问题中，考虑到3阶微量的精度时，所得出的平衡微分方程都相同。其区别只是在3阶微量（即更高阶微量）上，可以略去不计。 2－7 应用的基本假定是：平衡微分方程和几何方程─连续性和小变形，物理方程─理想弹性体。 2－8 在大边界上，应分别列出两个精确的边界条件；在小边界（即次要边界）上，按照圣维南原理可列出3个积分的近似边界条件来代替。 2－9 在小边界OA边上，对于图2－15（a）、（b）问题的三个积分边界条件相同，因此，这两个问题为静力等效。 2－10 参见本章小结。 2－11 参见本章小结。 2－12 参见本章小结。 2－13 注意按应力求解时，在单连体中应力分量必须满足 （1）平衡微分方程， （2）相容方程， （3）应力边界条件（假设)。 2－14 见教科书。 2－15 见教科书。 2－16 见教科书。 2－17 取 它们均满足平衡微分方程，相容方程及x=0和的应力边界条件，因此，它们是该问题的正确解答。 2－18 见教科书。 2－19 提示：求出任一点的位移分量和，及转动量，再令,便可得出。 第三章习题的提示与答案 3－1 本题属于逆解法，已经给出了应力函数，可按逆解法步骤求解： （1）校核相容条件是否满足， （2）求应力， （3）推求出每一边上的面力从而得出这个应力函数所能解决的问题。
3－2 用逆解法求解。由于本题中l>>h, x=0,l 属于次要边界（小边界），可将小边界上的面力化为主矢量和主矩表示。 3－3 见3-1例题。 3－4 本题也属于逆解法的问题。首先校核是否满足相容方程。再由求出应力后，并求对应的面力。本题的应力解答如习题3-10所示。应力对应的面力是： 主要边界： 所以在边界上无剪切面力作用。下边界无法向面力；上边界有向下的法向面力q。 次要边界： x=0面上无剪切面力作用；但其主矢量和主矩在x=0 面上均为零。 因此，本题可解决如习题3-10所示的问题。 3－5 按半逆解法步骤求解。 （1）可假设 （2）可推出 （3）代入相容方程可解出f、，得到 （4）由求应力。 （5）主要边界x=0,b上的条件为 次要边界y=0上，可应用圣维南原理，三个积分边界条件为 读者也可以按或的假设进行计算。 3－6 本题已给出了应力函数，应首先校核相容方程是否满足，然后再求应力，并考察边界条件。在各有两个应精确满足的边界条件，即 而在次要边界y=0 上，已满足，而的条件不可能精确满足（否则只有A=B=0， 使本题无解），可用积分条件代替： 3－7 见例题2。 3－8 同样，在的边界上，应考虑应用一般的应力边界条件(2-15)。
Mathematica 软件快速上手指南(11) ———数学课件制作简介 梁肇军 (华中师范大学数学系,湖北 武汉 430079) 中图分类号:G 434 文献标识码:B 文章编号:0488-7395(2001)09-0005-02 收稿日期:2001-02-20 作者简介:梁肇军(1938— ),男,广西柳州人,华中师范大学数学系教授.1 课件制作的基本程序 一篇优秀的数学课件,必须具备下述三个条件: 1.能体现现代、先进的教育思想,符合教育的科学 原理.2.能充分反映计算机在教学中的独特作用,利用其超强的计算能力,精确、快速的图形效果,能实时调控以及具有动画功能等.3.符合数学的基本原理,内容科学.因此,我们分下述三步来进行一篇数学课件的制作. 1.1 选取数学课件的脚本 数学课件的脚本,可以是一本教材,也可以是一个讲稿,或者是课本中的某一章某一节内容,按照讲义的要求处理,力求简练. 内容选取以后,我们需要把讲授的文稿设计成框图,框图里的信息力求简明扼要,框图里的内容要有先后顺序安排.同时,要设计一定的师生交流以及人机交流的内容.当然,巩固练习的安排也是必不可少的.由于计算机的内存有限,在这里我们还要考虑课件所需计算机存贮空间的分配. 1.2 按照教育教学原理,把脚本进行适当的编排、 增补 在排列学习项目时,应注意分析问题、逻辑推理的合理性和思维过程的流畅性.同时,适当地安排一些具有意外性的内容,对学生具有挑战性,能引起学生学习上的兴趣.为了避免学生在长期的紧张、连续的学习中产生疲劳,把学习过程分为若干个阶段是适宜而且应该的.在学习结束时,针对学生在练习中反映的问题,课件的最后应安排一定的小结内容. 在利用计算机进行辅助教学中,对计算机的实时调控,图形的逼真,动画的直观,以及教师与学生的对话、人机互动等等,要能深刻体现教育教学基本规律,要能充分调动学生学习的积极性. 1.3 数学式的软件语言处理 我们可以把课件内容大致分为两部分:一类是静态的,如中西文文字说明,图形,表格数据等,这部分内容制作者只需将其按Mathematica 系统的对其的输入要求录入到课件中即可;另一类是动态的,如代数式的实时运算,作图,动画,比较(包括代数式中某此数据修改后的结果对比,图形叠加,加色,数学模型中的适时调控等).下面将从中西文文字、作图、动画等三个方面来说明如何用Mathematica 软件语言处理上述数学式.至于Mathematica 系统在建立数学模型、开放实验等方面的有关问题,我们将另文说明. 中西文文字、图形、表格、数据(含代数式)在数学课件中主要起说明作用,除图形外其它的都可以直接利用键盘输入,只不过代数式的输入要按照Mathem 2 atica 系统规定的格式对它进行处理.对于图形,我们 可以通过函数作图或者是图元作图,即先输入一行正确的命令,然后运行它,输出结果即为需要的图形,存盘后即可保留在文稿上.对某些需要输出结果的代数式,也可如上办理.需要运算的代数式要放在独立的单元Cell 里,一个完整代数式建立一个独立的单元.由于计算机的计算速度快,一般几秒钟即可完成,为了教学的需要,我们需要把代数式分成几个部分,一步一步地让计算机执行.如化简一个代数式,我们可 5 2001年第9期 数学通讯
弹性力学简明教程（第四版）课后习题解答 第一章绪论 【1-1】试举例说明什么是均匀的各向异性体，什么是非均匀的各向同性体？ 【分析】均匀的各项异形体就是满足均匀性假定，但不满足各向同性假定；非均匀的各向异性体，就是不满足均匀性假定，但满足各向同性假定。 【解答】均匀的各项异形体如：竹材，木材。 非均匀的各向同性体如：混凝土。 【1-2】一般的混凝土构件和钢筋混凝土构件能否作为理想弹性体？一般的岩质地基和土质地基能否作为理想弹性体？ 【分析】能否作为理想弹性体，要判定能否满足四个假定：连续性，完全弹性，均匀性，各向同性假定。 【解答】一般的混凝土构件和土质地基可以作为理想弹性体；一般的钢筋混凝土构件和岩质地基不可以作为理想弹性体。 【1-3】五个基本假定在建立弹性力学基本方程时有什么作用？ 【解答】（1）连续性假定：假定物体是连续的，也就是假定整个物体的体积都被组成这个物体的介质所填满，不留下任何空隙。引用这一假定后，物体的应力、形变和位移等物理量就可以看成是连续的。因此，建立弹性力学的基本方程时就可以用坐标的连续函数来表示他们的变化规律。 完全弹性假定：假定物体是完全弹性的，即物体在对应形变的外力被去除后，能够完全恢复原型而无任何形变。这一假定，还包含形变与引起形变的应力成正比的涵义，亦即两者之间是成线性关系的，即引用这一假定后，应力与形变服从胡克定律，从而使物理方程成为线性的方程，其弹性常数不随应力或形变的大小而变。 均匀性假定：假定物体是均匀的，即整个物体是由同一材料组成的，引用这一假定后整个物体的所有各部分才具有相同的弹性，所研究物体的内部各质点的物理性质都是相同的，因而物体的弹性常数不随位置坐标而变化。 各向同性假定：假定物体是各向同性的，即物体的弹性在所有各个方向都相同，引用此假定后，物体的弹性常数不随方向而变。 小变形假定：假定位移和变形是微小的。亦即，假定物体受力以后整个物体所有各点的位移都远远小于物体原来的尺寸，而且应变和转角都远小于1。这样在建立物体变形以后的平衡方程时，就可以方便的用变形以前的尺寸来代替变形以后的尺寸。在考察物体的位移与形变的关系时，它们的二次幂或乘积相对于其本身都可以略去不计，使得弹性力学中的微分方程都简化为线性的微分方程。
Mathematic 9.0教程 1．基本介绍 1.1．Mathematica 的结构 Mathematica 是一个模块化的软件系统，其执行运算的内核与处理用户交互的前端是互相分离的。 表1 Mathematica 系统的基本部分 这样的设计比整体结构有许多优势。 例如，Mathematica 前端可运行在具有增强图形处理能力的本地计算机上，而Mathematica 内核可运行在更快地远程计算机上。或运行多个内核只需一个前端。 最常见的 Mathematica 工作方式是使用交互式文档称为笔记本。 笔记本把具有文字、图形、面板和其它资料的输入和输出放在一起。 用户使用笔记本既可进行运算，也可作为表达或发布自己的结果的工具。 其它常见的 Mathematica 界面包括基于文本的界面和 MathLink 接口。 表2 Mathematica 的常见界面种类 Mathematica 的一个重要特点是它不仅能与人交互，还能和其它程序交互. 这个功能是通过 MathLink 来实现的.。它是外部程序和 Mathematica 内核之间的标准双向通讯协议。 在众多可用的 MathLink 兼容的程序中，一些被用来作为 Mathematica 的前端。这些前端常常提供自己特有的用户界面，并把 Mathematica 内核纯粹作为嵌入的计算引擎。 Mathematica 内核 实际执行运算的部分 Mathematica 前端 处理与用户交互的部分 笔记本界面 交互式文档 基于文本的界面 MathLink 接口 由键盘输入的文本 与其它程序通讯
1.2．使用笔记本界面 如果使用图形界面，启动Mathematica的典型方法是双击Mathematica图标。如果使用文本操作系统，启动Mathematica的典型方法是键入mathematica。 表3 运行笔记本界面的Mathematica 使用图标或开始菜单图形界面中启动Mathematica 的典型操作mathematica 启动Mathematica 的操作系统命令 Mathematica 输入（在某些键盘上输入 文本输入以Shift+Enter 结束 Shift+Return） 选择退出菜单项退出Mathematica（在某些系统使用Quit）在"笔记本" 界面中，通过创建交互式文档与Mathematica 交互. 笔记本前端包含了许多菜单和图形化工具以便于创建和读取笔记本文档，发送和接收来自于Mathematica 内核的资料. 混合有文本，图形和Mathematica 的输入和输出的笔记本.
题提示和答案 《弹性力学简明教程》 习题提示和参考答案 第二章习题的提示与答案 2－1 是 2－2 是 2－3 按习题2－1分析。 2－4 按习题2－2分析。 2－5 在的条件中，将出现2、3阶微量。当略去3阶微量后，得出的切 应力互等定理完全相同。 2－6 同上题。在平面问题中，考虑到3阶微量的精度时，所得出的平衡微分方程都相同。其区别只是在3阶微量（即更高阶微量）上，可以略去不计。 2－7 应用的基本假定是：平衡微分方程和几何方程─连续性和小变形，物理方程─理想弹性体。 2－8 在大边界上，应分别列出两个精确的边界条件；在小边界（即次要边界）上，按照圣维南原理可列出3个积分的近似边界条件来代替。 2－9 在小边界OA边上，对于图2－15（a）、（b）问题的三个积分边界条件相同，因此，这两个问题为静力等效。 2－10 参见本章小结。 2－11 参见本章小结。 2－12 参见本章小结。 2－13 注意按应力求解时，在单连体中应力分量必须满足 （1）平衡微分方程， （2）相容方程， （3）应力边界条件（假设)。 2－14 见教科书。 2－15 见教科书。 2－16 见教科书。 2－17 取
它们均满足平衡微分方程，相容方程及x=0和的应力边界条件，因此，它们是该问题的正确解答。 2－18 见教科书。 2－19 提示：求出任一点的位移分量和，及转动量，再令,便可得出。 第三章习题的提示与答案 3－1 本题属于逆解法，已经给出了应力函数，可按逆解法步骤求解： （1）校核相容条件是否满足， （2）求应力， （3）推求出每一边上的面力从而得出这个应力函数所能解决的问题。 3－2 用逆解法求解。由于本题中 l>>h, x=0,l 属于次要边界（小边界），可将小边界上的面力化为主矢量和主矩表示。 3－3 见3-1例题。 3－4 本题也属于逆解法的问题。首先校核是否满足相容方程。再由求出 应力后，并求对应的面力。本题的应力解答如习题3-10所示。应力对应的面力是：主要边界： 所以在边界上无剪切面力作用。下边界无法向面力；上边 界有向下的法向面力q。 次要边界： x=0面上无剪切面力作用；但其主矢量和主矩在 x=0 面上均为零。 因此，本题可解决如习题3-10所示的问题。 3－5 按半逆解法步骤求解。 （1）可假设 （2）可推出 （3）代入相容方程可解出f、，得到