【溫馨提示】====【1】設(shè)計(jì)包含CAD圖紙 和 DOC文檔，均可以在線預(yù)覽，所見即所得，，dwg后綴的文件為CAD圖，超高清，可編輯，無任何水印，，充值下載得到【資源目錄】里展示的所有文件======【2】若題目上備注三維，則表示文件里包含三維源文件，由于三維組成零件數(shù)量較多，為保證預(yù)覽的簡潔性，店家將三維文件夾進(jìn)行了打包。三維預(yù)覽圖，均為店主電腦打開軟件進(jìn)行截圖的，保證能夠打開，下載后解壓即可。======【3】特價(jià)促銷，，拼團(tuán)購買，，均有不同程度的打折優(yōu)惠，，詳情可咨詢QQ：1304139763 或者 414951605======【4】 題目最后的備注【JA系列】為店主整理分類的代號，與課題內(nèi)容無關(guān)，請忽視
立體光照成型的注塑模具工藝的綜合模擬
摘要 功能性零部件都需要設(shè)計(jì)驗(yàn)證測試，車間試驗(yàn)，客戶評價(jià)，以及生產(chǎn)計(jì)劃。在小批量生產(chǎn)零件的時(shí)候，通過消除多重步驟，建立了有快速成型形成的注塑模具，這種方法可以保證縮短時(shí)間和節(jié)約成本。這種潛在的一體化由快速成型形成注塑模具的方法已經(jīng)被多次證明是可行的。無論是模具設(shè)計(jì)還是注塑成型的過程中，缺少的是對如何修改這個模具材料和快速成型制造過程的影響有最根本的認(rèn)識。此外，數(shù)字模擬技術(shù)現(xiàn)在已經(jīng)成為模具設(shè)計(jì)工程師和工藝工程師開注塑模具的有用的工具。但目前所有的做常規(guī)注塑模具的模擬包已經(jīng)不再適合這種新型的注塑模具，這主要是因?yàn)槟＞卟牧系某杀咀兓艽?。在本文中，以完成特定的?shù)字模擬注塑液塑造成快速成型模具的綜合方法已經(jīng)發(fā)明出來了，而且還建立了相應(yīng)的模擬系統(tǒng)。通過實(shí)驗(yàn)結(jié)果表明，目前這個方法非常適合處理快速成型模具中的問題。
關(guān)鍵詞 注塑成型,數(shù)字模擬,快速成型
１引言
在注塑成型中，聚合物熔體在高溫和高壓下進(jìn)入模具中。因此，模具的材料需要有足夠的熱性能和機(jī)械性能來經(jīng)受高溫和高壓的塑造循環(huán)。許多研究的焦點(diǎn)都是直接有快速成型形成注塑模具的過程。在生產(chǎn)小批量零件的時(shí)候，通過消除多重步驟，直接由快速成型形成的注塑模具可以保證縮短時(shí)間和節(jié)約成本。這種潛在的有快速成型形成注塑模具的方法已經(jīng)被證明成功了?？焖俪尚湍＞咴谛阅苌鲜怯袆e與傳統(tǒng)的金屬模具。主要差異是導(dǎo)熱性能和彈性模量（剛性）。舉例來說，在立體光照成型模具中的聚合物的導(dǎo)熱率小于鋁制的工具的千分之一。在用快速成型技術(shù)來制造鑄模時(shí)，整個模具設(shè)計(jì)和注塑成型工藝參數(shù)都需要修改和優(yōu)化，傳統(tǒng)的方法是改變徹底的刀具材料．不過，目前還沒有對如何修改這個模具材料的方法有根本的了解．在當(dāng)前的模具中，僅僅改變一些材料的性能是不能得到一個合理的結(jié)果的。同樣，使用傳統(tǒng)方法的時(shí)候，實(shí)際生產(chǎn)的零件也會有出先次品。因此，研究出一個快速成型過程，材料和注塑模具之間的互動關(guān)系是非?；鸺钡?。這樣就可以確定模具設(shè)計(jì)標(biāo)準(zhǔn)和快速模具的注塑的技術(shù)。
此外，計(jì)算機(jī)模擬是一種預(yù)測模塑件的質(zhì)量的有效的方法。目前，商用仿真軟件包已經(jīng)成為模具設(shè)計(jì)師和工藝工程師在注塑過程中例行性的工具。不幸的是，目前常規(guī)注塑成型的模擬程序已經(jīng)不再適用于這個快速成型模具，因?yàn)樗鼧O大的需要不同的刀具材料。例如，利用現(xiàn)在的仿真軟件在鋁和立體光照模具之間做個實(shí)驗(yàn)比較一下，雖然鋁模具模擬植的部分失真是合理的，但是結(jié)果是不可以接受的，因?yàn)檎`差超過了百分之五十。在注塑成型中，失真主要是由于塑料零件的收縮和翹曲，模具也是一樣的。對于通常模具，失真的主要因素是塑料件的收縮和翹曲，這個在目前的模擬中能測試準(zhǔn)確。但是對于快速成型模具，潛在的失真會更多，在當(dāng)前的測試中，其中就會有些失真會被忽視。例如，用一個簡單的三步驟模擬分析模具變形的時(shí)候，就會出現(xiàn)很多偏差。
在本文中，基于以上分析，一個新的快速成型模具的仿真系統(tǒng)已經(jīng)開發(fā)出來了。擬議制度著重于預(yù)測部分失真，主要是用與預(yù)測快速成型模具的缺陷。先進(jìn)的仿真系統(tǒng)可以用于預(yù)測快速成型模具設(shè)計(jì)和工藝是否最合理。我們的仿真系統(tǒng)已經(jīng)被我們的實(shí)驗(yàn)證明是沒有錯誤的。
雖然有很多材料可以用于快速成型技術(shù)，但是我們還是專注于利用立體光照模具的技術(shù)來制造聚合物模具．立體光照成型的過程是利用激光能量一層一層建立零件的部分。使用立體光照則可以體現(xiàn)出雙方在快速成型工業(yè)的商業(yè)優(yōu)勢，而且在以后也可以生產(chǎn)出準(zhǔn)確的，高品質(zhì)的零部件。直到最近，立體光照主要是用于建立物理模型，為了檢查視覺效果，僅僅只利用了它的一點(diǎn)點(diǎn)功能。不過，新一代的立體光照的光改善了立體化，機(jī)械性能，熱學(xué)性能，所以它可以更好的應(yīng)用于實(shí)際的模具中。
2 綜合仿真的成型過程
2.1 方法
為了在注塑成型過程中模擬立體光照模具的功能，反復(fù)的試驗(yàn)中得到了一個方法。不同的軟件組已經(jīng)開發(fā)出來了，而且也已經(jīng)做到了這一點(diǎn)。主要的假設(shè)是，溫度和負(fù)載邊界條件造成立體光照模具的扭曲，仿真步驟如下：
?。辈糠謳缀文Ｐ蛣t作為一個實(shí)體模型，這將通過流量分析軟件包被翻譯到一個文件中。
２模擬光聚合物模具中熔融體填充的過程，然后輸出溫度和壓力的資料。
?。吃谇耙徊将@得了熱負(fù)荷和邊界條件，然后對光模具進(jìn)行結(jié)構(gòu)分析，其中失真的計(jì)算是在該注塑過程中進(jìn)行的。
４如果模具的扭曲收斂了，那么直接進(jìn)行下一步．否則，扭曲的型腔（改動扭曲后的型腔的尺寸）返回第二個步驟，以熔體形式模擬注入扭曲的模具中。
?。等缓笞⑸涑尚土慵氖湛s和翹曲模擬就開始應(yīng)用了，算出該成型零件最終的扭曲部分．
上述的模擬流動中，基本上是三個仿真模塊。
2.2充型模擬的熔體
2.2.1數(shù)字建模
計(jì)算機(jī)仿真技術(shù)已經(jīng)能成功的預(yù)測到在極其復(fù)雜的幾何形狀下的填充情況。然而，目前大多數(shù)字模擬是基于一種混合有限元和有限差的中性平面上的。模擬軟件包的應(yīng)用過程基于這一模型說明圖１。然而，不同與ＣＡＤ系統(tǒng)中模具設(shè)計(jì)中的表面／實(shí)體模型，這里所謂的中性平面（如圖所示，圖１Ｂ）是一個假想的在中間型腔中有距離和方向的一個平面，這個平面可能會在應(yīng)用的過程中帶來很大的不便。舉例來說，模具表面常用于目前的快速成型系統(tǒng)中（通常是ＳＴＬ格式），所以當(dāng)用模擬軟件包的時(shí)候，第二次建模是不可避免的。那是因?yàn)槟Ｐ驮诳焖俪尚拖到y(tǒng)和仿真系統(tǒng)中是不一樣的?？紤]到這些缺點(diǎn)，在模擬系統(tǒng)中，型腔的表面將以基準(zhǔn)面來引入，而不是中性平面。
根據(jù)以往的調(diào)查，流量和溫度場的方程式可以寫為：
X,Y是中性平面坐標(biāo)系中的兩個平面，Ｚ是高度坐標(biāo)，Ｕ，Ｖ，Ｗ是Ｘ，Ｙ，Ｚ方向上的速度．Ｕ，Ｖ是整體的平均厚度，η, ρ,CP (T), K(T)分別表示聚合物的粘性，密度，周期熱，熱導(dǎo)率。
圖１ Ａ－Ｄ是中性平面的模擬程序．Ａ是３維表面模型，Ｂ是中性平面模型，Ｃ是網(wǎng)狀的平面模型，Ｄ是最后的模擬結(jié)果
此外，在高度方向上的邊界條件的誤差可以表示為：
正如圖２中的Ａ中表示，TW 是恒壁溫度.結(jié)合方程１－４和方程５－６，表明了u, v, T, P在Ｚ坐標(biāo)上面應(yīng)該是對稱的，因此在上半個高度中的平均u, v應(yīng)該和整個高度中的平均u, v是一樣的。根據(jù)這個特點(diǎn)，我們可以把整個型腔在上下高度上分為兩個部分，正如圖２Ｂ中的第一部分和第二部分。同時(shí)，型腔（如圖２Ｂ）表面產(chǎn)生的三角有限元將替代了中性平面（如圖２Ａ）。因此，在高度方向上的有限元誤差僅僅限于型腔表面，正如圖２Ｂ所示，高度上的誤差將從０到Ｂ。這是中性平面上的單一性。此外，從圖２Ａ到圖２Ｂ，坐標(biāo)也隨之改變了。為了配合上述調(diào)整，方程仍是用方程１－４。然而，原來的邊界條件高度方向則改寫為：
與此同時(shí)，為了保持在同一坐標(biāo)（７）上的兩部分能夠流動，那么更多的邊界條件必須滿足Ｚ＝Ｂ。
下標(biāo)I和II則分別代表第一部分和第二部分的參數(shù)．Cm-I 和Cm-II 則表示在填充階段中分開的兩個表面上的自由移動的熔融線。
應(yīng)該指出的是，方程９與１０和方程７與８不同，９和１０在數(shù)字模擬過程中將變的更難，主要原因是以下幾點(diǎn)：
?。蓖粋€斷層的表面都已經(jīng)都已經(jīng)有著特殊的網(wǎng)格，這將導(dǎo)致同一層上的獨(dú)特的格局．因此，在比較兩個熔接口的時(shí)候，應(yīng)該計(jì)算出各自的u, v, T, P。
?。惨?yàn)閮蓚€部分都有各自的流道通向節(jié)點(diǎn)Ａ和節(jié)點(diǎn)Ｃ（如圖２Ｂ所示）．在同一段中，有可能兩個都充滿，也有可能一個滿，一個空．這兩個情況應(yīng)該分開處理，應(yīng)該平均流動，使后者也分配到流動。
?。尺@意味著在前線熔合處出現(xiàn)一點(diǎn)點(diǎn)小的誤差是可以允許的．通過控制時(shí)間和選擇更好的位置來控制前線熔合節(jié)點(diǎn)。
?。疵總€流場的邊界都擴(kuò)張到熔線前線，所以核查方程１０是否準(zhǔn)確是相當(dāng)重要的。
?。佃b于上述分析，在同一個節(jié)點(diǎn)處的物理參數(shù)應(yīng)該加以比較和調(diào)整。所以在進(jìn)行模擬之前，描述同一節(jié)點(diǎn)有限元的信息應(yīng)該準(zhǔn)備好，也就是說，匹配的原理應(yīng)該先預(yù)備好。
圖２ Ａ－Ｂ表明表面模型中的中性平面Ｂ的高度方向Ａ上的邊界條件
2.2.2數(shù)字模擬
壓力場．在建模中，粘度 η是由于熔提的剪切速率，溫度和壓力引起的性能．剪切變稀后，這就代表一個跨越式的模式，例如：
其中對應(yīng)于冪律指數(shù)，τ的特點(diǎn)是在在牛頓和冪律漸近極限之間的剪應(yīng)力過渡區(qū)。無論在溫度還是壓力指數(shù)上，η0(T, P)都可以有合理的表示，詳情如下：
方程11和12構(gòu)成了一個五個常數(shù)，可以代表粘度，而且通過粘度的剪切速率的計(jì)算可以得到：
根據(jù)上述情況，通過方程1—4，我們可以推斷出一下充氣壓力方程：
其中S是由計(jì)算出來的。運(yùn)用伽遼金方法，對壓力的有限元方程推導(dǎo)為：
其中l(wèi)是所有要素的的導(dǎo)線，包括節(jié)點(diǎn)N，而且其中i和j代表此處的N節(jié)點(diǎn)的數(shù)目，的計(jì)算方法如下：
其中代表三角有限元，而代表有限元中的壓力。
溫度場中，為了確定高度方向上的誤差，應(yīng)該在模具表面上分為一層一層的三角有限元的網(wǎng)格。左邊的能量方程4可以表示為：
其中代表每一層N節(jié)點(diǎn)上的溫度。熱傳導(dǎo)的計(jì)算方法是：
其中l(wèi)是所有要素，包括節(jié)點(diǎn)N，而且i和j分別代表此處的N節(jié)點(diǎn)個數(shù)。
對流項(xiàng)的計(jì)算方法是：
當(dāng)是粘性熱時(shí)，計(jì)算方法是：
把方程17—20帶入方程4，溫度方程變?yōu)椋?
2.3 模具結(jié)構(gòu)分析
結(jié)構(gòu)分析的目的是預(yù)測在填充過程中，模具由于熱和機(jī)械壓力而產(chǎn)生的變形。這個模型是基于一個三維熱邊界元法。邊界元法是比較適合這個應(yīng)用的，因?yàn)橹挥凶冃蔚哪＞弑砻娌庞羞@樣的信息。此外，邊界元法有一個優(yōu)點(diǎn)，那就是在計(jì)算變形的模具的時(shí)候，它的計(jì)算是不會白費(fèi)的。
模具在所受載荷超過彈性范圍的時(shí)候會產(chǎn)生應(yīng)力。因此，在決定模具變形的時(shí)候，模具材料是一個基準(zhǔn)。模具的熱性能和力學(xué)性能是各向同性的，而且溫度也是獨(dú)立的。
盡管這個過程是循環(huán)的，但是相同時(shí)間的溫度和熱流都是可以用于計(jì)算模具變形的．通常情況下，在模具里面每個瞬間溫度都局限于型腔的表面和噴嘴的頂端。在觀察距離的時(shí)候，瞬間的衰減變化是很微笑的，小于２．５毫米．這說明在模具的噴嘴處的變形是很小的，因此，忽略這個影響也是合理的．穩(wěn)態(tài)溫度場滿足拉普拉斯方程?2T = 0的邊界條件。至于機(jī)械邊界條件，型腔表面受到熔體的壓力，模具的表面會連接到工作臺上的，而其他的外部表面將會假設(shè)是自由的.熱邊界的推導(dǎo)方程１０是大家都知道的，這是由于：
其中uk, pk和Ｔ分別是位移，牽引力和溫度。α, ν是代表材料的膨脹系數(shù)和泊松比。Ulk是在ＸＹ方向上基本的位移。在一個三維空間中，各向同性彈性區(qū)域中，由一個單元產(chǎn)生的負(fù)荷主要集中在xl方向上，它是以下面的形式產(chǎn)生的：
其中δlk是Kronecker三角函數(shù)，μ是該模具材料的剪切模量。Plk的基本收縮都是在模具表面的每個Ｎ節(jié)點(diǎn)處測量的，可以表示為：
整個Ｎ將分散在模具的表面上，轉(zhuǎn)變?yōu)榉匠蹋玻玻?
其中Γn是指在這個區(qū)域上的表面成分。
把恰當(dāng)?shù)木€性函數(shù)代入方程２５，得到的線性邊界方程就是模具的方程．這個方程適用于每個離散的模具表面，從而組合成線性方程組，其中Ｎ是節(jié)點(diǎn)的總數(shù)。每個節(jié)點(diǎn)有八個相關(guān)數(shù)量，三個位移組成部分，三個牽引組成部分，還有溫度和熱流量。在穩(wěn)態(tài)熱模型中，每個節(jié)點(diǎn)處的溫度和磁場是已知的，余下的６個量中，三個必須是已知的。此外，在若干個節(jié)點(diǎn)處的位移值的方程必須消除剛體運(yùn)動和剛體自轉(zhuǎn)的奇異系統(tǒng)。由此產(chǎn)生的系統(tǒng)方程式是一個集合起來的綜合矩陣，它可以為有限元方法求解。
基于方程１２的注塑假設(shè)，下面將給出元件的應(yīng)力和應(yīng)變：
該偏元件的應(yīng)力和應(yīng)變分別是：
用類似的方法可以預(yù)測在回火玻璃中的殘余應(yīng)力了。以積分的形式在平面上分析粘性和彈性結(jié)構(gòu)關(guān)系時(shí)，可以表示為以下公式：
其中G1是材料的的剪切模量。擴(kuò)張的應(yīng)變的情況如下：
其中Ｋ是材料體積的彈性模量，α和θ的定義是：
如果α(t) = α0，那么方程２７到方程２９的結(jié)果則為：
同樣的，利用方程３１到方程２８消除應(yīng)變εxx(z, t)，得到：
利用拉普拉斯變化方程３２，輔助系數(shù)R(ξ)由下面的方程得出：
利用上述方程３３，并簡化在模具中的應(yīng)力和應(yīng)變的形式，那么注塑中殘余的應(yīng)力在冷卻階段中，由下面的方程獲得：
方程３４可以通過梯形正交被解決。由于材料的時(shí)間在快速的變化，所以需要一個準(zhǔn)數(shù)控程序來檢測。輔助模量是檢測數(shù)控梯形的規(guī)則。
關(guān)于翹曲分析，節(jié)點(diǎn)位移和曲率將以殼單元表達(dá)為：
其中[ k ]單元剛度矩陣，[Be]是衍生算子矩陣，ttt9nbx是位移，{re}是 負(fù)載單元，可以由下面的方程得出：
使用完整的三維有限元分析法的好處就是可以準(zhǔn)確知道翹曲的結(jié)果。但是，當(dāng)零件的形狀很復(fù)雜的時(shí)候，它也是相當(dāng)麻煩的。在本文中，在殼體理論基礎(chǔ)上介紹了一種二維有限元分析方法。這種方法被大量使用是因?yàn)榇蠖鄶?shù)注塑模具的零件都有一些部分幾何的厚度遠(yuǎn)遠(yuǎn)小于其他部分。因此，那些部分則可以被作為一個集會的單元來預(yù)測翹曲。每三個節(jié)點(diǎn)殼單元組合成一個恒應(yīng)變?nèi)菃卧鸵粋€離散克?；舴蛉窃?，如圖３所示，因此翹曲可以分為平面伸展變形ＣＳＴ和板彎曲變形ＤＫＴ。并相應(yīng)的以單元剛度矩陣來描述翹曲的拉伸剛度矩陣和彎曲剛度矩陣。
圖３ a-c是殼單元在局部坐標(biāo)系統(tǒng)里的變形分解．a(chǎn)是平面伸展元素，b是平面彎曲元素，c是殼單元
三 實(shí)驗(yàn)驗(yàn)證
對提出的模型進(jìn)行了評定和發(fā)展，最后核查是非常重要的。從模型模擬中得到的扭曲數(shù)據(jù)將和文獻(xiàn)８中的立體光照模具數(shù)據(jù)比較。如圖４所示，有一個注塑尺寸36 × 36 × 6毫米和實(shí)驗(yàn)數(shù)據(jù)中是相同的。薄壁和加強(qiáng)筋的厚度都是1.5毫米，這個注塑材料是聚丙烯。注塑機(jī)的型號是ARGURY Hydronica320-210-750，它的工藝參數(shù)是，熔解溫度是２５０度，模具溫度是３０度，注塑壓力是１３.７９帕，保壓時(shí)間是３秒，冷卻時(shí)間是４８秒。立體光照模具材料使用杜邦SOMOSTM６１１０樹脂，能抵御高達(dá)３００度的高溫。如上所述，熱傳導(dǎo)是區(qū)分立體光照模具和傳統(tǒng)模具的一個重要因素。模具中的熱量轉(zhuǎn)移會產(chǎn)生溫度的不均勻分布，所以導(dǎo)致了成型零件的翹曲．立體光照成型模具的周期是可以預(yù)測的。以高的熱傳導(dǎo)率金屬為背面做的薄殼立體光照模具將會增加自身的熱傳導(dǎo)率。
圖４ 模型腔
圖５ 不同的熱傳導(dǎo)率下，在Ｘ方向上的扭曲失真比較．實(shí)驗(yàn)值，三步走和常規(guī)都是指最后的實(shí)驗(yàn)結(jié)果．常規(guī)是指實(shí)驗(yàn)中最好的結(jié)果．三步走步驟的模擬過程分別與傳統(tǒng)的注塑成型相似
圖６ 在不同的熱傳導(dǎo)率下，在Ｙ方向上的扭曲失真比較
圖７ 在不同熱傳導(dǎo)率下，在Ｚ方向上扭曲失真比較
圖８ 不同熱傳導(dǎo)率下各個捻度變量的比較
對于這個部分，扭曲包括三個方向上的位移和捻度（兩個最初的平行邊的夾角的誤差）．如圖５到圖８，實(shí)驗(yàn)結(jié)果表明，這些數(shù)值也包括通過傳統(tǒng)注塑模具模擬系統(tǒng)預(yù)測的扭曲值和報(bào)道［３］中的三步驟。
４結(jié)論
本文介紹了一個綜合模擬的快速成型模具的方法，并且建立了相應(yīng)的仿真系統(tǒng)。為了驗(yàn)證這個系統(tǒng)，實(shí)驗(yàn)還進(jìn)行了快速焊接立體光照成型模具。
很明顯，立體光照模具也會出現(xiàn)傳統(tǒng)的注塑模具模擬軟件一樣的故障．假設(shè)由于注射中的溫度和負(fù)載荷引起了扭曲．那么用三步驟完成的話，結(jié)果也會出現(xiàn)比較多的誤差。不過更先進(jìn)的模型會使結(jié)果更接近與實(shí)驗(yàn)。
立體光照模具改進(jìn)了熱傳導(dǎo)率極大的增加了零件質(zhì)量．由于溫度比壓力（負(fù)載）對模具的影響更大，所以改進(jìn)立體光照模具的熱傳導(dǎo)率可以更顯著的提高零件質(zhì)量。
無論零件多么復(fù)雜，快速成型技術(shù)可以使人們造型更快，更便捷，更便宜．在快速成型穩(wěn)步發(fā)展的基礎(chǔ)上，快速制造也將隨之而來，并且需要更多的精確工具來確定工藝過程的參數(shù)．現(xiàn)行的模擬工具不能滿足研究者研究模具相對的變化。正如本文中所述，對于一個綜合模型來說，要預(yù)測最后零件質(zhì)量是相當(dāng)重要的。在不久的將來，我們期待看到通過快速成型擴(kuò)展到快速模具制造的模擬程序。
參考文獻(xiàn)
[1] Wang KK (1980) System approach to injection molding process. Polym-Plast Technol Eng 14(1):75–93.
[2] Shelesh-Nezhad K, Siores E (1997) Intelligent system for plastic injection molding process design. J Mater Process Technol 63(1–3):458–462.
[3] Aluru R, Keefe M, Advani S (2001) Simulation of injection molding into rapid-prototyped molds. Rapid Prototyping J 7(1):42–51.
[4] Shen SF (1984) Simulation of polymeric flows in the injection molding process. Int J Numer Methods Fluids 4(2):171–184.
[5] Agassant JF, Alles H, Philipon S, Vincent M (1988) Experimental and theoretical study of the injection molding of thermoplastic materials.Polym Eng Sci 28(7):460–468.
[6] Chiang HH, Hieber CA, Wang KK (1991) A unified simulation of the filling and post-filling stages in injection molding. Part I: formulation.Polym Eng Sci 31(2):116–124.
[7] Zhou H, Li D (2001) A numerical simulation of the filling stage in injection molding based on a surface model. Adv Polym Technol 20(2):125–131.
[8] Himasekhar K, Lottey J, Wang KK (1992) CAE of mold cooling in injection molding using a three-dimensional numerical simulation. J EngInd Trans ASME 114(2):213–221.
[9] Tang LQ, Pochiraju K, Chassapis C, Manoochehri S (1998) Computeraided optimization approach for the design of injection mold cooling systems. J Mech Des, Trans ASME 120(2):165–174.
[10] Rizzo FJ, Shippy DJ (1977) An advanced boundary integral equation method for three-dimensional thermoelasticity. Int J Numer MethodsEng 11:1753–1768.
[11] Hartmann F (1980) Computing the C-matrix in non-smooth boundary points. In: New developments in boundary element methods, CML Publications,Southampton, pp 367–379.
[12] Chen X, Lama YC, Li DQ (2000) Analysis of thermal residual stress in plastic injection molding. J Mater Process Technol 101(1):275–280.
[13] Lee EH, Rogers TG (1960) Solution of viscoelastic stress analysis problems using measured creep or relaxation function. J Appl Mech 30(1):127–134.
[14] Li Y (1997) Studies in direct tooling using stereolithography. Dissertation, University of Delaware, Newark, DE.