三、（15分）如图所示刚架，已知所有杆件截面尺寸相同，弹性模量E=20Pa ，截面面积A=0.2m 2，截面惯性矩I z =1/25m 4，荷载q =12kN·m ，F =10kN ，结点4的位移为δ4= [0.78 22.308 5.618]，考虑杆件轴向变形，试求局部坐标下单元①的内力。
解： （1）局部坐标系下单元①刚度矩阵：
1
0010
000.150.300.150.300.30.800.30.4=10010000.150.3
00.150.30
0.3
0.4
0.30.8k -??
??-????
-?
?-????
---??
-??
①
（2）局部坐标系下单元①结点位移列阵：
[]0
00
0.78422.308 5.618f δ=①
（3）局部坐标系下单元①杆端力列阵：
[]0
5505
5
f =---①
F （4）局部坐标下单元①的内力 []T
0.
784 6.6619.4450.784
3.339
2.802
f
=+=----①
①①
①
F
k F δ
四、（15分）已知图示结构，整体坐标下两杆件的单元刚度矩阵如下，总体坐标系和局部坐标系如图所示，试按照先处理法计算结构的总体刚度矩阵。
20020012
121201201212012120200200
1216012812016120820020012
0121201201212012120200200
12
8
12
1612
8
12
16k k ----????????-????????--==????--??????
??----????--????????
①
②
， 解：
2
0200016
01282
014012012014120
8
12
12
32K -????-????=-??--????-??
五、（10分）（1）对图a 所示结构，按照先处理法进行离散（含单元、结点、结点位移编码，忽略结构轴向变形，在图b 上标注）；（2）按照先处理法离散后，结构的刚度矩阵K 的大小为（ ）行（ ）列（请将结果直接写入括号）。 解：
六、（35分）设有一等厚矩形薄板，边长l =1m ，划分为两个个单元，其边上受如图所示均布载荷=50N q 作用，弹性模量E=100，泊松比μ=0，厚度t =1，已知单元①刚度矩阵，按照图示结构离散后：（1）计算单元①的几何矩阵B ①
；（2）按照先处理法求其整体刚度矩阵K 和等效结点荷载；（3）求解结点位移。
50
500000252525250502575252500252575255002525252500
50
50k
-????--????---=??---????--??-????
①
解： （1）0.5;A =
1; 1; 0
0; 1; 1
bi bj bm ci cj cm =-====-= ? 10
10000001010
1111
0-????=-????--??
①
B （2）
750075K ??
=??
??
[]T
250P =
（3）
[]T
T 100.33303??
?==????
m
(a )
(b)
装订线(答题不得超过此线)
. 一、题目 如图1所示，一个厚度均匀的三角形薄板，在顶点作用沿板厚方向均匀分布的竖向载荷。已知：P=150N/m ，E=200GPa ，=0.25，t=0.1m ，忽略自重。试计算薄板的位移及应力分布。 要求： 1. 编写有限元计算机程序，计算节点位移及单元应力。（划分三角形 单元，单元数不得少于30个）； 2. 采用有限元软件分析该问题（有限元软件网格与程序设计网格必 须一致），详细给出有限元软件每一步的操作过程，并将结果与程序计算结果进行对比（任选取三个点，对比位移值）； 3. 提交程序编写过程的详细报告及计算机程序； 4. 所有同学参加答辩，并演示有限元计算程序。 有限元法中三节点三角形分析结构的步骤如下： 1）整理原始数据，如材料性质、荷载条件、约束条件等，离散结构并进行单元编码、结点编码、结点位移编码、选取坐标系。 2）单元分析，建立单元刚度矩阵。 3）整体分析，建立总刚矩阵。 4）建立整体结构的等效节点荷载和总荷载矩阵 5）边界条件处理。 6）解方程，求出节点位移。 7）求出各单元的单元应力。 8）计算结果整理。 一、程序设计 网格划分 如图，将薄板如图划分为6行，并建立坐标系，则
刚度矩阵的集成 建立与总刚度矩阵等维数的空矩阵，已变单元刚度矩阵的集成。 由单元分析已知节点、单元的排布规律，继而通过循环计算求得每个单元对应的节点序号。 通过循环逐个计算：（1）每个单元对应2种单元刚度矩阵中的哪一种； （2）该单元对应总刚度矩阵的那几行哪几列 （3）将该单元的单元刚度矩阵加入总刚度矩阵的对应行列 循环又分为3层循环：（1）最外层：逐行计算 （2）中间层：该行逐个计算 （3）最里层：区分为第 奇/偶 数个计算 单元刚度的集成：[ ][][][][][]' '''''215656665656266256561661e Z e e e Z e Z e e e e k k k K k k k k k k +?++=? =?==?==?=?????? 边界约束的处理：划0置1法 X Y P X Y P
结构有限元分析程序设计 绪论 §0.1 开设“有限元程序设计”课程的意义和目的 §0.2 课程特点 §0.3 课程安排 §0.4 课程要求 §0.5 基本方法复习 $0.1 意义和目的 1.有限元数值分析技术本身要求工程设计研究人员掌握 1). 有限元数值分析技术的完善标志着现代计算力学的真正成熟和实用化,已在各种 力学中得到了广泛的应用。比如：，已杨为工程结构分析中最得以收敛的技术手段，现代功用大致有： a). 现代结构论证。对结构设计从内力，位移等方面进行优劣评定，从而进 行结构优化设计。 b)可取代部份实验，局部实验+有限元分析，是现代工程设计研究方法的一大 特点。 c)结构的各种功能分析（疲劳断裂，可靠性分析等）都以有限元分析工具作为 核心的计算工具。 2). 有限元数值分析本身包括着理论+技术实现（本身功用所绝定的） 有限元数值分析本身包括着泛函理论+分片插值函数+程序设计 2. 有限元分析的技术实现（近十佘年的事）更依赖于计算机程序设计 有限元分析的技术取得的巨大的成就，从某种意义上说，得益于计算机硬件技术的发展和程序设计技术的发展，这两者的依赖性在当代表现得更加突出。（如可视化技术） 3.从学习的角度，不仅要学习理论，而且要从程序设计设计角度对这些理论的技术实现有 一个深入的了解，应当致力于掌握这些技术实现能力，从而开发它，发展它。（理论本身还有待于进一步完美相应的程序设计必须去开发） 4.程序设计不仅是实现有限元数值分析的工具和桥梁，而且在以下诸方面也有意义： 1). 精通基本概念，深化理论认识； 2). 锻炼实际工程分析，实际动手的能力； 3). 获得以后工作中必备的工具。（作业+老师给元素库） 目的：通过讲述有限元程序设计的技术与技巧，便能达到自编自读的能力。 §0.2 课程特点 总描述：理论+算法+数据结构（程序设计的意义） 理论：有限元算法，构造，步骤，解的等外性，收敛性，稳定性，误差分析 算法；指求解过程的技术方法，含两方面的含义；a. 有限元数值分析算法，b, 与数据结构有关的算法（总刚稀疏存贮，提取，节点优化编号等） 数据结构：指各向量矩阵存贮管理与实现，辅助管理结构（指针，数据记录等） 具体特点： 理论性强：能量泛函理论+有限元构造算法+数据结构构造算法 内容繁杂：理论方法+技术方法+技术技巧 技巧性强：排序，管理结构（指针生成，整型运算等）
湖南大学 《有限单元法》编程大作业 专业：土木工程 姓名： 学号： 2013年12月
目录 程序作业题目： (3) 1、程序编制总说明 (3) 2、Matlab程序编制流程图 (3) 3、程序主要标示符及变量说明 (4) 4、理论基础和求解过程 (5) 4.1、构造插值函数 (5) 4.2位移插值函数及应变应力求解 (5) 5.程序的验证 (6) 附录：程序代码 (15)
程序作业题目： 完成一个包含以下所列部分的完整的有限元程序( Project) 须提供如下内容的文字材料（1500字以上）： ①程序编制说明； ②方法的基本理论和基本公式； ③程序功能说明； ④程序所用主要标识符说明及主要流程框图； ⑤ 1～3 个考题：考题来源、输出结果、与他人成果的对比结果（误差百分比）； ⑥对程序的评价和结论（包括正确性、适用范围、优缺点及其他心得等）。 须提供源程序、可执行程序和算例的电子文档或文字材料。选题可根据各自的论文选题等决定。 1、程序编制总说明 a.该程序采用平面三角形等参单元，能解决弹性力学的平面应力、平面应变问题。 b.能计算单元受集中力的作用。 c.能计算结点的位移和单元应力。 d.考题计算结果与理论计算结果比较，并给出误差分析。 e.程序采用MATLAB R2008a编制而成。 2、Matlab程序编制流程图
图1 整个程序流程图 3、程序主要标示符及变量说明 1、变量说明： Node ------- 节点定义 gElement ---- 单元定义 gMaterial --- 材料定义，包括弹性模量，泊松比和厚度 gBC1 -------- 约束条件 gNF --------- 集中力 gk------------总刚 gDelta-------结点位移 输入结构控制参数 输入其它数据 形成整体刚度阵 引入支承条件 解方程，输出位移 求应力，输出应力 形成节点荷载向量 开始 结束 1 单元面积 求弹性矩阵 单元刚度矩阵 位移-应变矩阵 6 7 8 9 10 2 3 4 5
“Java程序设计基础”课程习题 一、填空 1.Java程序分两类___Applet___和application，Java Application 类型的程序，程序从 ___main方法___开始执行。 2.定义一个Java类时，通过关键字__extends____指明该类的父类。一个类可以有___1___ 个父类。 3.用public修饰的类称为_公有类或公用类__。用public修饰的类成员称为公有成员。被 说明为public的内容可以被__所有其他类___ 使用。如果public类文件与使用它的类文件不在同一目录中，需要通过__import____语句引入。 4.用___private___ 修饰的类成员称为私有成员。私有成员只能在__本类__ 中使用。 5.如果子类定义的成员变量与父类的成员变量同名，称为___方法覆盖___ ，要表明使用 子类的成员变量，可以在成员变量前加上关键字__super___ 。 6.____Object__ 类是Java类库中所有类的父类。 7.Java字符使用__16位的字符集,该字符集成为__Unicode____ 。 8.当子类中定义的方法与父类方法同名时，称子类方法___覆盖___ 父类方法，子类默认 使用自己的方法。使用父类的同名方法，必须用关键字__super__ 说明。 9.Java源程序文件名的后缀是___.java___，Java字节码文件名的后缀是_.class_____。 10.Java类名的第一个字母通常要求___大写___。 11.Java程序由____类__组成，每个程序有一个主类，Java程序文件名应与____主__类的 名称相同。 12.Java__Application_类型的程序需要main()方法，程序从__main____开始执行。 13.布尔型数据类型的关键字是_boolean__ ，占用位数是___1位___ ，有__true__ 和_false_ 两种值。整型数可以采用_十_ 、__八_ 和__十六_三种进制表示。 14.八进制整数以数字__0_开头。十六进制整数以_0x或0X_ 开头。 15.int整型数占用__32位内存。long整型数占用__64 位内存。 16.127L表示__长整型____ 常量。 17.根据占用内存长度将浮点常量分为_double_____ 和__float____ 两种。 18.单精度浮点常量占用__32_ 位内存，双精度浮点常量占用__64 位内存。 19.在Java语言中，字符串“ABC\tD\b\n”中包括__7个字符。 20.数学关系4
Java程序设计题库 第一部分绪论 1、下列关于Java语言的特点，描述错误的是（C） A．Java是跨平台的编程语言B．Java支持分布式计算 C．Java是面向过程的编程语言D．Java支持多线程 2、Java语言具有许多优点和特点，下列选项中，哪个反映了Java程序并行机制的特点？（B） A、安全性 B、多线性 C、跨平台 D、可移植 3、Java JDK中调试器的命令是（C）。 A、javac B、java C、jdb D、avah 4、运行jar文件中class文件需要在java命令后面加的参数为（A）。 A、-cp B-g C-d D-verbose 5、下面哪项在java中是不合法的标识符？（C） A、$user B、point C、You&me D、_endline 6、下面关于Java语言说法正确的是（ABCD）。 A、Java语言是面向对象的、解释执行的网络编程语言。 B、Java语言具有可移植性，是与平台无关的编程语言。 C、Java语言可对内存垃圾自动收集。 D、Java语言编写的程序虽然是“一次编译，到处运行”，但必须要有Java的运行环境。 7、在Java程序设计中，程序员创建（）文件，然后编译器把它们转化为（）文件。( B) A、源, HTML B、源, 字节代码 C、字节代码, 源 D、HTML, 字节代码 8、Java的JVM是指（B）。 A、由Java操作的家用设备（通常是一个漂亮的小烤箱） B、Java字节代码的解释程序 C、Java源代码的编译器 D、运行Java 所需的硬件设备 9、在Java中，关于CLASSPA TH环境变量的说法不正确的是（A）。 A、CLASSPATH一旦设置之后不可修改，但可以将目录添加到该环境变量中。 B、编译器用它来搜索各自的类文件。 C、CLASSPATH是一个目录列表。 D、解释器用它来搜索各自的类文件。 10、一个可以独立运行的Java应用程序（D）。
有限元分析软件 有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司，其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析： 1．在世界范围内的知名度：两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS，并且在五年前ANSYS 的界面是当时最好的界面之一，所以在中国，ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域：ANSYS 软件注重应用领域的拓展，目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究，致力于解决该领域的深层次实际问题。 3．性价比：ANSYS 软件由于价格政策灵活，具有多种销售方案，在解决常规的
J a v a程序设计总复习题 1、编写一个Java程序，用if-else语句判断某年份是否为闰年。(分支) // Programme Name LeapYear.java public class LeapYear{ public static void main(String args[]){ int year=2010; if(args.length!=0) year=Integer.parseInt(args[0]); if((year%4==0 && year%100!=0)||(year%400==0)) 年是闰年。"); else 年不是闰年。"); } }//if-else语句 2、编写一个Java程序在屏幕上输出1！+2！+3！+……+10！的和。（循 环） // programme name ForTest.java public class ForTest { public static void main( String args[] ) { int i,j,mul,sum=0; for(i=1;i<=10;i++) { mul=1; for(j=1,j<=i;j++) { mul=mul*j; } sum=sum+mul; } “1！+2！+3！+……+10！= ”+sum); } } 3、依次输入10个学生成绩，判断学生（优秀、良好、中等、及格、不及格） 并计算人数（switch） 4、使用冒泡排序（数组） public class BubbleSort { public static void main(String[] args) {
重庆大学本科学生课程设计任务书 课程设计题目有限元程序设计 学院资源及环境科学学院专业工程力学年级2010级 已知参数和设计要求： 1.独立完成有限元程序设计。 2.独立选择计算算例，并能通过算例判断程序的正确性。 3.独立完成程序设计报告，报告内容包括理论公式、程序框图、程序本 体、计算算例，算例结果分析、结论等。 学生应完成的工作： 1.复习掌握有限单元法的基本原理。 2.掌握弹性力学平面问题3节点三角形单元或4节点等参单元有限元方法 的计算流程，以及单元刚度矩阵、等效节点载荷、节点应变、节点应力 和高斯积分等的计算公式。 3.用Fortran语言编写弹性力学平面问题3节点三角形单元或4节点等参 单元的有限元程序。 4.在Visual Fortran 程序集成开发环境中完成有限元程序的编辑和调试 工作。 5.利用编写的有限元程序，计算算例，分析计算结果。 6.撰写课程设计报告。 目前资料收集情况（含指定参考资料）： 1.王勖成，有限单元法，北京：高等教育出版社，2002。 2.O.C. Zienkiewicz, R. L. Taylor, Finite Element Method, 5th Eition, McGraw-Hall Book Company Limited, 2000。 3.张汝清，董明，结构计算程序设计，重庆：重庆大学出版社，1988。 课程设计的工作计划： 1.第1周星期一上午：教师讲解程序设计方法，程序设计要求和任务安 排。 2.第1周星期一至星期二完成程序框图设计。 3.第1周星期三至第2周星期四完成程序设计。 4.第2周星期五完成课程设计报告。 任务下达日期 2013 年 6 月 6 日完成日期 2013 年 07 月 03 日 指导教师（签名） 学生（签名）
有限元分析软件的比较（购买必看）－转贴 随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element A nalysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面： 增加设计功能，减少设计成本； 缩短设计和分析的循环周期； 增加产品和工程的可靠性； 采用优化设计，降低材料的消耗或成本； 在产品制造或工程施工前预先发现潜在的问题； 模拟各种试验方案，减少试验时间和经费； 进行机械事故分析，查找事故原因。 在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件
Ｊavａ程序设计基础课后习题参考答案 第2章 1、关于Jａｖa Applicａｔion得入口方法maｉn（）得检验: ｍａin（）方法得参数名就是否可以改变? maｉｎ（）方法得参数个数就是否可以改变？ 该方法名就是否可以改变？ 参考答案:（1)maｉn（）方法得参数名可以改变.(2）maiｎ()方法得参数个数不可以改变。（3）该方法名不可以改变。 2、当一个程序没有main（)方法时，能编译吗?如果能编译，能运行吗? 参考答案:当一个程序没有ｍaｉn(）方法就是，就是可以编译通过得,但就是不能给运行,因为找不到一个主函数入口。 3、下列语句能否编译通过？ bytｅi =127； ｂｙｔｅj = １28； longｌ1 = ９99９９９； ｌoｎg l２= ９9９9999999； 参考答案：ｂyｔe i 与long l１可以编译通过。而bｙte j 与ｌongｌ2 超出自身数据类型范围,所以编译失败。 4、下列语句能否编译通过？ float f1 =3、５； floａt f2 = 3.5f; 参考答案：java中浮点型得数据在不声明得情况下都就是doublｅ型得，如果要表示一个数据就是ｆloaｔ型得,必须在数据后面加上“Ｆ”或“ｆ”；因此,floatｆ1 无法编译通过。 5、验证ｉnt 与char，ｉnt与doｕbｌｅ等类型就是否可以相互转换。 参考答案:(１）cｈar类型可以转换为iｎt 类型得，但就是inｔ类型无法转换为chａｒ类型得;（２）iｎt 可以转换为ｄoubｌe类型得，但就是ｄｏuble类型无法转换为int 类型得。 6、计算下列表达式,注意观察运算符优先级规则。若有表达式就是非法表达式,则指出不合法之处且进行解释。 （1)4+5 ＝= 6＊2 ?（2) （4=5）/6?? (3）９%2*7/3＞17（4）(4＋5）<=6/3 ? (5) 4＋5％３！=7－2????（６）４+5/6〉=10％2 参考答案：表达式（２）为不合法表达式，只能将值赋值给一个变量,因此其中（4=5）将5赋值给4就是不合法得. 7、下列（）就是合法得Ｊava标识符。 (１）Counter1 ??（2）$indeｘ， (3) nａmｅ-7 ??(4）_bytｅ
参考答案 第1章 3、计算机系统由哪两部分组成的？计算机硬件结构由哪几部分组成？它们各自有什么作用？ 一台计算机由硬件和软件组成。一台典型的计算机由五大部分组成。这五大部分是：运算器，控制器，存储器，输入设备和输出设备。 运算器是执行算术运算和逻辑运算的部件。 控制器是向计算机其他部分发送命令的部件。 存储器是计算机用来存储数据、信息的部件。 输入设备就是外界向计算机输入信息设备。 输出设备恰好与输入设备的作用相反，它将处理过后的信息输出呈现给用户。 9、将以下十进制数转换为对应的二进制数 （1）32 （2）97 （3）256 （4）500 （1）100000 （2）1100001 （3）100000000 （4）111110100 第2章 1．Java语言有哪些主要特点。 平台独立性 安全性 多线程 网络化 面向对象 3．Java Application的开发步骤有哪些。 Java Application的开发步骤： （1）下载JDK软件并安装； （2）配置相应的环境变量（path和classpath）； （3）编写Java源程序（文本编辑器或集成开发环境IDE）； （4）编译Java源程序，得到字节码文件（javac *.java）； （5）执行字节码文件（java 字节码文件名）。
4．什么是环境变量，设置环境变量的主要目的是什么。 环境变量的配置主要是为了进行“寻径”，也即让程序能找到它需要的文件，所以设置的内容就是一些路径。 第3章 1．Java语言对于合法标识符的规定是什么？指出以下哪些为合法标识符。 a a2 3a *a _a $a int a% 在Java语言中，标识符必须以字母、美元符号或者下划线打头，后接字母、数字、下划线或美元符号串。另外，Java语言对标识符的有效字符个数不做限定。 合法的标识符： a a2 _a $a 5．数据类型强制转换的原则是什么？如何转换？ 对于变窄转换，如long到short、double到float，或者不兼容转换：float到short、char 到short等，则需要进行强制转换。 float f = 11.5; short b ; b = （short）f; （强制转换） 第4章 5．用穷举法求出3位数中百、十、个位数的立方和就是该数的数。 public class Test { public static void main(String[] args) { int a,b,c,x=100; while(x<1000){ a=x%10; b=(x%100-a)/10; c=(x-x%100)/100; if(a*a*a+b*b*b+c*c*c==x) System.out.println(x); x+=1;
《JA V A程序设计》练习题 一、单选择题 1、编译Java Application 源程序文件将产生相应的字节码文件，这些字节码文件的扩展名为( B )。 A. java B. .class C. html D. .exe 2、设x = 1 , y = 2 , z = 3，则表达式y＋＝z－－/＋＋x 的值是( A)。 A. 3 B. 3. 5 C. 4 D. 5 3、不允许作为类及类成员的访问控制符的是( C )。 A. public B. private C. static D. protected 4、为AB类的一个无形式参数无返回值的方法method书写方法头，使得使用类名AB作为前缀就可以调用它，该方法头的形式为( A)。 A. static void method( ) B. public void method( ) C. final void method( ) D. abstract void method( ) 5、关于选择结构下列哪个说法正确？( B ) A．if语句和else语句必须成对出现 B．if语句可以没有else语句对应 C．switch结构中每个case语句中必须用break语句 D．switch结构中必须有default语句 6、while循环和do…while循环的区别是：( D) A．没有区别，这两个结构任何情况下效果一样 B．while循环比do…while循环执行效率高 C．while循环是先循环后判断，所以循环体至少被执行一次 D．do…while循环是先循环后判断，所以循环体至少被执行一次 7、关于for循环和while循环的说法哪个正确？( B) A．while循环先判断后执行，for循环先执行后判断。 B．while循环判断条件一般是程序结果，for循环的判断条件一般是非程序结果C．两种循环任何时候都不可以替换 D．两种循环结构中都必须有循环体，循环体不能为空
【程序1】 题目：古典问题：有一对兔子，从出生后第3个月起每个月都生一对兔子，小兔子长到第三个月后每个月又生一对兔子，假如兔子都不死，问每个月的兔子总数为多少 //这是一个菲波拉契数列问题 public class lianxi01 { public static void main(String[] args) { System.out.println("第1个月的兔子对数: 1"); System.out.println("第2个月的兔子对数: 1"); int f1 = 1, f2 = 1, f, M=24; for(int i=3; i<=M; i++) { f = f2; f2 = f1 + f2; f1 = f; System.out.println("第" + i +"个月的兔子对数: "+f2); } } } 【程序2】 题目：判断101-200之间有多少个素数，并输出所有素数。 程序分析：判断素数的方法：用一个数分别去除2到sqrt(这个数)，如果能被整除，则表明此数不是素数，反之是素数。 public class lianxi02 { public static void main(String[] args) { int count = 0; for(int i=101; i<200; i+=2) { boolean b = false; for(int j=2; j<=Math.sqrt(i); j++) { if(i % j == 0) { b = false; break; } else { b = true; } } if(b == true) {count ++;System.out.println(i );} } System.out.println( "素数个数是: " + count); } } 【程序3】 题目：打印出所有的"水仙花数"，所谓"水仙花数"是指一个三位数，其各位数字立方和等于该数本身。例如：153是一个"水仙花数"，因为153=1的三次方＋5的三次方＋3的三次方。 public class lianxi03 { public static void main(String[] args) { int b1, b2, b3; for(int m=101; m<1000; m++) {
有限元编程的c++实现算例 1. #include<> 2. #include<> 3. 4. 5. #define ne 3 #define nj 4 #define nz 6 #define npj 0 #define npf 1 #define nj3 1 2 #define dd 6 #define e0 #define a0 #define i0 #define pi 16. 17. 18. int jm[ne+1][3]={{0,0,0},{0,1,2},{0,2,3},{0,4,3}}; /*gghjghg*/ 19. double gc[ne+1]={,,,}; 20. double gj[ne+1]={,,,}; 21. double mj[ne+1]={,a0,a0,a0}; 22. double gx[ne+1]={,i0,i0,i0}; 23. int zc[nz+1]={0,1,2,3,10,11,12}; 24. double pj[npj+1][3]={{,,}}; 25. double pf[npf+1][5]={{0,0,0,0,0},{0,-20,,,}}; 26. double kz[nj3+1][dd+1],p[nj3+1]; 27. double pe[7],f[7],f0[7],t[7][7]; 28. double ke[7][7],kd[7][7]; 29. 30. 31. 36. void jdugd(int); 38. void zb(int); 39. void gdnl(int); 40. void dugd(int);
若有定义int a＝2,b＝2; 则表达式(a++)+(++b) 的值是( )。 A. 5 监听器接口的方法返回值的类型是（）。C. void 在浏览器中执行applet 程序，以下选项中的哪个方法将被最先执行（）。 D. init() 以下哪个不是Java的原始数据类型（） B. Boolean 向Applet传递参数的正确描述是（）。C. C. PARAM name＝"age,value＝20" 阅读下列代码; public class Test{ String s＝"One World One Dream"; public static void main(String args[]){ System.out.println(s); } } 其运行结果是＿ C. 编译时出错//除非创建Text类的对象 如果容器组件p的布局是BorderLayout，则在p的下边中添加一个按钮b，应该使用的语句是（）。 A. p.add(b,"South"); 编译和运行下面的应用程序，并在命令行界面输入12345 ，则回车后屏幕输出的结果是( )。public class A { public static void main(String args[]) throws IOException{ BufferedReader buf＝new BufferedReader( new I nputStreamReader(System.in)); String str＝buf.readLine(); int x＝Integ er.parseInt(str); System.out.println(x/100); } } B. 123 构造函数何时被调用？（） B. 创建对象时 以下哪项可能包含菜单条（）。 D. Frame
有限元编程的c++实现算例 1.#include
有限元及程序设计 1. 下列关于高精度单元描述正确的是（ ）。 A.等参元的位移模式和坐标变换采用不同的形函数 B.矩形单元形状规则，因而使用范围较 广 C.6结点三角形单元、10结点三角形单元、8结点矩形单元和 12结点矩形单元的单元 刚度矩阵的建立过程是不一样的 D.6结点三角形单元较容易模拟物体的边界形状 【参考答案】：D 2. 0 =cxy 能解决矩形板（ ）问题。 A.左右均布拉压 B.上下均布拉压 C 纯剪切 D.纯弯曲 【参考答案】：C 3. 下列关于等参元的叙述不正确的是（ A.精度较高 B.能较好的模拟边界条件 【参考答案】：D 4. 薄板的边界不包括（ A.简支边界 B.固定边界 C.自由边界和荷载边界 D.非固定边界 【参考答案】：D ）。 C.输入的信息量较少 D.输入的信息量较多 5. 下列属于平面应力问题的是（ A.平板坝的平板支墩 B 挡土墙 【参考答案】：A 6. 在应力函数上任意增减一个（ A.线性项 B.二次项 C.三次项 【参考答案】：A 7. 下列不属于提高单元精度的方法是 A.增加单元结点数目 B 在单元内增设结点 【参考答案】：C 8. 空间问题的基本平衡微分方程有（ A.2 B.3 C.4 D.5 【参考答案】：C ）。 C.重力水坝 D.受内水压力作用的圆管 ），对应力分量无影响。 D.常数项 ）。 C.减少单元结点数目 D.设等参元 个。 ）问题。 C 纯剪切 D.纯弯曲 A.左右均布拉压 B.上下均布拉压 【参考答案】：B 10. 下列属于不规则单元的有（ A.正四面体单元 B.正三棱体单元 【参考答案】：C 11. 空间问题的基本未知位移分量有（ A.2 B.3 C.4 D.5 【参考答案】：B 1. 薄板小挠度弯曲理论的基本假定是（ ）。 A.直法线假定 B.法向位移假定 C.中面位移假定 D.板内无挤压假定 【参考答案】：A|C|D 2. 弹性力学平面问题 按应力求解具体可分为（ ）两种。 A.逆解法 B.顺解法 C.半逆解法 D.半顺解法 ）。 C.任意四面体单元 D.正六面体单元 ）个。
《java程序设计》复习题 第一部分选择题（单选） 1. javac 是 java的（B） A.解释器 B.编译器 C.小程序浏览器 D.调试器 2. appletviewer是java的（C） A.解释器 B.编译器 C.小程序浏览器 D.调试器 3. 编译一个定义了三个类和三个方法的 java源码文件，总共会产生几个字节码文件，其扩展名是什么？（C） A.5个字节码文件，其扩展名是“.java” B.2个字节码文件，其扩展名是“.java” C.3个字节码文件，其扩展名是“.class” D.2个字节码文件，其扩展名是“.class” 4. java源码文件的名称与主类的名称（B） A.必须相同但不区分大小写 B.必须相同且区分大小写 C.必须不同 D.任意 5. 下面那些选项不能用做变量名的首字符（C） A.字母 B.下划线（_） C.数字 D.美元符号（$） 6. 下面哪个不是java的原始数据类型？（B） A.short B.boolean C.Int D.Float 7.下面哪条语句不会出现编译警告或错误？（A） A.float f=1.3F; B.chat c =”a”; C.float f=1.3;
D.boolean d=null; 8． 以下声明常量不正确的是（B） A.final int PI=3; B.final PI=3; C.static final int PI=3: D.public static final int PI=3; 9. java定义了4个整型数据类型：字节性（byte）短整型(short) 整型（int）长整型（long），其中byte占用内存为（A） A.8bit B.16bit C.32 bit D.64bit 10. 5%2.1的值是（D） A.1 B.0 C.0.8 D.不能通过编译 11. -11&3的值是（A） A.1 B.3 C.9 D.–11 12. 8<<2的值是（A） A.2 B.4 C.8 D.32 13 若有int a = 1, b= 2 , c = 3, d=4; 则表达式a>b?a:c>d?c:d的值是（D） A.1 B.2 C.3 D.4 14. 以下为合法类型转换的为（D） A.( int )”A” B.(char)”A” C.(boolean)(2-2) D.(char) (651/12)
P 9 m 9 m 一、题目 如图1所示，一个厚度均匀的三角形薄板，在顶点作用沿板厚方向均匀分布的竖向载荷。已知：P=150N/m，E=200GPa，=，t=，忽略自重。试计算薄板的位移及应力分布。 要求： 1.编写有限元计算机程序，计算节点位移及单元应力。（划分三角形 单元，单元数不得少于30个）； 2.采用有限元软件分析该问题（有限元软件网格与程序设计网格必 须一致），详细给出有限元软件每一步的操作过程，并将结果与程序计算结果进行对比（任选取三个点，对比位移值）； 3.提交程序编写过程的详细报告及计算机程序； 4.所有同学参加答辩，并演示有限元计算程序。 有限元法中三节点三角形分析结构的步骤如下： 1）整理原始数据，如材料性质、荷载条件、约束条件等，离散结构并进行单元编码、结点编码、结点位移编码、选取坐标系。 2）单元分析，建立单元刚度矩阵。 3）整体分析，建立总刚矩阵。 4）建立整体结构的等效节点荷载和总荷载矩阵 5）边界条件处理。 6）解方程，求出节点位移。 7）求出各单元的单元应力。 8）计算结果整理。 一、程序设计
网格划分 如图，将薄板如图划分为6行，并建立坐标系，则 X Y P X Y P 节点编号 单元编号
刚度矩阵的集成 建立与总刚度矩阵等维数的空矩阵，已变单元刚度矩阵的集成。 由单元分析已知节点、单元的排布规律，继而通过循环计算求得每个单元对应的节点序号。 通过循环逐个计算：（1）每个单元对应2种单元刚度矩阵中的哪一种； （2）该单元对应总刚度矩阵的那几行哪几列 （3）将该单元的单元刚度矩阵加入总刚度矩阵的对应行列 循环又分为3层循环：（1）最外层：逐行计算 （2）中间层：该行逐个计算 （3）最里层：区分为第奇/偶数个计算 单元刚度的集成： [][] [][] [][] ' ' ' ' ' ' 2 1 56 56 6 6 56 56 2 6 6 2 56 56 1 6 6 1 e Z e e e Z e Z e e e e k k k K k k k k k k + ? + + = ? = ? = = ? = = ? = ? ? ? ? ? ? 边界约束的处理：划0置1法 适用：这种方法适用于边界节点位移分量为已知(含为0)的各种约束。 做法： （1）将总刚矩阵〔K〕中相应于已知位移行主对角线元素置1，其他元素改为零；同 时将载荷列阵｛R｝中相应元素用已知位移置换。 ◎这样，由该方程求得的此位移值一定等于已知量。 （2）将〔K〕中已知位移相应的列的非主对角成元素也置0，以保持〔K〕的对称性。 ◎当然，在已知位移分量不为零的情况下，这样做就改变了方程左端的数值，为 保证方程成立，须在方程右端减去已知位移对该方程的贡献——已知位移和相应总刚元素的乘积。◎若约束为零位移约束时，此步则可省去。 特点： （1）经以上处理同样可以消除刚性位移(约束足够的前提下)，去掉未知约束反力。 （2）但这种方法不改变方程阶数，利于存贮。