不得不说，Mathematica真是个好东西，以前学习有限元的时候，对于书中的方程推导，看到了就看过去了，从没有想过要自己推导一遍，原因是手工推导太复杂。有了MM，原来很复杂的东西突然变得简单了。 1.单元几何描述 上图是纯弯梁单元，长度l，弹模E，面积A，惯性矩I。两个节点1和2的位移列阵为 \[...
分类：
其他好文 时间：
2015-08-11 21:02:31
阅读次数：
277
不得不说，Mathematica真是个好东西，以前学习有限元的时候，对于书中的方程推导，看到了就看过去了，从没有想过要自己推导一遍，原因是手工推导太复杂。有了MM，原来很复杂的东西突然变得简单了。 1.单元几何描述 上图是纯弯梁单元，长度l，弹模E，面积A，惯性矩I。两个节点1和2的位移列阵为 \[...
分类：
其他好文 时间：
2015-08-09 22:03:30
阅读次数：
123
表达式是变量之间的运算关系，表达式求值就是对变量赋值并运算出结果的过程。针对于Mathematica中的表达式，有两种方法获得其值。
一是对相应的变量直接赋值，这样对应的表达式在调用时便会直接运算得到结果。改变变量的值，表达式本身的值也会相应改变。不过这样做有一个缺点，那就是表达式本的关系被本身的值覆盖了，通过?才能获得表达式的关系
第二种方法是通过变量替换临时给表达式的变量赋值 他的语法规则如下：...
分类：
其他好文 时间：
2015-07-22 16:10:59
阅读次数：
145
使用%号做标记，获得文本界面之前的运算结果：...
分类：
其他好文 时间：
2015-07-22 13:22:58
阅读次数：
140
Mathematica 中获取帮助非常简单，因为应用程序本身就自带了帮助文档，在documentation中搜索一下就可以得到非常详细的结果。不过相比documentation来说还有一种简洁的方法，可以获得函数的功能说明。
那就是在函数前加一个问号并运行。如下：...
分类：
其他好文 时间：
2015-07-21 20:44:41
阅读次数：
154
mathematica 中为变量的赋值时永久性的。
为避免程序的程序错误和节省资源，mathematica 中用完一个变量后就应该马上清除他。有两种快速清除变量值的方法:
var=.
Clear[var]...
分类：
其他好文 时间：
2015-07-21 18:35:37
阅读次数：
194
问题碰到这样的问题，感觉很神奇。
定子方程，短幅内摆线方程：
{x1y1=(R?r)sinτ+esin(z2τ)?resinθ=(R?r)cosτ?ecos(z2τ)+recosθ\left\{\quad
\begin{array}{rl}
x_1&=(R-r) \sin \tau+e \sin \left(z_2 \tau\right)-r_e \sin \theta \ y_1&=(...
分类：
其他好文 时间：
2015-07-18 18:38:21
阅读次数：
204
问题看到有人问：
已知拱高10.4米的圆弧弧长36.7米，求圆的半径和圆心角（最好附带Matlab程序）。怎么限制x是锐角？
解答先根据大致情况推导一番
因为rr跟θ\theta有如下关系，只须求一个rr即可：
s=r?θs=r\cdot \theta可以利用的关系从图上看出还有：
x2=h(2r?h)x^2=h(2r-h)
sinθ2=xr\sin{\dfrac{\theta}{2}}...
分类：
其他好文 时间：
2015-07-15 11:12:35
阅读次数：
193
说画就画ClearAll["Global`*"]; x = Pi;
rotateMeHarder1[g_, vertical_, viewpoint0_, center_List: {0, 0, 0},nframes_Integer: 15, opts : OptionsPattern[]] :=
Module[{grlist}, grlist = Table[Show[g, ViewVer...
分类：
其他好文 时间：
2015-07-14 13:37:14
阅读次数：
115
多年来，我都对戈特弗里德·莱布尼茨很感兴趣，尤其是因为早在3个世纪以前，他就似乎想要制造一种类似Mathematica和Wolfram Alpha的工具，没准还可能会写本《一种新科学》。所以，在最近一次德国之旅中，我对能够拜访坐落在汉诺威的莱布尼茨文献馆感到兴奋不已。翻阅着他发黄的手稿（仍旧够挺，经...
分类：
其他好文 时间：
2015-07-13 06:32:58
阅读次数：
175
