最新ANSYS地震响应分析讨论.doc

《最新ANSYS地震响应分析讨论.doc》由会员分享，可在线阅读，更多相关《最新ANSYS地震响应分析讨论.doc（146页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-dateANSYS地震响应分析讨论地震响应分析地震响应分析1模态组合就是根据模态分析中的几阶振型（也可以少于这几阶，看你要求的精度）进行组合（类似于结构最不利组合），从而求出地震响应的最大值。2组合各振型反应的最大值，求得结构地震响应的最大值。这个问题在论坛上已经有很多人问过，也有各种各样的回答，但是至今没有令人满意的解答。我自己试过很多种方法，加上论坛上其他人提到的方法，大

2、致归类如下： 1.先做静力恒载工况分析，打开预应力pstres开关；然后转到时程分析。 结果：恒载对后面的时程计算不起作用，时程计算依然从0开始。 2.直接在antype，trans中考虑恒载：先把timint，off加acel,9.81,打开应力刚化，sstif，on，lswrite，1，然后timint，on开始时程计算。 结果：恒载9.81起作用了，但结果是错的，它被积分了。 3.不用什么前处理，直接把9.81加在地震波上acel，9.81+ac（i）。 结果，同2，9.81带入了积分，这个9.81相当于阶跃荷载，而不是产生恒载。 4.ansys帮助中施加初始加速度的方法（篇幅限制请自己

3、看帮助）。 结果，同2、3，9.81还是带进时间积分。 5.这种是我受到别人的启发，通过结构受ramp荷载的特点施加的，可以近似的解决问题。 即1）求出结构的自振一阶频率w 2）令tr1/w 3) 定义ramp荷载为从0到tr加到9.81，然后在整个时间积分中保持不变 4）antype,trans中分几个荷载步将荷载从0加到9.81 5) 在随后的荷载步中acel，9.81+ac（i） 这种做法虽然也是将9.81+加到地震波中，但是因为满足TR的要求，所以这个动力效应被削弱到了静力效应，它作用在结构上就像静载一样。对于单自由度结构理论上跟静载是完全一样的，但是多自由度会子静力效应上下很小的范围

4、内波动，所以可以认为相当于静载的作用，这样我们就可以达到考虑恒载的目的了。 第5种是我至今为止考虑恒载的做法，我也很想知道还有没有更简单精确的方法，或者在前4种方法中就有只是我使用不正确，希望大家能一起来讨论，彻底解决这个问题。谢谢！地震反应怎么考虑重力SOLU ANTYPE, TRANS TRNOPT,FULL TIMINT,OFF !*先关闭时间积分效应 TIME,1E-8 !*设一个极短的积分时间 acel,9.8 NSUBST,2 !有时候子步数要增大 KBC,1 LSWR,1 !*把这个写入第一步 TIMINT,ON !*然后再时间积分效应开关，以后就正常写载荷步了 这种方法应该是对

5、的,ANSYS帮助文件中也有提到, 可是,有一个问题:由于是阶跃荷载,就会产生动力效应,整个结构的变形大于实际的情况吧?这样与实际结构在重力下受到的变形就不一样了! 我想,之所以用KBC,1的原因是产生阶跃荷载,第一个子步产生变形,然后后几个子步使其稳定,不产生初始速度吧 自己作了一个例子, 终于彻底搞明白啦, KBC,1是产生阶跃荷载,不过,关键就在于此时关掉了时间积分效应,TIMINT,OFF那么,此时就相当于是静力分析了,当然就不会产生动力效应了,整个结构的变形就是实际变形.另外KBC,1的作用是不产生初始速度. 总结下:非线性时程分析时,结构在地震波来之前,已经存在初始位移和初始内力,

6、但是没有初始速度和加速度.因此,一定要注意KBC,1,同时TIMINT,OFF ,设置2个以上的子步,很小的时间. 一楼的总结有很多不妥的地方, 你那个方法不知道正确不,而且太麻烦了. tmh_cn在第12楼总结的方法是对的,其他人说的方法都有问题. 另外,不用打开应力刚化 .fini/config,nres,20000!/config,nproc,2lstep=1500 !el*DIM,aei,ARRAY,lstep,1,1, , ,*CREATE,ansuitmp*VREAD,aei(1),el,txt,1,(e16.6) *END/INPUT,ansuitmp csys,0/solu a

7、ntype,transtrnopt,fullTIMINT,OFF !*先关闭时间积分效应 SSTIF,OnKBC,1TIME,1e-3!*设一个极短的积分时间 acel,8.7 NSUBST,2!有时候子步数要增大 solvetimint,onkbc,0 !TINTP, ,0.25,0.5,-1,0.5,-1,ki=0.05kj=0.05oi=2*3.14159*0.09oj=2*3.14159*0.16ma=2*(ki*oj-kj*oi)*oi*oj/(oj*oj-oi*oi)mb=2*(kj*oj-ki*oi)/(oj*oj-oi*oi)ALPHAD,ma,BETAD,mb, allsel

8、,all*do,tm,1,1500,1time,(tm)*2/100NSUBST,1, , ,1 AUTOTS,1acel,aei(tm)/314.7 !EI波，只输入X向加速度solve*enddofini*关于动力学分析模拟自重对结构的影响，这个问题，现在偶已经解决这个问题啦，几行命令流就可以了：/PREP7ACEL,0,9800,0, !*先把重力加速度在预处理设好，以后就不在设了FINISH/SOLUANTYPE, TRANSTRNOPT,FULLTIMINT,OFF !*先关闭时间积分效应TIME,1E-8 !*设一个极短的积分时间NSUBST,2KBC,1LSWR,1 !*把这个写

9、入第一步!*TIMINT,ON !*然后再时间积分效应开关，以后就正常写载荷步了.虽然，这只是个很简单的问题，但如果不注意，会使计算的动力响应值受重力影响而偏大，我当时就很郁闷，怎么也不能把重力的影响去掉，现在终于搞定了！现与大家分享偶的一点心得，同时，也感谢各位大侠的关注：-） 模态叠加法（时间可能取的太短了）求地震反应! 步骤1：进行模态分析/SOLU! 1:MODE ANALYSISAnType,ModalModOpt,Lanb,20! EXTRACT 20 MODESMxpand,20,1SOLVEFINISH! 时程分析命令流! 采用模态叠加法*set,st,0.01*set,end

10、point,100*dim,xforce,array,endpoint*vread,xforce(1),Seismic_x,txt(7F10.5) *VOPER,xforce(1),xforce(1),MULT,9.8, , ,/SOLU !* ANTYPE,4TRNOPT,MSUP,20,1 ! ! ! 施加位移约束NSEL,S,LOC,Z,0,0! 选取底层所有节点D,ALL, , , , , ,ALL, , , , ,! 约束所有自由度 ALLSEL,ALL!* *dim,shijian,array,endpoint*vfill,shijian(1),ramp,st,stNt_star=

11、1Nt_end=100Nt_incr=1*DO,Nt,Nt_star,Nt_end,Nt_incrTIME,shijian(Nt)ACEL,xforce(Nt),NSEL,ALL!LSWRITE,NtSOLVE*ENDDOFINISH!模态扩展/SOLU EXPASS,ONNUMEXP,all,0,1OUTRES,all,allSOLVEFINISH*谱分析我想，作谱分析之前应该先弄明白谱分析的工程意义。实际上，谱分析是一种线性化的动力分析思路，简单点说，就是在获得结构自振特性以后，用一条符合本地场地土特性和设防烈度（地震响应谱分析）的地震谱或者一条符合本地风场特性的风速谱，作为外部作用施加与

12、结构上，以求得结构在各个频率对应的谱值作用下的结构响应。他实际也是一种时程分析，只不过时间轴为频率值，所以也称为频域分析。由于反应普的获得时基于结构的线性化假定，所以谱分析只能对线性结构进行分析。窗体顶端先试贴，如有人需要给留言发送给你 Q为问，A为回复。-Q:模态分析得到的结果是不是某个方向上的各阶频率啊？我要得到各个方向的一阶频率能做到吗？A：模态分析得到的结果是你所选择的自由度内的振型。如果需要得到某一个方向内的振型，用RUDECED法，选择你所希望的主自由度，如UX。将所求的频率值设置为一即可。Q：不过大型复杂结构的振型一般都不是一个方向的，应该是一个方向为主，其它方向为次，如果把所要

13、求的方向设为主自由度，会不会跟真是实际结果有出入？Q：对于复杂的大型实体，其主自由度无法给出，此时就无法应用reduce法。可不可以改约束？还是就认为一次得到前3阶频率就是x,y,z方向的一阶频率啊？A：正是因为复杂结构的主自由度不好确定，所以我还是倾向于用SUBSPACE和BLOCK LANCZOS法。约束不要随便改，它应该描述结构真实的状态。不能说前3阶频率即为X、Y、Z方向的第一阶频率。你可以好好想想所用方程中K和M的意义。-Q：轴承是用来支撑要研究的对象，把轴承等价为弹簧，将其刚度加在对象的节点上。这样就要定义刚度的大小在加。是否还要定义单元类型，然后才能定义其刚度？A：我曾经用LIN

14、K8单元模拟过，根据K=EA/L确定刚度，我当时取E=210e9,L=0.1,然后A=KL/E。经过检验，是可行的。A：1、用link单元模拟其刚度 ，k=EA/L，通过假定E,A,L的值使之和k相等，其密度要0；2、用combin14模拟其刚度，阻尼0-Q：ANSYS处理动力问题（如土层地震响应、或基础振动反应），可否从边界施加加速度荷载？若可以，如何？瞬态动力反应分析，时间子步为>200?A：如果你所希望求的最高频率的周期为f，则ANSYS取时间步长为1/20f为默认值。Q：板单元的后处理的显示在单元坐标还是在总体坐标? 比如位移,应力？A：应力位移在总体坐标，弯矩不一定。A：应力、

15、位移是显示在结果坐标系的，顺便说明一下，节点信息显示在显示坐标系。不过结果坐标系默认为全局坐标系，而显示坐标系也是默认为全局坐标系。你可以这样修改结果坐标系：Main Menu:General Postprocessor>Options for Output.你可以这样修改显示坐标系：Utility Menu:Workplane>Change Display CS toA：你可以用ETABLE命令，先存储对应于各单元的位移，等等，再用ANSYS的柱状图，或是别的画图软件PLOT出图形，你认为如何-Q：我现在需分析一个板梁结构，板已用SHELL63单元划分好；梁我是用板上的一条线划分

16、单元并添加截面而生成的。但现在运算时，发现板和梁是分开的，它们之间互不影响。请教各位高手，怎样将板和梁合并为一个整体？A：Preprocessor>Numbering Ctrls>Merge Items里element and node 合并-Q: 用SHELL63单元划分的板结构，想加垂直于板面的匀布载荷?A: 找到如下位置：ANSYS HELP->ANSYS ELEMENTS REFERENCE->ELEMENT LIBRARY->SHELL63 ELASTIC SHELL。在INPUT SUMMARY 中，有SURFACE LOADS PRESSURES。其

17、中，face 1 (I-J-K-L)(bottom,in +Z direction)的意思为：当Load Key=1时,pressure加在由节点IJKL确定的面上，在单元坐标中以+Z方向为正，也可以说垂直指向面IJKL时为正。至于均布，则将几个值都填成一样的就可以了，或者只填主要参数如value at I。自己多看看帮助文件，做一个例子验证一下。-Q: ANSYS可以自己定义输出格式吗?A: 可以，先打开一个文件*cfopen，再用*vwrite,label通过打开一个文件，读出数据，或以追加的方式写入数据，并可控制输入和输出的格式*cfopen,*vread,(f12.6)*vwrite,

18、(f12.6)*cfclos,-Q:作模态分析时，就是集中质量单元和刚度加法？A：在选用求解方法时，例如SUBSPACE和BLOCK LANCZOS，会有一个选项：LUMPM，默认为选用一致质量矩阵，你可以选为集中质量矩阵。A：用mass21，通过mesh keypoint来把点划分成mass21单元即可-Q：我在原有的静态分析的基础上，想继续作一下结构的模态分析。应该采取那些步骤啊？我用solution>new analysis-modal,then current LS,出现如下的提示：1.the eigenvalue extraction method must be specif

19、ied on MODOPT command. SUBSP or LANB are usually recommanded.2.please expand modes if any downstream MODAL analysis based calculations need to be performed.这是怎么一回事？A：你最好先看看帮助文件中的 ANSYS STRUCTURE REFERENCE中的BUCKLING ANALYSIS。如果你对理论不清楚，我劝你不要盲目计算。不然结果也许并不可信。这是我的建议。A：系统默认为0阶!-Q：1.在ANSYS中如何设定小数的有效位数？比如我输

20、入的应力值为0.0004495568而且我的一组数据都是最后几位差别，但是ANSYSY自做聪明把它四舍五入了，结果弄的都一样了，如何更改这个有效位数呢？2.ANSYS在输出时如何修改输出的数字显示模式？比如用科学记数法，或者将数字缩放多少倍那样显示？在什么地方能够修改呢？3.ANSYS5.7.1中什么地方能够出现那种材料库表？就是说不同ANSYS的什么material models对话框生成，而是通过table输入材料属性诸如应力应变曲线之类的东西?A：用a.bEn来表示，如2.45678e-9。材料特性用tb-Q：我是个ANSYS初学者，现在碰到下面几个问题，请各位高手帮忙看一下：1：是什么

21、原因使得相邻ELEMENT的共用节点解得的力不一样，应该是方法步骤的原因而不是误差吧？请简单介绍一下。2：假如我想用ANSYS模拟陶瓷的烧结过程，是不是有实际*作的可能性？如果有，其基本的原理和模型怎样？3：ANSYS可以用来模拟热循环下的应力分配不均造成的焊点（63%Sn&37%Pb)失效的过程，我看到有一篇文献介绍用在Sn球里随机分布Pb相和其它共晶相的方法建模的，那么，这个随机过程在ANSYS中怎么实现的？而且如果要生成足够多的晶相小颗粒，是不是只能手动一个一个的建立？A：可能是你提取数据的时候思路有问题，比如说，当你提取的是有多个连接单元的节点数据，而实际上你提取的是各个单元的。A：相

22、邻ELEMENT的共用节点你MERGE了吗？ANSYS是通过单元的节点传递数值的，从而得出结果的。-Q：在进行模拟计算结构多阶段施工过程计算时,在ansys里面说可以用死活单元的方法，但是死活单元每次只允许选择一个,我的模型比较大,如果采用这个方法不太现实,各位同志有什么比较好的方法没有？请赐教！A：涉及到两个方面的问题：一、以TABLE形式定义参数。二、DO-LOOP循环。你可以参考ANSYS GENERAL GUIDE中SOLUTION中对多重载荷步的叙述，很简单。-Q：比如在beam4和shell63中显示的输出结果有应力和变形，以及节点等效荷载。请问如何才能显示截面弯矩，扭矩和剪力呢？

23、A：在ELEMENT TABLE 中，首先定义变量。注意一：你要知道BEAM设置了几个积分点。注意二：要知道SEQUENCE NUMBER值。所以你需要参考ELEMENT REFERENCE 中的单元介绍。-Q：最近在计算一个模型，管状，在地下，分析中要施加地层给予的支撑，我想简化为弹性支撑。可是发现在ansys中找不到这样的支撑情况！A：首先，你要知道弹性支撑的弹性系数，具体的计算方法你应该知道的 。然后，用LINK8单元模拟，要注意的是刚度k=EA/L，由此确定LINK单元的面积A来体现单元的弹性刚度。我对此作过验证，结论是可行的。Q：用弹簧，可是问题很多，因为是一个管体，整个管壁都有支撑

24、。简化为弹簧出现的问题很多了。比如，用多少个合适？隔多长合适，线性支撑和非线性支撑能等同吗？我最近看校验手册发现V135这个例子可以借鉴，可是还是有许多麻烦。哪位大侠如果有这方面经验请指点指点！-Q：如何梁单元同时施加两个方向分布荷载，（竖向，侧向）？A：同时施加两个方向的均布荷载相当与分别施加不同的均布荷载，即分别加。主要是要选择好均布荷载的ID，参考单元参考手册中的说明。例如：BEAM4 INPUT SUMMARY 中：Surface Loads Pressure: 有 face2(i-j)(-Y normal direction) 这里，2就是LOAD ID，以单元坐标的-Y轴为正。你试

25、一试看。-Q：请问如何打印图形结果？A：用ANSYA中的功能就可以拉，GUI菜单命令：PLOTCTRLSCAPTURE IMAGES或RESTORE IMAGESA：我是用plotctrl>hard copy>to file,然后在剪贴板中就可以看见了。A：我都是用抓图工具抓下来，然后再打印的。A：/ui,copy,save是一个截取屏幕图形的命令,现在总算又实现了gui肯定可以通过相应的命令方式来实现的观点。A：Ansys 5.6已经做到了把图形窗口COPY到文件如：JPG/BMP.等格式，PLOTCTRL/HARD COPY/TO FILE。很好用！-Q：其实APDL语言是AN

26、SYS的二次开发语言，是在FORTRAN编程语言的基础上形成的，你只要会FORTRAN语言，这就很容易理解APDL语言了，也有循环，选择等编程结构，把枯燥的*作用计算机自动实现。然后形成一个宏*作拉。Q：是的，APDL语言确实很方便，而且和FORTRAN语言类似，很容易读懂。在计算电机电磁场中电机的反电势更加方便，前几天刚算了一个需要循环运算288次的程序。使用APDL语言使的我可以让电脑自己运行程序，保存结果数据，而我的任务就是睡一觉醒来，对结果数据进行处理。Q: 我看见好多人在讨论用命令流处理程序，很不明白命令流和gui的形式有很大区别吗？在什么情况下更适用命令呢？A：很容易控制Ansys

27、，且结果文件小；功能很强大；可以进行循环控制，分支等*作；运算速度快，充分利用APDL语言，对ansys进一步开发。-Q：用ansys5.6.2计算一个压力容器,划分网格9万余个.运行环境win2000,cpu:1G,内存640M,硬盘40G,模型建立,网格划分均无问题.但是运行到占用硬盘大约9.2G的时候,就出错.提示内存无法read.请问该如何解决?是内存的问题?还是ansys本身问题?还是系统问题?硬件问题?模型建立实在无法简化.A：我的硬盘是30G的，内存1G，OS windowsNT,Total Workingspace是1400M,再高就控制不了，Database 是56,我能解的

28、有效节点数是140000,只是因为磁盘空间不够，工作的分区是0G,当ansys探测到磁盘空间已满，就自动退出。尽量大的增大节点数，并不是好办法，因为运行速度太慢了。利用对称性减小求解区域，利用Lesize,esize等控制网格，都是好的办法。A：选用对称性的时候要注意边界条件。但的确是个好办法。A：以前在使用ALGOR时也遇到过类似情况，当时有10万个结点。对WIN2000来说，已不存在WIN9X中文件最大只能2G的瓶颈，但在处理超大虚拟内存需求时仍然会出问题。处理方法有二：1：使用PQMAGIC划分出一个20G的逻辑分区（一定要跟WIN2000系统在不同的分区），再在控制面板-系统中将虚拟内

29、存文件指定到这个大分区上，以尽量避免发生系统堆栈错误（虚拟文件大小仍交WINDOWS控制）。2：将至少两个分区的硬盘升级至NTFS，由于NTFS文件可以跨分区存储，且要小10%左右，比FAT32要好用得多。-此帖已经发过多遍，请注意发帖规则日月星辰 积分 28帖子 38#22004-6-17 20:42 Q：我在求解是出现如下的提示：1.real constant 1 referenced by at least element types 1 and 2.2.coefficient ratio exceeds 1.0e8.我全部忽略上述提示，得出了结果。结果能可信吗？我有错在哪里呢？模型我都

30、检查了数遍了，还是和初始文件名有关？我是在一个初始名下面建立若干相关的模型。由于我的单元只是板梁的结合，所以，我重新检查了一下。删除了一个影响不大的扁梁，（很扁，尺寸不大，系统可能认为是板了）运行就通过了。上述的提示可能是因为系统对某些单元有默认的尺寸限制。A：有些单元类型在默认值状态下不需要输入REAL CONSTANT。但是当你改变了某项选项，你就需要重新输入REAL CONSTANT。你可以用UTILITY>SELECT>ENTITIES命令中的SELECT ELEMENT BY ATTRIBUTE来查看单元材料，REAL CONSTANT等等是否定义明确，然后用APPLY

31、， REPLOT命令来显示每个不同材料的单元。A：我试了一下，按照你所说的步骤。出现了如下提示：a default command range is not provided for the ESELcommand.the ESEL command is ignored.所以，没有任何结果可见。是不是需要对某些default的命令做些改动呢？A：.real constant 1 referenced by at least element types 1 and 2.此句意思好象是说，你为1号单元和2号单元或是别的单元定义的real constant 均为1，你需要检查你的模型的单元类型和实常

32、数部分。A：第二个问题经常出现.第一个问题有点问题,你看看单元定义对没有如果单元都没有划分对,结果当然不可靠啦。-Q：slab on girder highway bridge 桥板为SOLID单元，由于girder上部flange嵌贴在板的下表面，所以建模时，我将桥板定义为SOLID单元，桥板下表面AREA定义为SHELL63单元。但SOLID和SHELL单元的自由度不同，这样建模肯定有不对的地方。运算结果与实验结果也相差了40%。有什么好的解决办法。A：你可以根据需要为不同的单元增加自由度（DOF），顺便问一句，你的模型的两个部分的节点是否完全重合，提高计算精度要用更细的网格划分，另外，要

33、保证模型中每个单元的长宽比和理，是否有使用大应变的需要，在ANSYS非线形手册中有影响精度的说明。A：本来SOLID体的下表面被我定义为SHELL，这时面单元的节点同时也是体单元下表面的节点。但由于两种不同的单元自由度不同，所以我放弃了这个模型。另建模型的SHELL单元和SOLID单元相贴面节点坐标完全相同，节点号不同，这样才可能把它们耦合，但是结果仍不理想。可有什么方法?A：比较麻烦,可以通过编程输入命令流!QA：我已经把SHELL和SOLID对应点的UX，UY，UZ全部耦合起来，让这些NODE的X，Y，Z三个方向的位移相同。但是结果仍然不理想。我想约束方程的功能也进乎如此，有没有更进一步的

34、建议，十分感谢。-Q：ansys有没有评价一个网格划分好坏的参数或命令？请告知菜单*作。A：PREPROCESSOR>MESHING>CHECK MESH>两个命令中一个用来CHECK网格划分，另一个用来CHECK连贯性。只是CHECK完了别哭就好。A：PLOT BAD ELEMENT!-Q：(1)在gui下运行宏时,重新定义原来已有的数组,ansys就会有个提示,想不要这个提示,怎么办?(/uis,msgpop,3这个行不通,我试过了)(2)ansys中数组的定义为a*b*c a<10040 b<255 c<7 而我实际上要定义更大的有什么好办法,难道非要

35、写到文件里再读进来?A：q1:有两种办法（1）在ultility menu / menuctrls/ message controls下设 (2)将 /uis,msgpop,3 加入到start.ans文件中去，以后每一次运行ansys时都会默认这个设置-Q：我现在做一个板壳的屈曲分析，想要得到该结构在荷载作用下的极限承载力，利用非线性屈曲分析，请问极限承载力的荷载比例因子在什么命令下可以查到。A：非线性屈曲分析与特征屈曲的不同之处是，它一开始就挠屈变形，随着荷载的增大，变形日益增大。我们认为结构最危险点的变形时间历程曲线中变形若从某点开始发生较大的转折，与先前完全不一样，就发生了屈曲，此时的

36、荷载就是屈曲荷载。故不存在荷载比例因子，即特征值问题。A：要进行非线性屈曲,必须加到可以使结构发生屈曲的荷载,即荷载要比结构的实际承载力要大,在进行非线性屈曲分析(极限承载力分析)时荷载逐步增加到破坏为止.计算时一定要输出每一步的计算结果,这样在post1中就会有一系列的荷载系数,而最后收敛的那个系数就是相对于极限承载力的系数,极限承载力=系数和总荷载的乘积。-Q：如果材料一定为各向异性，我用的材料为钢板，该如何选用单元类型？A：只要材料参数中有Ex,Ey,Ez单独输入项，就能考虑各向异性。A：波纹板的EX、EY可以输入，但是EZ该如何确定呢？我自己练习了一下，如果不输入EZ，则计算就不能进行下去，提示我材料的弹性模量要输入正值。这是为什么？A：如果是钢板，一般认为是各向同性的。厚板因轧制次数少，晶体不密实存在着Z向（厚度方向）弱化。对厚板，一般地Ex=Ey，Ez较小，通过试验确定。对薄板<16mm，认为各向同性。A：不管是什么板，材料的本构关系是微小单元之间应力与应变