第5章力学准则法PPT讲稿.ppt

上传人:石*** 文档编号:43657843 上传时间:2022-09-18 格式:PPT 页数:40 大小:2.61MB
返回 下载 相关 举报
第5章力学准则法PPT讲稿.ppt_第1页
第1页 / 共40页
第5章力学准则法PPT讲稿.ppt_第2页
第2页 / 共40页
点击查看更多>>
资源描述

《第5章力学准则法PPT讲稿.ppt》由会员分享,可在线阅读,更多相关《第5章力学准则法PPT讲稿.ppt(40页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第5章力学准则法第1页,共40页,编辑于2022年,星期一其本思路其本思路充分发挥材料的强度潜力,刚度潜力和贮能能力,使结构材料最省。充分发挥材料的强度潜力,刚度潜力和贮能能力,使结构材料最省。力学准则法只能寻求结构的最小体积或近似最小体积设计。力学准则法只能寻求结构的最小体积或近似最小体积设计。当结构各部分的容重当结构各部分的容重 大体一致时,最小体积设计就等于或接近最轻设计。大体一致时,最小体积设计就等于或接近最轻设计。分类分类1 1)等强度准则)等强度准则满应力设计满应力设计 要求尽可能使结构在使用过程中各构件的要求尽可能使结构在使用过程中各构件的 最大应力都能达到其允许值。最大应力都能

2、达到其允许值。2 2)同步失效准则)同步失效准则满约束准则满约束准则满约束设计满约束设计 要求在结构设计时使尽可能多的不等式约束同时达到临界。要求在结构设计时使尽可能多的不等式约束同时达到临界。当不等式约束主要是应力约束时,满约束准则就退化为满当不等式约束主要是应力约束时,满约束准则就退化为满 应力准则,后者是前者的特殊情况。应力准则,后者是前者的特殊情况。3 3)能量准则)能量准则 尽可能充分发挥材料的贮能(应变能)的能力。尽可能充分发挥材料的贮能(应变能)的能力。材料的贮能能力是与其强度和刚度相关联的,材料的贮能能力是与其强度和刚度相关联的,因此能量准则常常与满约束准则相一致或接近。因此能

3、量准则常常与满约束准则相一致或接近。第2页,共40页,编辑于2022年,星期一优点优点物理概念清楚,与过去的设计思想相衔接,容易为工程设计人员所接受;物理概念清楚,与过去的设计思想相衔接,容易为工程设计人员所接受;算法简单,迭代收敛较快,且结构重分析的次数与设计变量的数目无关,算法简单,迭代收敛较快,且结构重分析的次数与设计变量的数目无关,比较适合中小型和大型结构的优化设计。比较适合中小型和大型结构的优化设计。缺点缺点1 1)适用范围较窄,只能用于最小体积设计或最轻设计;)适用范围较窄,只能用于最小体积设计或最轻设计;2 2)在某些情况下有失效的可能性;)在某些情况下有失效的可能性;3 3)没

4、有直接建立与目标函数的关系,并不能保证使目标函数最小。)没有直接建立与目标函数的关系,并不能保证使目标函数最小。第3页,共40页,编辑于2022年,星期一5.2 5.2 满应力设计法满应力设计法一、满应力设计法的基本概念一、满应力设计法的基本概念 1)满应力设计的对象:)满应力设计的对象:满应力设计的对象一般是结构布局已订并具有应力约束和尺寸约束的结构。满应力设计的对象一般是结构布局已订并具有应力约束和尺寸约束的结构。2)满应力设计准则:)满应力设计准则:严格满足应力设计:严格满足应力设计:当只有应力约束时,满应力设计要求结构的每一个构件至少在一个工况当只有应力约束时,满应力设计要求结构的每一

5、个构件至少在一个工况 下达到满应力,即至少在一种荷载状态下应力等于其允许值。下达到满应力,即至少在一种荷载状态下应力等于其允许值。广义应力设计:广义应力设计:当还有截面最小尺寸的几何约束时,则要求每一构件在强度约束和尺寸当还有截面最小尺寸的几何约束时,则要求每一构件在强度约束和尺寸 约束中至少使其中一个达到临界。约束中至少使其中一个达到临界。3)静定结构的满应力设计:)静定结构的满应力设计:设某静定结构有设某静定结构有 根杆件,根杆件,种工况。种工况。第第 杆在第杆在第 工况中最不利内力的绝对值为工况中最不利内力的绝对值为 。设其中最不利者为设其中最不利者为 ,即,即第4页,共40页,编辑于2

6、022年,星期一(一)严格满应力设计(一)严格满应力设计 选择设计方案选择设计方案,使,使则则,可以证明:静定结构的严格满应力解,可以证明:静定结构的严格满应力解=最轻解。最轻解。(二)广义满应力设计(二)广义满应力设计(三)受弯构件(三)受弯构件第5页,共40页,编辑于2022年,星期一4)超静定结构的满应力设计)超静定结构的满应力设计 1.静定满应力解总存在,且不止一个;静定满应力解总存在,且不止一个;2.2.在单工况下,一般地说,不可能使全部构件达到满应力;在单工况下,一般地说,不可能使全部构件达到满应力;3.3.能否实现满应力设计,不仅取决于结构,而且与荷载和工况数有关;能否实现满应力

7、设计,不仅取决于结构,而且与荷载和工况数有关;4.4.保持原结构布局的超静定严格满应力解不一定存在。保持原结构布局的超静定严格满应力解不一定存在。5.5.超静定结构满应力解存在性的判别超静定结构满应力解存在性的判别 存在原结构布局严格满应力解的必要条件为:存在原结构布局严格满应力解的必要条件为:式中:式中:工况数;工况数;杆件数;杆件数;满应力度:满应力度:当允许对结构施加任意荷载时,在一个工况下可以达到满应力的当允许对结构施加任意荷载时,在一个工况下可以达到满应力的杆数的最大值。也就是在结构位移可行域边界各个顶点中,最多汇交的满应杆数的最大值。也就是在结构位移可行域边界各个顶点中,最多汇交的

8、满应力约束面的数目。力约束面的数目。满应力度:满应力度:6.由于超静定结构的内力与各杆件的截面性质有关,由于超静定结构的内力与各杆件的截面性质有关,没有一个简单的方法可以直接求出满应力解。一般可用迭代的方法逐步逼近满应力解。没有一个简单的方法可以直接求出满应力解。一般可用迭代的方法逐步逼近满应力解。第6页,共40页,编辑于2022年,星期一二、应力比法二、应力比法1)估选初始设计方案(选定初始点)估选初始设计方案(选定初始点)初始点:初始点:进行力学分析进行力学分析,求出各构件在各工况中的最不利应力求出各构件在各工况中的最不利应力 式中式中 为第为第 杆在第杆在第 工况的应力。工况的应力。为第

9、为第 杆在各工况中的最不利应力。杆在各工况中的最不利应力。2)修改截面)修改截面将将与其允许值与其允许值相比较。相比较。(1)(1)若若,表示材料未充分发挥作用,可减小其,表示材料未充分发挥作用,可减小其截面截面(2 2)若)若,表示材料超载,应增加其截面,表示材料超载,应增加其截面 第7页,共40页,编辑于2022年,星期一假设:假设:修改截面后各杆件内力不变。修改截面后各杆件内力不变。且且 达到满应力,即达到满应力,即 ,因此,因此式中式中 称为应力比。称为应力比。可用可用作为下一轮迭代循环的初始点作为下一轮迭代循环的初始点 迭代公式:迭代公式:若还规定有截面尺寸约束若还规定有截面尺寸约束

10、 迭代公式:迭代公式:这种逐次逼近的方法称为这种逐次逼近的方法称为应力比法应力比法,或,或比例满应力法比例满应力法。第8页,共40页,编辑于2022年,星期一3)迭代终止条件)迭代终止条件改进措施改进措施采用应力比法求超静定结构满应力解得迭代次数较多,一般需要十几次,采用应力比法求超静定结构满应力解得迭代次数较多,一般需要十几次,改进措施之一是在迭代公式中引入一个改进措施之一是在迭代公式中引入一个“超松弛系数超松弛系数”,以加快收敛速度:,以加快收敛速度:超松弛系数,或应力比指数超松弛系数,或应力比指数 对于拉杆,对于拉杆,可取可取 值大于值大于1(如(如1.3左右)。然后逐步减小到使之趋于左

11、右)。然后逐步减小到使之趋于1。对于压杆,可取对于压杆,可取,式中,式中为杆件的长细比。为杆件的长细比。第9页,共40页,编辑于2022年,星期一迭代收敛性迭代收敛性满应力设计的迭代收敛性,用数学方法证明很困难,实践中主要靠观察收敛过满应力设计的迭代收敛性,用数学方法证明很困难,实践中主要靠观察收敛过程和利用力学判断来解决。程和利用力学判断来解决。正常型结构:性质上接近静定结构,收敛的较快;正常型结构:性质上接近静定结构,收敛的较快;交感型结构:收敛较慢,甚至不能收敛。交感型结构:收敛较慢,甚至不能收敛。第10页,共40页,编辑于2022年,星期一三、分布优化法三、分布优化法 1.基本思路基本

12、思路对现行结构设计方案在各种工况下进行对现行结构设计方案在各种工况下进行(1 1)结构整体分析,得到它的内力分布,)结构整体分析,得到它的内力分布,(2 2)然后在暂时冻结内力的假定下,然后在暂时冻结内力的假定下,(3 3)修改各部分的设计变量,)修改各部分的设计变量,(4 4)再将各部分重新集合起来,得到一个新方案,)再将各部分重新集合起来,得到一个新方案,这样就完成了一次循环。接下去再进行下一次循环,直到收敛。这样就完成了一次循环。接下去再进行下一次循环,直到收敛。结构整体分析结构整体分析冻结内力冻结内力修改各部分的设计变量修改各部分的设计变量重新集合起来重新集合起来第11页,共40页,编

13、辑于2022年,星期一2.2.实现问题实现问题主要取决于约束条件的性质主要取决于约束条件的性质(1 1)局部性约束:)局部性约束:只涉及结构的各个部分。只涉及结构的各个部分。如截面尺寸约束、应力约束、局部稳定约束、局部相对变形约束等。如截面尺寸约束、应力约束、局部稳定约束、局部相对变形约束等。约束条件:约束条件:(2 2)整体性约束:)整体性约束:牵涉到整体结构。如结点位移约束、结构整体稳定性约束、频率约束。牵涉到整体结构。如结点位移约束、结构整体稳定性约束、频率约束。分部优化只能对局部性约束条件适用,如果在约束条件中包含了整体性约束,分部优化只能对局部性约束条件适用,如果在约束条件中包含了整

14、体性约束,则只能进行整体优化,不过可以利用一些分部优化的技巧来简化计算。则只能进行整体优化,不过可以利用一些分部优化的技巧来简化计算。第12页,共40页,编辑于2022年,星期一3.迭代步骤迭代步骤1)给出初始点)给出初始点;2)2)按按进行结构整体分析,求出各部分受力状态;进行结构整体分析,求出各部分受力状态;3)假定内力分布不变,对各部分进行优化,求出新的设计方案假定内力分布不变,对各部分进行优化,求出新的设计方案;若,则停止运算。否则将,则停止运算。否则将作为作为转入步骤转入步骤2 2)4)若)若第13页,共40页,编辑于2022年,星期一5.3 5.3 能量准则法能量准则法 一、基本思

15、路一、基本思路充分发挥材料的贮能(应变能)作用,使所用材料最省或结构最轻。充分发挥材料的贮能(应变能)作用,使所用材料最省或结构最轻。满应变能准则,满应变能准则,应变能密度准则,应变能密度准则,最大总应变能准则。最大总应变能准则。能量准则法能量准则法二、应变能密度准则二、应变能密度准则考虑在一组荷载考虑在一组荷载 作用下(单工况),具有能量约束的结构最轻设计问题。作用下(单工况),具有能量约束的结构最轻设计问题。式中式中 为结构的总应变能,为结构的总应变能,是总应变能的最大允许值,是总应变能的最大允许值 。第14页,共40页,编辑于2022年,星期一能量约束实际上代表对结构的一个广义的刚度要求

16、,是一个综合性的位移约束。能量约束实际上代表对结构的一个广义的刚度要求,是一个综合性的位移约束。最轻解要求此约束达到临界,亦即使最轻解要求此约束达到临界,亦即使 ,这时,这时Langrange函数为函数为 根据根据L-TL-T条件:条件:因为因为可以证明,可以证明,为单元的应变能。为单元的应变能。证明:证明:第15页,共40页,编辑于2022年,星期一证明:证明:证明:证明:第16页,共40页,编辑于2022年,星期一所有各构件的应变能密度相同,等于结构总的应变能密度。所有各构件的应变能密度相同,等于结构总的应变能密度。在结构最轻情况下,在结构最轻情况下,结构总应变能结构总应变能 这样就得到单

17、工况下,满足广义刚度要求这样就得到单工况下,满足广义刚度要求的结构最轻设计的能量准则为:的结构最轻设计的能量准则为:,因而,因而即在最轻情况下,结构所接受的应变能密度最大,即在最轻情况下,结构所接受的应变能密度最大,并且应变能在整个结构上按重量均匀分布。并且应变能在整个结构上按重量均匀分布。第17页,共40页,编辑于2022年,星期一三、应变能密度准则迭代公式三、应变能密度准则迭代公式 根据根据可得可得开方后等式两边均乘以开方后等式两边均乘以,得,得这样便得到这样便得到VenkayyaVenkayya建议的递推公式建议的递推公式第18页,共40页,编辑于2022年,星期一5.3 5.3 齿行法

18、齿行法一、概述一、概述分部优化设计分部优化设计 假定:假定:“在调整时各部分所受内力不变在调整时各部分所受内力不变”这个假定人为地这个假定人为地割断割断了个了个部分部分和和整体整体之间的联系,之间的联系,这样便带来三个后果:这样便带来三个后果:1 1)各部分分别局部优化后,整个结构是不协调的,因而必须多次迭代进行调整;)各部分分别局部优化后,整个结构是不协调的,因而必须多次迭代进行调整;2 2)不排除结构整体的)不排除结构整体的“不可行性不可行性”3)最后得到的结果并不能保证就是结构整体最优的方案。)最后得到的结果并不能保证就是结构整体最优的方案。第19页,共40页,编辑于2022年,星期一分

19、部优化法的优越性分部优化法的优越性 步骤简单步骤简单;收敛迅速;收敛迅速。结构设计中结构设计中局部性约束是大量的局部性约束是大量的;整体性约束是少量的;整体性约束是少量的。(1)齿行法的基本思路:)齿行法的基本思路:分别对待分别对待局部性约束局部性约束和和整体性约束整体性约束,通过,通过分部优化分部优化和和整体调整整体调整两种手段两种手段对结构体进行优化。对结构体进行优化。分部优化分部优化利用满应力(满位移)条件,决定搜索方向。利用满应力(满位移)条件,决定搜索方向。整体调整整体调整及时把设计点拉回到可行域的边界上来。及时把设计点拉回到可行域的边界上来。(2 2)齿行法的基本步骤:)齿行法的基

20、本步骤:分部优化分部优化满应力准则,满位移准则满应力准则,满位移准则;整体调整整体调整利用利用射线步射线步进行可行性调整。进行可行性调整。第20页,共40页,编辑于2022年,星期一二、射线步二、射线步最优解最优解落在可行域边界约束面的交点落在可行域边界约束面的交点 或边界面与目标函数等值或边界面与目标函数等值面面得公切点上。得公切点上。应该沿可行域的边界面搜索寻找最优解。应该沿可行域的边界面搜索寻找最优解。应把当前设计点拉回到当前最严约束面上。应把当前设计点拉回到当前最严约束面上。最简单的一种可行性调整办法最简单的一种可行性调整办法射线步或比例步射线步或比例步 第21页,共40页,编辑于20

21、22年,星期一1)结构响应与方案调整的关系)结构响应与方案调整的关系假设:假设:结构处于线弹性范围内;结构处于线弹性范围内;构件刚度为设计变量的线性齐次函数;构件刚度为设计变量的线性齐次函数;荷载向量荷载向量与设计向量与设计向量无关。无关。若按若按 调整设计方案,则有下述关系成立:调整设计方案,则有下述关系成立:2)3)4)1)第22页,共40页,编辑于2022年,星期一证明如下:证明如下:显然成立,例如绗架单元显然成立,例如绗架单元2.3.4.第23页,共40页,编辑于2022年,星期一3.射线步公式射线步公式当设计点当设计点求出后,需要把它拉回到最严的应力约束面上,求出后,需要把它拉回到最

22、严的应力约束面上,的各应力比的各应力比,令其最大者为,令其最大者为为此,首先求出为此,首先求出可行性调整比例系数可行性调整比例系数,最新设计点最新设计点 则则将沿射线方向落到最严约束面将沿射线方向落到最严约束面上,证明如下:上,证明如下:第24页,共40页,编辑于2022年,星期一把设计点拉到最严的位移约束面上把设计点拉到最严的位移约束面上 第25页,共40页,编辑于2022年,星期一三、只有应力约束时的齿行法三、只有应力约束时的齿行法1)基本思路:)基本思路:只有应力约束时的结构,可以引入如下的优化方法:只有应力约束时的结构,可以引入如下的优化方法:在每产生一个方案之后在每产生一个方案之后射

23、线步射线步使设计点落到当前最严的约束面上使设计点落到当前最严的约束面上满应力步满应力步直至连续二个可行设计点足够接近直至连续二个可行设计点足够接近或它们的目标函数值开始增加为止。或它们的目标函数值开始增加为止。由于这个迭代过程在设计空间的轨迹呈锯齿行,故称为齿行法由于这个迭代过程在设计空间的轨迹呈锯齿行,故称为齿行法 第26页,共40页,编辑于2022年,星期一第27页,共40页,编辑于2022年,星期一2.2.迭代程序迭代程序(1)(1)选取初始点选取初始点(2 2)射线步。)射线步。移到当前最严约束面上移到当前最严约束面上计算计算得目标函数值得目标函数值,并求出新的应力比,并求出新的应力比

24、(3)(3)调整步(满应力步)调整步(满应力步)按满应力准则求下一个设计点按满应力准则求下一个设计点(4 4)迭代收敛条件)迭代收敛条件以上完成一次循环(包括一次射线步,一次调整步),将所得以上完成一次循环(包括一次射线步,一次调整步),将所得作为下一次循环的初始点,如此反复迭代,直到作为下一次循环的初始点,如此反复迭代,直到进行力学分析,求出应力比进行力学分析,求出应力比 ,;或或为止。为止。(5 5)输出最优方案)输出最优方案近似局优解。近似局优解。第28页,共40页,编辑于2022年,星期一四、带有位移约束时的齿行法四、带有位移约束时的齿行法 1.基本思路基本思路 最严约束最严约束 位移

25、约束位移约束 应力约束应力约束 在每产生一个方案之后在每产生一个方案之后射线步射线步使设计点落到当前最严的约束面上使设计点落到当前最严的约束面上(同时考虑应力约束和位移约束同时考虑应力约束和位移约束)满应力步满应力步直至连续二个可行设计点足够接近直至连续二个可行设计点足够接近或它们的目标函数值开始增加为止。或它们的目标函数值开始增加为止。满位移步满位移步 第29页,共40页,编辑于2022年,星期一求设计变量求设计变量 使结构重量使结构重量 并满足约束并满足约束 优化设计问题描述:优化设计问题描述:每一个迭代循环中都用优化准则步得到新的设计点后,每一个迭代循环中都用优化准则步得到新的设计点后,

26、再检验尺寸约束,当某个再检验尺寸约束,当某个方案方案:在处理位移约束和应力约束时,先不管尺寸约束。:在处理位移约束和应力约束时,先不管尺寸约束。时,令时,令。第30页,共40页,编辑于2022年,星期一2.射线步射线步完成了射线步完成了射线步 并有总刚度矩阵并有总刚度矩阵 第31页,共40页,编辑于2022年,星期一3.3.满应力步满应力步 若若射线步射线步把设计点射到把设计点射到应力最严约束面上应力最严约束面上,则用满应力准则来修改设计变量:则用满应力准则来修改设计变量:式中,式中,即,即 ,是控制步长的阻尼系数。是控制步长的阻尼系数。无位移约束无位移约束有位移约束有位移约束 ,取与位移步相

27、同的值。,取与位移步相同的值。第32页,共40页,编辑于2022年,星期一4.4.满位移步满位移步若若射线步射线步把设计点射到把设计点射到位移最严约束面上位移最严约束面上,则用满位移准则来修改设计变量。,则用满位移准则来修改设计变量。这里需要求解一个子规划问题:这里需要求解一个子规划问题:这里沿最严约束搜索,只需考虑这一个约束,所以约束取等号这里沿最严约束搜索,只需考虑这一个约束,所以约束取等号 第33页,共40页,编辑于2022年,星期一构造构造Langrange函数:函数:根据根据K-T条件得条件得(1)(2)将式(将式(1)乘)乘 以并求和,可得以并求和,可得(3)上式中上式中(4)第3

28、4页,共40页,编辑于2022年,星期一结构分析的基本方程为结构分析的基本方程为 左右两端分别对左右两端分别对求偏导:求偏导:左右两端乘以左右两端乘以并求和:并求和:根据杆件刚度与根据杆件刚度与间的间的线性关系线性关系,有,有 证法证法 第35页,共40页,编辑于2022年,星期一将式(将式(4)代入式()代入式(3),并考虑到式(),并考虑到式(2),便可求出),便可求出Langrange乘子乘子(2)(3)(4)(6)(1)(7)第36页,共40页,编辑于2022年,星期一 式中式中为相应于最严约束的工况为相应于最严约束的工况作用下的作用下的号位移。其一阶导数为:号位移。其一阶导数为:=0

29、虚位移原理虚位移原理满位移步公式满位移步公式 满位移步公式满位移步公式 式中式中为工况为工况中杆件的内力,中杆件的内力,为沿位移为沿位移方向的单位力在杆方向的单位力在杆中引起的内力。它们都与设计点中引起的内力。它们都与设计点相应。相应。为控制步长的阻尼系数。为控制步长的阻尼系数。效果较好。效果较好。第37页,共40页,编辑于2022年,星期一当满位移迭代公式中圆括号内为负值时,该式无意义,说明第当满位移迭代公式中圆括号内为负值时,该式无意义,说明第杆截面不能用满位移步调整。这时该杆截面可保持不变,杆截面不能用满位移步调整。这时该杆截面可保持不变,或根据当前最严的应力约束,用满应力准则步修改截面

30、。或根据当前最严的应力约束,用满应力准则步修改截面。第38页,共40页,编辑于2022年,星期一5.5.迭代程序迭代程序1)选定初始方案)选定初始方案进行力学分析,求其应力比和挠度比:进行力学分析,求其应力比和挠度比:,2)射线步)射线步求出所以应力比和挠度比中的最大者求出所以应力比和挠度比中的最大者,将设计点引到当前最严的约束面上,将设计点引到当前最严的约束面上,按下式求应力比按下式求应力比 3 3)调整步)调整步若最严约束是应力约束,走满应力步,则其调整截面为若最严约束是应力约束,走满应力步,则其调整截面为 式中式中 若最严约束是位移约束,并且若最严约束是位移约束,并且,走满位移步,则其调

31、整截面为,走满位移步,则其调整截面为 检查尺寸约束,使各杆检查尺寸约束,使各杆,确定,确定 当当 时,可令时,可令 ,或根据当前最严应力约束确定,或根据当前最严应力约束确定 。第39页,共40页,编辑于2022年,星期一4 4)迭代收敛条件)迭代收敛条件 以上完成一次循环,将所得以上完成一次循环,将所得 作为下一次循环的初始点,作为下一次循环的初始点,如此反复迭代,直至射线步后所得目标函数不再减小,如此反复迭代,直至射线步后所得目标函数不再减小,即即 或满足其他收敛条件为止。或满足其他收敛条件为止。5)输出最优方案,以最后二次射线步后所得点之一作为近似最轻点。即)输出最优方案,以最后二次射线步后所得点之一作为近似最轻点。即或或 第40页,共40页,编辑于2022年,星期一

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 大学资料

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁