hand_up的博客

原创 瞬态动力学分析

显式算法则采用前一时刻的物理量与当前物理量对比，只要时间步足够小，一般不存在收敛问题。一般情况下，显式算法多用于作用时间非常短暂的分析，如爆炸、跌落冲击等， LS - DYNA 采用的即为显式算法，而隐式算法多用于作用时间相对较长的分析，如机械传动、地震作用等， Mechanical 的 TransientStructural 模块采用的即为隐式算法。可以求解上述方程，但模态叠加法的前提条件是结构需要满足线性条件，当结构存在接触行为、材料为非线性或结构存在大变形可能时，只能通过直接积分法求解动力学方程。

原创 ANSYS Workbench注意事项

应力集中和应力奇异在结构分析过程中其应力值的表现形式并不相同，两个概念不能混为一谈。应力奇异是由于网格不断细分其应力值呈现无穷大的趋势，应力值是发散的，而应力集中是由于结构的构造特性所引起的局部应力远大于其他部位的应力，例如本例中的 L 形梁的圆角过渡处及工程中常见的平面薄板的中间圆孔处等，随着网格的细化其应力值最终是收敛的。一般情况下，应力奇异的部位一定存在应力集中，而应力集中部位不一定有应力奇异。应力奇异现象一般发生在模型以直线作为突变的位置。

原创 结构屈曲分析

非线性屈曲分析的目的是获取结构可以承受载荷的最大值及结构失稳的变形状态，一般用户分析前需要获取一个载荷值用于非线性屈曲的加载，故在进行非线性屈曲分析前，一般以同样的结构和边界分析其线性屈曲的载荷值，将其载荷值放大1.5倍用于非线性屈曲分析的载荷值。线性屈曲分析需要具备较多的前提条件，如载荷无偏心、材料无缺陷等，在实际工程应用中结构制作过程和加载方式很难达到线性屈曲要求的状态，故线性屈曲得到的载荷值不能直接作为判定结构失稳的依据。(1）屈曲分析针对的结构处于受压状态，当结构处于受拉状态时，结构更易强度失效。

原创 强度理论介绍和惯性矩推导

强度理论选定后，当结构的某点应力大于许用应力时，则应根据结构形式和边界条件综合考虑改善结构以防止结构失效，一般而言：对于仅承受拉压载荷或者剪切载荷的构件，当强度不满足时可以增加构件的受力面积以增强结构的强度；当结构形式已经固定且无法改变结构形状或面积时，则可以考虑更换材料，例如以高强度钢代替的方式提高结构抵抗失效的能力。材料力学给出4个强度理论，分别为第一强度理论（最大拉应力理论）、第二强度理论（最大拉应变理论）、第三强度理论（最大切应力理论）、第四强度理论（形状改变比能理论）。选用第一、第二强度理论；

原创 PTA 币值转换

输入一个整数（位数不超过9位）代表一个人民币值（单位为元），请转换成财务要求的大写中文格式。如23108元，转换后变成“贰万叁仟壹百零捌”元。为了简化输出，用小写英文字母a-j顺序代表大写数字0-9，用S、B、Q、W、Y分别代表拾、百、仟、万、亿。于是23108元应被转换输出为“cWdQbBai”元。输入格式：输入在一行中给出一个不超过9位的非负整数。输出格式：在一行中输出转换后的结果。注意“零”的用法必须符合中文习惯。输入样例1：813227345输出样例1：iYbQ

原创 PDF转换出矢量图EMF

matlab打开图像--编辑-复制--粘贴。solidworks截图打印为PDF；CAD打印输出PDF。

原创 Visio使用

点击确定—确定（CAD转Visio会产生失真，字体等都需要重新调整填充需在Visio中重新绘制，然后再填充，记得按照设计要求修改字体），如下图所示。选中要成为上标的文字,ctrl+shift+"="最后一步复制—粘贴即可，WORD中插入后如图所示。AutoCAD 版本 2007 到 2017。选中要成为下标的文字,ctrl+"="导入 .dwg 或 .dxf 文件。右键图形—CAD绘图对象—转换。将默认选项改为如下图所示。

原创 Recurdyn导入，导出

装配层次结构：选上表示导入的装配体包括子装配体，再勾选body或者subsystem表示子装配体以。显示在recurdyn中；不选，表示导入的装配体全部以零件显示。

原创 matlab常用特殊字符表

英文中文。

原创 matlab绘制指定球心坐标及半径球面

【代码】matlab绘制指定球心坐标及半径球面。

原创 matlab函数图像最大值

【代码】matlab函数图像最大值。

原创 recurdyn实用操作及注意事项

方法二：在安装路径下找到文件：RecurDyn V * R * \ Toolkits \ MF tire \ Road data Files，然后将其拷贝到存储模型的同一个文件夹中，接下来单击：Ground>Road Data GRImport，选择刚刚保存的rdf文件即可。面过多不易选取，可以右击，在弹出的快捷命令中通过利用"Select Box”或者"Select List"命令来进行对象的选择。方法一：选择在建立路面的时候，就将路面曲线Curve绘制成自己需要的形状，然后再变成Road。

原创 recurdyn接触特征参数含义

所以，在滑动摩擦计算模型(Sliding Type)中，RecurDyn 利用图8-3 所示速度与摩擦系数之间的关系曲线来解决这个问题，用户可以指定静摩擦系数与动摩擦系数发生作用时的相应速度值，图 8-3 中us,ud分别表示静摩擦系数、动摩擦系数，vs、vd分别表示静滑动速度和动滑动速度。压痕指数):用于产生压痕衰减效果，当接触穿透值非常小时，计算的接触力有可能为负值，为避免这种情况，可以通过制定一个大于 的压痕指数来避免这种情况的发生，其所对应的指数即为法向接触力f的公式中的 m3。

原创 SW质量属性

Ix = (1.00, 0.00, 0.00)，Iy = (0.00, 0.00, -1.00)， Iz = (0.00, 1.00, 0.00)是3个特征向量，指出了惯性主轴的方向。(Moment of Inertia)，是刚体绕轴转动时惯性（回转物体保持其匀速圆周运动或静止的特性）的量度，用字母I或J表示。中的角色相当于线性动力学中的质量，可形式地理解为一个物体对于旋转运动的惯性，用于建立角动量、角速度、力矩和角加速度等数个量之间的关系。：指派质量、质量中心和惯性张量的值以覆写所计算的值。

原创 RecurDyn的分析速度过慢或者结果异常时

原创 recurdyn履带问题

用于搜索与链轮接触的链节的工具。有两种方式,Full Search，搜索所有的链节并与之进行接触，和Partial Search，搜索链轮附近的履带链节并与之进行接触，选用Partial Search可以减少仿真时间，但是可能会出现链轮与部分履带链节之间不发生接触导致穿透的情况，而Partial Search又可以让用户自定义搜索范围，就是User boundary」，因此如果没接触上建议先变为Full search。履带链节之间是通过力来连接的，有力的限制，一般情况下尺寸参数合适不会考虑干涉。

原创 recurdyn同步带装配问题

计算一段带长，平均分配到上述中。

原创 SW simulation结果打包功能

将本次开启文件的所有关联文件(零件、装配体、工程图、参考、 设计表、设计活页夹内容、贴图、外观及布景以及SOLIDWORKS Simulation结果)收集到一个文件夹或者zip(压缩)文件中。Pack and Go功能非常重要。模型文件信息繁多，且可能分布在不不同的文件夹中，只有通过Packand Go功能才能将所有的文件收集到一起，方便整体移动或者传输，否则容易出现文件缺失。

原创 SW simulation解算器选取

通过表5-4的数据对比可以看出，FFEPlus的求解性能非常优秀，但是考虑到解算器稳定性等多方因素，本书建议优先考虑Direct sparse解算器，当求解内存不足时再根据需要首先考虑FFEPlus解算器，若算例在FFEPlus解算器下报错，最后选择Large Problem Direct Sparse解算器。它可以用于大多数计算模型，堪称万能解算器，但是求解方式比较低效，内存利用率一般，所以会导致一些问题，比如速度较慢。4)算例4:大变形问题，第5章的弹簧，强制10mm位移，整体网格密度为0.4mm。

原创 ABAQUS振动切削

大家好，切削仿真由于涉及到大变形、高非线性，因此变的非常复杂，再加上振动切削也许就更复杂啦，现在又出来个椭圆振动，可能真的让有些同学有些吃不消了。大家不要怕，今天我们就通过一步一步的操作和大家共同完成椭圆振动切削的全部过程，帮助大家突破椭圆振动切削仿真的瓶颈，为大家的科研助上一臂之力！！！1、工件建模给定工件长为2.4mm，高为1.2mm，考虑到计算的速度与相关的切削参数，本文对工件进行了模型缩减，并进行了切分。决定取长为0.6mm，高为0.2mm，在工件高1/2处进行切分。建模如下所示：

原创 特定范围所有整数平方和（包含判断输入是否正确为整数）

#include<stdio.h>#include<stdbool.h>/*验证输入为一整数*/long get_long(void){ long input; char ch; while (scanf("%ld",&input) != 1)/*输入不是整数*/ { while ((ch = getchar()) != '\n') { putchar(ch);/*输出错误的输入*/ } printf(" is not an int.

原创 税收方案类别

编写一程序，让用户指定缴纳的税金种类及应纳税收入，自动计算税金。类别 税金 单身 17850美元按15%计，超出部分按28%计 户主 23900美元按15%计，超出部分按28%计 已婚，共有 29750美元按15%计，超出部分按28%计 已婚，离异 14875美元按15%计，超出部分按28%计 如：一位工资20000美元的单身人士，应缴纳的税费为0.15*7850+0.28*（20000-17850）美元。代码：#include<std.

原创 求一个数的约数+ 统计字符数、单词数、行数

#include<stdio.h>#include<stdbool.h>/*用true和false分别代替1 0*/int main(void){ unsigned long num;/*待测试的数*/ unsigned long div;/*可能的约数*/ bool isprime;/*素数标记*/ printf("请输入一个待分析的数："); printf("退出按q\n"); while (scanf("%lu", &num) == 1) {.

原创 PTA 质因子分解 （连续因子 相似）

本题目要求读入一个大于1的整数，编程将其分解成若干个质因子（素数因子）积的形式。输入格式:大于1的整数一个。输出格式:将输入的正整数分解成若干个质因子积的形式，质因子的出现顺序按从小到大排列。如：40=2*2*2*5；如果整数本身为质数或素数，直接输出，如：13=13。输入样例:12480输出样例:12480=2*2*2*2*2*2*3*5*13...

原创 PTA 连续因子（素数）

一个正整数N的因子中可能存在若干连续的数字。例如 630 可以分解为 3×5×6×7，其中 5、6、7 就是 3 个连续的数字。给定任一正整数N，要求编写程序求出最长连续因子的个数，并输出最小的连续因子序列。输入格式：输入在一行中给出一个正整数 N（1<N<231）。输出格式：首先在第 1 行输出最长连续因子的个数；然后在第 2 行中按 因子1*因子2*……*因子k 的格式输出最小的连续因子序列，其中因子按递增顺序输出，1 不算在内。输入样例：630输出...

原创 Matlab 入门

原创 PTA 谁先倒

划拳是古老中国酒文化的一个有趣的组成部分。酒桌上两人划拳的方法为：每人口中喊出一个数字，同时用手比划出一个数字。如果谁比划出的数字正好等于两人喊出的数字之和，谁就输了，输家罚一杯酒。两人同赢或两人同输则继续下一轮，直到唯一的赢家出现。下面给出甲、乙两人的酒量（最多能喝多少杯不倒）和划拳记录，请你判断两个人谁先倒。输入格式：输入第一行先后给出甲、乙两人的酒量（不超过100的非负整数），以空格分隔。下一行给出一个正整数N（≤100），随后N行，每行给出一轮划拳的记录，格式为：甲喊 甲划 乙喊

原创 PTA 兔子繁衍问题

一对兔子，从出生后第3个月起每个月都生一对兔子。小兔子长到第3个月后每个月又生一对兔子。假如兔子都不死，请问第1个月出生的一对兔子，至少需要繁衍到第几个月时兔子总数才可以达到N对？输入格式:输入在一行中给出一个不超过10000的正整数N。输出格式:在一行中输出兔子总数达到N最少需要的月数。输入样例:30输出样例:9#include<stdio.h> int main(){ int n; scanf("%d", &n); if(n==

原创 Average

You are going to read a serial of none-negative integers, which ends with a negative number that does not count as one of the data. Your program calculates the average of all the none-negative integers, and prints the average which rounds to two decimal pl

原创 PTA 上三角数字三角形 及 又来一个上三角数字三角形

输入一个正整数n，输出具有n层的上三角数字三角形。输入格式:只有一个正整数n，1<=n<=100。输出格式:一个上三角数字三角形，每个数字占四个字符位置。输入样例:5输出样例: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15代码：#include<stdio.h>int main(){ int n; scanf("%d",&n);//n层

原创 PTA 作品评分 (冒泡排序法）

全国中小学生Scratch作品大赛拉开了序幕。每个参赛选手可以通过网络直接上传作品。本次比赛人人可做评委。每个网络评委可以通过网络对每一件作品进行打分。评分系统也是请程序高手设计的，能自动去掉一个最高分和一个最低分，求出平均分。输入格式:输入数据包括两行: 第一行为n，表示n个评委，n>2。 第二行是n个评委的打分，分数之间有一个空格。打分是可以带有小数部分的。输出格式:输出平均分，结果保留两位小数。输入样例:610 9 8 7.9 9 9.5输出样例:8.88

原创 PTA 空心字母金字塔+数字金字塔+字母金字塔

输入一个大写的英文字母，输出空心的字母金字塔。输入格式:一个大写英文字母。输出格式:一个空心的大写英文字母金字塔，其中第1层的“A”在第1行的第40列，列从1开始计数。输入样例:E输出样例: A B B C C

原创 PTA 打印沙漏

本题要求你写个程序把给定的符号打印成沙漏的形状。例如给定17个“*”，要求按下列格式打印***** *** * ********所谓“沙漏形状”，是指每行输出奇数个符号；各行符号中心对齐；相邻两行符号数差2；符号数先从大到小顺序递减到1，再从小到大顺序递增；首尾符号数相等。给定任意N个符号，不一定能正好组成一个沙漏。要求打印出的沙漏能用掉尽可能多的符号。输入格式:输入在一行给出1个正整数N（≤1000）和一个符号，中间以空格分隔。输出格式:首先打印出由给定符号组成的最

原创 PTA 等腰直角三角形* + 字母等腰三角形

等腰直角三角形是指一个角是直角，且两条直角边相等的三角形。这里我们输出直角边长为n的等腰直角三角形的格式如下所示: 比如n=1，则输出:*n=2，输出:*** n=3，输出:*** *** 那么，你能用程序来实现么?输入格式:输入一个数n，表示三角形的边长。1<n<1000。输出格式:输出对应的用*表示的等腰直角三角形。输入样例:4输出样例:**********代码#include<stdio.h

原创 PTA Shuffling Machine-洗牌

Shuffling is a procedure used to randomize a deck of playing cards. Because standard shuffling techniques are seen as weak, and in order to avoid "inside jobs" where employees collaborate with gamblers by performing inadequate shuffles, many casinos employ

原创 PTA 整除光棍

这里所谓的“光棍”，并不是指单身汪啦~ 说的是全部由1组成的数字，比如1、11、111、1111等。传说任何一个光棍都能被一个不以5结尾的奇数整除。比如，111111就可以被13整除。 现在，你的程序要读入一个整数x，这个整数一定是奇数并且不以5结尾。然后，经过计算，输出两个数字：第一个数字s，表示x乘以s是一个光棍，第二个数字n是这个光棍的位数。这样的解当然不是唯一的,题目要求你输出最小的解。提示：一个显然的办法是逐渐增加光棍的位数，直到可以整除x为止。但难点在于，s可能是个非常大的数 —— 比如，程

原创 PTA N个数求和（最大公约数、最小公倍数）

本题的要求很简单，就是求N个数字的和。麻烦的是，这些数字是以有理数分子/分母的形式给出的，你输出的和也必须是有理数的形式。输入格式：输入第一行给出一个正整数N（≤100）。随后一行按格式a1/b1 a2/b2 ...给出N个有理数。题目保证所有分子和分母都在长整型范围内。另外，负数的符号一定出现在分子前面。输出格式：输出上述数字和的最简形式 —— 即将结果写成整数部分 分数部分，其中分数部分写成分子/分母，要求分子小于分母，且它们没有公因子。如果结果的整数部分为0，则只输出分数部分。输入

原创 PTA 求e的近似值

自然常数e可以用级数1+1/1!+1/2!+⋯+1/n!+⋯来近似计算。本题要求对给定的非负整数n，求该级数的前n+1项和。输入格式:输入第一行中给出非负整数n（≤1000）。输出格式:在一行中输出部分和的值，保留小数点后八位。输入样例:10输出样例:2.71828180代码：#include<stdio.h>double fact(int n)//阶乘n!函数定义 { double sum=1; while(n>0...

原创 PTA 人民币兑换

1元5角钱人民币兑换5分、2分和1分的硬币（每一种都要有）共100枚，会有很多种兑换方案。请编写程序给出各种兑换方案。输入格式:输入为一个正整数n，表示要求输出前n种可能的方案。方案的顺序，是按照5分硬币从少到多排列的。输出格式:显示前n种方案中5分、2分、1分硬币各多少枚。每行显示一种方案，数字之间空一格，最后一个数字后没有空格。注意：如果全部方案不到n种，就顺序输出全部可能的方案。输入样例:5输出样例:1 46 53 2 42 56 3 38 59 4

原创 PTA 换硬币

将一笔零钱换成5分、2分和1分的硬币，要求每种硬币至少有一枚，有几种不同的换法？输入格式:输入在一行中给出待换的零钱数额x∈(8,100)。输出格式:要求按5分、2分和1分硬币的数量依次从大到小的顺序，输出各种换法。每行输出一种换法，格式为：“fen5:5分硬币数量, fen2:2分硬币数量, fen1:1分硬币数量, total:硬币总数量”。最后一行输出“count = 换法个数”。输入样例:13输出样例:fen5:2, fen2:1, fen1:1, total:4